# Initial conditions-3D HTM model

Dear Users and Developers,

I am running a 3 model for HTM processes. I have completed the HT models and seem to be working well; however, when I implemented the third effect (M); I am not getting any solution at all. I have assigned the initial conditions as they are showed in the benchmarks, for example:

GINA - Initial Condition

#INITIAL_CONDITION

\$PCS_TYPE

HEAT_TRANSPORT

\$PRIMARY_VARIABLE

TEMPERATURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

CONSTANT 453.15

#INITIAL_CONDITION

\$PCS_TYPE

LIQUID_FLOW

\$PRIMARY_VARIABLE

PRESSURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.854e06+14.513e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.176e06+22.491e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION 0.3e06+25.60e03*z

#STOP

Since my calculations are at 4km depth, I was wondering what is the correct sign (+/-) for stresses that OGS recognizes. Is it (-) for compression and (+) for tension liek in convetional mechanics of materials. Or is it (+) for compression and negative for tension. Furthermore, what are the implications when using equations; is the depth(z) taken with the sign I modeled the geometry e.g. -4000m? I assume that all the units need to be transformed into Pascals like the pressures for the Hydraulic simulation.

1 Like

* The sign of stress or strain takes the conventional solid mechanics' one, i.e. + for tension.
* z takes the value of z coordinate of an node defined in the mesh data.

The initial stresses can be input as (in a benchmark directory of THM/init)
GeoSys-IC: Initial Conditions

···

------------------------------------------------
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
3
0 -32.1e6+0.55e5*y
1 -32.1e6+0.55e5*y
2 -32.1e6+0.55e5*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
3
0 -13243500+26487*y
1 -13243500+26487*y
2 -13243500+26487*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
SUB_DOMAIN
3
0 -10.6e6+0.2e5*y
1 -10.6e6+0.2e5*y
2 -10.6e6+0.2e5*y
#STOP

Since one one domain is shown in your data, you input can be, for example,
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
1
0 -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:

Dear Users and Developers,

I am running a 3 model for HTM processes. I have completed the HT models and seem to be working well; however, when I implemented the third effect (M); I am not getting any solution at all. I have assigned the initial conditions as they are showed in the benchmarks, for example:

GINA - Initial Condition
#INITIAL_CONDITION
\$PCS_TYPE
HEAT_TRANSPORT
\$PRIMARY_VARIABLE
TEMPERATURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
CONSTANT 453.15
#INITIAL_CONDITION
\$PCS_TYPE
LIQUID_FLOW
\$PRIMARY_VARIABLE
PRESSURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.854e06+14.513e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.176e06+22.491e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION 0.3e06+25.60e03*z
#STOP

Since my calculations are at 4km depth, I was wondering what is the correct sign (+/-) for stresses that OGS recognizes. Is it (-) for compression and (+) for tension liek in convetional mechanics of materials. Or is it (+) for compression and negative for tension. Furthermore, what are the implications when using equations; is the depth(z) taken with the sign I modeled the geometry e.g. -4000m? I assume that all the units need to be transformed into Pascals like the pressures for the Hydraulic simulation.

--
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>.

Dear Wenqing Wang

Regarding to your explanation about the sub-domain. I have 5 material groups actually; however I decided to assign the initial conditions to the WHOLE DOMAIN.

My material groups are numbered from 0,1,… to 4 and four materials are discrete fractures. The material No4 is the rock matrix. Let me know if I am mistaken please.

Thanks again!

Jhonny

···

On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:

``````  * The sign of stress or strain takes the conventional solid mechanics' one, i.e. + for tension.

* z takes the value of z coordinate of  an node defined in the mesh data.

The initial stresses can be input as (in a benchmark directory of THM/init)

GeoSys-IC: Initial Conditions ------------------------------------------------

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

SUB_DOMAIN

3

0   -32.1e6+0.55e5*y

1     -32.1e6+0.55e5*y

2    -32.1e6+0.55e5*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

3

0   -13243500+26487*y

1   -13243500+26487*y

2   -13243500+26487*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

SUB_DOMAIN

3

0   -10.6e6+0.2e5*y

1   -10.6e6+0.2e5*y

2   -10.6e6+0.2e5*y

#STOP

Since one one domain is shown in your data, you input can be, for example,

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

1

0  -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
``````

Dear Users and Developers,

``````      I am running a 3 model for HTM processes. I have completed the HT models and seem to be working well; however, when I implemented the third effect (M); I am not getting any solution at all. I have assigned the initial conditions as they are showed in the benchmarks, for example:
``````

GINA - Initial Condition

#INITIAL_CONDITION

\$PCS_TYPE

HEAT_TRANSPORT

\$PRIMARY_VARIABLE

TEMPERATURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

CONSTANT 453.15

#INITIAL_CONDITION

\$PCS_TYPE

LIQUID_FLOW

\$PRIMARY_VARIABLE

PRESSURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.854e06+14.513e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.176e06+22.491e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION 0.3e06+25.60e03*z

#STOP

``````      Since my calculations are at 4km depth, I was wondering what is the correct sign (+/-) for stresses that OGS recognizes. Is it (-) for compression and (+) for tension liek in convetional mechanics of materials. Or is it (+) for compression and negative for tension. Furthermore, what are the implications when using equations; is the depth(z) taken with the sign I modeled the geometry e.g. -4000m?  I assume that all the units need to be transformed into Pascals like the pressures for the Hydraulic simulation.
``````
``````      Please let me know your opinions and advice. I will be very glad.
``````

``````  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+...@googlegroups.com.

``````

Then you have to input the initial stress data as
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
5
0 [expression of stress]
1 [expression of stress]
2 [expression of stress]
3 [expression of stress]
4 [expression of stress]

···

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:

Dear Wenqing Wang

Regarding to your explanation about the sub-domain. I have 5 material groups actually; however I decided to assign the initial conditions to the WHOLE DOMAIN.
My material groups are numbered from 0,1,.. to 4 and four materials are discrete fractures. The material No4 is the rock matrix. Let me know if I am mistaken please.
Thanks again!

Jhonny

On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:

* The sign of stress or strain takes the conventional solid
mechanics' one, i.e. + for tension.
* z takes the value of z coordinate of an node defined in the
mesh data.

The initial stresses can be input as (in a benchmark directory of
THM/init)
GeoSys-IC: Initial Conditions
------------------------------------------------
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
3
0 -32.1e6+0.55e5*y
1 -32.1e6+0.55e5*y
2 -32.1e6+0.55e5*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
3
0 -13243500+26487*y
1 -13243500+26487*y
2 -13243500+26487*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
SUB_DOMAIN
3
0 -10.6e6+0.2e5*y
1 -10.6e6+0.2e5*y
2 -10.6e6+0.2e5*y
#STOP

Since one one domain is shown in your data, you input can be, for
example,
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
1
0 -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:

Dear Users and Developers,

I am running a 3 model for HTM processes. I have completed the HT
models and seem to be working well; however, when I implemented
the third effect (M); I am not getting any solution at all. I
have assigned the initial conditions as they are showed in the
benchmarks, for example:

