← Back to branch summary

~fluidity-core/fluidity/bug_1183080

« back to all changes in this revision

Viewing changes to tests/Helmholtz-anisotropic-scalar-mms-p1p1cg-noadv/MMS_E.flml

  • Committer: Tim Greaves
  • Date: 2013-04-19 11:53:43 UTC
  • mfrom: (3463.2.94 turbulence-clean)
  • Revision ID: tim.greaves@imperial.ac.uk-20130419115343-qq0m21v905aachuy
This merge is a replay of a previous merge which was reviewed in detail by
Christian Jacobs. There have been no changes since this was approved, hence I'm
fast-tracking this through as the concerns that led to the previous uncommit
have been resolved.

This merge brings in general updates from Jonathan Bull to his turbulence codes
and updates the BFS examples and various related tests to reflect this.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/Helmholtz-anisotropic-scalar-mms-p1p1cg-noadv/MMS_E.flml
1
 
<?xml version='1.0' encoding='utf-8'?>
2
 
<fluidity_options>
3
 
  <simulation_name>
4
 
    <string_value lines="1">MMS_E</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_mesh>
15
 
        <mesh name="InputMesh"/>
16
 
        <stat>
17
 
          <exclude_from_stat/>
18
 
        </stat>
19
 
      </from_mesh>
20
 
    </mesh>
21
 
    <mesh name="VelocityMesh">
22
 
      <from_mesh>
23
 
        <mesh name="InputMesh"/>
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
 
    <mesh name="PressureMesh">
35
 
      <from_mesh>
36
 
        <mesh name="InputMesh"/>
37
 
        <stat>
38
 
          <exclude_from_stat/>
39
 
        </stat>
40
 
      </from_mesh>
41
 
    </mesh>
42
 
    <mesh name="InputMesh">
43
 
      <from_file file_name="MMS_E">
44
 
        <format name="triangle"/>
45
 
        <stat>
46
 
          <include_in_stat/>
47
 
        </stat>
48
 
      </from_file>
49
 
    </mesh>
50
 
    <quadrature>
51
 
      <degree>
52
 
        <integer_value rank="0">4</integer_value>
53
 
      </degree>
54
 
    </quadrature>
55
 
  </geometry>
56
 
  <io>
57
 
    <dump_format>
58
 
      <string_value>vtk</string_value>
59
 
    </dump_format>
60
 
    <dump_period>
61
 
      <constant>
62
 
        <real_value rank="0">25</real_value>
63
 
      </constant>
64
 
    </dump_period>
65
 
    <output_mesh name="VelocityMesh"/>
66
 
    <stat/>
67
 
  </io>
68
 
  <timestepping>
69
 
    <current_time>
70
 
      <real_value rank="0">0.0</real_value>
71
 
    </current_time>
72
 
    <timestep>
73
 
      <real_value rank="0">0.01</real_value>
74
 
    </timestep>
75
 
    <finish_time>
76
 
      <real_value rank="0">1</real_value>
77
 
    </finish_time>
78
 
    <nonlinear_iterations>
79
 
      <integer_value rank="0">2</integer_value>
80
 
    </nonlinear_iterations>
81
 
    <steady_state>
82
 
      <tolerance>
83
 
        <real_value rank="0">1.e-8</real_value>
84
 
        <infinity_norm/>
85
 
      </tolerance>
86
 
    </steady_state>
87
 
  </timestepping>
88
 
  <material_phase name="NS">
89
 
    <vector_field name="Velocity" rank="1">
90
 
      <prescribed>
91
 
        <mesh name="VelocityMesh"/>
92
 
        <value name="WholeMesh">
93
 
          <constant>
94
 
            <real_value shape="2" dim1="dim" rank="1">0. 0.</real_value>
95
 
          </constant>
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="Tracer" rank="0">
107
 
      <prognostic>
108
 
        <mesh name="VelocityMesh"/>
109
 
        <equation name="AdvectionDiffusion"/>
110
 
        <spatial_discretisation>
111
 
          <continuous_galerkin>
112
 
            <stabilisation>
113
 
              <no_stabilisation/>
114
 
            </stabilisation>
115
 
            <advection_terms/>
116
 
            <mass_terms/>
117
 
          </continuous_galerkin>
118
 
          <conservative_advection>
119
 
            <real_value rank="0">0.</real_value>
120
 
          </conservative_advection>
121
 
        </spatial_discretisation>
122
 
        <temporal_discretisation>
123
 
          <theta>
