← Back to branch summary

~fluidity-core/fluidity/exorcised

« back to all changes in this revision

Viewing changes to tests/rad_1d_2grp_hom_zerobc_eig/rad_1d_2grp_hom_zerobc_eig_A.flml

Latest trunk changes, and modifications to the backward_facing_step_2d example, which are experimental and work is ongoing. It's better than the current trunk example in that at least the k_epsilon code is correct.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/rad_1d_2grp_hom_zerobc_eig/rad_1d_2grp_hom_zerobc_eig_A.flml
1
 
<?xml version='1.0' encoding='utf-8'?>
2
 
<fluidity_options>
3
 
  <simulation_name>
4
 
    <string_value lines="1">rad_1d_2grp_hom_zerobc_eig_A</string_value>
5
 
    <comment>Radiation problem:
6
 
- 1 dim
7
 
- 2 group
8
 
- 1 material (homogeneous)
9
 
- all zero BC
10
 
- eigenvalue run</comment>
11
 
  </simulation_name>
12
 
  <problem_type>
13
 
    <string_value lines="1">fluids</string_value>
14
 
  </problem_type>
15
 
  <geometry>
16
 
    <dimension>
17
 
      <integer_value rank="0">1</integer_value>
18
 
    </dimension>
19
 
    <mesh name="CoordinateMesh">
20
 
      <from_file file_name="rad_1d_2grp_hom_zerobc_eig_A">
21
 
        <format name="triangle"/>
22
 
        <stat>
23
 
          <include_in_stat/>
24
 
        </stat>
25
 
      </from_file>
26
 
    </mesh>
27
 
    <mesh name="VelocityMesh">
28
 
      <from_mesh>
29
 
        <mesh name="CoordinateMesh"/>
30
 
        <stat>
31
 
          <exclude_from_stat/>
32
 
        </stat>
33
 
      </from_mesh>
34
 
    </mesh>
35
 
    <mesh name="ParticleMesh">
36
 
      <from_mesh>
37
 
        <mesh name="CoordinateMesh"/>
38
 
        <stat>
39
 
          <exclude_from_stat/>
40
 
        </stat>
41
 
      </from_mesh>
42
 
    </mesh>
43
 
    <mesh name="ParticleMaterialMesh">
44
 
      <from_mesh>
45
 
        <mesh name="CoordinateMesh"/>
46
 
        <mesh_shape>
47
 
          <polynomial_degree>
48
 
            <integer_value rank="0">0</integer_value>
49
 
          </polynomial_degree>
50
 
        </mesh_shape>
51
 
        <mesh_continuity>
52
 
          <string_value>discontinuous</string_value>
53
 
        </mesh_continuity>
54
 
        <stat>
55
 
          <exclude_from_stat/>
56
 
        </stat>
57
 
      </from_mesh>
58
 
    </mesh>
59
 
    <quadrature>
60
 
      <degree>
61
 
        <integer_value rank="0">2</integer_value>
62
 
      </degree>
63
 
    </quadrature>
64
 
  </geometry>
65
 
  <io>
66
 
    <dump_format>
67
 
      <string_value>vtk</string_value>
68
 
    </dump_format>
69
 
    <dump_period>
70
 
      <constant>
71
 
        <real_value rank="0">1.0</real_value>
72
 
      </constant>
73
 
    </dump_period>
74
 
    <output_mesh name="CoordinateMesh"/>
75
 
    <stat/>
76
 
  </io>
77
 
  <timestepping>
78
 
    <current_time>
79
 
      <real_value rank="0">0.0</real_value>
80
 
    </current_time>
81
 
    <timestep>
82
 
      <real_value rank="0">1.0</real_value>
83
 
    </timestep>
84
 
    <finish_time>
85
 
      <real_value rank="0">1.0</real_value>
86
 
    </finish_time>
87
 
    <final_timestep>
88
 
      <integer_value rank="0">0</integer_value>
89
 
    </final_timestep>
90
 
  </timestepping>
91
 
  <material_phase name="dummy">
