/
CaseThwaites01_Def.xml
260 lines (260 loc) · 11.9 KB
/
CaseThwaites01_Def.xml
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
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
<?xml version="1.0" encoding="UTF-8" ?>
<!-- Case name: CaseThwaites01 -->
<case app="DesignSPHysics v0.5.1806-06" date="02-04-2019 11:56:59">
<casedef>
<constantsdef>
<lattice bound="1" fluid="1" />
<gravity x="0" y="0" z="-9.81" comment="Gravitational acceleration" units_comment="m/s^2" />
<rhop0 value="1000" comment="Reference density of the fluid" units_comment="kg/m^3" />
<hswl value="0" auto="true" comment="Maximum still water level to calculate speedofsound using coefsound" units_comment="metres (m)" />
<gamma value="7" comment="Polytropic constant for water used in the state equation" />
<speedsystem value="0" auto="true" comment="Maximum system speed (by default the dam-break propagation is used)" />
<coefsound value="20" comment="Coefficient to multiply speedsystem" />
<speedsound value="0" auto="true" comment="Speed of sound to use in the simulation (by default speedofsound=coefsound*speedsystem)" />
<coefh value="0.91924" comment="Coefficient to calculate the smoothing length (h=coefh*sqrt(3*dp^2) in 3D)" />
<cflnumber value="0.1" comment="Coefficient to multiply dt" />
</constantsdef>
<mkconfig boundcount="240" fluidcount="10">
</mkconfig>
<geometry>
<definition dp="0.007" comment="Initial inter-particle distance" units_comment="metres (m)">
<pointmin x="1e-04" y="0.0" z="-1.2" />
<pointmax x="10.7999" y="0.0" z="1.05" />
</definition>
<commands>
<mainlist>
<setshapemode>dp | bound | fluid | mk</setshapemode>
<matrixreset />
<setmkbound mk="0"/>
<setdrawmode mode="face"/>
<drawfilestl file="z_ice_bed.stl" objname="Bed">
<drawscale x="0.001" y="0.001" z="0.001" />
</drawfilestl>
<!-- Note: the following commented command was used to convert large STL files to small, quantized VTKs. -->
<!--<shapeout file="z_ice_bed" />-->
<matrixreset />
<setmkfluid mk="0"/>
<setdrawmode mode="face"/>
<move x="-0.01" y="0" z="-0.105" />
<drawfilestl file="z_mw.stl" objname="Meltwater">
<drawscale x="0.001" y="0.001" z="0.001" />
</drawfilestl>
<!-- Note: the following commented command was used to convert large STL files to small, quantized VTKs. -->
<!--<shapeout file="z_mw" />-->
<matrixreset />
<setmkbound mk="1"/>
<setdrawmode mode="face"/>
<drawfilestl file="z_ice_sfc.stl" objname="Ice surface">
<drawscale x="0.001" y="0.001" z="0.001" />
</drawfilestl>
<!-- Note: the following commented command was used to convert large STL files to small, quantized VTKs. -->
<!--<shapeout file="z_ice_sfc" />-->
<matrixreset />
<setmkbound mk="1"/>
<setdrawmode mode="face"/>
<drawfilestl file="z_ice_bott.stl" objname="Ice-bottom">
<drawscale x="0.001" y="0.001" z="0.001" />
</drawfilestl>
<!-- Note: the following commented command was used to convert large STL files to small, quantized VTKs. -->
<!--<shapeout file="z_ice_bott" />-->
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="0.0" y="-0.05" z="0.4" />
<drawbox objname="wind000">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="0.01" y="0.1" z="0.8" />
</drawbox>
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="2.15" y="-0.05" z="0.4" />
<drawbox objname="wind001">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="0.01" y="0.1" z="0.8" />
</drawbox>
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="4.3" y="-0.05" z="0.4" />
<drawbox objname="wind002">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="0.01" y="0.1" z="0.8" />
</drawbox>
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="6.45" y="-0.05" z="0.4" />
<drawbox objname="wind003">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="0.01" y="0.1" z="0.8" />
</drawbox>
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="8.6" y="-0.05" z="0.4" />
<drawbox objname="wind004">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="0.01" y="0.1" z="0.8" />
</drawbox>
<matrixreset />
<setmkbound mk="2"/>
<setdrawmode mode="full"/>
<move x="-0.05" y="-0.05" z="1.04" />
<drawbox objname="top">
<boxfill>solid</boxfill>
<point x="0" y="0" z="0" />
<size x="10.85" y="0.1" z="0.01" />
</drawbox>
<matrixreset />
<setmkbound mk="0"/>
<setdrawmode mode="solid"/>
<move x="0.0" y="-0.05" z="-1.3" />
<fillbox x="3.0" y="0.05" z="0.3" objname="FillBox-bed">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="10.801" y="0.1" z="1.0" />
<matrixreset />
</fillbox>
<matrixreset />
<setmkbound mk="1"/>
<setdrawmode mode="solid"/>
<move x="-5e-04" y="-0.05" z="-0.9" />
<fillbox x="3.0005" y="0.05" z="0.8" objname="FillBox-ice">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="4.5" y="0.1" z="1.23" />
<matrixreset />
</fillbox>
<matrixreset />
<setmkfluid mk="0"/>
<setdrawmode mode="solid"/>
<move x="0.9" y="-0.05" z="-0.9" />
<fillbox x="1.4" y="0.05" z="0.58" objname="FillBox-meltwater">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="4.15" y="0.1" z="0.9001" />
<matrixreset />
</fillbox>
<matrixreset />
<setmkfluid mk="1"/>
<setdrawmode mode="solid"/>
<move x="0.9" y="-0.05" z="-1.2" />
<fillbox x="6.1" y="0.05" z="0.88" objname="FillBox-cdw">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="9.901" y="0.1" z="1.2001" />
<matrixreset />
</fillbox>
<matrixreset />
<setmkfluid mk="2"/>
<setdrawmode mode="solid"/>
<move x="-5e-04" y="-0.05" z="-0.007" />
<fillbox x="5.0005" y="0.05" z="0.502" objname="FillBox-air">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="10.801" y="0.1" z="1.203" />
<matrixreset />
</fillbox>
<shapeout file="" />
</mainlist>
</commands>
</geometry>
<properties> %Material properties for liquid-gas flows are set here.
<links>
<link mkbound="0" property="basalt"/>
<link mkbound="1" property="ice"/>
<link mkbound="2" property="air"/>
<link mkfluid="0" property="meltwater" />
<link mkfluid="1" property="cdw"/>
<link mkfluid="2" property="air"/>
<link mkfluid="3" property="air"/>
</links>
<property name="basalt">
<rhop0 value="3300"/>
<gamma value="7"/>
<cs0 value="150"/>
<visco value="0.01"/>
<phasetype value="0"/>
</property>
<property name="ice">
<rhop0 value="900"/>
<gamma value="7"/>
<cs0 value="150"/>
<visco value="0.01"/>
<phasetype value="0"/>
</property>
<property name="air">
<rhop0 value="1.44"/>
<gamma value="1.4"/>
<cs0 value="150"/>
<visco value="0.01"/>
<phasetype value="1"/>
</property>
<property name="meltwater">
<rhop0 value="1020"/>
<gamma value="7"/>
<cs0 value="20"/>
<visco value="0.1"/>
<phasetype value="0"/>
</property>
<property name="cdw">
<!-- CDW T/sal: 0.1–1.4C (but warming overall and locally warmer), 34.62–34.73 per mil -->
<!-- Source: Cook et al 2016 - http://science.sciencemag.org.prxy4.ursus.maine.edu/content/353/6296/283 -->
<!-- Fig. 2 https://www.nature.com/scitable/knowledge/library/key-physical-variables-in-the-ocean-temperature-102805293 -->
<rhop0 value="1028"/>
<gamma value="7"/>
<cs0 value="20"/>
<visco value="0.1"/>
<phasetype value="0"/>
</property>
</properties>
<motion>
<objreal ref="2">
<!-- Movement Name: Catabatic wind -->
<begin mov="1" start="0"/>
<mvrect id="1" duration="1000.0">
<vel x="3.0" y="0.0" z="0.0"/>
</mvrect>
</objreal>
</motion>
</casedef>
<execution>
<special>
</special>
<parameters>
<parameter key="PosDouble" value="1" comment="Precision in particle interaction 0:Simple, 1:Double, 2:Uses and saves double (default=0)" />
<parameter key="FlowType" value="2" comment="Select Single-phase (1) or multi-phase flow(2)" />
<parameter key="StepAlgorithm" value="2" comment="Step Algorithm 1:Verlet, 2:Symplectic (default=1)" />
<parameter key="VerletSteps" value="40" comment="Verlet only: Number of steps to apply Euler timestepping (default=40)" />
<parameter key="Kernel" value="2" comment="Interaction Kernel 1:Cubic Spline, 2:Wendland (default=2)" />
<parameter key="ViscoTreatment" value="1" comment="Viscosity formulation 1:Artificial, 2:Laminar+SPS (default=1)" />
<parameter key="Visco" value="0.0" comment="Viscosity value" /> % Note it is not used for multi-phase flows. The "visco" variable in the properties section will be used instead
<parameter key="ViscoBoundFactor" value="0" comment="Multiply viscosity value with boundary (default=1)" />
<parameter key="DeltaSPH" value="0.1" comment="DeltaSPH value, 0.1 is the typical value, with 0 disabled (default=0)" />
<parameter key="Shifting" value="1" comment="Shifting mode 0:None, 1:Ignore bound, 2:Ignore fixed, 3:Full (default=0)" />
<parameter key="ShiftCoef" value="-2" comment="Coefficient for shifting computation (default=-2)" />
<parameter key="ShiftTFS" value="1.5" comment="Threshold to detect free surface. Typically 1.5 for 2D and 2.75 for 3D (default=0)" />
<parameter key="MPCoefficient" value="0.1" comment="Colagrossi and Landrini multiphase coefficient." />
<parameter key="RigidAlgorithm" value="1" comment="Rigid Algorithm 1:SPH, 2:DEM (default=1)" />
<parameter key="#FtPause" value="0.0" comment="Time to freeze the floatings at simulation start (warmup) (default=0)" units_comment="seconds" />
<parameter key="#CoefDtMin" value="0.05" comment="Coefficient to calculate minimum time step dtmin=coefdtmin*h/speedsound (default=0.05)" />
<parameter key="DtIni" value="0.00001" comment="Initial time step (default=h/speedsound)" units_comment="seconds" />
<parameter key="DtMin" value="0.0000001" comment="Minimum time step (default=coefdtmin*h/speedsound)" units_comment="seconds" />
<parameter key="#DtFixed" value="DtFixed.dat" comment="Dt values are loaded from file (default=disabled)" />
<parameter key="DtAllParticles" value="0" comment="Velocity of particles used to calculate DT. 1:All, 0:Only fluid/floating (default=0)" />
<parameter key="TimeMax" value="30" comment="Time of simulation" units_comment="seconds" />
<parameter key="TimeOut" value="0.01" comment="Time out data" units_comment="seconds" />
<parameter key="IncZ" value="0" comment="Increase of Z+" units_comment="decimal" />
<parameter key="PartsOutMax" value="1" comment="%/100 of fluid particles allowed to be excluded from domain (default=1)" units_comment="decimal" />
<parameter key="RhopOutMin" value="0" comment="Minimum rhop valid (default=700)" units_comment="kg/m^3" />
<parameter key="RhopOutMax" value="1300" comment="Maximum rhop valid (default=1300)" units_comment="kg/m^3" />
<parameter key="XPeriodicIncZ" value="0" comment="Increase of Z with periodic BC" />
</parameters>
</execution>
</case>