GINA - Initial Condition
#INITIAL_CONDITION
\$PCS_TYPE
HEAT_TRANSPORT
\$PRIMARY_VARIABLE
TEMPERATURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
CONSTANT 453.15
#INITIAL_CONDITION
\$PCS_TYPE
LIQUID_FLOW
\$PRIMARY_VARIABLE
PRESSURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.854e06+14.513e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.176e06+22.491e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION 0.3e06+25.60e03*z
#STOP

Since my calculations are at 4km depth, I was wondering what is
the correct sign (+/-) for stresses that OGS recognizes. Is it
(-) for compression and (+) for tension liek in convetional
mechanics of materials. Or is it (+) for compression and negative
for tension. Furthermore, what are the implications when using
equations; is the depth(z) taken with the sign I modeled the
geometry e.g. -4000m? I assume that all the units need to be
transformed into Pascals like the pressures for the Hydraulic
simulation.

-- You received this message because you are subscribed to the
To unsubscribe from this group and stop receiving emails from it,
send an email to ogs-users+...@googlegroups.com <javascript:>.

--
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>.

Thanks for the help

The initial stresses are ok now. But I am getting result of deformation = 0.00 for every node; I am wondering what is going on?? In the VTK file, there is a column full of zeros. OGS is taking more than an hour to solve just the deformation process. Do you imagine what the problem would be? I apologize for the many inquiries. I am doing a coursework and there is not much time left for the submission.

Thanks,

Jhonny

···

On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:

``````  Then you have to input the initial stress data as

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

SUB_DOMAIN

5

0  [expression of stress]

1  [expression of stress]

2  [expression of stress]

3  [expression of stress]

4  [expression of stress]

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
``````

Dear Wenqing Wang

``````        Regarding to your explanation about the sub-domain. I have 5 material groups actually; however I decided to assign the initial conditions to the WHOLE DOMAIN.
``````
``````        My material groups are numbered from 0,1,.. to 4  and four materials are discrete fractures. The material No4 is the rock matrix. Let me know if I am mistaken please.
``````

Thanks again!

Jhonny

``````      On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
``````
``````            * The sign of stress or strain takes the conventional solid mechanics' one, i.e. + for tension.

* z takes the value of z coordinate of  an node defined in the mesh data.

The initial stresses can be input as (in a benchmark directory of THM/init)

GeoSys-IC: Initial Conditions ------------------------------------------------

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

SUB_DOMAIN

3

0   -32.1e6+0.55e5*y

1     -32.1e6+0.55e5*y

2    -32.1e6+0.55e5*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

3

0   -13243500+26487*y

1   -13243500+26487*y

2   -13243500+26487*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

SUB_DOMAIN

3

0   -10.6e6+0.2e5*y

1   -10.6e6+0.2e5*y

2   -10.6e6+0.2e5*y

#STOP

Since one one domain is shown in your data, you input can be, for example,

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

1

0  -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
``````

Dear Users and Developers,

``````                I am running a 3 model for HTM processes. I have completed the HT models and seem to be working well; however, when I implemented the third effect (M); I am not getting any solution at all. I have assigned the initial conditions as they are showed in the benchmarks, for example:
``````

GINA - Initial Condition

#INITIAL_CONDITION

\$PCS_TYPE

HEAT_TRANSPORT

\$PRIMARY_VARIABLE

TEMPERATURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

CONSTANT 453.15

#INITIAL_CONDITION

\$PCS_TYPE

LIQUID_FLOW

\$PRIMARY_VARIABLE

PRESSURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.854e06+14.513e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.176e06+22.491e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION 0.3e06+25.60e03*z

#STOP

``````                Since my calculations are at 4km depth, I was wondering what is the correct sign (+/-) for stresses that OGS recognizes. Is it (-) for compression and (+) for tension liek in convetional mechanics of materials. Or is it (+) for compression and negative for tension. Furthermore, what are the implications when using equations; is the depth(z) taken with the sign I modeled the geometry e.g. -4000m?  I assume that all the units need to be transformed into Pascals like the pressures for the Hydraulic simulation.
``````
``````                Please let me know your opinions and advice. I will be very glad.
``````

``````            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+...@googlegroups.com.

``````

``````  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+...@googlegroups.com.

``````

By the way,

I was thinking that the problem could be associated to the MSP. I have assigned the Elasticity properties as follows.

GINA - Material Solid Properties

#SOLID_PROPERTIES

\$NAME

MSPFAULT

\$DENSITY

1 2.85e+03

\$ELASTICITY

POISSION 0.25

YOUNGS_MODULUS

1 6.00e+09

\$THERMAL

EXPANSION

1.000e-05

CAPACITY

1 1.017e+003

CONDUCTIVITY

1 3.00e+000

\$BIOT_CONSTANT

1.0

BISHOP_COEFFICIENT

3 0.0 ; model_3: bishop=0.0 or 1.0

\$GRAVITY_CONSTANT

0

\$SOLID_BULK_MODULUS

1.e-9

#SOLID_PROPERTIES

\$NAME

MSPFRACT2

\$DENSITY

1 2.85e+003

\$ELASTICITY

POISSION 0.25

YOUNGS_MODULUS

1 6.00e+09

\$THERMAL

EXPANSION

1.00e-005

CAPACITY

1 1.017e+003

CONDUCTIVITY

1 3.00e+000

; BISHOP-BIOT Model:

; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w

\$BIOT_CONSTANT

1.0

BISHOP_COEFFICIENT

3 0.0 ; model_3: bishop=0.0 or 1.0

\$GRAVITY_CONSTANT

0

\$SOLID_BULK_MODULUS

1.e-9

#SOLID_PROPERTIES

\$NAME

MSPFRACT3

\$DENSITY

1 2.85e+03

\$ELASTICITY

POISSION 0.25

YOUNGS_MODULUS

1 6.00e+09

\$THERMAL

EXPANSION

1.000e-05

CAPACITY

1 1.017e+003

CONDUCTIVITY

1 3.00e+000

; BISHOP-BIOT Model:

; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w

\$BIOT_CONSTANT

1.0

BISHOP_COEFFICIENT

3 0.0 ; model_3: bishop=0.0 or 1.0

\$GRAVITY_CONSTANT

0

\$SOLID_BULK_MODULUS

1.e-9

#SOLID_PROPERTIES

\$NAME

MSPFRACT4

\$DENSITY

1 2.85e+03

\$ELASTICITY

POISSION 0.25

YOUNGS_MODULUS

1 6.00e+09

\$THERMAL

EXPANSION

1.00e-05

CAPACITY

1 1.017e+003

CONDUCTIVITY

1 3.00e+000

; BISHOP-BIOT Model:

; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w

\$BIOT_CONSTANT

1.0

BISHOP_COEFFICIENT

3 0.0 ; model_3: bishop=0.0 or 1.0

\$GRAVITY_CONSTANT

0

\$SOLID_BULK_MODULUS

1.e-9

#SOLID_PROPERTIES

\$NAME

MSPROCKMATRIX

\$DENSITY

1 2.850e+003

