TH2M Test case

Mailing List Archives ogs-users
1

Hi All,

I am trying to run a TH2M case. But, coupling is not working. I checked TH2M example problem, and used similar *.num file. After Process 1: HEAT_TRANSPORT

"Inner coupling loop 1/3 complete

Max coupling error (relative error tolerance):0"

After Process 2: DEFORMATION

Starting loading step 1/1. Load factor :1

Assembling equation system…

Then terminated. Any idea? Most likely, problem originated from the *.num file.

Any cooperation will be appreciable. In addition, is there any way to take complete error message and log in a text file?

Regards,

Mohammad

···

*.bc

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

POLYLINE TOP

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_Y1

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_X1

$GEO_TYPE

POLYLINE left

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_X1

$GEO_TYPE

POLYLINE right

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

HEAT_TRANSPORT

$PRIMARY_VARIABLE

TEMPERATURE1

$GEO_TYPE

POLYLINE TOP_C

$DIS_TYPE

CONSTANT 70

#STOP

=============================

*.ic

#INITIAL_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 0

#INITIAL_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 0

#INITIAL_CONDITION

$PCS_TYPE

HEAT_TRANSPORT

$PRIMARY_VARIABLE

TEMPERATURE1

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 30.0

#STOP

*.mfp

#FLUID_PROPERTIES

$FLUID_TYPE

LIQUID

$PCS_TYPE

PRESSURE1

$DENSITY

1 1.000000e+003

$VISCOSITY

1 1.000000e-003

$SPECIFIC_HEAT_CAPACITY

1 4.280000e+003

$HEAT_CONDUCTIVITY

1 6.000000e-001

#FLUID_PROPERTIES

$FLUID_TYPE

GAS

$PCS_TYPE

PRESSURE2

$DENSITY

7

$VISCOSITY

1 1.8e-005

$SPECIFIC_HEAT_CAPACITY

1 1.01e+3

$HEAT_CONDUCTIVITY

1 0.026

$PHASE_DIFFUSION

1 2.13e-6

#STOP

*.mmp

#MEDIUM_PROPERTIES

$GEOMETRY_DIMENSION

2

$GEOMETRY_AREA

1.000000e+000

$POROSITY

1 2.975000e-001

$TORTUOSITY

1 0.8

$DIFFUSION

273

$PERMEABILITY_TENSOR

ISOTROPIC 4.500000e-013

$PERMEABILITY_SATURATION

0 2

0 2

$CAPILLARY_PRESSURE

0 1

#STOP

*.msp

#SOLID_PROPERTIES

$DENSITY

1 -2000

$THERMAL

EXPANSION: 1.0e-5

CAPACITY:

1 1.091000e+003

CONDUCTIVITY:

1 0.4200000e+000

$ELASTICITY

POISSION 0.4

YOUNGS_MODULUS:

1 1.300000e+006

$BIOT_CONSTANT

1

#STOP

*.num

GeoSys-NUM: Numerical Parameter ----------------------------------------

$OVERALL_COUPLING

;min_iter – max_iter

1 1

#NUMERICS

$PCS_TYPE

MULTI_PHASE_FLOW

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 2 1.e-10 2000 1.0 100 4

$ELE_MASS_LUMPING

0

$ELE_GAUSS_POINTS

2

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

PICARD 1.0e-4 30 0.0

$COUPLING_CONTROL

;error method – tolerances

LMAX 1.e-3 1.e-3

$COUPLED_PROCESS

;process or variable – min_iter – max_iter

HEAT_TRANSPORT 1 3

#NUMERICS

$PCS_TYPE

HEAT_TRANSPORT

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 1 1.e-012 1000 0.5 100 4

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

PICARD 1.0e-6 100 0.0

$COUPLING_CONTROL

;error method – tolerances

LMAX 1.e-6

#NUMERICS

$PCS_TYPE

DEFORMATION

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

NEWTON 1e-2 1e-10 100 0.0

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 0 1.e-011 2000 1.0 100 4

$COUPLING_CONTROL

;error method – tolerances

LMAX 1.e-3

#STOP

2

Hi,

> In addition, is there any way to take complete error message and log in a text file?

you can do it as

ogs &> log.txt

···

On 06/22/2017 03:31 AM, Mohammad Islam wrote:

Hi All,
I am trying to run a TH2M case. But, coupling is not working. I checked TH2M example problem, and used similar *.num file. After Process 1: HEAT_TRANSPORT

"Inner coupling loop 1/3 complete
Max coupling error (relative error tolerance):0"

After Process 2: DEFORMATION
Starting loading step 1/1. Load factor :1
Assembling equation system...

Then terminated. Any idea? Most likely, problem originated from the *.num file.

Any cooperation will be appreciable. In addition, is there any way to take complete error message and log in a text file?

