#*******************************************************************
#**
#**    v e m b l d e x m 0 4
#**
#**  time-dependent velocity driven diffusion on the 3-dimensional
#**  unit cube. 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 'vembldexm04.equation') and the second part
#**  defines the control parameters  (please copy it to
#**  'vembldexm04.resource'). The FORTRAN code for the solution
#**  of the problem is generated by entering
#**  'vembuild vembldexm04' into your shell.
#**
#*******************************************************************
#>>>>>>> cut here to get vembldexm04.equation <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem is the velocity driven, 3-D diffusion problem
#**  with Dirichlet and Neuman boundary conditions.
#**
#**  The domain is [-.5,5]^3 unit cube, where the mesh is
#**  generated by I-DEAS. The mesh uses tetrahedron elements of
#**  order 2 and triangle elements of order 2 for the boundaries.
#**  One Dirichlet condition is set. cube.unv is the I-DEAS
#**  universal file.
#**
#*******************************************************************
#**
#**  (w1,w2,w3) specifies the driving velocity profile:
#**
     w1=0
     w2=0
     w3=(x1-.5)*(x1+.5)*(x2-.5)*(x2+.5)*16.
#**
#*******************************************************************
#**
#**  The functions u01, b1, r1, g1, g2 and g2 are selected, so that
#**
#**     u1 = x3 * exp(sin(t))
#**
#**  gets the exact solution of the problem.
#**
     fac= exp(sin(t))
     u01=x3 * fac
     b1=u01
     r1=w3 * fac + x3 * fac * cos(t)
     g1=0
     g2=0
     g3=fac
#**
#*******************************************************************
#**
#** this is the outer normal direction:
#**
     n1 = tau21*tau32-tau31*tau22
     n2 = tau31*tau12-tau11*tau32
     n3 = tau11*tau22-tau21*tau12
     nn = - sqrt( n1^2 + n2^2 + n3^2 )
#**         |
#**         |- because of the orientation of the area elements in
#**            I-DEAS
#**
#*******************************************************************
#**
#**    The Dirichlet conditions:
#**
     u1=b1
#**
#*******************************************************************
#**
#**  the functional equation :
#**
    volume{v1x1 * u1x1 + v1x2 * u1x2 + v1x3 * u1x3 +
           v1*( w1 * u1x1 + w2 * u1x2 + w3 * u1x3 + ut1 - r1)}
  + area{-v1*(n1*g1+n2*g2+n3*g3)/nn} =0
#**
#*******************************************************************
>>>>>>>> cut here to get vembldexm04.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 1000
#**  nodes and 1000 elements are allowed:
#**
    PROCESS_STORAGE=50
    PROCESS_MAXNN=8000
    PROCESS_MAXNE=2000
#**
#*******************************************************************
#**
#**  the is read from the I-DEAS universal file cube.unv.
#**
    MESH_PREP=i-deas
    MESH_POSTP=i-deas
    MESH_FILEIN=cube.unv
#**
    SOLVER_TOL=1.E-2
    SOLVER_ERRSTP=0
#**
#*******************************************************************
#**
#**  activate the nonsteady solver :
#**
    SOLVER_STEADY=0
    SOLVER_T0=0
    SOLVER_TEND=10
    SOLVER_DT=1
    SOLVER_INTERP=0
#**
#*******************************************************************
#**
#**  The solution components are written into the file solution.unv
#**  and the indicator is written into file error.unv.
#**
    OUTPUT_FILE=solution.unv
    OUTPUT_ERRFILE=error.unv
#**
#*******************************************************************