← Back to branch summary

~reducedmodelling/fluidity/ROM_Non-intrusive-ann

« back to all changes in this revision

Viewing changes to tests/darcy_p1dgp2_velBCinlet_1d/darcy_p1dgp2_velBCinlet_1d.xml

  • Committer: Brendan Tollit
  • Date: 2011-11-20 13:22:14 UTC
  • mto: This revision was merged to the branch mainline in revision 3860.
  • Revision ID: brendan.tollit05@imperial.ac.uk-20111120132214-no8ljc9argzddwon
Remove the two darcy test cases that were there in 1d and 
replace them with versions that do 1,2 and 3d.

Change the sign of the ct_rhs add to in Divergence_CV.F90 
for incompressible flow.

Adjust the darcy test that was there to test this so that 
it runs and the 1,2 and 3d version use weak velocity BC.

Add a similar test that uses two phase but one with a large 
absorption such as to be the solid porous media. This tests
the multiphase ct_rhs add to sign change.

The mms_ns case still has issues.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/darcy_p1dgp2_velBCinlet_1d/darcy_p1dgp2_velBCinlet_1d.xml
1
 
<?xml version="1.0" encoding="UTF-8" ?>
2
 
 
3
 
<testproblem>
4
 
  <name>darcy_p1dgp2_velBCinlet</name>
5
 
  <owner userid="btollit"/>
6
 
  <tags>flml</tags>
7
 
  <problem_definition length="short" nprocs="1">
8
 
    <command_line>
9
 
       fluidity darcy_p1dgp2_velBCinlet_1d.flml 
10
 
    </command_line>
11
 
    <!-- One dimensional problem for darcy flow with one region and one material using p1dgp2 element type testing the pressure gradient against analytic and the interstitial velocity. This tests using a diagnostic velocity absorption field associated with a mesh and momentum dg for darcy flow-->
12
 
  </problem_definition>
13
 
  <variables>
14
 
    <variable name="pressure" language="python">
15
 
import vtktools
16
 
v = vtktools.vtu("darcy_p1dgp2_velBCinlet_1d_1.vtu")
17
 
pressure = v.GetScalarRange("Pressure")
18
 
    </variable>
19
 
    <variable name="inter_vel" language="python">
20
 
import vtktools
21
 
v = vtktools.vtu("darcy_p1dgp2_velBCinlet_1d_1.vtu")
22
 
inter_vel = max(v.GetScalarRange("interstitial_velocity"))
23
 
    </variable>
24
 
  </variables>
25
 
  <pass_tests>
26
 
    <test name="change_P should equal domain_length*visc*darcy_vel_BC/perm, check relative difference to be under tolerance 1.0e-05" language="python">
27
 
change_P = abs(max(pressure) - min(pressure))
28
 
visc = 1.0e-04
29
 
darcy_vel_BC = 5.0
30
 
perm = 1.0e-10
31
 
domain_length = 300.0
32
 
assert abs((change_P - domain_length*visc*darcy_vel_BC/perm)/(domain_length*visc*darcy_vel_BC/perm)) &lt; 1.0e-05
33
 
    </test>
34
 
    <test name="interstitial velocity should equal darcy_vel/porosity), check relative difference to be under tolerance 1.0e-05" language="python">
35
 
analytic_vel = 10.0
36
 
assert abs((inter_vel - analytic_vel)/analytic_vel) &lt; 1.0e-05
37
 
    </test>
38
 
  </pass_tests>
39
 
  <warn_tests>
40
 
  </warn_tests>
41
 
</testproblem>