Regards,
Mohammad

*******************
*.bc
  ***************

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

POLYLINE TOP

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_Y1

$GEO_TYPE

POLYLINE BOTTOM

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_X1

$GEO_TYPE

POLYLINE left

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

DISPLACEMENT_X1

$GEO_TYPE

POLYLINE right

$DIS_TYPE

CONSTANT 0

#BOUNDARY_CONDITION

$PCS_TYPE

HEAT_TRANSPORT

$PRIMARY_VARIABLE

TEMPERATURE1

$GEO_TYPE

POLYLINE TOP_C

$DIS_TYPE

CONSTANT 70

#STOP

=============================

*.ic

*******************************************

#INITIAL_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 0

#INITIAL_CONDITION

$PCS_TYPE

MULTI_PHASE_FLOW

$PRIMARY_VARIABLE

PRESSURE2

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 0

#INITIAL_CONDITION

$PCS_TYPE

HEAT_TRANSPORT

$PRIMARY_VARIABLE

TEMPERATURE1

$GEO_TYPE

DOMAIN

$DIS_TYPE

CONSTANT 30.0

#STOP

*************************

*.mfp

************************

#FLUID_PROPERTIES

$FLUID_TYPE

LIQUID

$PCS_TYPE

PRESSURE1

$DENSITY

1 1.000000e+003

$VISCOSITY

1 1.000000e-003

$SPECIFIC_HEAT_CAPACITY

1 4.280000e+003

$HEAT_CONDUCTIVITY

1 6.000000e-001

#FLUID_PROPERTIES

$FLUID_TYPE

GAS

$PCS_TYPE

PRESSURE2

$DENSITY

7

$VISCOSITY

1 1.8e-005

$SPECIFIC_HEAT_CAPACITY

1 1.01e+3

$HEAT_CONDUCTIVITY

1 0.026

$PHASE_DIFFUSION

1 2.13e-6

#STOP

***************************************

*.mmp

*************************************

#MEDIUM_PROPERTIES

$GEOMETRY_DIMENSION

2

$GEOMETRY_AREA

1.000000e+000

$POROSITY

1 2.975000e-001

$TORTUOSITY

1 0.8

$DIFFUSION

273

$PERMEABILITY_TENSOR

ISOTROPIC 4.500000e-013

$PERMEABILITY_SATURATION

0 2

0 2

$CAPILLARY_PRESSURE

0 1

#STOP

**************************************

*.msp

*************************************

#SOLID_PROPERTIES

$DENSITY

1 -2000

$THERMAL

EXPANSION: 1.0e-5

CAPACITY:

1 1.091000e+003

CONDUCTIVITY:

1 0.4200000e+000

$ELASTICITY

POISSION 0.4

YOUNGS_MODULUS:

1 1.300000e+006

$BIOT_CONSTANT

1

#STOP

********************************

*.num

*****************************

GeoSys-NUM: Numerical Parameter ----------------------------------------

$OVERALL_COUPLING

;min_iter -- max_iter

1 1

#NUMERICS

$PCS_TYPE

MULTI_PHASE_FLOW

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 2 1.e-10 2000 1.0 100 4

$ELE_MASS_LUMPING

0

$ELE_GAUSS_POINTS

2

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

PICARD 1.0e-4 30 0.0

$COUPLING_CONTROL

;error method -- tolerances

LMAX 1.e-3 1.e-3

$COUPLED_PROCESS

;process or variable -- min_iter -- max_iter

HEAT_TRANSPORT 1 3

#NUMERICS

$PCS_TYPE

HEAT_TRANSPORT

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 1 1.e-012 1000 0.5 100 4

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

PICARD 1.0e-6 100 0.0

$COUPLING_CONTROL

;error method -- tolerances

LMAX 1.e-6

#NUMERICS

$PCS_TYPE

DEFORMATION

$NON_LINEAR_SOLVER

; method error_tolerance max_iterations relaxation

NEWTON 1e-2 1e-10 100 0.0

$LINEAR_SOLVER

; method error_tolerance max_iterations theta precond storage

2 0 1.e-011 2000 1.0 100 4

$COUPLING_CONTROL

;error method -- tolerances

LMAX 1.e-3

#STOP

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ogs-users+unsubscribe@googlegroups.com <mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout\.

3

Dear Nori,

Thank you.

Regards,

Mohammad

···
4

Hi All,

Still I am having trouble with TH2M coupling. H2M case individually working fine, but when I added Thermal part, OGS stopped working with error messages as follows

SolidProperties: HeatConductivity Tensor (): conductivity mode is not supported

Error: in CMediumProperties:: Heat Capacity: no valid material model

After Process 2: DEFORMATION

OGS stopped

Any idea?

