pdd_ns_supg/include/pdh_ns_supg.h

Go to the documentation of this file.

/************************************************************************
File pdh_ns_supg.h - problem module's types and functions
 BY ASSUMPTION NOT CALLED BY OTHER MODULES' FUNCTIONS
 (functions called by other modules are in pdh_ns_supg_weakform.h)

Contains problem module defines (see below)


Contains declarations of routines:
 BY ASSUMPTION NOT CALLED BY OTHER MODULES' FUNCTIONS
 (functions called by other modules are in pdh_ns_supg_problem.h)
  pdr_ns_supg_time_integration - time integration driver (in time_integration)
  pdr_ns_supg_ZZ_error - to compute estimated norm of error based on 
        recovered first derivatives - the notorious ZZ error estimate (in adapt) 
  pdr_ns_supg_err_indi_ZZ - to return error indicator for an element,
        based on ZZ first derivative recovery 
  pdr_ns_supg_err_indi_explicit - to return error indicator for an element,
        based on explicit formulae 
  pdr_ns_supg_adapt - to enforce adaptation strategy for ns problem (in adapt)
  pdr_ns_supg_dump_data - dump data to files (in input_output)
  pdr_ns_supg_write_paraview - to write graphics data to file (in input_output)

utilities:
  pdr_ns_supg_post_process
  pdr_ns_supg_write_profile
  pdr_ns_supg_initial_condition

------------------------------
History:
    initial version - Krzysztof Banas
    2011    - Przemyslaw Plaszewski (pplaszew@agh.edu.pl)
    2011    - Aleksander Siwek (Aleksander.Siwek@agh.edu.pl)
        2012    - Krzysztof Banas (pobanas@cyf-kr.edu.pl)
*************************************************************************/

#ifndef PDH_NS_SUPG
#define PDH_NS_SUPG

#include <stdio.h>

/* problem dependent interface with the PDEs  */
#include "pdh_intf.h"

/* types and functions related to problem structures */
#include "pdh_ns_supg_problem.h" 
// bc and material header files are included in problem header files

//#define TURBULENTFLOW
#ifdef TURBULENTFLOW
#include "../../pdd_turbulent/supg.h"
void GetLocalViscosity(int elementId,int GaussPointId, double* viscLocal,int problem_id);
#endif //TURBULENTFLOW

#ifdef __cplusplus
extern "C" 
{
#endif

/**************************************/
/* DEFINES                            */
/**************************************/
/* Rules:
/* - always uppercase */
/* - name starts with PDC_ */

#define PDC_NS_SUPG_ID  1

/**************************************/
/* TYPES                              */
/**************************************/
/* Rules:
/* - type name starts always witn pdt_ */


/**************************************/
/* GLOBAL VARIABLES                   */
/**************************************/
/* Rules:
/* - name always begins with pdv_ */
/* - constants always uppercase and start with PDC_ */

extern double pdv_ns_supg_timer_all;
extern double pdv_ns_supg_timer_pdr_comp_el_stiff_mat;
extern double pdv_ns_supg_timer_pdr_comp_fa_stiff_mat;


// ID of the current problem
extern int pdv_ns_supg_current_problem_id;  /* ID of the current problem */
// problem structure for ns_supg module
extern pdt_ns_supg_problem pdv_ns_supg_problem;

#ifdef PARALLEL
  extern int pdv_exchange_table_index; // for a single problem
#endif

/**************************************/
/* INTERNAL PROCEDURES                */
/**************************************/
/* Rules:
/* - name always begins with pdr_ */
/* - argument names start uppercase */


extern void pdr_ns_supg_time_integration(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output);


extern int pdr_ns_supg_adapt(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output);

double pdr_ns_supg_ZZ_error(        
       /* returns  - Zienkiewicz-Zhu error for the whole mesh */
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output
);

extern double pdr_ns_supg_err_indi_ZZ(      
                        /* returns error indicator for an element */
  int Problem_id,   /* in: data structure to be used  */
  int El    /* in: element number */
    );

extern double pdr_ns_supg_err_indi_explicit(    
                        /* returns error indicator for an element */
  int Problem_id,   /* in: data structure to be used  */
  int El    /* in: element number */
    );

extern int pdr_ns_supg_dump_data(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

extern int pdr_ns_supg_write_paraview(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output
  );


  // IN: MAIN/PDS_NS_SUPG_UTIL.C

double pdr_ns_supg_post_process(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

int pdr_ns_supg_write_profile(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

double pdr_ns_supg_initial_condition(
  int Field_id, // field_id - each problem should know its field id
  double *Coor,   // point coordinates
  int Sol_comp_id // solution component
);


#ifdef __cplusplus
}
#endif

#endif

---