\$ELASTICITY

POISSION 0.25

YOUNGS_MODULUS

1 6.00e+010

\$THERMAL

EXPANSION

1.00e-05

CAPACITY

1 1.017e+003

CONDUCTIVITY

1 3.00e+000

; BISHOP-BIOT Model:

; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w

\$BIOT_CONSTANT

1.0

BISHOP_COEFFICIENT

3 0.0 ; model_3: bishop=0.0 or 1.0

\$GRAVITY_CONSTANT

0

\$SOLID_BULK_MODULUS

1.e-9

#STOP

\$PLASTICITY

MOHR-COULOMB

4903500

45 ;friction angle

0

2451750

0

0

···

I was wondering what characteristic I am missing. Are the friction angles in radiant? I was also thinking that the problem might be in the numerical NUM. Do I have to consider the overall coupling so that the thermal effects cause strain in my material?

Thanks!

On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:

``````  Then you have to input the initial stress data as

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

SUB_DOMAIN

5

0  [expression of stress]

1  [expression of stress]

2  [expression of stress]

3  [expression of stress]

4  [expression of stress]

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
``````

Dear Wenqing Wang

``````        Regarding to your explanation about the sub-domain. I have 5 material groups actually; however I decided to assign the initial conditions to the WHOLE DOMAIN.
``````
``````        My material groups are numbered from 0,1,.. to 4  and four materials are discrete fractures. The material No4 is the rock matrix. Let me know if I am mistaken please.
``````

Thanks again!

Jhonny

``````      On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
``````
``````            * The sign of stress or strain takes the conventional solid mechanics' one, i.e. + for tension.

* z takes the value of z coordinate of  an node defined in the mesh data.

The initial stresses can be input as (in a benchmark directory of THM/init)

GeoSys-IC: Initial Conditions ------------------------------------------------

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

SUB_DOMAIN

3

0   -32.1e6+0.55e5*y

1     -32.1e6+0.55e5*y

2    -32.1e6+0.55e5*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

3

0   -13243500+26487*y

1   -13243500+26487*y

2   -13243500+26487*y

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

SUB_DOMAIN

3

0   -10.6e6+0.2e5*y

1   -10.6e6+0.2e5*y

2   -10.6e6+0.2e5*y

#STOP

Since one one domain is shown in your data, you input can be, for example,

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

SUB_DOMAIN

1

0  -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
``````

Dear Users and Developers,

``````                I am running a 3 model for HTM processes. I have completed the HT models and seem to be working well; however, when I implemented the third effect (M); I am not getting any solution at all. I have assigned the initial conditions as they are showed in the benchmarks, for example:
``````

GINA - Initial Condition

#INITIAL_CONDITION

\$PCS_TYPE

HEAT_TRANSPORT

\$PRIMARY_VARIABLE

TEMPERATURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

CONSTANT 453.15

#INITIAL_CONDITION

\$PCS_TYPE

LIQUID_FLOW

\$PRIMARY_VARIABLE

PRESSURE1

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_XX

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.854e06+14.513e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_YY

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION -1.176e06+22.491e03*z

#INITIAL_CONDITION

\$PCS_TYPE

DEFORMATION

\$PRIMARY_VARIABLE

STRESS_ZZ

\$GEO_TYPE

DOMAIN

\$DIS_TYPE

FUNCTION 0.3e06+25.60e03*z

#STOP

``````                Since my calculations are at 4km depth, I was wondering what is the correct sign (+/-) for stresses that OGS recognizes. Is it (-) for compression and (+) for tension liek in convetional mechanics of materials. Or is it (+) for compression and negative for tension. Furthermore, what are the implications when using equations; is the depth(z) taken with the sign I modeled the geometry e.g. -4000m?  I assume that all the units need to be transformed into Pascals like the pressures for the Hydraulic simulation.
``````
``````                Please let me know your opinions and advice. I will be very glad.
``````

``````            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+...@googlegroups.com.

``````

``````  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+...@googlegroups.com.

``````

The MSP looks not bad, though I don't think the following lines are necessary (why the gravity constant to zero?).

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

Have you checked whether a linear solver successfully converged?

Best,
Nori

···

On 07/13/2016 11:40 AM, Jhonny Encalada wrote:

By the way,

I was thinking that the problem could be associated to the MSP. I have
assigned the Elasticity properties as follows.

GINA - Material Solid Properties
#SOLID_PROPERTIES
\$NAME
MSPFAULT
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT2
\$DENSITY
1 2.85e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-005
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT3
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT4
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPROCKMATRIX
\$DENSITY
1 2.850e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+010
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9
#STOP

\$PLASTICITY
MOHR-COULOMB
4903500
45 ;friction angle
0
2451750
0
----------------------------------------------------------------------------------------------------------
I was wondering what characteristic I am missing. Are the friction
angles in radiant? I was also thinking that the problem might be in the
numerical NUM. Do I have to consider the overall coupling so that the
thermal effects cause strain in my material?
Thanks!

On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:

Then you have to input the initial stress data as
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
5
0 [expression of stress]
1 [expression of stress]
2 [expression of stress]
3 [expression of stress]
4 [expression of stress]

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
> Dear Wenqing Wang
>
>
> material groups actually; however I decided to assign the initial
> conditions to the WHOLE DOMAIN.
> My material groups are numbered from 0,1,.. to 4 and four
> materials are discrete fractures. The material No4 is the rock
> matrix. Let me know if I am mistaken please.
> Thanks again!
>
> Jhonny
>
> On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
>
> * The sign of stress or strain takes the conventional solid
> mechanics' one, i.e. + for tension.
> * z takes the value of z coordinate of an node defined in the
> mesh data.
>
> The initial stresses can be input as (in a benchmark directory
> of THM/init)
> GeoSys-IC: Initial Conditions
> ------------------------------------------------
> #INITIAL_CONDITION
> \$PCS_TYPE
> DEFORMATION
> \$PRIMARY_VARIABLE
> STRESS_XX
> \$GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -32.1e6+0.55e5*y
> 1 -32.1e6+0.55e5*y
> 2 -32.1e6+0.55e5*y
> #INITIAL_CONDITION
> \$PCS_TYPE
> DEFORMATION
> \$PRIMARY_VARIABLE
> STRESS_YY
> \$GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -13243500+26487*y
> 1 -13243500+26487*y
> 2 -13243500+26487*y
> #INITIAL_CONDITION
> \$PCS_TYPE
> DEFORMATION
> \$PRIMARY_VARIABLE
> STRESS_ZZ
> \$GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -10.6e6+0.2e5*y
> 1 -10.6e6+0.2e5*y
> 2 -10.6e6+0.2e5*y
> #STOP
>
> Since one one domain is shown in your data, you input can be,
> for example,
> #INITIAL_CONDITION
> \$PCS_TYPE
> DEFORMATION
> \$PRIMARY_VARIABLE
> STRESS_YY
> \$GEO_TYPE
> SUB_DOMAIN
> 1
> 0 -1.176e06+22.491e03*z
>
> On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
>> Dear Users and Developers,
>>
>> I am running a 3 model for HTM processes. I have completed
>> the HT models and seem to be working well; however, when I
>> implemented the third effect (M); I am not getting any
>> solution at all. I have assigned the initial conditions as
>> they are showed in the benchmarks, for example:
>>
>> GINA - Initial Condition
>> #INITIAL_CONDITION
>> \$PCS_TYPE
>> HEAT_TRANSPORT
>> \$PRIMARY_VARIABLE
>> TEMPERATURE1
>> \$GEO_TYPE
>> DOMAIN
>> \$DIS_TYPE
>> CONSTANT 453.15
>> #INITIAL_CONDITION
>> \$PCS_TYPE
>> LIQUID_FLOW
>> \$PRIMARY_VARIABLE
>> PRESSURE1
>> \$GEO_TYPE
>> DOMAIN
>> \$DIS_TYPE
>> #INITIAL_CONDITION
>> \$PCS_TYPE
>> DEFORMATION
>> \$PRIMARY_VARIABLE
>> STRESS_XX
>> \$GEO_TYPE
>> DOMAIN
>> \$DIS_TYPE
>> FUNCTION -1.854e06+14.513e03*z
>> #INITIAL_CONDITION
>> \$PCS_TYPE
>> DEFORMATION
>> \$PRIMARY_VARIABLE
>> STRESS_YY
>> \$GEO_TYPE
>> DOMAIN
>> \$DIS_TYPE
>> FUNCTION -1.176e06+22.491e03*z
>> #INITIAL_CONDITION
>> \$PCS_TYPE
>> DEFORMATION
>> \$PRIMARY_VARIABLE
>> STRESS_ZZ
>> \$GEO_TYPE
>> DOMAIN
>> \$DIS_TYPE
>> FUNCTION 0.3e06+25.60e03*z
>> #STOP
>>
>> Since my calculations are at 4km depth, I was wondering what
>> is the correct sign (+/-) for stresses that OGS recognizes.
>> Is it (-) for compression and (+) for tension liek in
>> convetional mechanics of materials. Or is it (+) for
>> compression and negative for tension. Furthermore, what are
>> the implications when using equations; is the depth(z) taken
>> with the sign I modeled the geometry e.g. -4000m? I assume
>> that all the units need to be transformed into Pascals like
>> the pressures for the Hydraulic simulation.
>>
>> --
>> You received this message because you are subscribed to the
>> To unsubscribe from this group and stop receiving emails from
>> it, send an email to ogs-users+...@googlegroups.com.
>> For more options, visit https://groups.google.com/d/optout
>
> --
> 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+...@googlegroups.com <javascript:>.
> For more options, visit https://groups.google.com/d/optout

--
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch

Dear Norihiro,

Thanks for the reply. I have modified the MSP file and eliminated the unnecessary information. I also modified to NUM file so the calculation can converge fast. I dont have problems with the H and T processes; however the Deformation PCS is still with defect. In the results, I am getting “zero effective strain” and the stress invariants “p” and “q” have reasonable values for some sones “sinks”. The diagram of stresses in the VTK files is giving me a correct initial stress distribution; the stress in the fracture zones is decreasing with time; however, there are !no strains to show!. The deformation process takes about 30 min per step. and H and T take 2-3 min per step. When I use overall coupling, the deformations process does not converge and I have been running 1 step for almost 4 hours without any result.

About the gravity constant, I assumed it was zero, because I assumed an initial stress field, and I thought the thermal expansion was the unique effect to consider and not gravitational loads. Is this assumption correct? I am a little confused about the “gravity constant”. I was wondering whether it is “gravity aceleratio” or any other constant. I don’t understand its function very well.

Regarding to the PERMEABILITY FUCTION STRAIN. I found some functions, but are for PERMEABILITY_SATURATION. In the OGS, is the relative permeability in the range [0, 1]?. I want to know this, because I am assuming any function to run the model.

I will be glad that you give me some clues about the defects of my model. Thanks in advance!

···

On Monday, 18 July 2016 06:14:10 UTC+1, Norihiro Watanabe wrote:

The MSP looks not bad, though I don’t think the following lines are necessary (why the gravity constant to zero?).

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

Have you checked whether a linear solver successfully converged?

Best,
Nori

On 07/13/2016 11:40 AM, Jhonny Encalada wrote:

By the way,

I was thinking that the problem could be associated to the MSP. I have
assigned the Elasticity properties as follows.

GINA - Material Solid Properties
#SOLID_PROPERTIES
\$NAME
MSPFAULT
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT2
\$DENSITY
1 2.85e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-005
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT3
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT4
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPROCKMATRIX
\$DENSITY
1 2.850e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+010
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9
#STOP

## \$PLASTICITY MOHR-COULOMB 4903500 45 ;friction angle 0 2451750 0 0

I was wondering what characteristic I am missing. Are the friction
angles in radiant? I was also thinking that the problem might be in the
numerical NUM. Do I have to consider the overall coupling so that the
thermal effects cause strain in my material?
Thanks!

On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:

``````Then you have to input the initial stress data as
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
5
0  [expression of stress]
1  [expression of stress]
2  [expression of stress]
3  [expression of stress]
4  [expression of stress]

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
``````
``````Dear Wenqing Wang

material groups actually; however I decided to assign the initial
conditions to the WHOLE DOMAIN.
My material groups are numbered from 0,1,.. to 4  and four
materials are discrete fractures. The material No4 is the rock
matrix. Let me know if I am mistaken please.
Thanks again!

Jhonny

On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:

* The sign of stress or strain takes the conventional solid
mechanics' one, i.e. + for tension.
* z takes the value of z coordinate of  an node defined in the
mesh data.

The initial stresses can be input as (in a benchmark directory
of THM/init)
GeoSys-IC: Initial Conditions
------------------------------------------------
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
3
0   -32.1e6+0.55e5*y
1     -32.1e6+0.55e5*y
2    -32.1e6+0.55e5*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
3
0   -13243500+26487*y
1   -13243500+26487*y
2   -13243500+26487*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
SUB_DOMAIN
3
0   -10.6e6+0.2e5*y
1   -10.6e6+0.2e5*y
2   -10.6e6+0.2e5*y
#STOP

Since one one domain is shown in your data, you input can be,
for example,
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
1
0  -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
``````
``````    Dear Users and Developers,

I am running a 3 model for HTM processes. I have completed
the HT models and seem to be working well; however, when I
implemented the third effect (M); I am not getting any
solution at all. I have assigned the initial conditions as
they are showed in the benchmarks, for example:

GINA - Initial Condition
#INITIAL_CONDITION
\$PCS_TYPE
HEAT_TRANSPORT
\$PRIMARY_VARIABLE
TEMPERATURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
CONSTANT   453.15
#INITIAL_CONDITION
\$PCS_TYPE
LIQUID_FLOW
\$PRIMARY_VARIABLE
PRESSURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.854e06+14.513e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION  -1.176e06+22.491e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION 0.3e06+25.60e03*z
#STOP