92
 
    <vector_field name="Velocity" rank="1">
93
 
      <prescribed>
94
 
        <mesh name="VelocityMesh"/>
95
 
        <value name="WholeMesh">
96
 
          <constant>
97
 
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
98
 
          </constant>
99
 
        </value>
100
 
        <output/>
101
 
        <stat>
102
 
          <include_in_stat/>
103
 
        </stat>
104
 
        <detectors>
105
 
          <exclude_from_detectors/>
106
 
        </detectors>
107
 
      </prescribed>
108
 
    </vector_field>
109
 
  </material_phase>
110
 
  <embedded_models>
111
 
    <radiation>
112
 
      <particle_type name="Neutron">
113
 
        <equation name="Eigenvalue">
114
 
          <power_iteration>
115
 
            <maximum>
116
 
              <integer_value rank="0">100</integer_value>
117
 
            </maximum>
118
 
            <flux_tolerance_absolute>
119
 
              <real_value rank="0">1.0e-6</real_value>
120
 
            </flux_tolerance_absolute>
121
 
            <keff_tolerance_relative>
122
 
              <real_value rank="0">1.0e-6</real_value>
123
 
            </keff_tolerance_relative>
124
 
          </power_iteration>
125
 
          <flux_normalisation name="TotalFlux">
126
 
            <value>
127
 
              <real_value rank="0">1.0</real_value>
128
 
            </value>
129
 
          </flux_normalisation>
130
 
        </equation>
131
 
        <energy_discretisation name="WholeSpatialMesh">
132
 
          <continuity>
133
 
            <string_value>discontinuous</string_value>
134
 
          </continuity>
135
 
          <degree>
136
 
            <integer_value rank="0">0</integer_value>
137
 
          </degree>
138
 
          <energy_group_set name="all">
139
 
            <number_of_energy_groups>
140
 
              <integer_value rank="0">2</integer_value>
141
 
            </number_of_energy_groups>
142
 
            <angular_discretisation name="WholeSpatialMesh">
143
 
              <method name="SphericalHarmonic">
144
 
                <parity name="Even">
145
 
                  <exclude_odd_parity_time_source/>
146
 
                  <angular_moment_set name="AllMoments">
147
 
                    <number_of_angular_moments>
148
 
                      <integer_value rank="0">1</integer_value>
149
 
                    </number_of_angular_moments>
150
 
                    <mesh name="ParticleMaterialMesh"/>
151
 
                    <tensor_field name="Diffusivity">
152
 
                      <diagnostic>
153
 
                        <algorithm name="Internal"/>
154
 
                        <output>
155
 
                          <exclude_from_vtu/>
156
 
                        </output>
157
 
                        <stat>
158
 
                          <exclude_from_stat/>
159
 
                        </stat>
160
 
                        <detectors>
161
 
                          <exclude_from_detectors/>
162
 
                        </detectors>
163
 
                      </diagnostic>
164
 
                    </tensor_field>
165
 
                    <scalar_field name="Absorption">
166
 
                      <diagnostic>
167
 
                        <algorithm name="Internal"/>
168
 
                        <output>
169
 
                          <exclude_from_vtu/>
170
 
                        </output>
171
 
                        <stat>
172
 
                          <exclude_from_stat/>
173
 
                        </stat>
174
 
                        <detectors>
175
 
                          <exclude_from_detectors/>
176
 
                        </detectors>
177
 
                      </diagnostic>
178
 
                    </scalar_field>
179
 
                    <scalar_field name="ParticleFlux">
180
 
                      <prognostic>
181
 
                        <mesh name="ParticleMesh"/>
182
 
                        <equation name="AdvectionDiffusion"/>
183
 
                        <spatial_discretisation>
184
 
                          <continuous_galerkin>
185
 
                            <advection_terms>
186
 
                              <exclude_advection_terms/>
187
 
                            </advection_terms>
188
 
                            <mass_terms>
189
 
                              <exclude_mass_terms/>
190
 
                            </mass_terms>
191
 
                          </continuous_galerkin>
192
 
                          <conservative_advection>
