~reducedmodelling/fluidity/ROM_Non-intrusive-ann

« back to all changes in this revision

Viewing changes to tests/periodic_2plus1_adaptivity/periodic_2plus1_adaptivity.flml

  • Committer: pfarrell
  • Date: 2010-03-19 13:13:14 UTC
  • Revision ID: svn-v4:5bf5533e-7014-46e3-b1bb-cce4b9d03719:trunk:729
Add a 2+1 periodic adaptivity test. Not working yet

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">periodic_adaptivity</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">fluids</string_value>
 
8
  </problem_type>
 
9
  <geometry>
 
10
    <dimension replaces="NDIM">
 
11
      <integer_value rank="0">3</integer_value>
 
12
    </dimension>
 
13
    <mesh name="CoordinateMesh">
 
14
      <from_mesh>
 
15
        <mesh name="ExtrudedMesh"/>
 
16
        <periodic_boundary_conditions name="x">
 
17
          <physical_boundary_ids>
 
18
            <integer_value shape="1" rank="1">10</integer_value>
 
19
          </physical_boundary_ids>
 
20
          <aliased_boundary_ids>
 
21
            <integer_value shape="1" rank="1">8</integer_value>
 
22
          </aliased_boundary_ids>
 
23
          <coordinate_map>
 
24
            <string_value lines="20" type="python">def val(X,t):
 
25
   from math import pi
 
26
   return(X[0]-1, X[1], X[2])</string_value>
 
27
          </coordinate_map>
 
28
          <inverse_coordinate_map>
 
29
            <string_value lines="20" type="python">def val(X,t):
 
30
   from math import pi
 
31
   return(X[0]+1, X[1], X[2])</string_value>
 
32
          </inverse_coordinate_map>
 
33
          <remove_periodicity/>
 
34
        </periodic_boundary_conditions>
 
35
        <stat>
 
36
          <exclude_from_stat/>
 
37
        </stat>
 
38
      </from_mesh>
 
39
    </mesh>
 
40
    <mesh name="VelocityMesh">
 
41
      <from_mesh>
 
42
        <mesh name="ExtrudedMesh"/>
 
43
        <stat>
 
44
          <exclude_from_stat/>
 
45
        </stat>
 
46
      </from_mesh>
 
47
    </mesh>
 
48
    <mesh name="PressureMesh">
 
49
      <from_mesh>
 
50
        <mesh name="ExtrudedMesh"/>
 
51
        <stat>
 
52
          <exclude_from_stat/>
 
53
        </stat>
 
54
      </from_mesh>
 
55
    </mesh>
 
56
    <mesh name="PeriodicMesh">
 
57
      <from_mesh>
 
58
        <mesh name="ExternalMesh"/>
 
59
        <periodic_boundary_conditions name="x">
 
60
          <physical_boundary_ids>
 
61
            <integer_value shape="1" rank="1">10</integer_value>
 
62
          </physical_boundary_ids>
 
63
          <aliased_boundary_ids>
 
64
            <integer_value shape="1" rank="1">8</integer_value>
 
65
          </aliased_boundary_ids>
 
66
          <coordinate_map>
 
67
            <string_value lines="20" type="python">def val(X,t):
 
68
   from math import pi
 
69
   return(X[0]-1, X[1])</string_value>
 
70
          </coordinate_map>
 
71
          <inverse_coordinate_map>
 
72
            <string_value lines="20" type="python">def val(X,t):
 
73
   from math import pi
 
74
   return(X[0]+1, X[1])</string_value>
 
75
          </inverse_coordinate_map>
 
76
        </periodic_boundary_conditions>
 
77
        <stat>
 
78
          <exclude_from_stat/>
 
79
        </stat>
 
80
      </from_mesh>
 
81
    </mesh>
 
82
    <mesh name="ExternalMesh">
 
83
      <from_file file_name="MMS_B">
 
84
        <format name="triangle"/>
 
85
        <stat>
 