Since my calculations are at 4km depth, I was wondering what
is the correct sign (+/-) for stresses that OGS recognizes.
Is it (-) for compression and (+) for tension liek in
convetional mechanics of materials. Or is it (+) for
compression and negative for tension. Furthermore, what are
the implications when using equations; is the depth(z) taken
with the sign I modeled the geometry e.g. -4000m?  I assume
that all the units need to be transformed into Pascals like
the pressures for the Hydraulic simulation.

--
You received this message because you are subscribed to the
To unsubscribe from this group and stop receiving emails from
it, send an email to ogs-users+...@googlegroups.com.
``````
``````    <[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>.
``````
``````--
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it,
send an email to ogs-users+...@googlegroups.com <javascript:>.
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>.
``````

Groups “ogs-users” group.
To unsubscribe from this group and stop receiving emails from it, send

Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch

By the way,

I enclose a compilation of the files I wrote. I hope you can check it.

Thanks!

model.txt (6.59 KB)

···

On Monday, 18 July 2016 06:14:10 UTC+1, Norihiro Watanabe wrote:

The MSP looks not bad, though I don’t think the following lines are necessary (why the gravity constant to zero?).

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

Have you checked whether a linear solver successfully converged?

Best,
Nori

On 07/13/2016 11:40 AM, Jhonny Encalada wrote:

By the way,

I was thinking that the problem could be associated to the MSP. I have
assigned the Elasticity properties as follows.

GINA - Material Solid Properties
#SOLID_PROPERTIES
\$NAME
MSPFAULT
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT2
\$DENSITY
1 2.85e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-005
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT3
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPFRACT4
\$DENSITY
1 2.85e+03
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

#SOLID_PROPERTIES
\$NAME
MSPROCKMATRIX
\$DENSITY
1 2.850e+003
\$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+010
\$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000

; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9
#STOP

## \$PLASTICITY MOHR-COULOMB 4903500 45 ;friction angle 0 2451750 0 0

I was wondering what characteristic I am missing. Are the friction
angles in radiant? I was also thinking that the problem might be in the
numerical NUM. Do I have to consider the overall coupling so that the
thermal effects cause strain in my material?
Thanks!

On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:

``````Then you have to input the initial stress data as
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
5
0  [expression of stress]
1  [expression of stress]
2  [expression of stress]
3  [expression of stress]
4  [expression of stress]

On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
``````
``````Dear Wenqing Wang

material groups actually; however I decided to assign the initial
conditions to the WHOLE DOMAIN.
My material groups are numbered from 0,1,.. to 4  and four
materials are discrete fractures. The material No4 is the rock
matrix. Let me know if I am mistaken please.
Thanks again!

Jhonny

On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:

* The sign of stress or strain takes the conventional solid
mechanics' one, i.e. + for tension.
* z takes the value of z coordinate of  an node defined in the
mesh data.

The initial stresses can be input as (in a benchmark directory
of THM/init)
GeoSys-IC: Initial Conditions
------------------------------------------------
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
SUB_DOMAIN
3
0   -32.1e6+0.55e5*y
1     -32.1e6+0.55e5*y
2    -32.1e6+0.55e5*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
3
0   -13243500+26487*y
1   -13243500+26487*y
2   -13243500+26487*y
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
SUB_DOMAIN
3
0   -10.6e6+0.2e5*y
1   -10.6e6+0.2e5*y
2   -10.6e6+0.2e5*y
#STOP

Since one one domain is shown in your data, you input can be,
for example,
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
SUB_DOMAIN
1
0  -1.176e06+22.491e03*z

On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
``````
``````    Dear Users and Developers,

I am running a 3 model for HTM processes. I have completed
the HT models and seem to be working well; however, when I
implemented the third effect (M); I am not getting any
solution at all. I have assigned the initial conditions as
they are showed in the benchmarks, for example:

GINA - Initial Condition
#INITIAL_CONDITION
\$PCS_TYPE
HEAT_TRANSPORT
\$PRIMARY_VARIABLE
TEMPERATURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
CONSTANT   453.15
#INITIAL_CONDITION
\$PCS_TYPE
LIQUID_FLOW
\$PRIMARY_VARIABLE
PRESSURE1
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_XX
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION -1.854e06+14.513e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_YY
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION  -1.176e06+22.491e03*z
#INITIAL_CONDITION
\$PCS_TYPE
DEFORMATION
\$PRIMARY_VARIABLE
STRESS_ZZ
\$GEO_TYPE
DOMAIN
\$DIS_TYPE
FUNCTION 0.3e06+25.60e03*z
#STOP

Since my calculations are at 4km depth, I was wondering what
is the correct sign (+/-) for stresses that OGS recognizes.
Is it (-) for compression and (+) for tension liek in
convetional mechanics of materials. Or is it (+) for
compression and negative for tension. Furthermore, what are
the implications when using equations; is the depth(z) taken
with the sign I modeled the geometry e.g. -4000m?  I assume
that all the units need to be transformed into Pascals like
the pressures for the Hydraulic simulation.

--
You received this message because you are subscribed to the
To unsubscribe from this group and stop receiving emails from
it, send an email to ogs-users+...@googlegroups.com.
``````
``````    <[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>.
``````
``````--
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it,
send an email to ogs-users+...@googlegroups.com <javascript:>.
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>.
``````

Groups “ogs-users” group.
To unsubscribe from this group and stop receiving emails from it, send

Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch

Thanks for the file. It helps

The problem in your setting is that you are trying to use fracture elements for DEFORMATION process. I'm afraid that current OGS does not support it yet. To include mechanical effects into your simulation, you need to either

1) exclude them from mechanic simulation with \$DEACTIVATED_SUBDOMAIN keyword in PCS (i.e. no opening and slip effect of fractures)

https://svn.ufz.de/ogs/wiki/public/doc-auto/by_ext/pcs/S_deactivated_subdomain

or 2) use thin volumetric elements for the fractures

Regards,
Nori

···

On 07/18/2016 10:02 AM, Jhonny Encalada wrote:

By the way,

I enclose a compilation of the files I wrote. I hope you can check it.

Thanks!

On Monday, 18 July 2016 06:14:10 UTC+1, Norihiro Watanabe wrote:

The MSP looks not bad, though I don't think the following lines are
necessary (why the gravity constant to zero?).

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

Have you checked whether a linear solver successfully converged?

Best,
Nori

