VECFEM3 Reference Manual: vembldexm10.html

    #*******************************************************************
    #**
    #**    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
    #**
    #*******************************************************************