#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <assert.h>
#include "pdh_intf.h"
#include "pdh_control_intf.h"
#include "mmh_intf.h"
#include "aph_intf.h"
#include "uth_intf.h"
#include "../include/pdh_heat_problem.h"
#include "../include/pdh_heat_weakform.h"
#include "../include/pdh_heat.h"
Defines | |
#define | PDC_HEAT_NUM_DOFS_MAX 18 |
Functions | |
int | pdr_heat_get_velocity_at_point (int Problem_id, int El_id, double *X_loc, double *Base_phi, double *Base_dphix, double *Base_dphiy, double *Base_dphiz, double *Velocity, double *DVel_dx, double *DVel_dy, double *DVel_dz) |
int | pdr_heat_select_el_coeff_vect (int Problem_id, int *Coeff_vect_ind) |
double * | pdr_heat_select_el_coeff (int Problem_id, double **Mval, double **Axx, double **Axy, double **Axz, double **Ayx, double **Ayy, double **Ayz, double **Azx, double **Azy, double **Azz, double **Bx, double **By, double **Bz, double **Tx, double **Ty, double **Tz, double **Cval, double **Lval, double **Qx, double **Qy, double **Qz, double **Sval) |
int | pdr_heat_el_coeff (int Problem_id, int El_id, int Mat_num, double Hsize, int Pdeg, double *X_loc, double *Base_phi, double *Base_dphix, double *Base_dphiy, double *Base_dphiz, double *Xcoor, double *Uk_val, double *Uk_x, double *Uk_y, double *Uk_z, double *Un_val, double *Un_x, double *Un_y, double *Un_z, double *Mval, double *Axx, double *Axy, double *Axz, double *Ayx, double *Ayy, double *Ayz, double *Azx, double *Azy, double *Azz, double *Bx, double *By, double *Bz, double *Tx, double *Ty, double *Tz, double *Cval, double *Lval, double *Qx, double *Qy, double *Qz, double *Sval, double *Velocity, double Thermal_conductivity, double Density_times_specific_heat, double Delta_t, double Implicitness_coeff) |
int | pdr_heat_comp_el_stiff_mat (int Problem_id, int El_id, int Comp_sm, int Pdeg_in, int *Nrdofs_loc, double *Stiff_mat, double *Rhs_vect, char *Rewr_dofs) |
int | pdr_heat_comp_fa_stiff_mat (int Problem_id, int Fa_id, int Comp_sm, int Pdeg_in, int *Nrdofs_loc, double *Stiff_mat, double *Rhs_vect, char *Rewr_dofs) |
int | pdr_heat_get_temperature_at_point (int Problem_id, int El_id, double *X_loc, double *Base_phi, double *Base_dphix, double *Base_dphiy, double *Base_dphiz, double *Temp, double *DTemp_dx, double *DTemp_dy, double *DTemp_dz) |
int | pdr_heat_get_dtdt_at_point (int Problem_id, int El_id, double *X_loc, double *Base_phi, double *Base_dphix, double *Base_dphiy, double *Base_dphiz, double *Temp, double *Ddtdt_dx, double *Ddtdt_dy, double *Ddtdt_dz) |
int | pdr_heat_compute_CFL (int Problem_id, double *CFL_min_p, double *CFL_max_p, double *CFL_ave_p) |
#define PDC_HEAT_NUM_DOFS_MAX 18 |
int pdr_heat_comp_el_stiff_mat | ( | int | Problem_id, | |
int | El_id, | |||
int | Comp_sm, | |||
int | Pdeg_in, | |||
int * | Nrdofs_loc, | |||
double * | Stiff_mat, | |||
double * | Rhs_vect, | |||
char * | Rewr_dofs | |||
) |
----------------------------------------------------------- pdr_heat_comp_el_stiff_mat - to construct a stiffness matrix and a load vector for an element ------------------------------------------------------------
!! coarse element dofs should be supplied by calling routine !!!
int pdr_heat_comp_fa_stiff_mat | ( | int | Problem_id, | |
int | Fa_id, | |||
int | Comp_sm, | |||
int | Pdeg_in, | |||
int * | Nrdofs_loc, | |||
double * | Stiff_mat, | |||
double * | Rhs_vect, | |||
char * | Rewr_dofs | |||
) |
----------------------------------------------------------- pdr_heat_comp_fa_stiff_mat - to construct a stiffness matrix and a load vector for a face ------------------------------------------------------------
DIAGNOSTICS - uncomment !
-------------------------------------------------------------------!
--------------- THERMAL BOUNDARY CONDITIONS -----------------------!
-------------------------------------------------------------------!
------------------------------------------------------------------!
-------------------------- ISOTHERMAL ---------------------------!
DIAGNOSTICS - uncomment !
DIAGNOSTICS - uncomment !
DIAGNOSTICS - uncomment !
DIAGNOSTICS - uncomment !
------------------------------------------------------------------!
-------------------------- HEAT FLUX (NEUMANN)-------------------!
------------------------------------------------------------------!
-------------------------- OUTFLOW -------------------!
------------------------------------------------------------------!
-------------------- RADCONV (NEUMANN)----------------------------!
------------------------------------------------------------------!
----------------- GOLDAK HEAT SOURCE (NEUMANN)--------------------!
note: it's position is updated elsewhere - in time integration !
------------------------------------------------------------------!
----------------- CONTACT BC --------------------!
note: it's position is updated elsewhere - in time integration !
------------------------------------------------------------------!
----------------- END OF CONTACT BC --------------------!
int pdr_heat_compute_CFL | ( | int | Problem_id, | |
double * | CFL_min_p, | |||
double * | CFL_max_p, | |||
double * | CFL_ave_p | |||
) |
int pdr_heat_el_coeff | ( | int | Problem_id, | |
int | Elem, | |||
int | Mat_num, | |||
double | Hsize, | |||
int | Pdeg, | |||
double * | X_loc, | |||
double * | Base_phi, | |||
double * | Base_dphix, | |||
double * | Base_dphiy, | |||
double * | Base_dphiz, | |||
double * | Xcoor, | |||
double * | Uk_val, | |||
double * | Uk_x, | |||
double * | Uk_y, | |||
double * | Uk_z, | |||
double * | Un_val, | |||
double * | Un_x, | |||
double * | Un_y, | |||
double * | Un_z, | |||
double * | Mval, | |||
double * | Axx, | |||
double * | Axy, | |||
double * | Axz, | |||
double * | Ayx, | |||
double * | Ayy, | |||
double * | Ayz, | |||
double * | Azx, | |||
double * | Azy, | |||
double * | Azz, | |||
double * | Bx, | |||
double * | By, | |||
double * | Bz, | |||
double * | Tx, | |||
double * | Ty, | |||
double * | Tz, | |||
double * | Cval, | |||
double * | Lval, | |||
double * | Qx, | |||
double * | Qy, | |||
double * | Qz, | |||
double * | Sval, | |||
double * | Velocity, | |||
double | Thermal_conductivity, | |||
double | Density_times_specific_heat, | |||
double | Delta_t, | |||
double | Implicitness_coeff | |||
) |
-------------------------------------------------------- pdr_heat_el_coeff - to return coefficients for element integrals ----------------------------------------------------------
----------------------------------------------------------------------!
------------ CALCULATE ELEM. "SIZE" (h_k) AT GAUSS POINT----------- --!
----------------------------------------------------------------------!
----------------------------------------------------------------------!
------------------ CALCULATE STABILIZATION COEFFS ----------------- --!
----------------------------------------------------------------------!
int pdr_heat_get_dtdt_at_point | ( | int | Problem_id, | |
int | El_id, | |||
double * | X_loc, | |||
double * | Base_phi, | |||
double * | Base_dphix, | |||
double * | Base_dphiy, | |||
double * | Base_dphiz, | |||
double * | dtdt, | |||
double * | Ddtdt_dx, | |||
double * | Ddtdt_dy, | |||
double * | Ddtdt_dz | |||
) |
----------------------------------------------------------- pdr_heat_get_dtdt_at_point - to provide the eqivalent specific heat and its gradient at a particular point with local coordinates within an element MODULE PROVIDES IMPLEMENTATION FOR ALL OTHER MODULES (in pds_heat_weakform.c) ------------------------------------------------------------
int pdr_heat_get_temperature_at_point | ( | int | Problem_id, | |
int | El_id, | |||
double * | X_loc, | |||
double * | Base_phi, | |||
double * | Base_dphix, | |||
double * | Base_dphiy, | |||
double * | Base_dphiz, | |||
double * | Temp, | |||
double * | DTemp_dx, | |||
double * | DTemp_dy, | |||
double * | DTemp_dz | |||
) |
----------------------------------------------------------- pdr_heat_get_temperature_at_point - to provide the temperature and its gradient at a particular point with local coordinates within an element MODULE PROVIDES IMPLEMENTATION FOR ALL OTHER MODULES (in pds_heat_weakform.c) ------------------------------------------------------------
int pdr_heat_get_velocity_at_point | ( | int | Problem_id, | |
int | El_id, | |||
double * | X_loc, | |||
double * | Base_phi, | |||
double * | Base_dphix, | |||
double * | Base_dphiy, | |||
double * | Base_dphiz, | |||
double * | Velocity, | |||
double * | DVel_dx, | |||
double * | DVel_dy, | |||
double * | DVel_dz | |||
) |
double* pdr_heat_select_el_coeff | ( | int | Problem_id, | |
double ** | Mval, | |||
double ** | Axx, | |||
double ** | Axy, | |||
double ** | Axz, | |||
double ** | Ayx, | |||
double ** | Ayy, | |||
double ** | Ayz, | |||
double ** | Azx, | |||
double ** | Azy, | |||
double ** | Azz, | |||
double ** | Bx, | |||
double ** | By, | |||
double ** | Bz, | |||
double ** | Tx, | |||
double ** | Ty, | |||
double ** | Tz, | |||
double ** | Cval, | |||
double ** | Lval, | |||
double ** | Qx, | |||
double ** | Qy, | |||
double ** | Qz, | |||
double ** | Sval | |||
) |
!!!!! OLD OBSOLETE VERSION !!!!!!
int pdr_heat_select_el_coeff_vect | ( | int | Problem_id, | |
int * | Coeff_vect_ind | |||
) |
----------------------------------------------------------- pdr_heat_select_el_coeff_vect - to select coefficients returned to approximation routines for element integrals ------------------------------------------------------------