Skip to content
Snippets Groups Projects
Commit 56510ed4 authored by Christoph Lehmann's avatar Christoph Lehmann
Browse files

[T/CT] Added advection diffusion reaction test case

parent f3b0b092
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 577 additions and 0 deletions
ProcessLib/ComponentTransport/Tests.cmake
+ 1
− 0
View file @ 56510ed4
......@@ -930,3 +930,4 @@ endif()
OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/only_grad_c.prj RUNTIME 1)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p.prj RUNTIME 3)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/only_grad_p.prj RUNTIME 1)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p_and_r.prj RUNTIME 3)
Tests/Data/Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p_and_r.prj 0 → 100644
+ 338
− 0
View file @ 56510ed4
<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<meshes>
<mesh>grad_c_and_grad_p_and_r.vtu</mesh>
<mesh>grad_c_and_grad_p_and_r_left.vtu</mesh>
<mesh>grad_c_and_grad_p_and_r_right.vtu</mesh>
</meshes>
<processes>
<process>
<name>HC</name>
<type>ComponentTransport</type>
<integration_order>4</integration_order>
<process_variables>
<concentration>C</concentration>
<pressure>pressure</pressure>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity"/>
<secondary_variable internal_name="CFlux" output_name="CFlux"/>
</secondary_variables>
<specific_body_force>0 0 0</specific_body_force>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>C</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>1.e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
</components>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1.e-9</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.15</value>
</property>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
</properties>
</medium>
<medium id="1">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>C</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>1.e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
</components>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1.e-9</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.65</value>
</property>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
</properties>
</medium>
</media>
<parameters>
<parameter>
<name>decay</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>C_right</name>
<type>Constant</type>
<value>1</value>
</parameter>
<parameter>
<name>C_left</name>
<type>Constant</type>
<value>0.5</value>
</parameter>
<parameter>
<name>C_ini</name>
<type>MeshNode</type>
<field_name>C_ini</field_name>
</parameter>
<parameter>
<name>p_left</name>
<type>Constant</type>
<value>0.3</value>
</parameter>
<parameter>
<name>p_right</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p_ini</name>
<type>MeshNode</type>
<field_name>p_ini</field_name>
</parameter>
<parameter>
<name>r</name>
<type>Constant</type>
<value>-5e-10</value>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>p_ini</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>grad_c_and_grad_p_and_r_left</mesh>
<type>Dirichlet</type>
<parameter>p_left</parameter>
</boundary_condition>
<boundary_condition>
<mesh>grad_c_and_grad_p_and_r_right</mesh>
<type>Dirichlet</type>
<parameter>p_right</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>C</name>
<components>1</components>
<order>1</order>
<initial_condition>C_ini</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>grad_c_and_grad_p_and_r_right</mesh>
<type>Dirichlet</type>
<parameter>C_right</parameter>
</boundary_condition>
<boundary_condition>
<mesh>grad_c_and_grad_p_and_r_left</mesh>
<type>Dirichlet</type>
<parameter>C_left</parameter>
</boundary_condition>
</boundary_conditions>
<source_terms>
<source_term>
<mesh>grad_c_and_grad_p_and_r</mesh>
<type>Volumetric</type>
<parameter>r</parameter>
</source_term>
</source_terms>
</process_variable>
</process_variables>
<time_loop>
<processes>
<process ref="HC">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>INFINITY_N</norm_type>
<!--<reltols>1e-12 1e-10</reltols>-->
<abstols>1e-15 1e-15</abstols>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 4e8 </t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>1e8</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>grad_c_and_grad_p_and_r</prefix>
<timesteps>
<pair>
<repeat>1000</repeat>
<each_steps>1</each_steps>
</pair>
</timesteps>
<variables>
</variables>
<output_extrapolation_residuals>true</output_extrapolation_residuals>
</output>
</time_loop>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>50</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>SparseLU</solver_type>
<scaling>true</scaling>
</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>
<test_definition>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>pressure</field>
<absolute_tolerance>4.5e-9</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>C</field>
<absolute_tolerance>2.8e-8</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>CFlux</field>
<absolute_tolerance>3.4e-13</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>CFlux_residual</field>
<absolute_tolerance>3.1e-13</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>darcy_velocity</field>
<absolute_tolerance>5.6e-15</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<file>grad_c_and_grad_p_and_r_ts_4_t_400000000.000000.vtu</file>
<field>darcy_velocity_residual</field>
<absolute_tolerance>1.7e-15</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
</test_definition>
</OpenGeoSysProject>
Tests/Data/Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p_and_r_left.vtu 0 → 100644
+ 53
− 0
View file @ 56510ed4
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<Piece NumberOfPoints="4" NumberOfCells="1">
<PointData Scalars="C_ini">
<DataArray type="Float64" Name="C_ini" format="binary" RangeMin="0.5" RangeMax="0.5">
AQAAAACAAAAgAAAAEAAAAA==eJxjYACBB/YMOGgAPewEfQ==
</DataArray>
<DataArray type="Float64" Name="p_ini" format="binary" RangeMin="0.3" RangeMax="0.3">
AQAAAACAAAAgAAAAEAAAAA==eJwzNgaBy/bGOGgAjsAJEQ==
</DataArray>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="3003">
AQAAAACAAAAgAAAAFgAAAA==eJxjYICAS+wQ+iUzhN7NDaEBKVMCjA==
</DataArray>
<DataArray type="Int64" IdType="1" Name="vtkOriginalPointIds" format="binary" RangeMin="0" RangeMax="6">
AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAZijNAqXZoDQAANgADg==
</DataArray>
</PointData>
<CellData Scalars="MaterialIDs">
<DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
AQAAAACAAAAEAAAADAAAAA==eJxjYGBgAAAABAAB
</DataArray>
<DataArray type="Int64" IdType="1" Name="vtkOriginalCellIds" format="binary" RangeMin="2" RangeMax="2">
AQAAAACAAAAIAAAACwAAAA==eJxjYoAAAAAYAAM=
</DataArray>
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0.0011313708213176444">
AQAAAACAAABgAAAAFgAAAA==eJxjYCAaPHhk5mWPh48uDqYBAjsKBQ==
<InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
<Value index="0">
0
</Value>
<Value index="1">
0.0011313708213
</Value>
</InformationKey>
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="3">
AQAAAACAAAAgAAAAEwAAAA==eJxjYIAARijNDKWZoDQAAHgABw==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="4">
AQAAAACAAAAIAAAACwAAAA==eJxjYYAAAAAoAAU=
</DataArray>
<DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
AQAAAACAAAABAAAACQAAAA==eJzjBAAACgAK
</DataArray>
</Cells>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p_and_r_right.vtu 0 → 100644
+ 53
− 0
View file @ 56510ed4
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<Piece NumberOfPoints="4" NumberOfCells="1">
<PointData Scalars="C_ini">
<DataArray type="Float64" Name="C_ini" format="binary" RangeMin="1.0000000102426312" RangeMax="1.0000000102426312">
AQAAAACAAAAgAAAAEwAAAA==eJxzu7efiYHhg70bDhoAz3QMUQ==
</DataArray>
<DataArray type="Float64" Name="p_ini" format="binary" RangeMin="-4.470348341500951e-9" RangeMax="-4.470348341500951e-9">
AQAAAACAAAAgAAAAEwAAAA==eJz7//+/obGx8b7/OGgAOsMSFQ==
</DataArray>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="1000" RangeMax="4003">
AQAAAACAAAAgAAAAGAAAAA==eJx7wcwABru4IfRFdgi9mB9CAwBC5AM7
</DataArray>
<DataArray type="Int64" IdType="1" Name="vtkOriginalPointIds" format="binary" RangeMin="3001" RangeMax="3004">
AQAAAACAAAAgAAAAFQAAAA==eJzbyc0ABrug9G4ovQdKAwA9hAMX
</DataArray>
</PointData>
<CellData Scalars="MaterialIDs">
<DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="1" RangeMax="1">
AQAAAACAAAAEAAAADAAAAA==eJxjZGBgAAAACAAC
</DataArray>
<DataArray type="Int64" IdType="1" Name="vtkOriginalCellIds" format="binary" RangeMin="2001" RangeMax="2001">
AQAAAACAAAAIAAAADAAAAA==eJy7yM4ABgAGwQDZ
</DataArray>
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0.800000011920929" RangeMax="0.8000008119204767">
AQAAAACAAABgAAAAIwAAAA==eJxjYGBYMHPmS3sGTIAu/uCRmZc9kjiMj64eJg6mAZL5Fe0=
<InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
<Value index="0">
0.80000001192
</Value>
<Value index="1">
0.80000081192
</Value>
</InformationKey>
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="3">
AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAJijNDKUZoTQAAIgABw==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="4">
AQAAAACAAAAIAAAACwAAAA==eJxjYYAAAAAoAAU=
</DataArray>
<DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
AQAAAACAAAABAAAACQAAAA==eJzjBAAACgAK
</DataArray>
</Cells>
</Piece>
</UnstructuredGrid>
</VTKFile>
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