86
          <include_in_stat/>
 
87
        </stat>
 
88
      </from_file>
 
89
    </mesh>
 
90
    <mesh name="ExtrudedMesh">
 
91
      <from_mesh>
 
92
        <mesh name="PeriodicMesh"/>
 
93
        <extrude>
 
94
          <bottom_depth>
 
95
            <constant>
 
96
              <real_value rank="0">1.0</real_value>
 
97
            </constant>
 
98
          </bottom_depth>
 
99
          <sizing_function>
 
100
            <constant>
 
101
              <real_value rank="0">0.1</real_value>
 
102
            </constant>
 
103
          </sizing_function>
 
104
        </extrude>
 
105
        <stat>
 
106
          <exclude_from_stat/>
 
107
        </stat>
 
108
      </from_mesh>
 
109
    </mesh>
 
110
    <quadrature>
 
111
      <degree replaces="NGI">
 
112
        <integer_value rank="0">4</integer_value>
 
113
      </degree>
 
114
    </quadrature>
 
115
  </geometry>
 
116
  <io>
 
117
    <dump_format>
 
118
      <string_value>vtk</string_value>
 
119
    </dump_format>
 
120
    <dump_period_in_timesteps>
 
121
      <integer_value rank="0">1</integer_value>
 
122
    </dump_period_in_timesteps>
 
123
    <output_mesh name="CoordinateMesh"/>
 
124
    <stat/>
 
125
  </io>
 
126
  <timestepping>
 
127
    <current_time replaces="ACCTIM">
 
128
      <real_value rank="0">0.0</real_value>
 
129
    </current_time>
 
130
    <timestep replaces="DT">
 
131
      <real_value rank="0">1.0</real_value>
 
132
    </timestep>
 
133
    <finish_time replaces="LTIME">
 
134
      <real_value rank="0">20.0</real_value>
 
135
    </finish_time>
 
136
  </timestepping>
 
137
  <physical_parameters/>
 
138
  <material_phase name="Main">
 
139
    <vector_field name="Velocity" rank="1">
 
140
      <prescribed>
 
141
        <mesh name="VelocityMesh"/>
 
142
        <value name="WholeMesh">
 
143
          <constant>
 
144
            <real_value shape="2" dim1="dim" rank="1">0.0 0.0 0.0</real_value>
 
145
          </constant>
 
146
        </value>
 
147
        <output/>
 
148
        <stat>
 
149
          <include_in_stat/>
 
150
        </stat>
 
151
        <detectors>
 
152
          <exclude_from_detectors/>
 
153
        </detectors>
 
154
      </prescribed>
 
155
    </vector_field>
 
156
    <scalar_field name="Unity" rank="0">
 
157
      <prescribed>
 
158
        <mesh name="VelocityMesh"/>
 
159
        <value name="WholeMesh">
 
160
          <python>
 
161
            <string_value lines="20" type="python">def val(X, t):
 
162
  #return (X[0]**2 - X[0])/-0.25
 
163
  return 1.0</string_value>
 
164
          </python>
 
165
        </value>
 
166
        <output/>
 
167
        <stat/>
 
168
        <detectors>
 
169
          <exclude_from_detectors/>
 
170
        </detectors>
 
171
        <adaptivity_options>
 
172
          <absolute_measure replaces="ADOPTT = 0">
 
173
            <scalar_field replaces="ADWEIT" name="InterpolationErrorBound" rank="0">
 
174
              <prescribed>
 
175
                <value name="WholeMesh">
 
176
                  <constant>
 
177
                    <real_value rank="0">0.01</real_value>
 
178
                  </constant>
 
179
                </value>
 
180
                <output/>
 
181
                <stat/>
 
182
                <detectors>
 
183
                  <exclude_from_detectors/>
 
184
                </detectors>
 
185
              </prescribed>
 
186
            </scalar_field>
 
187
          </absolute_measure>
 
188
        </adaptivity_options>
 
189
      </prescribed>
 
190
    </scalar_field>
 
191
    <scalar_field name="PrognosticFieldConsistent" rank="0">
 
192
      <prognostic>
 
193
        <mesh name="VelocityMesh"/>
 
194
        <equation name="AdvectionDiffusion"/>
 
195
        <spatial_discretisation>
 
196
          <continuous_galerkin>
 
197
            <stabilisation>
 
198
              <no_stabilisation/>
 
199
            </stabilisation>
 
200
            <advection_terms/>
 
201
            <mass_terms/>
 
202
          </continuous_galerkin>
 
203
          <conservative_advection replaces="TBETA, DEFALT (TBETA = 0.0)">
 
204
            <real_value rank="0">1.0</real_value>
 
205
          </conservative_advection>
 
206
        </spatial_discretisation>
 
207
        <temporal_discretisation>
 
208
          <theta replaces="TTHETA, DEFALT (TTHETA = 0.5)">
 
209
            <real_value rank="0">0.5</real_value>
 
210
          </theta>
 
211
        </temporal_discretisation>
 
212
        <solver>
 
213
          <iterative_method name="gmres">
 
214
            <restart>
 
215
              <integer_value rank="0">30</integer_value>
 
216
            </restart>
 
217
          </iterative_method>
 
218
          <preconditioner name="sor"/>
 
219
          <relative_error>
 
220
            <real_value rank="0">1.0e-8</real_value>
 
221
          </relative_error>
 
222
          <max_iterations>
 
223
            <integer_value rank="0">10000</integer_value>
 
224
          </max_iterations>
 
225
          <never_ignore_solver_failures/>
 
226
          <diagnostics>
 
227
            <monitors/>
 
228
          </diagnostics>
 
229
        </solver>
 
230
        <initial_condition name="WholeMesh">
 
231
          <constant>
 
232
            <real_value rank="0">1.0</real_value>
 
233
          </constant>
 
234
        </initial_condition>
 
235
        <output/>
 
236
        <stat/>
 
237
        <convergence>
 
238
          <include_in_convergence/>
 
239
        </convergence>
 
240
        <detectors>
 
241
          <include_in_detectors/>
 
242
        </detectors>
 
243
        <steady_state>
 
244
          <include_in_steady_state/>
 
245
        </steady_state>
 
246
        <consistent_interpolation/>
 
247
      </prognostic>
 
248
    </scalar_field>
 
249
  </material_phase>
 
250
  <mesh_adaptivity>
 
251
    <hr_adaptivity replaces="ADMESH">
 
252
      <period replaces="TIMMES">
 
253
        <real_value rank="0">2.0</real_value>
 
254
      </period>
 
255
      <maximum_number_of_nodes replaces="MXNODS">
 
256
        <integer_value rank="0">10000</integer_value>
 
257
      </maximum_number_of_nodes>
 
258
      <disable_gradation/>
 
259
      <tensor_field name="MinimumEdgeLengths">
 
260
        <anisotropic_symmetric>
 
261
          <constant>
 
262
            <real_value symmetric="true" dim2="dim" shape="3 3" dim1="dim" rank="2">0.01 0.0 0.0 0.0 0.01 0.0 0.0 0.0 0.01</real_value>
 
263
          </constant>
 
264
        </anisotropic_symmetric>
 
265
      </tensor_field>
 
266
      <tensor_field name="MaximumEdgeLengths">
 
267
        <anisotropic_symmetric>
 
268
          <constant>
 
269
            <real_value symmetric="true" dim2="dim" shape="3 3" dim1="dim" rank="2">0.1 0.0 0.0 0.0 0.1 0.0 0.0 0.0 0.1</real_value>
 
270
          </constant>
 
271
        </anisotropic_symmetric>
 
272
      </tensor_field>
 
273
      <vertically_structured_adaptivity/>
 
274
    </hr_adaptivity>
 
275
  </mesh_adaptivity>
 
276
</fluidity_options>