Example 339: Solving 3d Navier-Stokes equations

Solving the 3d non-dimensional Navier-Stokes pipe flow.

The comfile run by this example is as follows:

# Example 339 - Solving 3d non-dimensional Navier Stokes flow in a pipe

fem def para;r;pipe3d;example # Define the array parameters to allow Voronoi mesh use

fem def coord 3,1             # Define 3d rectangular Cartesian coordinates
fem def node;r;pipe3d;example # Define the Delaunay/Voronoi nodes
fem def base;r;pipe3d;example # Define 3d simplex basis functions 
fem def elem;r;pipe3d;example # Define 3d Delaunay elements
fem up dela                   # Update the Delaunay triangulation to ensure Delaunayism
fem def mesh;r;pipe3d;example # Define a Voronoi mesh to be calculated
fem def equa;r;pipe3d;example # Define equation type. This will be a time integration,
                              # Navier Stokes equation, of a general type, solved 
                              # using finite volumes.       
fem def mate;r;pipe3d;example #  Define material. This is unimportant as the equations
                              # solved are non-dimensional, ie only depend on the 
                              # Reynolds number. These parameters are used only for
                              # conversion back to dimensional form. This has not
                              # been implemented yet.
fem def init;r;pipe3d;example # Define initial conditions. At one end of the pipe will
                              # be a specified inlet velocity (Type inlet), while at the
                              # other end will be an outlet boundary condition. All other
                              # Bc's will be a Wall with zero velocity. Be sure to use
                              # duct flow outlet boundary conditions as they are superior
fem def solv;r;pipe3d;example # Define the solution parameters. 
                              # Always use automatic stepping unless stupid.
                              # As the mesh is fixed, you want to use a sparse LU 
                              # solver (only needs to be factorised once). 
                              # Use QUICK advection (best for steady state cases).
                              # The Courant number scales the maximum allowable
                              # time step by a factor of safety, should be between
                              # 0.3 and 0.7.
fem solv                      # Solve to steady state.

Additional testing commands:

fem list node solution

Files used by this example are:

Name             Modified     Size

example_339.com  20-Nov-2001  2.2k
pipe3d.ipbase    10-Apr-2000  2.1k
pipe3d.ipelem    10-Apr-2000  265k
pipe3d.ipequa    26-May-2003  1.4k
pipe3d.ipinit    13-Mar-2001   31k
pipe3d.ipmate    10-Apr-2000  684 
pipe3d.ipmesh    06-Nov-2001  978 
pipe3d.ipnode    10-Apr-2000   33k
pipe3d.ippara    12-Nov-2002  5.6k
pipe3d.ipsolv    13-Apr-2007  2.7k
test_output.com  10-Apr-2000   23 

Download the entire example:

Name                      Modified     Size

examples_3_33_339.tar.gz  14-Apr-2007   21k

Testing status by version:

Testing status by file:

• cmiss_test.log
  

  Success	cm_aix_mp_1:	cmiss_test.log.retain.
  Success	cm_aix_mp_2:	cmiss_test.log.retain.
  Success	cm_aix_mp_4:	cmiss_test.log.retain.
  Success	cm_aix_mp_8:	cmiss_test.log.retain.
  Success	cm_n32_mp_1:	cmiss_test.log.retain.
  Success	cm_n32_mp_2:	cmiss_test.log.retain.
  Success	cm_n32_mp_4:	cmiss_test.log.retain.
  Success	cm_n32_mp_8:	cmiss_test.log.retain.
  i686-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  mips-irix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
  rs6000-aix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  x86_64-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.

  
  
• cmiss_test.out
  

  Success	cm_aix_mp_1:	ndiff test: no significant differences with generic answer.
  Success	cm_aix_mp_2:	ndiff test: no significant differences with generic answer.
  Success	cm_aix_mp_4:	ndiff test: no significant differences with generic answer.
  Success	cm_aix_mp_8:	ndiff test: no significant differences with generic answer.
  Success	cm_n32_mp_1:	ndiff test: no significant differences with generic answer.
  Success	cm_n32_mp_2:	ndiff test: no significant differences with generic answer.
  Success	cm_n32_mp_4:	ndiff test: no significant differences with generic answer.
  Success	cm_n32_mp_8:	ndiff test: no significant differences with generic answer.
  i686-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  mips-irix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  rs6000-aix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  x86_64-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.

  
  

Html last generated: Sun Mar 6 05:50:18 2016

Input last modified: Fri Apr 13 10:32:43 2007