|
#******************************************************************* #** #** v e m b l d e x m 0 1 #** #** two coupled Poisson equation on the 2-dimensional unit cube. #** The mesh is read from a vecfem input file (e.g. generated by #** vemgen2d). #** #** 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 'vembldexm01.equation') and the second part #** defines the control parameters (please copy it to #** 'vembldexm01.resource'). The FORTRAN code for the solution #** of the problem is generated by entering #** 'vembuild vembldexm01' into your shell. #** #******************************************************************* ####### cut here to get vembldexm01.equation <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The problem is a system of two coupled 2-dimensional Poisson #** equations with Dirichlet and Neuman boundary conditions. #** #** The domain is [0,1] x [0,1], where the mesh is givem by #** rec.vem or is generated by vemgen2d. The mesh uses #** triangular elements of order 2 and line elements of order 2. #** Dirichlet conditions are set on all boundary points for the #** first component, but there is no Dirichlet condition for #** the second component. #** #******************************************************************* #** #** The functions b1, b2, f1, f2, g1 and g2 are selected, so that #** #** u1 = 1 #** u2 = 2 #** #** gets the exact solution of the problem. #** b1=1 b2=2 f1=10*b1+b2 f2=b1+10*b2 g1=0 g2=0 #** #******************************************************************* #** #** The Dirichlet conditions: #** u1=b1 u2=b2 #** #** 'u2=b2' is not necessary, since no Dirichlet conditions #** for the second component are specified in the input mesh. #** #******************************************************************* #** #** the functional equation : #** area{v1x1*u1x1+v1x2*u1x2+(10.*u1+ u2)*v1 - f1*v1 + v2x1*u2x1+v2x2*u2x2+( u1+10.*u2)*v2 - f2*v2} + line{v1*g1+v2*g2} =0 #** #******************************************************************* >>>>>>> cut here to get vembldexm01.resource <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The problem has a two dimensional domain and two solution #** component: #** DIM=2 NK=2 #** #******************************************************************* #** #** One processor with maximal 20 Mbytes are used. Maximal 1000 #** nodes and 1000 elements are allowed: #** PROCESS_STORAGE=20 PROCESS_MAXNN=1000 PROCESS_MAXNE=1000 #** #******************************************************************* #** #** the is read from the file rec.vem. other meshes can be #** created vemgen2d. #** MESH_PREP=print MESH_FILEIN=rec.vem #** #******************************************************************* #** #** The solution components are printed into the file solution.out #** and the indicator, which considers the error of all components #** is printed into file error.out. #** MESH_POSTP=print OUTPUT_INDEX=11 OUTPUT_FILE=solution.out OUTPUT_ERRINDEX=11 OUTPUT_ERRFILE=error.out #** #******************************************************************* |