#*******************************************************************
#**
#**    v e m b l d e x m 0 5
#**
#**  time-dependent thermal diffusion with temperature-dependent
#**  material coefficients on a 3-dimensional body. The mesh is
#**  read from an I-DEAS universal file.
#**
#**   by L. Grosz                           Karlsruhe, Jan. 1995
#**
#*******************************************************************
#**
#**  The data set of this examples has two parts (search for
#**  'cut here'). The first part specifies the problem
#**  (please copy it to 'vembldexm05.equation') and the second part
#**  defines the control parameters  (please copy it to
#**  'vembldexm05.resource'). The FORTRAN code for the solution
#**  of the problem is generated by entering
#**  'vembuild vembldexm05' into your shell.
#**
#*******************************************************************
#**
#**  The searched temperature in a thermal diffusion problem is
#**  given by partial differential equation of the Poisson type.
#**  Here we assumes a 3-dimensional body. On special portions of
#**  the boundary of the body the temperature is prescribed
#**  (=> Dirichlet conditions) and on the remainder portions
#**  convection boundary conditions are assumed. So you have a
#**  configuration like this:
#**
#**      /---------------------------------------\
#**      |   /--------\              / ------/   |   environment
#**      |   |  hole  |             /  hole /    |        u1=20
#**      |   | u1=800 |   body     /  u1=  /     |
#**      |   \ ------ /           /  1000 /      |
#**      |                      / ------/        |
#**      \---------------------------------------/
#**
#*******************************************************************
#>>>>>> cut here for vembldexm05.equation <<<<<<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  u1 is the searched temperature distribution.
#**
#**
#**  temperature of the environment :
#**
     tenv=20
#**
#**  thermal conductivity and capacity:
#**
      k=0.034*(1.+u1/100.)
      c=0.0045
#**
#**  no heat generation :
#**
      qb=0
#**
#**  convection boundary condition :
#**
      qs=0.015 * (u1 - tenv)
#**
#**  at the begin the body has the environ temperature:
#**
      u01=tenv
#**
#**  The temperature on the hole surfaces is set by the
#**  preprocessor. The value is increased from the initial
#**  temperature tenv to the actual value prevalue (this
#**  avoids an oscillation in space direction):
#**
      u1=(prevalue-tenv) * (1-exp(-100*t)) + tenv
#**
#**  The actual temperature distribution u1 is given by the
#**  minimal energy:
#**
     volume { k * ( u1x1 * v1x1 + u1x2 * v1x2 + u1x3 * v1x3) +
             (qb + c*ut1 ) * v1 } +
     area { qs * v1 } = 0
#**
#*******************************************************************
>>>>>>>> cut here to vembldexm05.resource <<<<<<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem has a three dimensional domain and one solution
#**  component:
#**
     DIM=3
     NK=1
#**
#*******************************************************************
#**
#**  One processor with maximal 20 Mbytes are used. Maximal 5500
#**  nodes and 1500 elements are allowed:
#**
    PROCESS_STORAGE=20
    PROCESS_MAXNN=5500
    PROCESS_MAXNE=1500
#**
#*******************************************************************
#**
#**  The pre- and the postprocessor is I-DEAS:
#**
    MESH_PREP=i-deas
    MESH_POSTP=i-deas
#**
#**  The mesh data are read from the I-DEAS universal file
#**  cooler.unv.
#**
    MESH_FILEIN= cooler.unv
#**
#*******************************************************************
#**
#** these are parameters to control the solver:
#**
    SOLVER_STEADY=0
    SOLVER_TOL=1.E-2
    SOLVER_T0=0.
    SOLVER_H=.05
    SOLVER_TEND=10.
    SOLVER_DT=.5
    SOLVER_INTERP=1
#**
#*******************************************************************
#**
#**  The first solution component is written to file temp.unv
#**  with the title 'temperature'. The error indicator is written
#**  to file error.unv :
#**
    OUTPUT_INDEX=1
    OUTPUT_FILE=temp.unv
    OUTPUT_TITLE=temperature
    OUTPUT_ERRFILE=error.unv
#**
#*******************************************************************