Example 338: Solving 2d Navier-Stokes equations

Solving the 2d non-dimensional Navier-Stokes equations in a rotating cylinder.

The comfile run by this example is as follows:

# example_338.com - Solving Navier-Stokes flow in a rotating cylinder

fem
fem def para;r;circle;example # Define parameters to allow Voronoi arrays to be used.

fem def node;r;circle;example # Define the nodes of the Voronoi/Delaunay mesh.
fem def base;r;circle;example # Define 2d simplex basis functions.
fem def elem;r;circle;example # Define triangular simplex elements.
fem update delaunay           # Update the Delaunay triangulation to ensure Delaunayism
fem def mesh;r;circle;example # Define a Voronoi mesh. 
fem def equa;r;circle;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;circle;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;circle;example # Define initial (boundary) conditions. This involves
                              # defining the type (and velocity if using inlet or
                              # wall boundary conditions) at each boundary (B) Voronoi
                              # node.
fem def solv;r;circle;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 it to steady state.
 

Additional testing commands:

fem list node solution

Files used by this example are:

Name                  Modified     Size

example_338.com       20-Nov-2001  2.0k
circle.ipbase         10-Apr-2000  1.0k
circle.ipelem         10-Apr-2000   38k
circle.ipequa         26-May-2003  1.3k
circle.ipinit         10-Apr-2000  7.9k
circle.ipinit.8jun99  10-Apr-2000  7.9k
circle.ipmate         10-Apr-2000  686 
circle.ipmesh         06-Nov-2001  947 
circle.ipnode         10-Apr-2000   13k
circle.ippara         12-Nov-2002  5.8k
circle.ipsolv         13-Apr-2007  2.7k
test_output.com       10-Apr-2000   23 

Download the entire example:

Name                      Modified     Size

examples_3_33_338.tar.gz  14-Apr-2007   15k

Testing status by version:

Testing status by file:

• cmiss_test.log
  

  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
  

  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:41 2007