193
 
                            <real_value rank="0">1.0</real_value>
194
 
                          </conservative_advection>
195
 
                        </spatial_discretisation>
196
 
                        <temporal_discretisation>
197
 
                          <theta>
198
 
                            <real_value rank="0">1.0</real_value>
199
 
                          </theta>
200
 
                        </temporal_discretisation>
201
 
                        <solver>
202
 
                          <iterative_method name="cg"/>
203
 
                          <preconditioner name="sor"/>
204
 
                          <relative_error>
205
 
                            <real_value rank="0">1.0e-10</real_value>
206
 
                          </relative_error>
207
 
                          <absolute_error>
208
 
                            <real_value rank="0">1.0e-10</real_value>
209
 
                          </absolute_error>
210
 
                          <max_iterations>
211
 
                            <integer_value rank="0">1000</integer_value>
212
 
                          </max_iterations>
213
 
                          <never_ignore_solver_failures/>
214
 
                          <diagnostics>
215
 
                            <monitors/>
216
 
                          </diagnostics>
217
 
                        </solver>
218
 
                        <initial_condition name="WholeMesh">
219
 
                          <constant>
220
 
                            <real_value rank="0">1.0</real_value>
221
 
                          </constant>
222
 
                        </initial_condition>
223
 
                        <boundary_conditions name="left">
224
 
                          <surface_ids>
225
 
                            <integer_value shape="1" rank="1">1</integer_value>
226
 
                          </surface_ids>
227
 
                          <type name="dirichlet">
228
 
                            <constant>
229
 
                              <real_value rank="0">0.0</real_value>
230
 
                            </constant>
231
 
                          </type>
232
 
                        </boundary_conditions>
233
 
                        <boundary_conditions name="right">
234
 
                          <surface_ids>
235
 
                            <integer_value shape="1" rank="1">2</integer_value>
236
 
                          </surface_ids>
237
 
                          <type name="dirichlet">
238
 
                            <constant>
239
 
                              <real_value rank="0">0.0</real_value>
240
 
                            </constant>
241
 
                          </type>
242
 
                        </boundary_conditions>
243
 
                        <output/>
244
 
                        <stat/>
245
 
                        <detectors>
246
 
                          <include_in_detectors/>
247
 
                        </detectors>
248
 
                        <consistent_interpolation/>
249
 
                      </prognostic>
250
 
                    </scalar_field>
251
 
                  </angular_moment_set>
252
 
                </parity>
253
 
              </method>
254
 
            </angular_discretisation>
255
 
          </energy_group_set>
256
 
        </energy_discretisation>
257
 
        <material_data_set name="dataset1">
258
 
          <from_file file_name="rad_1d_2grp_hom_zerobc_eig.radmats">
259
 
            <number_of_scatter_moments>
260
 
              <integer_value rank="0">1</integer_value>
261
 
            </number_of_scatter_moments>
262
 
            <format_radmats>
263
 
              <maximum_record_length>
264
 
                <integer_value rank="0">132</integer_value>
265
 
              </maximum_record_length>
266
 
            </format_radmats>
267
 
            <physical_material name="fuel">
268
 
              <interpolation_dimension name="first">
269
 
                <interpolation_attribute_none/>
270
 
                <interpolation_values>
271
 
                  <real_value shape="1" rank="1">1.0</real_value>
272
 
                </interpolation_values>
273
 
              </interpolation_dimension>
274
 
            </physical_material>
275
 
          </from_file>
276
 
        </material_data_set>
277
 
        <region_id_material_mapping>
278
 
          <region_to_physical_radiation_material_map name="whole">
279
 
            <region_id>
280
 
              <integer_value shape="1" rank="1">1</integer_value>
281
 
            </region_id>
282
 
            <data_set name="dataset1"/>
283
 
            <physical_material name="fuel"/>
284
 
          </region_to_physical_radiation_material_map>
285
 
        </region_id_material_mapping>
286
 
      </particle_type>
287
 
    </radiation>
288
 
  </embedded_models>
289
 
</fluidity_options>