← Back to branch summary

~fluidity-core/fluidity/remove_tictoc

« back to all changes in this revision

Viewing changes to tests/mms_tracer_cv_3d_adv_excludemass/MMS_E_cv.flml

  • Committer: Jon Hill
  • Date: 2012-02-28 21:36:27 UTC
  • mfrom: (3880.1.59 fluidity)
  • Revision ID: jon.hill@imperial.ac.uk-20120228213627-75tqkd0za9nlcjrt
Merge from trunk

Show diffs side-by-side

added added

removed removed

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