Skip to content
Snippets Groups Projects
Commit f077de26 authored by renchao.lu's avatar renchao.lu
Browse files

[T] Benchmark of tracer test.

parent c6d7d8c8
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 347 additions and 0 deletions
ProcessLib/ComponentTransport/Tests.cmake
+ 21
− 0
View file @ f077de26
......@@ -451,6 +451,27 @@ AddTest(
cube_1x1x1_hex_1e3_right_pcs_1_ts_1_t_1.000000.vtu cube_1x1x1_hex_1e3_right_pcs_1_ts_1_t_1.000000.vtu Si Si 1e-10 1e-16
)
AddTest(
NAME 2D_StaggeredScheme_ComponentTransport_TracerSimulation
PATH Parabolic/ComponentTransport/TracerSimulation
EXECUTABLE ogs
EXECUTABLE_ARGS TracerSimulation.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
TracerSimulation_pcs_1_ts_20_t_20000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_20_t_20000.000000.vtu Cs Cs 4e-7 1e-10
TracerSimulation_pcs_1_ts_40_t_40000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_40_t_40000.000000.vtu Cs Cs 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_60_t_60000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_60_t_60000.000000.vtu Cs Cs 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_80_t_80000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_80_t_80000.000000.vtu Cs Cs 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_100_t_100000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_100_t_100000.000000.vtu Cs Cs 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_20_t_20000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_20_t_20000.000000.vtu pressure pressure 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_40_t_40000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_40_t_40000.000000.vtu pressure pressure 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_60_t_60000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_60_t_60000.000000.vtu pressure pressure 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_80_t_80000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_80_t_80000.000000.vtu pressure pressure 1.3e-7 1e-10
TracerSimulation_pcs_1_ts_100_t_100000.000000_expected_ogs5.vtu TracerSimulation_pcs_1_ts_100_t_100000.000000.vtu pressure pressure 1.3e-7 1e-10
)
AddTest(
NAME LARGE_2D_ComponentTransport_Goswami
PATH Parabolic/ComponentTransport/goswami
......
Tests/Data/Parabolic/ComponentTransport/TracerSimulation/TracerSimulation.gml 0 → 100644 LFS
+ 3
− 0
View file @ f077de26
Source diff could not be displayed: it is stored in LFS. Options to address this: view the blob.
Tests/Data/Parabolic/ComponentTransport/TracerSimulation/TracerSimulation.prj 0 → 100644
+ 290
− 0
View file @ f077de26
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh>TracerSimulation.vtu</mesh>
<geometry>TracerSimulation.gml</geometry>
<processes>
<process>
<name>hc</name>
<type>ComponentTransport</type>
<integration_order>2</integration_order>
<coupling_scheme>staggered</coupling_scheme>
<process_variables>
<concentration>Cs</concentration>
<pressure>pressure</pressure>
</process_variables>
<fluid>
<density>
<type>Constant</type>
<value>1e3</value>
</density>
<viscosity>
<type>Constant</type>
<value>1e-3</value>
</viscosity>
</fluid>
<porous_medium>
<porous_medium id="0">
<permeability>
<type>Constant</type>
<permeability_tensor_entries>kappa</permeability_tensor_entries>
</permeability>
<porosity>
<type>Constant</type>
<porosity_parameter>porosity</porosity_parameter>
</porosity>
<storage>
<type>Constant</type>
<value>0.0</value>
</storage>
</porous_medium>
</porous_medium>
<fluid_reference_density>rho_fluid</fluid_reference_density>
<retardation_factor>retardation</retardation_factor>
<decay_rate>decay</decay_rate>
<specific_body_force>0 0</specific_body_force>
<secondary_variables>
<secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
</secondary_variables>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>Cs</name>
<properties>
<property>
<name>molecular_diffusion</name>
<type>Constant</type>
<value>1e-9</value>
</property>
</properties>
</component>
</components>
</phase>
</phases>
<properties>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.2</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.2</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<global_process_coupling>
<max_iter>6</max_iter>
<convergence_criteria>
<!-- convergence criterion for the first process (p) -->
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<!-- convergence criterion for the second process (Cs) -->
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
</convergence_criteria>
</global_process_coupling>
<processes>
<!-- convergence criterion for hydraulic equation -->
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>100000</t_end>
<timesteps>
<pair>
<repeat>100</repeat>
<delta_t>1000</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<!-- convergence criterion for component transport equation (Cs) -->
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>100000</t_end>
<timesteps>
<pair>
<repeat>100</repeat>
<delta_t>1000</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>TracerSimulation</prefix>
<timesteps>
<pair>
<repeat>5</repeat>
<each_steps>20</each_steps>
</pair>
</timesteps>
<variables>
<variable>Cs</variable>
<variable>pressure</variable>
<variable>darcy_velocity</variable>
</variables>
</output>
</time_loop>
<parameters>
<parameter>
<name>kappa</name>
<type>Constant</type>
<values>1.019368e-11 0 0 1.019368e-11</values>
</parameter>
<parameter>
<name>porosity</name>
<type>Constant</type>
<value>0.2</value>
</parameter>
<parameter>
<name>rho_fluid</name>
<type>Constant</type>
<value>1000</value>
</parameter>
<parameter>
<name>retardation</name>
<type>Constant</type>
<value>1</value>
</parameter>
<parameter>
<name>decay</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p_Neumann_left</name>
<type>Constant</type>
<value>3e-2</value>
</parameter>
<parameter>
<name>p_Dirichlet_right</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>c0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>c_spatial</name>
<type>Constant</type>
<value>1</value>
</parameter>
<parameter>
<name>c_Dirichlet_left</name>
<type>CurveScaled</type>
<curve>c_temporal</curve>
<parameter>c_spatial</parameter>
</parameter>
</parameters>
<curves>
<curve>
<name>c_temporal</name>
<coords>0 10000 10000.001 100000</coords>
<values>1e-5 1e-5 0 0</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>p0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>geometry</geometrical_set>
<geometry>left</geometry>
<type>Neumann</type>
<parameter>p_Neumann_left</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>geometry</geometrical_set>
<geometry>right</geometry>
<type>Dirichlet</type>
<parameter>p_Dirichlet_right</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Cs</name>
<components>1</components>
<order>1</order>
<initial_condition>c0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>geometry</geometrical_set>
<geometry>left</geometry>
<type>Dirichlet</type>
<parameter>c_Dirichlet_left</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 bicgstab -p ilut -tol 1e-8 -maxiter 20000</lis>
<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>
<petsc>
<prefix>hc</prefix>
<parameters>-hc_ksp_type bcgs -hc_pc_type bjacobi -hc_ksp_rtol 1e-8 -hc_ksp_max_it 20000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
Tests/Data/Parabolic/ComponentTransport/TracerSimulation/TracerSimulation.vtu 0 → 100644 LFS
+ 3
− 0
View file @ f077de26
File added
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment