← Back to branch summary

~fluidity-core/fluidity/sediment

« back to all changes in this revision

Viewing changes to tests/Helmholtz-mms-p2p1cg-noadv/MMS_D.flml

  • Committer: Sam Parkinson
  • Date: 2013-05-22 12:32:03 UTC
  • mfrom: (3405.125.284 clean-trunk)
  • Revision ID: s.parkinson11@imperial.ac.uk-20130522123203-0c8jq0a41obsm6na
merge with trunk

Show diffs side-by-side

added added

removed removed

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