Page 209 - Process Modelling and Simulation With Finite Element Methods

196 Process Modelling and Simulation with Finite Element Methods

function a=pinit (x)
presgrad= [
183.59
183.471
...
4.01706
2.00851
0.000001 ;
xlist=[O:O.l:lOl ;
a=interpl (xlist, presgrad, x, I spline I ) ;
pinit( x - x) = jexp ( R (1 - ~gc(5 )))it
0
where X is the domain length, taken here as X=10.

Figure 5.8 Initial pressure profile for concentration given by I-erfc(x)
Many intermediate points are left out.

Now pull down the Mesh menu and select the Initialize Mesh option.
Select the time dependent solver and set output times to 0:5:200. Solve. The
final concentration and pressure profile should look as Figure 5.9.

Max 1 Max 1
Time=200 Surface concentratlon of c (c) Contour veloclty field (U-dl) 1
11
1 09
09 08
08 07
07
06
06
05 05
04 04
03 03
02 02
01
0 01
-0 1
0 05 1 1 5 2 2 5 3 3 5 4 4 5 5 5 5 6 65 7 75 8 85 9 9510 Mln 915Y0O96e-006
Figure 5.9 Concentration and pressure profiles for Darcy's Law and ConvectionDiffusion model.