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

GRADIENT 9810.00

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

Jhonny Encalada

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
  GRADIENT 9810.00
#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.

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

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

GRADIENT 9810.00

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

Jhonny Encalada

‚Äď

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

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
      GRADIENT 9810.00
    #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.

    Jhonny Encalada
    -- 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
    <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+unsubscribe@googlegroups.com <mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.

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

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

GRADIENT 9810.00

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

Jhonny Encalada

‚Äď

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

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

GRADIENT 9810.00

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

Jhonny Encalada

‚Äď

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

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
>> GRADIENT 9810.00
>> #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.
>>
>> Jhonny Encalada
>> --
>> 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 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
> <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+unsubscribe@googlegroups.com
<mailto:ogs-users+unsubscribe@googlegroups.com>.
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.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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

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!

Jhonny Encalada

···

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
      GRADIENT   9810.00
    #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.

    Jhonny Encalada
    --
    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)
    <[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:>.
For more options, visit [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

mailto:ogs-users+...@googlegroups.com.
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 / 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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

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

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
      GRADIENT   9810.00
    #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.

    Jhonny Encalada
    --
    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)
    <[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:>.
For more options, visit [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

mailto:ogs-users+...@googlegroups.com.
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 / 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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

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
     > >> GRADIENT 9810.00
     > >> #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.
     > >>
     > >> Jhonny Encalada
     > >> --
     > >> 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
    <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
    <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:>
     > <mailto:ogs-users+...@googlegroups.com <javascript:>>.
     > For more options, visit 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
    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
¬†¬†¬†¬†Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

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

--
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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

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
 > >>           GRADIENT   9810.00
 > >>         #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.
 > >>
 > >>         Jhonny Encalada
 > >>         --
 > >>         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) <[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>
 > >>         <[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:>.
 > >     For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout)
<[https://groups.google.com/d/optout](https://groups.google.com/d/optout)>
 > >     <[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:>
 > <mailto:ogs-users+...@googlegroups.com <javascript:>>.
 > For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout)
<[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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann

‚Äď
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

mailto:ogs-users+...@googlegroups.com.
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 / 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
Administrativer Gesch√§ftsf√ľhrer: Dr. Heike Gra√ümann