Example 325: Flow from a cavity into an unsaturated porous media

This example calculates the flow from a cavity into an unsaturated porous media using the Yukawa (modified Helmholtz) equation.


The comfile run by this example is as follows:

#Example_325 Flow from a cavity into an unsaturated porous media

#           The Yukawa equation (or modified Helmholtz equation) arises
#           in the study of infiltration into unsaturated soils.  The 
#           problem considered below is that of steady infiltration 
#           from a circular cylindrical cavity buried at a large depth 
#           in a porous medium.

fem
fem define node;r;cavityflow;example #the nodal coordinates for the cavity
#                                !  are stored in CAVITYFLOW.ipnode in 
#                                !  the example directory.
fem define window;c  
fem define base;r;;example           #defines a BEM Lagrange/Hermite basis
#                                !  function with quadratic
#                                !  interpolation in 1 Xi coordinate.
fem define element;r;;example        #defines 10 BEM elements.
fem draw node
fem draw element
fem draw line
fem define equation;r;;example       #defines a static analysis of the
#                                !  modified Helmholtz equation by 
#                                !  direct BEM.
fem list equation                #List the equation
fem define material;r;;example       #defines the parameter s
#                                !  (dimensionless size of the cavity) 
#                                !  as constant spatially with value 1.
fem define initial;r;;example        #defines the initial conditions around
#                                !  the cavity as H=e^-s*y.
fem list initial                 #List the initial conditions
fem define solve;r;;example          #defines the solution procedure.
fem solve
fem draw field                   #Draw the solution field
# To obtain the flux from the cavity use the "fem list output" command.
# The cavity flux should be 1.30 (analytic is 1.2902).

Additional testing commands:

FEM list node solution
FEM list node flux

Files used by this example are:

Name               Modified     Size

example_325.com 10-Apr-2000 1.8k cavityflow.ipbase 10-Apr-2000 859 cavityflow.ipdata 10-Apr-2000 277 cavityflow.ipelem 10-Apr-2000 6.5k cavityflow.ipequa 26-May-2003 1.2k cavityflow.ipinit 10-Apr-2000 5.2k cavityflow.ipmate 10-Apr-2000 310 cavityflow.ipnode 10-Apr-2000 3.4k cavityflow.ipsolv 13-Apr-2007 1.1k test_output.com 10-Apr-2000 42

Download the entire example:

Name                      Modified     Size

examples_3_32_325.tar.gz 14-Apr-2007 5.7k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:47 20160
cm-debugSuccessSat Mar 5 00:02:03 20160
mips-irix
cmSuccessSun Aug 19 01:24:40 20073
cm-debugSuccessWed Aug 15 01:20:28 20075
cm-debug-clear-mallocSuccessSat Aug 18 01:24:44 20079
cm-debug-clear-malloc7SuccessMon Aug 20 01:22:51 20079
cm64SuccessSun Aug 19 01:24:43 20073
cm64-debugSuccessTue Aug 21 01:19:34 20075
cm64-debug-clear-mallocSuccessThu Apr 1 10:13:04 20045
rs6000-aix
cmSuccessWed Mar 4 01:07:04 20091
cm-debugSuccessMon Mar 2 01:06:17 20091
cm64SuccessWed Mar 4 01:07:04 20091
cm64-debugSuccessTue Mar 3 01:11:39 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20160
cm-debugSuccessSat Mar 5 00:01:12 20160

Testing status by file:


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

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


CMISS Help / Examples / 3 / 32 / 325