124
 
            <real_value rank="0">0.5</real_value>
125
 
          </theta>
126
 
        </temporal_discretisation>
127
 
        <solver>
128
 
          <iterative_method name="gmres">
129
 
            <restart>
130
 
              <integer_value rank="0">30</integer_value>
131
 
            </restart>
132
 
          </iterative_method>
133
 
          <preconditioner name="sor"/>
134
 
          <relative_error>
135
 
            <real_value rank="0">1.e-7</real_value>
136
 
          </relative_error>
137
 
          <max_iterations>
138
 
            <integer_value rank="0">1000</integer_value>
139
 
          </max_iterations>
140
 
          <never_ignore_solver_failures/>
141
 
          <diagnostics>
142
 
            <monitors/>
143
 
          </diagnostics>
144
 
        </solver>
145
 
        <initial_condition name="WholeMesh">
146
 
          <python>
147
 
            <string_value lines="20" type="code" language="python">def val(X,t):
148
 
    from math import pi, exp 
149
 
    x = X[0]; y = X[1]
150
 
    # Gaussian
151
 
    r=(x**2+y**2)**.5
152
 
    u=exp(-pi*r**2)
153
 
    return u</string_value>
154
 
          </python>
155
 
        </initial_condition>
156
 
        <boundary_conditions name="sides">
157
 
          <surface_ids>
158
 
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
159
 
          </surface_ids>
160
 
          <type name="dirichlet">
161
 
            <python>
162
 
              <string_value lines="20" type="code" language="python">def val(X,t):
163
 
    from math import pi, exp 
164
 
    x = X[0]; y = X[1]
165
 
    # Gaussian
166
 
    r=(x**2+y**2)**.5
167
 
    u=exp(-pi*r**2)
168
 
    return u</string_value>
169
 
            </python>
170
 
          </type>
171
 
        </boundary_conditions>
172
 
        <output/>
173
 
        <stat/>
174
 
        <convergence>
175
 
          <include_in_convergence/>
176
 
        </convergence>
177
 
        <detectors>
178
 
          <include_in_detectors/>
179
 
        </detectors>
180
 
        <steady_state>
181
 
          <include_in_steady_state/>
182
 
        </steady_state>
183
 
        <consistent_interpolation/>
184
 
      </prognostic>
185
 
    </scalar_field>
186
 
    <scalar_field name="FilteredTracer" rank="0">
187
 
      <diagnostic>
188
 
        <algorithm source_field_type="scalar" material_phase_support="single" name="helmholtz_anisotropic_smoothed_scalar" source_field_name="Tracer">
189
 
          <smoothing_length_scale>
190
 
            <real_value rank="0">2.</real_value>
191
 
          </smoothing_length_scale>
192
 
          <solver>
193
 
            <iterative_method name="cg"/>
194
 
            <preconditioner name="sor"/>
195
 
            <relative_error>
196
 
              <real_value rank="0">1.e-7</real_value>
197
 
            </relative_error>
198
 
            <max_iterations>
199
 
              <integer_value rank="0">1000</integer_value>
200
 
            </max_iterations>
201
 
            <never_ignore_solver_failures/>
202
 
            <diagnostics>
203
 
              <monitors/>
204
 
            </diagnostics>
205
 
          </solver>
206
 
        </algorithm>
207
 
        <mesh name="VelocityMesh"/>
208
 
        <output/>
209
 
        <stat/>
210
 
        <convergence>
211
 
          <include_in_convergence/>
212
 
        </convergence>
213
 
        <detectors>
214
 
          <include_in_detectors/>
215
 
        </detectors>
216
 
        <steady_state>
217
 
          <include_in_steady_state/>
218
 
        </steady_state>
219
 
      </diagnostic>
220
 
    </scalar_field>
221
 
    <scalar_field name="Difference" rank="0">
222
 
      <diagnostic>
223
 
        <algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="Tracer" source_field_2_name="FilteredTracer" material_phase_support="single" source_field_1_type="scalar"/>
224
 
        <mesh name="VelocityMesh"/>
225
 
        <output/>
226
 
        <stat/>
227
 
        <convergence>
228
 
          <include_in_convergence/>
229
 
        </convergence>
230
 
        <detectors>
231
 
          <include_in_detectors/>
232
 
        </detectors>
233
 
        <steady_state>
234
 
          <include_in_steady_state/>
235
 
        </steady_state>
236
 
      </diagnostic>
237
 
    </scalar_field>
238
 
  </material_phase>
239
 
</fluidity_options>