On 07/13/2016 11:40 AM, Jhonny Encalada wrote:
> By the way,
>
> I was thinking that the problem could be associated to the MSP. I
have
> assigned the Elasticity properties as follows.
>
> GINA - Material Solid Properties
> #SOLID_PROPERTIES
> \$NAME
> MSPFAULT
> \$DENSITY
> 1 2.85e+03
> \$ELASTICITY
> POISSION 0.25
> YOUNGS_MODULUS
> 1 6.00e+09
> \$THERMAL
> EXPANSION
> 1.000e-05
> CAPACITY
> 1 1.017e+003
> CONDUCTIVITY
> 1 3.00e+000
>
> \$BIOT_CONSTANT
> 1.0
> BISHOP_COEFFICIENT
> 3 0.0 ; model_3: bishop=0.0 or 1.0
> \$GRAVITY_CONSTANT
> 0
> \$SOLID_BULK_MODULUS
> 1.e-9
>
> #SOLID_PROPERTIES
> \$NAME
> MSPFRACT2
> \$DENSITY
> 1 2.85e+003
> \$ELASTICITY
> POISSION 0.25
> YOUNGS_MODULUS
> 1 6.00e+09
> \$THERMAL
> EXPANSION
> 1.00e-005
> CAPACITY
> 1 1.017e+003
> CONDUCTIVITY
> 1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
> \$BIOT_CONSTANT
> 1.0
> BISHOP_COEFFICIENT
> 3 0.0 ; model_3: bishop=0.0 or 1.0
> \$GRAVITY_CONSTANT
> 0
> \$SOLID_BULK_MODULUS
> 1.e-9
>
> #SOLID_PROPERTIES
> \$NAME
> MSPFRACT3
> \$DENSITY
> 1 2.85e+03
> \$ELASTICITY
> POISSION 0.25
> YOUNGS_MODULUS
> 1 6.00e+09
> \$THERMAL
> EXPANSION
> 1.000e-05
> CAPACITY
> 1 1.017e+003
> CONDUCTIVITY
> 1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
> \$BIOT_CONSTANT
> 1.0
> BISHOP_COEFFICIENT
> 3 0.0 ; model_3: bishop=0.0 or 1.0
> \$GRAVITY_CONSTANT
> 0
> \$SOLID_BULK_MODULUS
> 1.e-9
>
> #SOLID_PROPERTIES
> \$NAME
> MSPFRACT4
> \$DENSITY
> 1 2.85e+03
> \$ELASTICITY
> POISSION 0.25
> YOUNGS_MODULUS
> 1 6.00e+09
> \$THERMAL
> EXPANSION
> 1.00e-05
> CAPACITY
> 1 1.017e+003
> CONDUCTIVITY
> 1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
> \$BIOT_CONSTANT
> 1.0
> BISHOP_COEFFICIENT
> 3 0.0 ; model_3: bishop=0.0 or 1.0
> \$GRAVITY_CONSTANT
> 0
> \$SOLID_BULK_MODULUS
> 1.e-9
>
> #SOLID_PROPERTIES
> \$NAME
> MSPROCKMATRIX
> \$DENSITY
> 1 2.850e+003
> \$ELASTICITY
> POISSION 0.25
> YOUNGS_MODULUS
> 1 6.00e+010
> \$THERMAL
> EXPANSION
> 1.00e-05
> CAPACITY
> 1 1.017e+003
> CONDUCTIVITY
> 1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
> \$BIOT_CONSTANT
> 1.0
> BISHOP_COEFFICIENT
> 3 0.0 ; model_3: bishop=0.0 or 1.0
> \$GRAVITY_CONSTANT
> 0
> \$SOLID_BULK_MODULUS
> 1.e-9
> #STOP
>
> \$PLASTICITY
> MOHR-COULOMB
> 4903500
> 45 ;friction angle
> 0
> 2451750
> 0
>
----------------------------------------------------------------------------------------------------------

> I was wondering what characteristic I am missing. Are the friction
> angles in radiant? I was also thinking that the problem might be
in the
> numerical NUM. Do I have to consider the overall coupling so that
the
> thermal effects cause strain in my material?
> Thanks!
>
> On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:
>
> Then you have to input the initial stress data as
> #INITIAL_CONDITION
> \$PCS_TYPE
> DEFORMATION
> \$PRIMARY_VARIABLE
> STRESS_XX
> \$GEO_TYPE
> SUB_DOMAIN
> 5
> 0 [expression of stress]
> 1 [expression of stress]
> 2 [expression of stress]
> 3 [expression of stress]
> 4 [expression of stress]
>
> On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
> > Dear Wenqing Wang
> >
> > Thanks for your help,
> >
> > Regarding to your explanation about the sub-domain. I have 5
> > material groups actually; however I decided to assign the
initial
> > conditions to the WHOLE DOMAIN.
> > My material groups are numbered from 0,1,.. to 4 and four
> > materials are discrete fractures. The material No4 is the rock
> > matrix. Let me know if I am mistaken please.
> > Thanks again!
> >
> > Jhonny
> >
> > On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
> >
> > * The sign of stress or strain takes the conventional
solid
> > mechanics' one, i.e. + for tension.
> > * z takes the value of z coordinate of an node defined
in the
> > mesh data.
> >
> > The initial stresses can be input as (in a benchmark
directory
> > of THM/init)
> > GeoSys-IC: Initial Conditions
> > ------------------------------------------------
> > #INITIAL_CONDITION
> > \$PCS_TYPE
> > DEFORMATION
> > \$PRIMARY_VARIABLE
> > STRESS_XX
> > \$GEO_TYPE
> > SUB_DOMAIN
> > 3
> > 0 -32.1e6+0.55e5*y
> > 1 -32.1e6+0.55e5*y
> > 2 -32.1e6+0.55e5*y
> > #INITIAL_CONDITION
> > \$PCS_TYPE
> > DEFORMATION
> > \$PRIMARY_VARIABLE
> > STRESS_YY
> > \$GEO_TYPE
> > SUB_DOMAIN
> > 3
> > 0 -13243500+26487*y
> > 1 -13243500+26487*y
> > 2 -13243500+26487*y
> > #INITIAL_CONDITION
> > \$PCS_TYPE
> > DEFORMATION
> > \$PRIMARY_VARIABLE
> > STRESS_ZZ
> > \$GEO_TYPE
> > SUB_DOMAIN
> > 3
> > 0 -10.6e6+0.2e5*y
> > 1 -10.6e6+0.2e5*y
> > 2 -10.6e6+0.2e5*y
> > #STOP
> >
> > Since one one domain is shown in your data, you input
can be,
> > for example,
> > #INITIAL_CONDITION
> > \$PCS_TYPE
> > DEFORMATION
> > \$PRIMARY_VARIABLE
> > STRESS_YY
> > \$GEO_TYPE
> > SUB_DOMAIN
> > 1
> > 0 -1.176e06+22.491e03*z
> >
> > On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
> >> Dear Users and Developers,
> >>
> >> I am running a 3 model for HTM processes. I have
completed
> >> the HT models and seem to be working well; however,
when I
> >> implemented the third effect (M); I am not getting any
> >> solution at all. I have assigned the initial
conditions as
> >> they are showed in the benchmarks, for example:
> >>
> >> GINA - Initial Condition
> >> #INITIAL_CONDITION
> >> \$PCS_TYPE
> >> HEAT_TRANSPORT
> >> \$PRIMARY_VARIABLE
> >> TEMPERATURE1
> >> \$GEO_TYPE
> >> DOMAIN
> >> \$DIS_TYPE
> >> CONSTANT 453.15
> >> #INITIAL_CONDITION
> >> \$PCS_TYPE
> >> LIQUID_FLOW
> >> \$PRIMARY_VARIABLE
> >> PRESSURE1
> >> \$GEO_TYPE
> >> DOMAIN
> >> \$DIS_TYPE
> >> #INITIAL_CONDITION
> >> \$PCS_TYPE
> >> DEFORMATION
> >> \$PRIMARY_VARIABLE
> >> STRESS_XX
> >> \$GEO_TYPE
> >> DOMAIN
> >> \$DIS_TYPE
> >> FUNCTION -1.854e06+14.513e03*z
> >> #INITIAL_CONDITION
> >> \$PCS_TYPE
> >> DEFORMATION
> >> \$PRIMARY_VARIABLE
> >> STRESS_YY
> >> \$GEO_TYPE
> >> DOMAIN
> >> \$DIS_TYPE
> >> FUNCTION -1.176e06+22.491e03*z
> >> #INITIAL_CONDITION
> >> \$PCS_TYPE
> >> DEFORMATION
> >> \$PRIMARY_VARIABLE
> >> STRESS_ZZ
> >> \$GEO_TYPE
> >> DOMAIN
> >> \$DIS_TYPE
> >> FUNCTION 0.3e06+25.60e03*z
> >> #STOP
> >>
> >> Since my calculations are at 4km depth, I was
wondering what
> >> is the correct sign (+/-) for stresses that OGS
recognizes.
> >> Is it (-) for compression and (+) for tension liek in
> >> convetional mechanics of materials. Or is it (+) for
> >> compression and negative for tension. Furthermore,
what are
> >> the implications when using equations; is the depth(z)
taken
> >> with the sign I modeled the geometry e.g. -4000m? I
assume
> >> that all the units need to be transformed into Pascals
like
> >> the pressures for the Hydraulic simulation.
> >>
> >> --
> >> 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+...@googlegroups.com.
> >> For more options, visit
> >
> > --
> > You received this message because you are subscribed to the
> > Groups "ogs-users" group.
> > To unsubscribe from this group and stop receiving emails
from it,
> > send an email to ogs-users+...@googlegroups.com <javascript:>.
> > For more options, visit https://groups.google.com/d/optout
>
> --
> 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+...@googlegroups.com <javascript:>
> For more options, visit https://groups.google.com/d/optout

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro...@ufz.de <javascript:> / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch

--
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch

Dear Norihiro,

The first option sounds good; for the second I would have to change everything I have done.

Thank you !

···

On Monday, 18 July 2016 09:11:41 UTC+1, Norihiro Watanabe wrote:

Thanks for the file. It helps

The problem in your setting is that you are trying to use fracture elements for DEFORMATION process. I’m afraid that current OGS does not support it yet. To include mechanical effects into your simulation, you need to either

1. exclude them from mechanic simulation with \$DEACTIVATED_SUBDOMAIN keyword in PCS (i.e. no opening and slip effect of fractures)

https://svn.ufz.de/ogs/wiki/public/doc-auto/by_ext/pcs/S_deactivated_subdomain

or 2) use thin volumetric elements for the fractures

Regards,
Nori

On 07/18/2016 10:02 AM, Jhonny Encalada wrote:

By the way,

I enclose a compilation of the files I wrote. I hope you can check it.

Thanks!

On Monday, 18 July 2016 06:14:10 UTC+1, Norihiro Watanabe wrote:

``````The MSP looks not bad, though I don't think the following lines are
necessary (why the gravity constant to zero?).

\$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3   0.0 ; model_3: bishop=0.0 or 1.0
\$GRAVITY_CONSTANT
0
\$SOLID_BULK_MODULUS
1.e-9

Have you checked whether a linear solver successfully converged?

Best,
Nori

On 07/13/2016 11:40 AM, Jhonny Encalada wrote:
> By the way,
>
> I was thinking that the problem could be associated to the MSP. I
have
> assigned the Elasticity properties as follows.
>
>
> GINA - Material Solid Properties
> #SOLID_PROPERTIES
>   \$NAME
> MSPFAULT
>   \$DENSITY
>    1 2.85e+03
>   \$ELASTICITY
>    POISSION 0.25
>    YOUNGS_MODULUS
>    1 6.00e+09
>   \$THERMAL
>    EXPANSION
>    1.000e-05
>    CAPACITY
>    1 1.017e+003
>    CONDUCTIVITY
>    1 3.00e+000
>
>   \$BIOT_CONSTANT
>    1.0
>   BISHOP_COEFFICIENT
>    3   0.0 ; model_3: bishop=0.0 or 1.0
>   \$GRAVITY_CONSTANT
>    0
>   \$SOLID_BULK_MODULUS
>    1.e-9
>
> #SOLID_PROPERTIES
>   \$NAME
> MSPFRACT2
>   \$DENSITY
>    1 2.85e+003
>   \$ELASTICITY
>    POISSION 0.25
>    YOUNGS_MODULUS
>    1 6.00e+09
>   \$THERMAL
>    EXPANSION
>    1.00e-005
>    CAPACITY
>    1 1.017e+003
>    CONDUCTIVITY
>    1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
>   \$BIOT_CONSTANT
>    1.0
>   BISHOP_COEFFICIENT
>    3   0.0 ; model_3: bishop=0.0 or 1.0
>   \$GRAVITY_CONSTANT
>    0
>   \$SOLID_BULK_MODULUS
>    1.e-9
>
> #SOLID_PROPERTIES
>   \$NAME
> MSPFRACT3
>   \$DENSITY
>    1 2.85e+03
>   \$ELASTICITY
>    POISSION 0.25
>    YOUNGS_MODULUS
>    1 6.00e+09
>   \$THERMAL
>    EXPANSION
>    1.000e-05
>    CAPACITY
>    1 1.017e+003
>    CONDUCTIVITY
>    1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
>   \$BIOT_CONSTANT
>    1.0
>   BISHOP_COEFFICIENT
>    3   0.0 ; model_3: bishop=0.0 or 1.0
>   \$GRAVITY_CONSTANT
>    0
>   \$SOLID_BULK_MODULUS
>    1.e-9
>
> #SOLID_PROPERTIES
>   \$NAME
> MSPFRACT4
>   \$DENSITY
>    1 2.85e+03
>   \$ELASTICITY
>    POISSION 0.25
>    YOUNGS_MODULUS
>    1 6.00e+09
>   \$THERMAL
>    EXPANSION
>    1.00e-05
>    CAPACITY
>    1 1.017e+003
>    CONDUCTIVITY
>    1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
>   \$BIOT_CONSTANT
>    1.0
>   BISHOP_COEFFICIENT
>    3   0.0 ; model_3: bishop=0.0 or 1.0
>   \$GRAVITY_CONSTANT
>    0
>   \$SOLID_BULK_MODULUS
>    1.e-9
>
> #SOLID_PROPERTIES
>   \$NAME
>   MSPROCKMATRIX
>   \$DENSITY
>    1 2.850e+003
>   \$ELASTICITY
>    POISSION 0.25
>    YOUNGS_MODULUS
>    1 6.00e+010
>   \$THERMAL
>    EXPANSION
>    1.00e-05
>    CAPACITY
>    1 1.017e+003
>    CONDUCTIVITY
>    1 3.00e+000
>
> ; BISHOP-BIOT Model:
> ; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
>   \$BIOT_CONSTANT
>    1.0
>   BISHOP_COEFFICIENT
>    3   0.0 ; model_3: bishop=0.0 or 1.0
>   \$GRAVITY_CONSTANT
>    0
>   \$SOLID_BULK_MODULUS
>    1.e-9
> #STOP
>
>
> \$PLASTICITY
>    MOHR-COULOMB
>    4903500
>    45                    ;friction angle
>    0
>    2451750
>    0
>    0
>
----------------------------------------------------------------------------------------------------------

