root/trunk/rheoplast/heatcond.h

/***************************************
  $Header$

  The typedefs and prototypes for heat conduction.  This is not meant to be
  included on its own, only when someone #includes "rheoplast.h".
  ***************************************/


#ifndef HEATCOND_H
#define HEATCOND_H    /*+ To stop multiple inclusions. +*/

/* These are to keep cxref happy */
#include <petsc.h>
#ifndef RHEOPLAST_H
typedef void AppCtx;
#endif

/*+ Structure typedef for heat conduction parameters. +*/
typedef struct {
  PetscScalar conductivity, thermdensity, Cp, /* k and rho*cp */
    latentheat,                         /* Latent heat of solidification */
    Tinit, Tcool, C;                  /* Initial temp and forced cooling parameters */
  int Tvar;                               /* Temperature field variable # */
} heatparm;

#define HEATCOND_HELP \
"Heat conduction is very simple.  To use it, add option:\n\
  -with_heat\n\
and control it with properties:\n\
  -conductivity <k>   thermal conductivity     [1.0]\n\
  -density      <rho> density                  [1.0]\n\
  -heatcap      <cp>  heat capacity            [1.0]\n\
  -latentheat   <L>   for vectorphase coupling [1.0]\n\
  -cooling   <C>   for vectorphase cooling     [0.0]\n\
  -T_cool <T_cool> for vectorphase cooling     [0.7]\n\
\n"
#endif /* HEATCOND_H */

#ifdef APPCTX_DEFINED
/* Function prototypes to make sure we're together */
int heatcond_first_setup
(PetscTruth threedee, int *vars, int *tempvars, int *stencilwid, AppCtx *data);
int heatcond_labels_initcond
(PetscScalar *globalarray, int nx,int ny,int nz, int xm,int ym,int zm, int xs,
 int ys,int zs, int vars, AppCtx *data, PetscScalar *max_explicit_deltat);
int heatcond_temp_parameters_line
(PetscScalar *x, PetscScalar *temp, int points, int gxm,int gym,
 PetscScalar xmin,PetscScalar xmax, PetscScalar ycoord,PetscScalar zcoord,
 PetscScalar time, AppCtx *data);
int heatcond_interior_line_function
(PetscScalar *x, PetscScalar *func, PetscScalar *temp,
 PetscTruth **mixed_constraints, int points, int gxm,int gym, PetscScalar xmin,
 PetscScalar xmax, PetscScalar ycoord, PetscScalar zcoord, PetscScalar time,
 AppCtx *data);
#endif /* APPCTX_DEFINED */