forked from FluidityProject/fluidity
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtemporal_discretisation.rng
120 lines (117 loc) · 4.83 KB
/
temporal_discretisation.rng
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
<?xml version="1.0" encoding="UTF-8"?>
<grammar xmlns:a="http://relaxng.org/ns/compatibility/annotations/1.0" xmlns="http://relaxng.org/ns/structure/1.0">
<define name="temporal_discretisation_options">
<element name="temporal_discretisation">
<element name="theta">
<a:documentation>Implicit/explicit control (TTHETA)
=0. -- explicit
=0.5 -- Crank-Nicolson
=1. -- implicit</a:documentation>
<ref name="real"/>
</element>
<optional>
<ref name="temporal_control_volume_options"/>
</optional>
<optional>
<ref name="temporal_discontinuous_galerkin_options"/>
</optional>
</element>
</define>
<define name="temporal_control_volume_options">
<element name="control_volumes">
<a:documentation>Temporal discretisation options that are only relevant if a control volume or coupled_cv spatial discretisation is selected for this field.</a:documentation>
<optional>
<element name="number_advection_iterations">
<a:documentation>Number of iterations within an advection solve.
This increases the accuracy of the face values and ensures that
the pivoted solution is cancelled out.
Defaults to 1 if unselected.</a:documentation>
<ref name="integer"/>
<optional>
<element name="tolerance">
<a:documentation>Cut short advection_iterations if the specified tolerance
is reached.
This only works for pure control volume discretisations.</a:documentation>
<ref name="real"/>
<choice>
<element name="infinity_norm">
<a:documentation>Select the norm with which you want the tolerance to be tested.
The infinity norm.</a:documentation>
<empty/>
</element>
<element name="l2_norm">
<a:documentation>Select the norm with which you want the tolerance to be tested.
The l2 norm.</a:documentation>
<empty/>
</element>
<element name="cv_l2_norm">
<a:documentation>Select the norm with which you want the tolerance to be tested.
The l2 norm evaluated on a control volume mesh.</a:documentation>
<empty/>
</element>
</choice>
</element>
</optional>
</element>
</optional>
<optional>
<choice>
<element name="maximum_courant_number_per_subcycle">
<a:documentation>Use timestep subcycling to solve this equation.
Specify the maximum courant number per subcycle.
This only works for pure control volume discretisations.</a:documentation>
<ref name="real"/>
<ref name="field_based_cfl_number_options"/>
</element>
<element name="number_advection_subcycles">
<a:documentation>Use timestep subcycling to solve this equation.
Specify the number of subcycles.
This only works for pure control volume discretisations.</a:documentation>
<ref name="integer"/>
</element>
</choice>
</optional>
<optional>
<element name="limit_theta">
<a:documentation>Only works if a control volume or coupled_cv spatial discretisation is selected.
If not active then the theta specified above will be used.
Otherwise use variable limited theta on individual faces.</a:documentation>
<empty/>
</element>
</optional>
<optional>
<element name="pivot_theta">
<a:documentation>Only works if a control volume or coupled_cv spatial discretisation is selected.
Time discretisation of upwind discretisation off which the
higher order solution is pivotted.
- pivot_theta = 1 - implicit pivot (default if not set and
best choice if not intentionally modifying
scheme to be explicit)
- pivot_theta = 0 - explicit pivot</a:documentation>
<ref name="real"/>
</element>
</optional>
</element>
</define>
<define name="temporal_discontinuous_galerkin_options">
<element name="discontinuous_galerkin">
<a:documentation>This enables DG-specific timestepping options, such as
explicit advection subcycling. </a:documentation>
<optional>
<choice>
<element name="maximum_courant_number_per_subcycle">
<a:documentation>Use timestep subcycling to solve this equation.
Specify the maximum courant number per subcycle.</a:documentation>
<ref name="real"/>
<ref name="dg_field_based_cfl_number_options"/>
</element>
<element name="number_advection_subcycles">
<a:documentation>Use timestep subcycling to solve this equation.
Specify the number of subcycles.</a:documentation>
<ref name="integer"/>
</element>
</choice>
</optional>
</element>
</define>
</grammar>