← Back to branch summary

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

« back to all changes in this revision

Viewing changes to examples/backward_facing_step_2d/backward_facing_step_2d_kim_reference.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:
examples/backward_facing_step_2d/backward_facing_step_2d_kim_reference.flml
57
57
      <real_value rank="0">0.0</real_value>
58
58
    </current_time>
59
59
    <timestep>
60
 
      <real_value rank="0">0.01</real_value>
 
60
      <real_value rank="0">0.02</real_value>
61
61
    </timestep>
62
62
    <finish_time>
63
 
      <real_value rank="0">50.0</real_value>
 
63
      <real_value rank="0">300.0</real_value>
64
64
    </finish_time>
65
65
    <nonlinear_iterations>
66
66
      <integer_value rank="0">2</integer_value>
67
67
    </nonlinear_iterations>
68
 
    <adaptive_timestep>
69
 
      <requested_cfl>
70
 
        <real_value rank="0">1.0</real_value>
71
 
      </requested_cfl>
72
 
      <courant_number name="CFLNumber">
73
 
        <mesh name="VelocityMesh"/>
74
 
      </courant_number>
75
 
    </adaptive_timestep>
76
68
    <steady_state>
77
69
      <tolerance>
78
 
        <real_value rank="0">1.e-5</real_value>
 
70
        <real_value rank="0">1.e-4</real_value>
79
71
        <infinity_norm/>
80
72
      </tolerance>
81
73
      <steady_state_file>
112
104
            <monitors/>
113
105
          </diagnostics>
114
106
        </solver>
115
 
        <boundary_conditions name="out">
116
 
          <surface_ids>
117
 
            <integer_value shape="1" rank="1">11</integer_value>
118
 
          </surface_ids>
119
 
          <type name="dirichlet">
120
 
            <constant>
121
 
              <real_value rank="0">0.0</real_value>
122
 
            </constant>
123
 
          </type>
124
 
        </boundary_conditions>
125
107
        <output/>
126
108
        <stat/>
127
109
        <convergence>
143
125
        <spatial_discretisation>
144
126
          <continuous_galerkin>
145
127
            <stabilisation>
146
 
              <no_stabilisation/>
 
128
              <streamline_upwind>
 
129
                <nu_bar_optimal/>
 
130
                <nu_scale name="0.5">
 
131
                  <real_value shape="1" rank="0">0.5</real_value>
 
132
                </nu_scale>
 
133
              </streamline_upwind>
147
134
            </stabilisation>
148
135
            <mass_terms>
149
136
              <lump_mass_matrix/>
222
209
          <surface_ids>
223
210
            <integer_value shape="2" rank="1">10 12</integer_value>
224
211
          </surface_ids>
 
212
          <type name="log_law_of_wall">
 
213
            <surface_roughness>
 
214
              <real_value rank="0">1.e-6</real_value>
 
215
            </surface_roughness>
 
216
          </type>
 
217
        </boundary_conditions>
 
218
        <boundary_conditions name="no_normal">
 
219
          <surface_ids>
 
220
            <integer_value shape="2" rank="1">10 12</integer_value>
 
221
          </surface_ids>
225
222
          <type name="dirichlet">
226
 
            <align_bc_with_cartesian>
227
 
              <x_component>
228
 
                <constant>
229
 
                  <real_value rank="0">0.0</real_value>
230
 
                </constant>
231
 
              </x_component>
232
 
              <y_component>
233
 
                <constant>
234
 
                  <real_value rank="0">0.0</real_value>
235
 
                </constant>
236
 
              </y_component>
237
 
            </align_bc_with_cartesian>
 
223
            <align_bc_with_surface>
 
224
              <normal_component>
 
225
                <constant>
 
226
                  <real_value rank="0">0.</real_value>
 
227
                </constant>
 
228
              </normal_component>
 
229
            </align_bc_with_surface>
238
230
          </type>
239
231
        </boundary_conditions>
240
232
        <tensor_field name="Viscosity" rank="2">
242
234
            <value name="WholeMesh">
243
235
              <isotropic>
244
236
                <constant>
245
 
                  <real_value rank="0">1.24e-4</real_value>
 
237
                  <real_value rank="0">1.74e-5</real_value>
246
238
                </constant>
247
239
              </isotropic>
248
240
            </value>
249
241
            <output/>
250
242
          </prescribed>
251
243
        </tensor_field>
252
 
        <vector_field name="Absorption" rank="1">
253
 
          <prescribed>
254
 
            <value name="WholeMesh">
255
 
              <python>
256
 
                <string_value lines="20" type="code" language="python">def val(X, t):
257
 
    #absorption coefficient
258
 
    k=0.5
259
 
    if(X[0]&gt;30.0):
260
 
        A=(X[0]-30.0)/5.0*k
261
 
    else:
262
 
        A=0.0
263
 
    return [A,0.0]</string_value>
264
 
              </python>
265
 
            </value>
266
 
            <output/>
267
 
            <stat>
268
 
              <include_in_stat/>
269
 
            </stat>
270
 
            <detectors>
271
 
              <exclude_from_detectors/>
272
 
            </detectors>
273
 
          </prescribed>
274
 
          <default_absorption/>
275
 
        </vector_field>
276
244
        <output/>
277
245
        <stat>
278
246
          <include_in_stat/>
312
280
        </steady_state>
313
281
      </diagnostic>
314
282
    </scalar_field>
 
283
    <scalar_field name="AveragePressure" rank="0">
 
284
      <diagnostic>
 
285
        <algorithm source_field_type="scalar" material_phase_support="single" name="time_averaged_scalar" source_field_name="Pressure">
 
286
          <spin_up_time>
 
287
            <real_value rank="0">50</real_value>
 
288
          </spin_up_time>
 
289
        </algorithm>
 
290
        <mesh name="VelocityMesh"/>
 
291
        <output/>
 
292
        <stat/>
 
293
        <convergence>
 
294
          <include_in_convergence/>
 
295
        </convergence>
 
296
        <detectors>
 
297
          <include_in_detectors/>
 
298
        </detectors>
 
299
        <steady_state>
 
300
          <include_in_steady_state/>
 
301
        </steady_state>
 
302
      </diagnostic>
 
303
    </scalar_field>
 
304
    <scalar_field name="FluctuatingPressure" rank="0">
 
305
      <diagnostic>
 
306
        <algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="Pressure" source_field_2_name="AveragePressure" material_phase_support="single" source_field_1_type="scalar"/>
 
307
        <mesh name="VelocityMesh"/>
 
308
        <output/>
 
309
        <stat/>
 
310
        <convergence>
 
311
          <include_in_convergence/>
 
312
        </convergence>
 
313
        <detectors>
 
314
          <include_in_detectors/>
 
315
        </detectors>
 
316
        <steady_state>
 
317
          <include_in_steady_state/>
 
318
        </steady_state>
 
319
      </diagnostic>
 
320
    </scalar_field>
 
321
    <vector_field name="AverageVelocity" rank="1">
 
322
      <diagnostic>
 
323
        <algorithm source_field_type="vector" material_phase_support="single" name="time_averaged_vector" source_field_name="Velocity">
 
324
          <spin_up_time>
 
325
            <real_value rank="0">50</real_value>
 
326
          </spin_up_time>
 
327
        </algorithm>
 
328
        <mesh name="VelocityMesh"/>
 
329
        <output/>
 
330
        <stat>
 
331
          <include_in_stat/>
 
332
        </stat>
 
333
        <convergence>
 
334
          <include_in_convergence/>
 
335
        </convergence>
 
336
        <detectors>
 
337
          <include_in_detectors/>
 
338
        </detectors>
 
339
        <steady_state>
 
340
          <include_in_steady_state/>
 
341
        </steady_state>
 
342
      </diagnostic>
 
343
    </vector_field>
 
344
    <vector_field name="FluctuatingVelocity" rank="1">
 
345
      <diagnostic>
 
346
        <algorithm source_field_2_type="vector" name="vector_difference" source_field_1_name="Velocity" source_field_2_name="AverageVelocity" material_phase_support="single" source_field_1_type="vector"/>
 
347
        <mesh name="VelocityMesh"/>
 
348
        <output/>
 
349
        <stat>
 
350
          <include_in_stat/>
 
351
        </stat>
 
352
        <convergence>
 
353
          <include_in_convergence/>
 
354
        </convergence>
 
355
        <detectors>
 
356
          <include_in_detectors/>
 
357
        </detectors>
 
358
        <steady_state>
 
359
          <include_in_steady_state/>
 
360
        </steady_state>
 
361
      </diagnostic>
 
362
    </vector_field>
315
363
  </material_phase>
316
364
</fluidity_options>