For convenience I attached H2M and TH2M case files. Source of the file: Benchmark Problem>> TH2M/H2M_TEM. In addition, it is to note that under JFNK folder, example problem is not running due to JFNK solver (it is no longer supported). Furthermore, under Newton folder, process is DEFORMATION_H2. Does it represent TH2M case? I attempted to modify for TH2M case, but it is not working. Finally, in same folder for Benchmark problem, th2m_quad.** , after initial state there is no evolution with respect to time. I checked it and found capillary pressure, gas pressure, saturation are same, which is not rational. Any co-operation will be appreciable.

In addition, I am also reading attached paper by Kolditz et al. (2012) " A systematic benchmarking approach for geologic CO2 injection and storage", Env. Earth Sc. 67: 613-632 . In Table 1: Non-isothermal two-phase flow consolidation TH2M case presented. I am wondering is there any way to get the benchmark example?

Thanks in advance.

Regards,

Mohammad

w_exp.bc (1.06 KB)

w_exp.gli (840 Bytes)

w_exp.ic (709 Bytes)

w_exp.mfp (285 Bytes)

w_exp.mmp (375 Bytes)

w_exp.msh (1.94 KB)

w_exp.num (764 Bytes)

w_exp.out (2.63 KB)

w_exp.pcs (410 Bytes)

w_exp.rfd (5.87 KB)

w_exp.st (402 Bytes)

w_exp.tim (264 Bytes)

th2m-w_exp.bc (1.2 KB)

th2m-w_exp.gli (996 Bytes)

th2m-w_exp.ic (842 Bytes)

th2m-w_exp.mfp (481 Bytes)

th2m-w_exp.mmp (417 Bytes)

th2m-w_exp.msh (1.96 KB)

th2m-w_exp.num (1.5 KB)

th2m-w_exp.out (3.75 KB)

th2m-w_exp.pcs (452 Bytes)

w_exp.rfd (5.87 KB)

w_exp.st (402 Bytes)

th2m-w_exp.tim (631 Bytes)

10.1007_s12665-012-1656-5.pdf (1.53 MB)

···

th2m-w_exp.msp

GeoSys-MAT-SP: Solid Properties ------------------------------------

#SOLID_PROPERTIES

$DENSITY

1 2010.0

$THERMAL

EXPANSION 1.0e-5

CAPACITY:

3

CONDUCTIVITY:

3

$ELASTICITY

POISSION 0.35

YOUNGS_MODULUS:

1 3.17e+8

$SWELLING_PRESSURE_TYPE

2 4.76e6

#STOP

5

Hi, can you also send us MSP files? I can't find them in the attachments.

n

···

On 06/28/2017 06:03 AM, Mohammad Islam wrote:

Hi All,
Still I am having trouble with TH2M coupling. H2M case individually working fine, but when I added Thermal part, OGS stopped working with error messages as follows

SolidProperties: HeatConductivity Tensor (): conductivity mode is not supported
Error: in CMediumProperties:: Heat Capacity: no valid material model

After Process 2: DEFORMATION
  OGS stopped

Any idea?

For convenience I attached H2M and TH2M case files. Source of the file: Benchmark Problem>> TH2M/H2M_TEM. In addition, it is to note that under JFNK folder, example problem is not running due to JFNK solver (it is no longer supported). Furthermore, under Newton folder, process is DEFORMATION_H2. Does it represent TH2M case? I attempted to modify for TH2M case, but it is not working. Finally, in same folder for Benchmark problem, th2m_quad.** , after initial state there is no evolution with respect to time. I checked it and found capillary pressure, gas pressure, saturation are same, which is not rational. Any co-operation will be appreciable.

In addition, I am also reading attached paper by Kolditz et al. (2012) " A systematic benchmarking approach for geologic CO2 injection and storage", Env. Earth Sc. 67: 613-632 . In Table 1: Non-isothermal two-phase flow consolidation TH2M case presented. I am wondering is there any way to get the benchmark example?

Thanks in advance.

Regards,
Mohammad

-----
th2m-w_exp.msp

GeoSys-MAT-SP: Solid Properties ------------------------------------

#SOLID_PROPERTIES

$DENSITY

1 2010.0

$THERMAL

EXPANSION 1.0e-5

CAPACITY:

3

CONDUCTIVITY:

3

$ELASTICITY

POISSION 0.35

YOUNGS_MODULUS:

1 3.17e+8

$SWELLING_PRESSURE_TYPE

2 4.76e6

#STOP

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ogs-users+unsubscribe@googlegroups.com <mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout\.

6

sorry, the MSP part was included in the message body. In the file, you need to add a function type before giving the value like

CAPACITY:
1 800.0

CONDUCTIVITY:
1 3.0

The first number 1 means you use a constant value for the solid phase conductivity. Depending on a parameter type, you may be able to choose other function type such as f(t) or f(T).

