← Back to branch summary

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

« back to all changes in this revision

Viewing changes to tests/mms_tracer_cv_fesweby_adv_steady/MMS_D_cv.flml

  • Committer: Simon Mouradian
  • Date: 2012-03-30 11:11:47 UTC
  • mfrom: (3520.32.239 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120330111147-y4n4y4th52x2dqcf
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mms_tracer_cv_fesweby_adv_steady/MMS_D_cv.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">MMS_D_cv</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="MMS_D">
 
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="CoordinateMesh"/>
 
24
        <stat>
 
25
          <exclude_from_stat/>
 
26
        </stat>
 
27
      </from_mesh>
 
28
    </mesh>
 
29
    <mesh name="PressureMesh">
 
30
      <from_mesh>
 
31
        <mesh name="CoordinateMesh"/>
 
32
        <stat>
 
33
          <exclude_from_stat/>
 
34
        </stat>
 
35
      </from_mesh>
 
36
    </mesh>
 
37
    <quadrature>
 
38
      <degree>
 
39
        <integer_value rank="0">5</integer_value>
 
40
      </degree>
 
41
      <controlvolume_surface_degree>
 
42
        <integer_value rank="0">5</integer_value>
 
43
      </controlvolume_surface_degree>
 
44
    </quadrature>
 
45
  </geometry>
 
46
  <io>
 
47
    <dump_format>
 
48
      <string_value>vtk</string_value>
 
49
    </dump_format>
 
50
    <dump_period>
 
51
      <constant>
 
52
        <real_value rank="0">1000.0</real_value>
 
53
      </constant>
 
54
    </dump_period>
 
55
    <output_mesh name="VelocityMesh"/>
 
56
    <stat/>
 
57
  </io>
 
58
  <timestepping>
 
59
    <current_time>
 
60
      <real_value rank="0">0.0</real_value>
 
61
    </current_time>
 
62
    <timestep>
 
63
      <real_value rank="0">0.0041263986301734437</real_value>
 
64
      <comment>0.0020631993150867218 gives a cv cfl number of 0.5</comment>
 
65
    </timestep>
 
66
    <finish_time>
 
67
      <real_value rank="0">20.0</real_value>
 
68
    </finish_time>
 
69
    <steady_state>
 
70
      <tolerance>
 
71
        <real_value rank="0">1.E-5</real_value>
 
72
        <infinity_norm/>
 
73
      </tolerance>
 
74
      <acceleration_form/>
 
75
      <steady_state_file>
 
76
        <plain_text_output/>
 
77
      </steady_state_file>
 
78
    </steady_state>
 
79
  </timestepping>
 
80
  <material_phase name="Burgers">
 
81
    <vector_field name="Velocity" rank="1">
 
82
      <prescribed>
 
83
        <mesh name="VelocityMesh"/>
 
84
        <value name="WholeMesh">
 
85
          <python>
 
86
            <string_value lines="20" type="code" language="python">def val(XX, t):
 
87
   from math import sin,cos
 
88
   x = XX[0];
 
89
   y = XX[1];
 
90
   x2 = x*x;
 
91
   y2 = y*y;
 
92
   u = sin(5*(x2+y2));
 
93
   v = cos(3*(x2-y2));
 
94
   return (u, v)</string_value>
 
95
          </python>
 
96
        </value>
 
97
        <output/>
 
98
        <stat>
 
99
          <include_in_stat/>
 
100
        </stat>
 
101
        <detectors>
 
102
          <exclude_from_detectors/>
 
103
        </detectors>
 
104
      </prescribed>
 
105
    </vector_field>
 
106
    <scalar_field name="NumericalSolution" rank="0">
 
107
      <prognostic>
 
108
        <mesh name="VelocityMesh"/>
 
109
        <equation name="AdvectionDiffusion"/>
 
110
        <spatial_discretisation>
 
111
          <control_volumes>
 
112
            <face_value name="FiniteElement">
 
113
              <limit_face_value>
 
114
                <limiter name="Sweby">
 
115
                  <project_upwind_value_from_point>
 
116
                    <reflect_off_domain_boundaries/>
 
117
                    <bound_projected_value_locally/>
 
118
                    <store_upwind_elements>
 
119
                      <store_upwind_quadrature/>
 
120
                    </store_upwind_elements>
 
121
                  </project_upwind_value_from_point>
 
122
                </limiter>
 
123
              </limit_face_value>
 
124
            </face_value>
 
125
            <diffusion_scheme name="BassiRebay"/>
 
126
          </control_volumes>
 
127
          <conservative_advection>
 
128
            <real_value rank="0">0.0</real_value>
 
129
          </conservative_advection>
 
130
        </spatial_discretisation>
 
131
        <temporal_discretisation>
 
132
          <theta>
 
133
            <real_value rank="0">1.0</real_value>
 
134
          </theta>
 
135
          <control_volumes>
 
136
            <number_advection_iterations>
 
137
              <integer_value rank="0">3</integer_value>
 
138
            </number_advection_iterations>
 
139
            <pivot_theta>
 
140
              <real_value rank="0">1.0</real_value>
 
141
            </pivot_theta>
 
142
          </control_volumes>
 
143
        </temporal_discretisation>
 
144
        <solver>
 
145
          <iterative_method name="gmres">
 
146
            <restart>
 
147
              <integer_value rank="0">30</integer_value>
 
148
            </restart>
 
149
          </iterative_method>
 
150
          <preconditioner name="sor"/>
 
151
          <relative_error>
 
152
            <real_value rank="0">1.0e-10</real_value>
 
153
          </relative_error>
 
154
          <max_iterations>
 
155
            <integer_value rank="0">1000</integer_value>
 
156
          </max_iterations>
 
157
          <never_ignore_solver_failures/>
 
158
          <diagnostics>
 
159
            <monitors/>
 
160
          </diagnostics>
 
161
        </solver>
 
162
        <initial_condition name="WholeMesh">
 
163
          <constant>
 
164
            <real_value rank="0">0.0</real_value>
 
165
          </constant>
 
166
        </initial_condition>
 
167
        <boundary_conditions name="sides">
 
168
          <surface_ids>
 
169
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
 
170
          </surface_ids>
 
171
          <type name="dirichlet">
 
172
            <apply_weakly/>
 
173
            <python>
 
174
              <string_value lines="20" type="code" language="python">def val(XX, t):
 
175
   from math import sin,cos,sqrt
 
176
   omega = 0.0;
 
177
   x = XX[0];
 
178
   y = XX[1];
 
179
   x4 = x*x*x*x;
 
180
   y2 = y*y;
 
181
   u = sin(25*x*y+omega*t)-2*y/(sqrt(x));
 
182
   return u</string_value>
 
183
            </python>
 
184
          </type>
 
185
        </boundary_conditions>
 
186
        <scalar_field name="Source" rank="0">
 
187
          <prescribed>
 
188
            <value name="WholeMesh">
 
189
              <python>
 
190
                <string_value lines="20" type="code" language="python">def val(XX, t):
 
191
   from math import sin,cos,sqrt
 
192
   omega = 0.0;
 
193
   x = XX[0];
 
194
   y = XX[1];
 
195
   x2 = x*x;
 
196
   y2 = y*y;
 
197
   xp5 = sqrt(x)
 
198
   x1p5 = xp5*x
 
199
   S = (25*y*cos(25*x*y) + 1.0*y/x1p5)*sin(5*(y2 + x2)) + (25*x*cos(25*x*y) - 2/xp5)*cos(3*(x2 - y2));
 
200
   return S</string_value>
 
201
              </python>
 
202
            </value>
 
203
            <output/>
 
204
            <stat/>
 
205
            <detectors>
 
206
              <exclude_from_detectors/>
 
207
            </detectors>
 
208
          </prescribed>
 
209
        </scalar_field>
 
210
        <output>
 
211
          <include_previous_time_step/>
 
212
        </output>
 
213
        <stat/>
 
214
        <convergence>
 
215
          <include_in_convergence/>
 
216
        </convergence>
 
217
        <detectors>
 
218
          <include_in_detectors/>
 
219
        </detectors>
 
220
        <steady_state>
 
221
          <include_in_steady_state/>
 
222
        </steady_state>
 
223
        <consistent_interpolation/>
 
224
      </prognostic>
 
225
    </scalar_field>
 
226
    <scalar_field name="AnalyticalSolution" rank="0">
 
227
      <prescribed>
 
228
        <mesh name="VelocityMesh"/>
 
229
        <value name="WholeMesh">
 
230
          <python>
 
231
            <string_value lines="20" type="code" language="python">def val(XX, t):
 
232
   from math import sin,cos,sqrt
 
233
   omega = 0.0;
 
234
   x = XX[0];
 
235
   y = XX[1];
 
236
   u = sin(25*x*y + omega*t) - 2*y/(sqrt(x));
 
237
   return u</string_value>
 
238
          </python>
 
239
        </value>
 
240
        <output/>
 
241
        <stat/>
 
242
        <detectors>
 
243
          <exclude_from_detectors/>
 
244
        </detectors>
 
245
      </prescribed>
 
246
    </scalar_field>
 
247
    <scalar_field name="ControlVolumeCFLNumber" rank="0">
 
248
      <diagnostic>
 
249
        <algorithm name="Internal" material_phase_support="multiple"/>
 
250
        <mesh name="VelocityMesh"/>
 
251
        <output/>
 
252
        <stat/>
 
253
        <convergence>
 
254
          <include_in_convergence/>
 
255
        </convergence>
 
256
        <detectors>
 
257
          <include_in_detectors/>
 
258
        </detectors>
 
259
        <steady_state>
 
260
          <include_in_steady_state/>
 
261
        </steady_state>
 
262
      </diagnostic>
 
263
    </scalar_field>
 
264
    <scalar_field name="AbsoluteDifference" rank="0">
 
265
      <diagnostic field_name_b="NumericalSolution" field_name_a="AnalyticalSolution">
 
266
        <algorithm name="Internal" material_phase_support="multiple"/>
 
267
        <mesh name="VelocityMesh"/>
 
268
        <output/>
 
269
        <stat>
 
270
          <include_cv_stats/>
 
271
        </stat>
 
272
        <convergence>
 
273
          <include_in_convergence/>
 
274
        </convergence>
 
275
        <detectors>
 
276
          <include_in_detectors/>
 
277
        </detectors>
 
278
        <steady_state>
 
279
          <include_in_steady_state/>
 
280
        </steady_state>
 
281
      </diagnostic>
 
282
    </scalar_field>
 
283
  </material_phase>
 
284
</fluidity_options>