Cannot run "recipitation/dissolution equilibrium reactions in a saturated column" benchmark

Usability
1

Hello,
When I try to run the benchmark (calcite.prj), I am getting an error shown in figure below.

calciteprj
calciteprj1848×923 176 KB

emre@emre-VirtualBox:~/Desktop/OGS1/EquilibriumPhase$ ~/Desktop/OGS1/build-release/bin/ogs -l debug calcite.prj

Am I missing something?

info: This is OpenGeoSys-6 version 6.3.3-573-g4e0cd0af9.
info: OGS started on 2021-03-11 11:59:19+0300.
debug: Project configuration from file ‘calcite.prj’ read.
debug: Reading multiple meshes.
debug: Reading mesh file ‘calcite.vtu’.
debug: Reading mesh file ‘calcite_upstream.vtu’.
debug: Reading mesh file ‘calcite_downstream.vtu’.
debug: Reading mesh file ‘calcite_ReactiveDomain.vtu’.
debug: Reading parameters:
info: ConstantParameter: kappa
debug: Using following values for the constant parameter:
debug: 1.157e-12
info: ConstantParameter: porosity
debug: Using value 0.32 for constant parameter.
info: ConstantParameter: decay
debug: Using value 0 for constant parameter.
info: ConstantParameter: p0
debug: Using value 1 for constant parameter.
info: ConstantParameter: p_upstream
debug: Using value 100000 for constant parameter.
info: ConstantParameter: p_downstream_Neumann
debug: Using value -0.00299769 for constant parameter.
info: ConstantParameter: c0_C(4)
debug: Using value 0.000123 for constant parameter.
info: ConstantParameter: c0_Ca
debug: Using value 0.000123 for constant parameter.
info: ConstantParameter: c0_Mg
debug: Using value 1e-12 for constant parameter.
info: ConstantParameter: c0_Cl
debug: Using value 1e-12 for constant parameter.
info: ConstantParameter: c0_H
debug: Using value 1.23027e-10 for constant parameter.
info: ConstantParameter: c_C(4)
debug: Using value 1e-10 for constant parameter.
info: ConstantParameter: c_Ca
debug: Using value 1e-10 for constant parameter.
info: ConstantParameter: c_Mg
debug: Using value 0.001 for constant parameter.
info: ConstantParameter: c_Cl
debug: Using value 0.002 for constant parameter.
info: ConstantParameter: c_H
debug: Using value 1e-07 for constant parameter.
debug: Parse process variables:
debug: Found parameter p0'. debug: Constructing process variable pressure debug: Found mesh 'calcite_upstream' with id 1. debug: Found mesh 'calcite_downstream' with id 2. info: No source terms for process variable 'pressure' found. debug: Found parameter c0_C(4)’.
debug: Constructing process variable C(4)
debug: Found mesh ‘calcite_upstream’ with id 1.
info: No source terms for process variable ‘C(4)’ found.
debug: Found parameter c0_Ca'. debug: Constructing process variable Ca debug: Found mesh 'calcite_upstream' with id 1. info: No source terms for process variable 'Ca' found. debug: Found parameter c0_Mg’.
debug: Constructing process variable Mg
debug: Found mesh ‘calcite_upstream’ with id 1.
info: No source terms for process variable ‘Mg’ found.
debug: Found parameter c0_Cl'. debug: Constructing process variable Cl debug: Found mesh 'calcite_upstream' with id 1. info: No source terms for process variable 'Cl' found. debug: Found parameter c0_H’.
debug: Constructing process variable H
debug: Found mesh ‘calcite_upstream’ with id 1.
info: No source terms for process variable ‘H’ found.
debug: Reading media:
debug: Create Constant property pore_diffusion.
critical: /home/emre/Desktop/OGS1/ogs/MaterialLib/MPL/PropertyType.h:319 convertStringToProperty()
warning: Property is not used: ‘property ‘pore_diffusion’ defined for medium’
error: The property name ‘pore_diffusion’ does not correspond to any known property
info: OGS terminated on 2021-03-11 11:59:19+0300.
error: OGS terminated with error.

2

Hi Emre,

The warning and error message indicate the termination may be raised by the setting of the material property “pore_diffusion”. I guess you are using the latest developing version. Could you please upload your project file of the calcite benchmark to facilitate me to have an insight into the error?

Best regards,
Renchao

3

Hi Renchao,

Here is the project file.

EquilibriumPhase.zip (105.8 KB)

Thank you for your help.

Emre

4

Hi Emre,

I have checked your uploaded project file for the calcite benchmark. It is up-to-date. I run this benchmark with the latest version of ogs6 on my local. It works well.

Therefore, I recommend you to update the source code. Alternatively, you can also replace the property “pore_diffusion” with “molecular_diffusion” throughout the project file. Your current source code shall work with this change.

Best regards,
Renchao

5

Hi Renchao,

Changing pore_diffusion with molecular_diffusion solved the problem.

Thank you.