Home LibAPMath Animath PWDB
Cavity Flow Poisson Linear Elasticity

animath logo


Version 1.6: Example Applications

Lid-Driven Cavity

This package provides a debugged flowsolver module for Animath along with the classic lid-driven cavity test case.
The SURC solvers are global-timestep, consistent mass-matrix CFD solvers whose formulation is based on the Charactersitic-Based split method, which is mostly the work of O.C. Zienkievicz and P. Nithiarasu.
The SURC solvers are not yet adapted to problems with inhomogenous pressure boundary conditions. The work is still heavily in progress, so please do not consider this production code. The findings regarding upwind stabilization that were made with the SURC4 solver, are sure to manifest in changes in the cfd_c class of Animath 1.7.

Some Images that were generated with the SURC2 solver
Driving velocity (0, 400, 0) on plane (0, y, z)
centerline x 400 centerline y 400
v-value on (x, 0.5, 0.5) cavity center line u-value on (0.5, y, 0.5) cavity center line
Flow-direction-indicating particles released from the driving wall

Poisson on a Cylinder

This example program is designed for verification of the quadratic interpolation. To this end the domain on which the Poisson problem is solved, is chosen to be circular, i.e. a cylinder. The relative error of the numerical solution is written to a file together the mesh resolution.
Computation with curved elements Computation with straight elements

Linear Elasticity

The static elasticity test problem in this program stems from the set of Dr. Fukunari. The relative error of the numerical solution is written to a file together the mesh resolution. Please note that the analytic solution is immediately obtained, there is hence no conventional convergence to the analytic solution.
The problem consists of a linear-elastic unit cube with a Young's modulus of 1MPa and compressibility of nu = 0.45, that is exposed to a homogeneous stress-field of (500, 0.0, 0.0, 0.0, 0.0, 0.0) Pa.
The deviatroic setup differs in that the RHS has no divergence, i.e the pressure of 500/3 is subtracted from the 3 first components of the stress.