|
#******************************************************************* #** #** v e m b l d e x m 0 8 #** #** Bingham fluid in 2-dimensional channel. 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 'vembldexm08.equation') and the second part #** defines the control parameters (please copy it to #** 'vembldexm08.resource'). The FORTRAN code for the solution #** of the problem is generated by entering #** 'vembuild vembldexm08' into your shell. #** #******************************************************************* #>>>>>>> cut here to get vembldexm08.equation <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The searched velocity (u1,u2) and pressure u2 in a fluid #** are the solution of the incompressible Navier-Stokes equation #** with stress dependent viscosity. It is a system of three #** nonlinear partial differential equations. The velocity is #** prescribed on all boundary portions (no slip) and the pressure #** is prescribed at a single point. #** #******************************************************************* #** #** material parameters: #** #** m - parameter #** etab - Bingham viscosity #** tt0 - yield point #** rho - density #** m=1 etab=0.1 tt0=1 rho=1 #** #** the boundary velocities are defined in the mesh data set : #** u1=prevalue u2=prevalue u3=1. #** #** viscosity : #** Ds11=(u1x1+u1x1)/2 Ds21=(u2x1+u1x2)/2 Ds12=(u1x2+u2x1)/2 Ds22=(u2x2+u2x2)/2 D=Ds11^2+Ds12^2+Ds21^2+Ds22^2 eta=2*(etab+tt0/sqrt(2*D)*(1-exp(-m*sqrt(2*D)))) #** #** stress tensor : #** Ts11=eta*Ds11 Ts21=eta*Ds21 Ts12=eta*Ds12 Ts22=eta*Ds22 #** #** the initial solution avoids that eta is defined in the first #** iteration step: #** u01=x1 u02=x2 u03=1. #** #** momentum equations : #** area { v1x1*Ts11+ v1x2*Ts12 + rho*v1*(u1*u1x1+u2*u1x2)+v1x1*u3 + v2x1*Ts21+ v2x2*Ts22 + rho*v2*(u1*u2x1+u2*u2x2)+v2x2*u3 #** #** equation of continuity #** + v3*(u1x1+u2x2)}=0 #** #******************************************************************* >>>>>>>> cut here to get vembldexm08.resource <<<<<<<<<<<<<<<<<<<<<<<<< #******************************************************************* #** #** The problem has a two dimensional domain and three solution #** component: #** DIM=2 NK=3 #** #******************************************************************* #** #** 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 pre- and the postprocessor is I-DEAS: #** MESH_PREP=i-deas MESH_POSTP=i-deas #** #******************************************************************* #** #** the is read from the I-DEAS universal file Lshape.unv. the #** domain is a L-shaped domain with prescribed input and output #** profile. at the walls the velocity is set to zero. the #** pressure is prescribed at a single node. #** MESH_FILEIN=Lshape.unv #** #** approximation: #** MESH_REDUCE=001 #** #******************************************************************* #** #** The solution component 1,2 are written to file velo.unv #** with the title 'velocity' and the third solution component #** is written to file pres.unv with the title 'pressure'. The #** considers only the error in the velocities and is written #** to file error.unv : #** OUTPUT_INDEX=110 001 OUTPUT_FILE=velo.unv, pres.unv OUTPUT_TITLE=velocity,pressure OUTPUT_ERRINDEX=110 OUTPUT_ERRFILE=error.unv OUTPUT_ERRSCAL=100. #** #******************************************************************* |