see for example, https://svn.ufz.de/ogs/browser/trunk/benchmarks/TH2M/th2m_quad.msp

nori

···

On 06/28/2017 06:12 AM, Norihiro Watanabe wrote:

Hi, can you also send us MSP files? I can't find them in the attachments.

n

On 06/28/2017 06:03 AM, Mohammad Islam wrote:

Hi All,
Still I am having trouble with TH2M coupling. H2M case individually working fine, but when I added Thermal part, OGS stopped working with error messages as follows

SolidProperties: HeatConductivity Tensor (): conductivity mode is not supported
Error: in CMediumProperties:: Heat Capacity: no valid material model

After Process 2: DEFORMATION
  OGS stopped

Any idea?

For convenience I attached H2M and TH2M case files. Source of the file: Benchmark Problem>> TH2M/H2M_TEM. In addition, it is to note that under JFNK folder, example problem is not running due to JFNK solver (it is no longer supported). Furthermore, under Newton folder, process is DEFORMATION_H2. Does it represent TH2M case? I attempted to modify for TH2M case, but it is not working. Finally, in same folder for Benchmark problem, th2m_quad.** , after initial state there is no evolution with respect to time. I checked it and found capillary pressure, gas pressure, saturation are same, which is not rational. Any co-operation will be appreciable.

In addition, I am also reading attached paper by Kolditz et al. (2012) " A systematic benchmarking approach for geologic CO2 injection and storage", Env. Earth Sc. 67: 613-632 . In Table 1: Non-isothermal two-phase flow consolidation TH2M case presented. I am wondering is there any way to get the benchmark example?

Thanks in advance.

Regards,
Mohammad

-----
th2m-w_exp.msp

GeoSys-MAT-SP: Solid Properties ------------------------------------

#SOLID_PROPERTIES

$DENSITY

1 2010.0

$THERMAL

EXPANSION 1.0e-5

CAPACITY:

3

CONDUCTIVITY:

3

$ELASTICITY

POISSION 0.35

YOUNGS_MODULUS:

1 3.17e+8

$SWELLING_PRESSURE_TYPE

2 4.76e6

#STOP

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ogs-users+unsubscribe@googlegroups.com <mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout\.

7

Dear Nori,
Thank you. Error message resolved. But, still having problem.

I prepared two sets of test: HM (Richards Flow) & H2M (Multiphase flow), then added Thermal part. HM, H2M, THM cases are running, but in TH2M case after initial state in results only .#IND and #QNAN (not a number?). On the other hand, in THM case results are identical to the HM case. I attached comparison.xls file (THM with & without heat injection and HM case).

HM.MSP

GeoSys-MAT-SP: Solid Properties ------------------------------------

#SOLID_PROPERTIES

$DENSITY

1 -2000

$ELASTICITY

POISSION 0.4

YOUNGS_MODULUS:

1 1.300000e+006

#STOP

THM.MSP

#SOLID_PROPERTIES

$DENSITY

1 -2000

$THERMAL

EXPANSION: 1.0e-5

CAPACITY:

1 1.091000e+003

CONDUCTIVITY:

1 0.4200000e+000

$ELASTICITY

POISSION 0.4

YOUNGS_MODULUS:

1 1.300000e+006

#STOP

Thanks in advance.

Regards,

Mohammad

comparison.xlsx (42.2 KB)

HM-Un-Lia-Q.bc (896 Bytes)

HM-Un-Lia-Q.gli (1.05 KB)

HM-Un-Lia-Q.ic (258 Bytes)

HM-Un-Lia-Q.mfp (405 Bytes)

HM-Un-Lia-Q.mmp (538 Bytes)

HM-Un-Lia-Q.msh (2.33 KB)

HM-Un-Lia-Q.num (1.59 KB)

HM-Un-Lia-Q.out (840 Bytes)

HM-Un-Lia-Q.pcs (171 Bytes)

HM-Un-Lia-Q.rfd (1.19 KB)

HM-Un-Lia-Q.st (1 KB)

HM-Un-Lia-Q.tim (763 Bytes)

THM-Un-Lia-Q.bc (1.06 KB)

HM-Un-Lia-Q.gli (1.05 KB)

THM-Un-Lia-Q.ic (435 Bytes)

THM-Un-Lia-Q.mfp (489 Bytes)

THM-Un-Lia-Q.mmp (580 Bytes)

HM-Un-Lia-Q.msh (2.33 KB)

THM-Un-Lia-Q.num (1.91 KB)

THM-Un-Lia-Q.out (885 Bytes)

THM-Un-Lia-Q.pcs (236 Bytes)

HM-Un-Lia-Q.rfd (1.19 KB)

THM-Un-Lia-Q.st (1 KB)

THM-Un-Lia-Q.tim (1.03 KB)