HemePure/cases/elasticWallsExample/inputAdvanced.xml at master · hemelb-codes/HemePure · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
<?xml version="1.0"?>
<hemelbsettings version="3">
<!-- Mesh Number = 0 -->
  <simulation>
    <step_length units="s" value="0.00033333333333333333333333333333333"/>
    <steps units="lattice" value="2400"/>
    <stresstype value="1"/>
    <voxel_size units="m" value="0.0002"/>
    <origin units="m" value="(0.0,0.0,0.0)"/>
    <elastic_wall_stiffness units="lattice" value="0.0444444444444444"/>
  </simulation>
 <geometry>
    <datafile path="pipe.gmy"/>
  </geometry>
  <initialconditions>
    <pressure>
      <uniform units="mmHg" value="0.0"/>
    </pressure>
  </initialconditions>
  <monitoring>
    <incompressibility/>
  </monitoring>

  <inlets>
    <inlet>
      <!-- index value=0-->
      <condition type="velocity" subtype="womersleyElastic">
        <pressure_gradient_amplitude value="-1e-1" units="mmHg/m"/>
        <period value="2.0" units="s"/>
        <womersley_number value="3.544907701811032" units="dimensionless"/>
        <radius value="0.004" units="m"/>
        <axial_position value="0.0" units="m"/>
        <poisson_ratio value="0.5" units="dimensionless"/>
        <youngs_modulus value="2400.0" units="Pa"/>
      </condition>
      <normal units="dimensionless" value="(0.0,0.0,1)"/>
      <position units="lattice" value="(23,23,3)"/>
    </inlet>
  </inlets>
  <outlets>
    <outlet>
      <!-- index value=0 -->
      <condition subtype="cosine" type="pressure">
        <amplitude units="mmHg" value="0.009930298683079167"/>
        <mean units="mmHg" value="0.0"/>
        <phase units="rad" value="-1.783329599538158"/>
	<period units="s" value="2.0"/>-->
      </condition>
      <normal units="dimensionless" value="(0.0,0.0,-1)"/>
      <position units="lattice" value="(23,23,203)"/>
    </outlet>
  </outlets>

  <properties>
   <propertyoutput file="inlet.dat" period="100000">
     <geometry type="inlet" />
     <field type="velocity" />
     <field type="pressure" />
   </propertyoutput>
   <propertyoutput file="outlet.dat" period="100000">
     <geometry type="outlet" />
     <field type="velocity" />
     <field type="pressure" />
   </propertyoutput>
   <propertyoutput file="planeX.dat" period="600">
     <geometry type="plane" >
	<normal units="dimensionless" value="(1.0,0.0,0.0)" />
       	<point units="m" value="(0.0046,0.0046,0.0206)" />
     </geometry>
     <field type="velocity" />
     <field type="pressure" />
   </propertyoutput>
    <propertyoutput file="surfaceshearstress.dat" period="600">
      <geometry type="surface" />
      <field type="shearstress"/>
    </propertyoutput>

   <propertyoutput file="whole.dat" period="100000">
     <geometry type="whole" />
     <field type="velocity" />
     <field type="pressure" />
   </propertyoutput>
  </properties>
</hemelbsettings>