> I was wondering what characteristic I am missing. Are the friction
> angles in radiant? I was also thinking that the problem might be
in the
> numerical NUM. Do I have to consider the overall coupling so that
the
> thermal effects cause strain in my material?
> Thanks!
>
> On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:
>
>     Then you have to input the initial stress data as
>     #INITIAL_CONDITION
>       \$PCS_TYPE
>        DEFORMATION
>       \$PRIMARY_VARIABLE
>        STRESS_XX
>       \$GEO_TYPE
>        SUB_DOMAIN
>         5
>         0  [expression of stress]
>         1  [expression of stress]
>         2  [expression of stress]
>         3  [expression of stress]
>         4  [expression of stress]
>
>
>     On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
> >     Dear Wenqing Wang
> >
> >     Thanks for your help,
> >
> >     Regarding to your explanation about the sub-domain. I have 5
> >     material groups actually; however I decided to assign the
initial
> >     conditions to the WHOLE DOMAIN.
> >     My material groups are numbered from 0,1,.. to 4  and four
> >     materials are discrete fractures. The material No4 is the rock
> >     matrix. Let me know if I am mistaken please.
> >     Thanks again!
> >
> >
> >     Jhonny
> >
> >     On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
> >
> >         * The sign of stress or strain takes the conventional
solid
> >         mechanics' one, i.e. + for tension.
> >         * z takes the value of z coordinate of  an node defined
in the
> >         mesh data.
> >
> >         The initial stresses can be input as (in a benchmark
directory
> >         of THM/init)
> >         GeoSys-IC: Initial Conditions
> >         ------------------------------------------------
> >         #INITIAL_CONDITION
> >          \$PCS_TYPE
> >           DEFORMATION
> >          \$PRIMARY_VARIABLE
> >           STRESS_XX
> >          \$GEO_TYPE
> >           SUB_DOMAIN
> >            3
> >            0   -32.1e6+0.55e5*y
> >            1     -32.1e6+0.55e5*y
> >            2    -32.1e6+0.55e5*y
> >         #INITIAL_CONDITION
> >          \$PCS_TYPE
> >           DEFORMATION
> >          \$PRIMARY_VARIABLE
> >           STRESS_YY
> >          \$GEO_TYPE
> >           SUB_DOMAIN
> >            3
> >            0   -13243500+26487*y
> >            1   -13243500+26487*y
> >            2   -13243500+26487*y
> >         #INITIAL_CONDITION
> >          \$PCS_TYPE
> >           DEFORMATION
> >          \$PRIMARY_VARIABLE
> >           STRESS_ZZ
> >          \$GEO_TYPE
> >           SUB_DOMAIN
> >            3
> >            0   -10.6e6+0.2e5*y
> >            1   -10.6e6+0.2e5*y
> >            2   -10.6e6+0.2e5*y
> >         #STOP
> >
> >         Since one one domain is shown in your data, you input
can be,
> >         for example,
> >         #INITIAL_CONDITION
> >          \$PCS_TYPE
> >           DEFORMATION
> >          \$PRIMARY_VARIABLE
> >           STRESS_YY
> >          \$GEO_TYPE
> >           SUB_DOMAIN
> >            1
> >            0  -1.176e06+22.491e03*z
> >
> >
> >         On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
> >>         Dear Users and Developers,
> >>
> >>         I am running a 3 model for HTM processes. I have
completed
> >>         the HT models and seem to be working well; however,
when I
> >>         implemented the third effect (M); I am not getting any
> >>         solution at all. I have assigned the initial
conditions as
> >>         they are showed in the benchmarks, for example:
> >>
> >>         GINA - Initial Condition
> >>         #INITIAL_CONDITION
> >>          \$PCS_TYPE
> >>           HEAT_TRANSPORT
> >>          \$PRIMARY_VARIABLE
> >>           TEMPERATURE1
> >>          \$GEO_TYPE
> >>           DOMAIN
> >>          \$DIS_TYPE
> >>           CONSTANT   453.15
> >>         #INITIAL_CONDITION
> >>          \$PCS_TYPE
> >>           LIQUID_FLOW
> >>          \$PRIMARY_VARIABLE
> >>           PRESSURE1
> >>          \$GEO_TYPE
> >>           DOMAIN
> >>          \$DIS_TYPE
> >>         #INITIAL_CONDITION
> >>          \$PCS_TYPE
> >>           DEFORMATION
> >>          \$PRIMARY_VARIABLE
> >>           STRESS_XX
> >>          \$GEO_TYPE
> >>           DOMAIN
> >>          \$DIS_TYPE
> >>           FUNCTION -1.854e06+14.513e03*z
> >>         #INITIAL_CONDITION
> >>          \$PCS_TYPE
> >>           DEFORMATION
> >>          \$PRIMARY_VARIABLE
> >>           STRESS_YY
> >>          \$GEO_TYPE
> >>           DOMAIN
> >>          \$DIS_TYPE
> >>           FUNCTION  -1.176e06+22.491e03*z
> >>         #INITIAL_CONDITION
> >>          \$PCS_TYPE
> >>           DEFORMATION
> >>          \$PRIMARY_VARIABLE
> >>           STRESS_ZZ
> >>          \$GEO_TYPE
> >>           DOMAIN
> >>          \$DIS_TYPE
> >>            FUNCTION 0.3e06+25.60e03*z
> >>         #STOP
> >>
> >>         Since my calculations are at 4km depth, I was
wondering what
> >>         is the correct sign (+/-) for stresses that OGS
recognizes.
> >>         Is it (-) for compression and (+) for tension liek in
> >>         convetional mechanics of materials. Or is it (+) for
> >>         compression and negative for tension. Furthermore,
what are
> >>         the implications when using equations; is the depth(z)
taken
> >>         with the sign I modeled the geometry e.g. -4000m?  I
assume
> >>         that all the units need to be transformed into Pascals
like
> >>         the pressures for the Hydraulic simulation.
> >>
> >>         --
> >>         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+...@googlegroups.com.
> >>         For more options, visit
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>>.
> >
> >     --
> >     You received this message because you are subscribed to the
> >     Groups "ogs-users" group.
> >     To unsubscribe from this group and stop receiving emails
from it,
> >     send an email to ogs-users+...@googlegroups.com <javascript:>.
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>>.
>
> --
> 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+...@googlegroups.com <javascript:>
``````
``````<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>.

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro...@ufz.de <javascript:> / [http://www.ufz.de](http://www.ufz.de)
``````
``````Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch
``````

Groups “ogs-users” group.
To unsubscribe from this group and stop receiving emails from it, send

Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch