← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to tests/sediment_critical_shear/sediment_critical_shear.flml

  • Committer: Simon Mouradian
  • Date: 2012-10-19 10:35:59 UTC
  • mfrom: (3520.32.371 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20121019103559-y36qa47phc69q8sc
merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/sediment_critical_shear/sediment_critical_shear.flml
1
 
<?xml version='1.0' encoding='utf-8'?>
2
 
<fluidity_options>
3
 
  <simulation_name>
4
 
    <string_value lines="1">Sediment_Reentrainment</string_value>
5
 
  </simulation_name>
6
 
  <problem_type>
7
 
    <string_value lines="1">fluids</string_value>
8
 
  </problem_type>
9
 
  <geometry>
10
 
    <dimension>
11
 
      <integer_value rank="0">2</integer_value>
12
 
    </dimension>
13
 
    <mesh name="CoordinateMesh">
14
 
      <from_file file_name="channel">
15
 
        <format name="gmsh"/>
16
 
        <stat>
17
 
          <include_in_stat/>
18
 
        </stat>
19
 
      </from_file>
20
 
    </mesh>
21
 
    <mesh name="VelocityMesh">
22
 
      <from_mesh>
23
 
        <mesh name="PeriodicMesh"/>
24
 
        <stat>
25
 
          <exclude_from_stat/>
26
 
        </stat>
27
 
      </from_mesh>
28
 
    </mesh>
29
 
    <mesh name="PressureMesh">
30
 
      <from_mesh>
31
 
        <mesh name="PeriodicMesh"/>
32
 
        <mesh_shape>
33
 
          <polynomial_degree>
34
 
            <integer_value rank="0">2</integer_value>
35
 
          </polynomial_degree>
36
 
        </mesh_shape>
37
 
        <stat>
38
 
          <exclude_from_stat/>
39
 
        </stat>
40
 
      </from_mesh>
41
 
    </mesh>
42
 
    <mesh name="PeriodicMesh">
43
 
      <from_mesh>
44
 
        <mesh name="CoordinateMesh"/>
45
 
        <periodic_boundary_conditions name="x">
46
 
          <physical_boundary_ids>
47
 
            <integer_value shape="1" rank="1">6</integer_value>
48
 
          </physical_boundary_ids>
49
 
          <aliased_boundary_ids>
50
 
            <integer_value shape="1" rank="1">7</integer_value>
51
 
          </aliased_boundary_ids>
52
 
          <coordinate_map>
53
 
            <string_value lines="20" type="code" language="python">def val(X,t):
54
 
   return(X[0]+200, X[1])</string_value>
55
 
          </coordinate_map>
56
 
        </periodic_boundary_conditions>
57
 
        <stat>
58
 
          <exclude_from_stat/>
59
 
        </stat>
60
 
      </from_mesh>
61
 
    </mesh>
62
 
    <quadrature>
63
 
      <degree>
64
 
        <integer_value rank="0">4</integer_value>
65
 
      </degree>
66
 
    </quadrature>
67
 
  </geometry>
68
 
  <io>
69
 
    <dump_format>
70
 
      <string_value>vtk</string_value>
71
 
    </dump_format>
72
 
    <dump_period>
73
 
      <constant>
74
 
        <real_value rank="0">10</real_value>
75
 
      </constant>
76
 
    </dump_period>
77
 
    <output_mesh name="CoordinateMesh"/>
78
 
    <stat/>
79
 
  </io>
80
 
  <timestepping>
81
 
    <current_time>
82
 
      <real_value rank="0">0</real_value>
83
 
    </current_time>
84
 
    <timestep>
85
 
      <real_value rank="0">1</real_value>
86
 
    </timestep>
87
 
    <finish_time>
88
 
      <real_value rank="0">400</real_value>
89
 
    </finish_time>
90
 
    <steady_state>
91
 
      <tolerance>
92
 
        <real_value rank="0">1e-10</real_value>
93
 
        <infinity_norm/>
94
 
      </tolerance>
95
 
    </steady_state>
96
 
  </timestepping>
97
 
  <physical_parameters>
98
 
    <gravity>
99
 
      <magnitude>
100
 
        <real_value rank="0">9.81</real_value>
101
 
      </magnitude>
102
 
      <vector_field name="GravityDirection" rank="1">
103
 
        <prescribed>
104
 
          <mesh name="CoordinateMesh"/>
105
 
          <value name="WholeMesh">
106
 
            <constant>
107
 
              <real_value shape="2" dim1="dim" rank="1">0 -1</real_value>
108
 
            </constant>
109
 
          </value>
110
 
          <output/>
111
 
          <stat>
112
 
            <include_in_stat/>
113
 
          </stat>
114
 
          <detectors>
115
 
            <exclude_from_detectors/>
116
 
          </detectors>
117
 
        </prescribed>
118
 
      </vector_field>
119
 
    </gravity>
120
 
    <coriolis>
121
 
      <f_plane>
122
 
        <f>
123
 
          <real_value rank="0">1</real_value>
124
 
        </f>
125
 
      </f_plane>
126
 
    </coriolis>
127
 
  </physical_parameters>
128
 
  <material_phase name="Fluid">
129
 
    <equation_of_state>
130
 
      <fluids>
131
 
        <linear>
132
 
          <reference_density>
133
 
            <real_value rank="0">1</real_value>
134
 
          </reference_density>
135
 
          <subtract_out_hydrostatic_level/>
136
 
        </linear>
137
 
      </fluids>
138
 
    </equation_of_state>
139
 
    <scalar_field name="Pressure" rank="0">
140
 
      <prognostic>
141
 
        <mesh name="PressureMesh"/>
142
 
        <spatial_discretisation>
143
 
          <continuous_galerkin>
144
 
            <remove_stabilisation_term/>
145
 
            <integrate_continuity_by_parts/>
146
 
          </continuous_galerkin>
147
 
        </spatial_discretisation>
148
 
        <reference_node>
149
 
          <integer_value rank="0">1</integer_value>
150
 
        </reference_node>
151
 
        <scheme>
152
 
          <poisson_pressure_solution>
153
 
            <string_value lines="1">never</string_value>
154
 
          </poisson_pressure_solution>
155
 
          <use_projection_method/>
156
 
        </scheme>
157
 
        <solver>
158
 
          <iterative_method name="cg"/>
159
 
          <preconditioner name="sor"/>
160
 
          <relative_error>
161
 
            <real_value rank="0">1e-10</real_value>
162
 
          </relative_error>
163
 
          <max_iterations>
164
 
            <integer_value rank="0">1000</integer_value>
165
 
          </max_iterations>
166
 
          <never_ignore_solver_failures/>
167
 
          <diagnostics>
168
 
            <monitors/>
169
 
          </diagnostics>
170
 
        </solver>
171
 
        <output>
172
 
          <include_previous_time_step/>
173
 
        </output>
174
 
        <stat/>
175
 
        <convergence>
176
 
          <include_in_convergence/>
177
 
        </convergence>
178
 
        <detectors>
179
 
          <exclude_from_detectors/>
180
 
        </detectors>
181
 
        <steady_state>
182
 
          <include_in_steady_state/>
183
 
        </steady_state>
184
 
        <consistent_interpolation/>
185
 
      </prognostic>
186
 
    </scalar_field>
187
 
    <scalar_field name="Density" rank="0">
188
 
      <diagnostic>
189
 
        <algorithm name="Internal" material_phase_support="multiple"/>
190
 
        <mesh name="VelocityMesh"/>
191
 
        <output/>
192
 
        <stat/>
193
 
        <convergence>
194
 
          <include_in_convergence/>
195
 
        </convergence>
196
 
        <detectors>
197
 
          <include_in_detectors/>
198
 
        </detectors>
199
 
        <steady_state>
200
 
          <include_in_steady_state/>
201
 
        </steady_state>
202
 
      </diagnostic>
203
 
    </scalar_field>
204
 
    <vector_field name="Velocity" rank="1">
205
 
      <prescribed>
206
 
        <mesh name="VelocityMesh"/>
207
 
        <value name="WholeMesh">
208
 
          <python>
209
 
            <string_value lines="20" type="code" language="python">def val(X,t):
210
 
        if (t &lt; 200):
211
 
                return [0,0]
212
 
        else:
213
 
                return [X[1]/50,0]</string_value>
214
 
          </python>
215
 
        </value>
216
 
        <output/>
217
 
        <stat>
218
 
          <include_in_stat/>
219
 
        </stat>
220
 
        <detectors>
221
 
          <exclude_from_detectors/>
222
 
        </detectors>
223
 
      </prescribed>
224
 
    </vector_field>
225
 
    <vector_field name="BedShearStress" rank="1">
226
 
      <diagnostic>
227
 
        <algorithm name="Internal" material_phase_support="multiple"/>
228
 
        <mesh name="VelocityMesh"/>
229
 
        <density>
230
 
          <real_value rank="0">1000</real_value>
231
 
        </density>
232
 
        <drag_coefficient>
233
 
          <real_value rank="0">0.001</real_value>
234
 
        </drag_coefficient>
235
 
        <output/>
236
 
        <stat>
237
 
          <include_in_stat/>
238
 
        </stat>
239
 
        <detectors>
240
 
          <include_in_detectors/>
241
 
        </detectors>
242
 
      </diagnostic>
243
 
    </vector_field>
244
 
    <sediment>
245
 
      <scalar_field name="SedimentTemplate" rank="0">
246
 
        <prognostic>
247
 
          <mesh name="VelocityMesh"/>
248
 
          <equation name="AdvectionDiffusion"/>
249
 
          <spatial_discretisation>
250
 
            <continuous_galerkin>
251
 
              <stabilisation>
252
 
                <streamline_upwind>
253
 
                  <nu_bar_optimal/>
254
 
                  <nu_scale name="0.5">
255
 
                    <real_value shape="1" rank="0">0.5</real_value>
256
 
                  </nu_scale>
257
 
                </streamline_upwind>
258
 
              </stabilisation>
259
 
              <advection_terms/>
260
 
              <mass_terms/>
261
 
            </continuous_galerkin>
262
 
            <conservative_advection>
263
 
              <real_value rank="0">1</real_value>
264
 
            </conservative_advection>
265
 
          </spatial_discretisation>
266
 
          <temporal_discretisation>
267
 
            <theta>
268
 
              <real_value rank="0">1</real_value>
269
 
            </theta>
270
 
          </temporal_discretisation>
271
 
          <solver>
272
 
            <iterative_method name="gmres">
273
 
              <restart>
274
 
                <integer_value rank="0">30</integer_value>
275
 
              </restart>
276
 
            </iterative_method>
277
 
            <preconditioner name="sor"/>
278
 
            <relative_error>
279
 
              <real_value rank="0">1e-7</real_value>
280
 
            </relative_error>
281
 
            <max_iterations>
282
 
              <integer_value rank="0">3000</integer_value>
283
 
            </max_iterations>
284
 
            <never_ignore_solver_failures/>
285
 
            <diagnostics>
286
 
              <monitors/>
287
 
            </diagnostics>
288
 
          </solver>
289
 
          <initial_condition name="WholeMesh">
290
 
            <constant>
291
 
              <real_value rank="0">0.1</real_value>
292
 
            </constant>
293
 
          </initial_condition>
294
 
          <boundary_conditions name="SolidLid">
295
 
            <surface_ids>
296
 
              <integer_value shape="1" rank="1">11</integer_value>
297
 
            </surface_ids>
298
 
            <type name="dirichlet">
299
 
              <constant>
300
 
                <real_value rank="0">0</real_value>
301
 
              </constant>
302
 
            </type>
303
 
          </boundary_conditions>
304
 
          <boundary_conditions name="Reentrainment">
305
 
            <surface_ids>
306
 
              <integer_value shape="1" rank="1">8</integer_value>
307
 
            </surface_ids>
308
 
            <type name="sediment_reentrainment"/>
309
 
          </boundary_conditions>
310
 
          <tensor_field name="Diffusivity" rank="2">
311
 
            <prescribed>
312
 
              <value name="WholeMesh">
313
 
                <isotropic>
314
 
                  <constant>
315
 
                    <real_value rank="0">0.01</real_value>
316
 
                  </constant>
317
 
                </isotropic>
318
 
              </value>
319
 
              <output/>
320
 
            </prescribed>
321
 
          </tensor_field>
322
 
          <scalar_field name="SinkingVelocity" rank="0">
323
 
            <prescribed>
324
 
              <value name="WholeMesh">
325
 
                <python>
326
 
                  <string_value lines="20" type="code" language="python">def val(X,t):
327
 
        if (t &lt; 200):
328
 
                return 0.1
329
 
        else:
330
 
                return 0</string_value>
331
 
                </python>
332
 
              </value>
333
 
              <output/>
334
 
              <stat/>
335
 
              <detectors>
336
 
                <exclude_from_detectors/>
337
 
              </detectors>
338
 
            </prescribed>
339
 
          </scalar_field>
340
 
          <output/>
341
 
          <stat/>
342
 
          <convergence>
343
 
            <include_in_convergence/>
344
 
          </convergence>
345
 
          <detectors>
346
 
            <include_in_detectors/>
347
 
          </detectors>
348
 
          <steady_state>
349
 
            <include_in_steady_state/>
350
 
          </steady_state>
351
 
          <consistent_interpolation/>
352
 
        </prognostic>
353
 
        <density>
354
 
          <real_value rank="0">2000</real_value>
355
 
        </density>
356
 
        <diameter>
357
 
          <real_value rank="0">0.5</real_value>
358
 
        </diameter>
359
 
        <porosity>
360
 
          <real_value rank="0">0.5</real_value>
361
 
        </porosity>
362
 
      </scalar_field>
363
 
      <scalar_field name="SedimentFluxTemplate" rank="0">
364
 
        <diagnostic>
365
 
          <surface_ids>
366
 
            <integer_value shape="1" rank="1">8</integer_value>
367
 
          </surface_ids>
368
 
          <mesh name="VelocityMesh"/>
369
 
          <output/>
370
 
          <stat/>
371
 
          <convergence>
372
 
            <include_in_convergence/>
373
 
          </convergence>
374
 
          <detectors>
375
 
            <include_in_detectors/>
376
 
          </detectors>
377
 
          <steady_state>
378
 
            <include_in_steady_state/>
379
 
          </steady_state>
380
 
        </diagnostic>
381
 
      </scalar_field>
382
 
      <sediment_class name="Single">
383
 
        <density>
384
 
          <real_value rank="0">2000</real_value>
385
 
        </density>
386
 
        <diameter>
387
 
          <real_value rank="0">0.5</real_value>
388
 
        </diameter>
389
 
        <erodability>
390
 
          <real_value rank="0">0.00518674459</real_value>
391
 
          <comment>Fixed to give a flux of 0.01 with a bed shear stress of 1Nm-2</comment>
392
 
        </erodability>
393
 
      </sediment_class>
394
 
    </sediment>
395
 
  </material_phase>
396
 
</fluidity_options>