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