transient_flow_in_inclined_2D_plane.prj 7.91 KB
Newer Older
1
2
3
4
5
6
7
8
9
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
    <mesh>inclined_2D_mesh.vtu</mesh>
    <geometry>plane_r.gml</geometry>
    <processes>
        <process>
            <name>LiquidFlow</name>
            <type>LIQUID_FLOW</type>
            <integration_order>2</integration_order>
10
            <process_variables>
11
12
13
14
15
16
                <process_variable>pressure</process_variable>
            </process_variables>
            <specific_body_force>0 0 -9.8</specific_body_force>
            <secondary_variables>
                <secondary_variable internal_name="darcy_velocity" output_name="v"/>
            </secondary_variables>
17
18
19
            <aperture_size>
                <parameter>fracture_thickness</parameter>
            </aperture_size>
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
        </process>
    </processes>
    <media>
        <medium>
            <phases>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value> 1.e-3 </value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000.0</reference_value>
                            <independent_variable>
                                <variable_name>phase_pressure
39
                                </variable_name>
40
                                <reference_condition>1e5
41
                                </reference_condition>
42
                                <slope>4.6511627906976743356e-10
43
                                </slope>
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
                            </independent_variable>
                        </property>
                    </properties>
                </phase>
            </phases>
            <properties>
                <property>
                    <name>permeability</name>
                    <type>Parameter</type>
                    <parameter_name>k</parameter_name>
                </property>
                <property>
                    <name>reference_temperature</name>
                    <type>Constant</type>
                    <value>293.15</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.2</value>
                </property>
                <property>
                    <name>storage</name>
                    <type>Constant</type>
                    <value> 1.0e-10 </value>
                </property>
            </properties>
        </medium>
    </media>
    <time_loop>
        <processes>
            <process ref="LiquidFlow">
                <nonlinear_solver>basic_picard</nonlinear_solver>
                <convergence_criterion>
                    <type>DeltaX</type>
                    <norm_type>NORM2</norm_type>
                    <abstol>1.e-6</abstol>
                </convergence_criterion>
                <time_discretization>
                    <type>BackwardEuler</type>
                </time_discretization>
                <time_stepping>
                    <type>FixedTimeStepping</type>
                    <t_initial> 0.0 </t_initial>
                    <t_end> 864000 </t_end>
                    <timesteps>
                        <pair>
                            <repeat>10</repeat>
                            <delta_t>86400</delta_t>
                        </pair>
                    </timesteps>
                </time_stepping>
            </process>
        </processes>
        <output>
            <type>VTK</type>
            <prefix>transient_flow_in_inclined_2D_plane</prefix>
            <timesteps>
                <pair>
                    <repeat> 5 </repeat>
                    <each_steps> 2 </each_steps>
                </pair>
            </timesteps>
            <variables>
                <variable> pressure </variable>
                <variable> v </variable>
            </variables>
            <suffix>_ts_t_{:time}</suffix>
        </output>
    </time_loop>
    <parameters>
115
116
117
118
119
        <parameter>
            <name>fracture_thickness</name>
            <type>MeshElement</type>
            <field_name>aperture_size</field_name>
        </parameter>
120
121
122
        <parameter>
            <name>k</name>
            <type>Constant</type>
123
124
125
126
127
128
129
130
131
132
133
134
            <!--Permeability tensor in the local coordinate system on the
            inclined plane. Its components associated with the off plane
            direction have to be set to zero. In other words, the local
            permeability tensor has to input as a diagonal block of the full
            size tensor in the global coordinate system. The bases of the
            local coordinate system has to be defined in order to rotate
            this local tensor to the global coordinate system.-->
            <values>1.e-12 0.0  0.0
                0.0  5.e-12  0.0
                0.0 0.0 0.0</values>
            <use_local_coordinate_system>true
            </use_local_coordinate_system>
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
        </parameter>
        <parameter>
            <name>p0</name>
            <type>Constant</type>
            <value>0</value>
        </parameter>
        <parameter>
            <name>p_top</name>
            <type>Constant</type>
            <value>1.0e+5</value>
        </parameter>
        <parameter>
            <name>p_bottom</name>
            <type>Constant</type>
            <value>1.0e+6</value>
        </parameter>
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
        <parameter>
            <name>b0</name>
            <type>Constant</type>
            <values> -0.5000001748438417 0.0 0.8660253028382761</values>
        </parameter>
        <parameter>
            <name>b1</name>
            <type>Constant</type>
            <values>0.0 -1.0 0.0</values>
        </parameter>
        <parameter>
            <name>b2</name>
            <!-- off plane direction -->
            <type>Constant</type>
            <values>0.8660253028382761 0.0 0.5000001748438417</values>
        </parameter>
167
    </parameters>
168
169
170
171
172
    <local_coordinate_system>
        <basis_vector_0>b0</basis_vector_0>
        <basis_vector_1>b1</basis_vector_1>
        <basis_vector_2>b2</basis_vector_2>
    </local_coordinate_system>
173
174
175
176
177
178
179
180
181
    <process_variables>
        <process_variable>
            <name>pressure</name>
            <components>1</components>
            <order>1</order>
            <initial_condition>p0</initial_condition>
            <boundary_conditions>
                <boundary_condition>
                    <type>Dirichlet</type>
182
                    <mesh>geometry_top</mesh>
183
184
185
186
                    <parameter>p_top</parameter>
                </boundary_condition>
                <boundary_condition>
                    <type>Dirichlet</type>
187
                    <mesh>geometry_bottom</mesh>
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
                    <parameter>p_bottom</parameter>
                </boundary_condition>
            </boundary_conditions>
        </process_variable>
    </process_variables>
    <nonlinear_solvers>
        <nonlinear_solver>
            <name>basic_picard</name>
            <type>Picard</type>
            <max_iter>10</max_iter>
            <linear_solver>general_linear_solver</linear_solver>
        </nonlinear_solver>
    </nonlinear_solvers>
    <linear_solvers>
        <linear_solver>
            <name>general_linear_solver</name>
            <lis>-i cg -p jacobi -tol 1e-20 -maxiter 10000</lis>
            <eigen>
                <solver_type>CG</solver_type>
                <precon_type>DIAGONAL</precon_type>
                <max_iteration_step>10000</max_iteration_step>
                <error_tolerance>1e-20</error_tolerance>
            </eigen>
            <petsc>
                <prefix>lf</prefix>
                <parameters>-lf_ksp_type cg -lf_pc_type bjacobi -lf_ksp_rtol 1e-16 -lf_ksp_max_it 10000</parameters>
            </petsc>
        </linear_solver>
    </linear_solvers>
</OpenGeoSysProject>