#*******************************************************************
#**
#**    v e m b l d e x m 1 0
#**
#**  the diffusion driven by a velocity field w on the
#**  2-dimensional unit cube [0,1]^2. The mesh is read from the
#**  vecfem input file (e.g. generated by vemgen2dq).
#**
#**   by L. Grosz                           Karlsruhe, June 1995
#**
#*******************************************************************
#**
#**  The data set of this examples has two parts (search for
#**  'cut here'). The first part specifies the problem
#**  (please copy it to 'vembldexm10.equation') and the second part
#**  defines the control parameters  (please copy it to
#**  'vembldexm10.resource'). The FORTRAN code for the solution
#**  of the problem is generated by entering
#**  'vembuild maple vembldexm10' into your shell.
#**
#*******************************************************************
#>>>>>>> cut here to get vembldexm10.equation <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem is the velocity driven, 2-D diffusion problem
#**  on the unit cube [0,1]^2 for a searched temperature
#**  distribution u1. On the lower boundary X2=0 the temperature is
#**  prescribed (=> Dirichlet condition) and at the upper boundary
#**  X2=1 a radiation boundary condition is set. At the remaining
#**  boundaries X1=0 and X1=1 the domain is isolated.
#**
#**  The domain is subdivided into quadrilateral elements and the
#**  boundary portion, where the radiation boundary condition is
#**  set, is subdivided into line elements. The mesh is given in
#**  square.vem or can be generated by vemgen2dq.
#**
#*******************************************************************
#**
#** these are the material constants:
#**
k=0.02          # thermal conductivity
c=10.           # heat capacity
alpha=7.6       # radiation coefficient
q=0             # there is no thermal source
#**
#*******************************************************************
#**
#** the driving velocity field:
#**
w1=x1/sqrt(x1^2+x2^2)*10
w2=x2/sqrt(x1^2+x2^2)*10
#**
#*******************************************************************
#**
#** the initial temperature at time 0:
#**
u01=1
#**
#*******************************************************************
#**
#** the temperature at the Dirichlet nodes is :
#**
u1=10
#**
#*******************************************************************
#**
#** the diffusion equation in the weak formulation:
#**
area{ k*(v1x1*u1x1+v1x2*u1x2) + v1*(w1*u1x1+w2*u1x2+c*ut1-q) }+
     line{v1*alpha*(u1-1) }=0
#**   /|
#**    ---- this is the radiation boundary condition
#**
#*******************************************************************
>>>>>>>> cut here to get vembldexm10.resource <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem has a two dimensional domain and one solution
#**  component:
#**
NK=1
DIM=2
#**
#*******************************************************************
#**
#**  the is read from the vecfem input file square.vem. Other
#**  meshes can be generated by vemgen2dq. the postprocessor is
#**  I-DEAS.
#**
MESH_PREP=print
MESH_FILEIN=square.vem
MESH_POSTP=i-deas
MESH_FILEOUT=mesh.unv
#**
#*******************************************************************
#**
#**  One processor with maximal 5 Mbytes are used. Maximal 700
#**  nodes and 200 elements are allowed:
#**
PROCESS_STORAGE=5
PROCESS_MAXNN=1000
PROCESS_MAXNE=300
#**
#*******************************************************************
#**
#**  activate the nonsteady solver :
#**
SOLVER_STEADY=0
SOLVER_H=0.01
SOLVER_T0=0
SOLVER_TEND=1
SOLVER_DT=0.1
SOLVER_INTERP=1 # the solution is computed at equidistant time steps
#**
#*******************************************************************
#**
#**  The solution component is written to the file temp.unv
#**  and the indicator is written into file error.unv.
#**
OUTPUT_ERRFILE=error.unv
OUTPUT_ERRSCAL=1.
OUTPUT_INDEX=1
OUTPUT_FILE=temp.unv
OUTPUT_TITLE=temperature
#**
#*******************************************************************