Skip to content

Commit

Permalink
[HT/Ctest] Add a ctest for fracture flow with HT.
Browse files Browse the repository at this point in the history
  • Loading branch information
ChaofanChen authored and endJunction committed Jul 19, 2023
1 parent 6d6532b commit 559ef76
Show file tree
Hide file tree
Showing 7 changed files with 470 additions and 0 deletions.
16 changes: 16 additions & 0 deletions ProcessLib/HT/Tests.cmake
Original file line number Diff line number Diff line change
Expand Up @@ -513,3 +513,19 @@ AddTest(
HeatTransportInStationaryFlow_ts_50_t_50000_000000_2.vtu HeatTransportInStationaryFlow_ts_50_t_50000_000000_2.vtu temperature temperature 1.e-9 1.0e-8
HeatTransportInStationaryFlow_ts_50_t_50000_000000_2.vtu HeatTransportInStationaryFlow_ts_50_t_50000_000000_2.vtu pressure pressure 1.e-9 1.0e-8
)

AddTest(
NAME 2D_single_fracture_HT
PATH Parabolic/HT/LowerDimensionalFracture
EXECUTABLE ogs
EXECUTABLE_ARGS 2D_single_fracture_HT.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
RUNTIME 55
DIFF_DATA
2D_single_fracture_HT_ts_3_t_1800.000000.vtu 2D_single_fracture_HT_ts_3_t_1800.000000.vtu temperature temperature 1.e-9 1.0e-11
2D_single_fracture_HT_ts_3_t_1800.000000.vtu 2D_single_fracture_HT_ts_3_t_1800.000000.vtu pressure pressure 1.e-9 1.0e-11
2D_single_fracture_HT_ts_6_t_3600.000000.vtu 2D_single_fracture_HT_ts_6_t_3600.000000.vtu temperature temperature 1.e-9 1.0e-11
2D_single_fracture_HT_ts_6_t_3600.000000.vtu 2D_single_fracture_HT_ts_6_t_3600.000000.vtu pressure pressure 1.e-9 1.0e-11
)

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
@@ -0,0 +1,327 @@
<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<meshes>
<mesh>2D_single_fracture.vtu</mesh>
<mesh>fracture_outflow.vtu</mesh>
<mesh>fracture_injection.vtu</mesh>
</meshes>
<processes>
<process>
<name>HeatTransportLowerDimensionalFracture</name>
<type>HT</type>
<integration_order>2</integration_order>
<process_variables>
<temperature>temperature</temperature>
<pressure>pressure</pressure>
</process_variables>
<specific_body_force>0.0 0.0</specific_body_force>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity"/>
</secondary_variables>
<numerical_stabilization>
<type>IsotropicDiffusion</type>
<tuning_parameter>0.15</tuning_parameter>
<cutoff_velocity>0</cutoff_velocity>
</numerical_stabilization>
<aperture_size>
<parameter>fracture_thickness</parameter>
</aperture_size>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4200</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0.6</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1000</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>storage</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>2700</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>3.0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>1000</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>thermal_longitudinal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1e-100</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>EffectiveThermalConductivityPorosityMixing</type>
</property>
</properties>
</medium>
<medium id="1">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4200</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0.6</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1000</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>storage</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>2850</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>0.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>thermal_longitudinal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>8.3333e-10</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>EffectiveThermalConductivityPorosityMixing</type>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="HeatTransportLowerDimensionalFracture">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-12</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 3600 </t_end>
<timesteps>
<pair>
<repeat>1</repeat>
<delta_t>600</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>2D_single_fracture_HT</prefix>
<timesteps>
<pair>
<repeat> 2 </repeat>
<each_steps> 3 </each_steps>
</pair>
</timesteps>
<variables>
<variable>temperature</variable>
<variable>pressure</variable>
<variable>darcy_velocity</variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>353.0</value>
</parameter>
<parameter>
<name>P0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p_left</name>
<type>Constant</type>
<value>972222</value>
</parameter>
<parameter>
<name>p_right</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>T_left</name>
<type>Constant</type>
<value>303.</value>
</parameter>
<parameter>
<name>fracture_thickness</name>
<type>Group</type>
<group_id_property>MaterialIDs</group_id_property>
<index_values>
<index>0</index>
<value>1</value>
</index_values>
<index_values>
<index>1</index>
<value>1e-4</value>
</index_values>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>temperature</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>fracture_injection</mesh>
<type>Dirichlet</type>
<parameter>T_left</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>P0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>fracture_injection</mesh>
<type>Dirichlet</type>
<parameter>p_left</parameter>
</boundary_condition>
<boundary_condition>
<mesh>fracture_outflow</mesh>
<type>Dirichlet</type>
<parameter>p_right</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>500</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>ILUT</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-12</error_tolerance>
</eigen>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>

Large diffs are not rendered by default.

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
@@ -0,0 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
<UnstructuredGrid>
<Piece NumberOfPoints="1" NumberOfCells="1" >
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="4" RangeMax="4" offset="0" />
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="4" RangeMax="4" offset="24" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="2" RangeMax="2" offset="48" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="92" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="116" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="140" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_CAAAAAAAAAAEAAAAAAAAAA==CAAAAAAAAAAEAAAAAAAAAA==GAAAAAAAAAAAAAAAAAAAAAAAAAAAAABAAAAAAAAAAAA=CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAABAAAAAAAAAA==AQAAAAAAAAAB
</AppendedData>
</VTKFile>
Original file line number Diff line number Diff line change
@@ -0,0 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
<UnstructuredGrid>
<Piece NumberOfPoints="1" NumberOfCells="1" >
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="5" RangeMax="5" offset="0" />
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="5" RangeMax="5" offset="24" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="7.2801098893" RangeMax="7.2801098893" offset="48" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="92" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="116" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="140" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_CAAAAAAAAAAFAAAAAAAAAA==CAAAAAAAAAAFAAAAAAAAAA==GAAAAAAAAAAAAAAAAAAcQAAAAAAAAABAAAAAAAAAAAA=CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAABAAAAAAAAAA==AQAAAAAAAAAB
</AppendedData>
</VTKFile>

0 comments on commit 559ef76

Please sign in to comment.