AW: Bugs fixed --> there are still some

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:
  1. 2D problem doesn't work if the problem in y-z coordinate system
  2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy
  3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today's commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

--
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.
For more options, visit https://groups.google.com/groups/opt_out.

I have forgot one to be ask:
  If I define q (m/s), I use constant_neumann in the .st file for the groundwater flow model, should I use constant in the .st file as FORSYTH case in the TRUNK.

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Shao, Hua
Gesendet: Donnerstag, 16. Januar 2014 11:57
An: ogs-users@googlegroups.com
Betreff: AW: [ogs-users] Bugs fixed --> there are still some

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:
  1. 2D problem doesn't work if the problem in y-z coordinate system
  2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy
  3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today's commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

--
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.
For more options, visit https://groups.google.com/groups/opt_out.

--
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.
For more options, visit https://groups.google.com/groups/opt_out.

Dear Hua,

as I know from Wenqing, output is in X and Y for H2M, even if you use X-Z-section. Maybe it's the same for Richards flow.

···

------------------------------------
There is another shortcoming in OGS5 for output 2D deformation result. Here are tips:
*Telplot: Whichever a coordinate system is chosen for 2D plane strain deformation analysis, i.e. (x, y, 0) or (x, 0, z), the output for this problem always uses x and y coordinates. In that way, the results of a 2D deformation problems are ux, uy; strain_xx, strain_yy, strain_xy; stress_xx, , stress_yy, , stress_zz (off-plane direction), stress_xy.
*PVD of Paraview: Displacements are stick to coordinates in the manner: (x, y, 0) makes output of ux and uy while (x, 0, z) makes output of ux and uz. Strain and stress are 'strain_xx, strain_yy, strain_xy; stress_xx, , stress_yy, , stress_zz (off-plane direction), stress_xy'.

Therefore, I suggest to use (x, 0, z) coordinates for your 2D example, and check the result of uz according to the above description.

--------------------------------

Cheers

Katharina

Katharina Benisch
Diplom-Geologin
Abteilung Geohydromodellierung
Institut f�r Geowissenschaften, Raum 351
Christian-Albrechts-Universit�t Kiel
Ludewig-Meyn-Strasse 10, 24118 Kiel
Tel.: 0431-880 1347
Fax.: 0431-880 6706
Web: www.gpi.uni-kiel.de

Am 16.01.2014 11:56, schrieb Shao, Hua:

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:
  1. 2D problem doesn't work if the problem in y-z coordinate system
  2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy
  3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Urspr�ngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today's commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

Thanks Katharina. I have accepted this already. Hua

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Katharina Benisch
Gesendet: Donnerstag, 16. Januar 2014 12:03
An: ogs-users@googlegroups.com
Betreff: Re: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

as I know from Wenqing, output is in X and Y for H2M, even if you use X-Z-section. Maybe it's the same for Richards flow.

------------------------------------
There is another shortcoming in OGS5 for output 2D deformation result.
Here are tips:
*Telplot: Whichever a coordinate system is chosen for 2D plane strain deformation analysis, i.e. (x, y, 0) or (x, 0, z), the output for this problem always uses x and y coordinates. In that way, the results of a 2D deformation problems are ux, uy; strain_xx, strain_yy, strain_xy; stress_xx, , stress_yy, , stress_zz (off-plane direction), stress_xy.
*PVD of Paraview: Displacements are stick to coordinates in the manner:
(x, y, 0) makes output of ux and uy while (x, 0, z) makes output of ux and uz. Strain and stress are 'strain_xx, strain_yy, strain_xy; stress_xx, , stress_yy, , stress_zz (off-plane direction), stress_xy'.

Therefore, I suggest to use (x, 0, z) coordinates for your 2D example, and check the result of uz according to the above description.

--------------------------------

Cheers

Katharina

Katharina Benisch
Diplom-Geologin
Abteilung Geohydromodellierung
Institut für Geowissenschaften, Raum 351 Christian-Albrechts-Universität Kiel Ludewig-Meyn-Strasse 10, 24118 Kiel
Tel.: 0431-880 1347
Fax.: 0431-880 6706
Web: www.gpi.uni-kiel.de

Am 16.01.2014 11:56, schrieb Shao, Hua:

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:
  1. 2D problem doesn't work if the problem in y-z coordinate system
  2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy
  3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im
Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today's commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

Dear Hua,

Could you upload the setup to your branch and give a link, that we can take a look?

I don’t understand what you mean by 1. “it doesn’t work” - can you give some details?

  1. For 2D, I think OGS will always write results in the first two components of the vectors. You can easily transform the vector in ParaView or TecPlot.

  2. The setup would help to understand why this happens :slight_smile:

Best,

Marc

···

Am 16.01.2014 11:56 schrieb “Shao, Hua” Hua.Shao@bgr.de:

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:

    1. 2D problem doesn't work if the problem in y-z coordinate system

    2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy

    3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the ‘RICHARD’s flow’?

Hua

-----Ursprüngliche Nachricht-----

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang

Gesendet: Dienstag, 14. Januar 2014 12:55

An: ogs-users@googlegroups.com

Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today’s commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

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.

For more options, visit https://groups.google.com/groups/opt_out.

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.

For more options, visit https://groups.google.com/groups/opt_out.

Hi all,
In the PETSC benchmark/Translay2d folder is another application with
reactive transport & Richards flow for X-Z coordinates. What I
experienced is, that output of "velocities" into VTK files may need some
rework. Josh and me did something in the routine, but I am not sure it
works for all different flow modules. And of course, as mentioned
earlier, for 2D problems only two velocities are defined.
AND: If I am not mistaken - at least for VTK - output are DARCY fluxes
and not velocities.
Greetings, Georg

···

On 01/16/2014 12:25 PM, Marc Walther wrote:

Dear Hua,

Could you upload the setup to your branch and give a link, that we can
take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

    Dear all

    We are currently using Richards-model for one project. Based on
    some Benchmark-example, e.g. Foryth, we got a very quick start.
    However, we have also problems as:
            1. 2D problem doesn't work if the problem in y-z
    coordinate system
            2. Definition of coordinate system is quite confuse, i.e.
    2D problem in x-z section, output of velocity is ux, uy
            3. The calculation of velocity is with great uncertainty:
    in some test cases, we got even vx=vy (???) in the output file.

    Does anyone have an idea or experience on the 'RICHARD's flow'?

    Hua

    -----Urspr�ngliche Nachricht-----
    Von: ogs-users@googlegroups.com
    <mailto:ogs-users@googlegroups.com>
    [mailto:ogs-users@googlegroups.com
    <mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
    Gesendet: Dienstag, 14. Januar 2014 12:55
    An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
    Betreff: [ogs-users] Bugs fixed

    Dear all,

    If you solve THM coupling problem with deactivated domains, please
    patch your source code with my today's commits , which contain the
    bug fixing for the boundary condition assignment.

    Wenqing

    --
    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%2Bunsubscribe@googlegroups.com>.
    For more options, visit https://groups.google.com/groups/opt_out\.

    --
    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%2Bunsubscribe@googlegroups.com>.
    For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email: georg.kosakowski@psi.ch h323:+41563104743
homepage: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

Thank you all for lots tips. I don’t know where is my branch, I attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as boundary conditions. The calculation of velocity, at least direction looks properly.

The xz coordinate system runs well, however, the velocity changes direction, especially around tunnel in the late time.

Hua

xz.7z (30.8 KB)

xy.7z (30.8 KB)

···

Am 16.01.2014 11:56 schrieb “Shao, Hua” Hua.Shao@bgr.de:

Dear all

We are currently using Richards-model for one project. Based on some Benchmark-example, e.g. Foryth, we got a very quick start. However, we have also problems as:

    1. 2D problem doesn't work if the problem in y-z coordinate system

    2. Definition of coordinate system is quite confuse, i.e. 2D problem in x-z section, output of velocity is ux, uy

    3. The calculation of velocity is with great uncertainty: in some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the ‘RICHARD’s flow’?

Hua

-----Ursprüngliche Nachricht-----

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang

Gesendet: Dienstag, 14. Januar 2014 12:55

An: ogs-users@googlegroups.com

Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please patch your source code with my today’s commits , which contain the bug fixing for the boundary condition assignment.

Wenqing

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.

For more options, visit
https://groups.google.com/groups/opt_out
.

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.

For more options, visit
https://groups.google.com/groups/opt_out
.

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.

For more options, visit https://groups.google.com/groups/opt_out.

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

···

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don't know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn't allow me to used gradient as boundary
conditions. The calculation of velocity, at least direction looks
properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we can
take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However, we
have also problems as:
        1. 2D problem doesn't work if the problem in y-z coordinate system
        2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
        3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Urspr�ngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the bug
fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email: georg.kosakowski@psi.ch h323:+41563104743
homepage: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

Hi Hua,

yes, I remember Georg having some problems with high gradients of permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE

LINEAR [number of points defined below]

[point_ID1] [PrimVar_value] [name of point with ID1]

[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z **$NAME** [name of point with ID1]

Maybe this helps you with applying a gradient.

Best,

Marc
···

mailto:ogs-users@googlegroups.com

ogs-users@googlegroups.com

Hua.Shao@bgr.de
mailto:Hua.Shao@bgr.de

ogs-users@googlegroups.commailto:ogs-users@googlegroups.com
mailto:ogs-users@googlegroups.com
mailto:ogs-users@googlegroups.com

ogs-users@googlegroups.commailto:ogs-users@googlegroups.com

ogs-users+unsubscribe@googlegroups.com
mailto:ogs-users%2Bunsubscribe@googlegroups.com
https://groups.google.com/groups/opt_out
https://groups.google.com/groups/opt_out

ogs-users+unsubscribe@googlegroups.com
mailto:ogs-users%2Bunsubscribe@googlegroups.com
https://groups.google.com/groups/opt_out
https://groups.google.com/groups/opt_out

ogs-users+unsubscribe@googlegroups.com
mailto:ogs-users+unsubscribe@googlegroups.com
https://groups.google.com/groups/opt_out

ogs-users+unsubscribe@googlegroups.com
https://groups.google.com/groups/opt_out
Email:georg.kosakowski@psi.ch
http://people.web.psi.ch/kosakowski/

ogs-users+unsubscribe@googlegroups.com
https://groups.google.com/groups/opt_out

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
$PCS_TYPE
  RICHARDS_FLOW
$PRIMARY_VARIABLE
  PRESSURE1
$GEO_TYPE
  POLYLINE LEFT_BOUNDARY
$DIS_TYPE
  GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

···

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email: georg.kosakowski@psi.ch h323:+41563104743
homepage: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

you can also try a distribution type FUNCTION, which f(x,y,z)=ax+by+cz. Example is

  $DIS_TYPE
   FUNCTION
   0+-11453*y

best,
nori

···

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
Dr.-Ing. Norihiro Watanabe
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

This is cool; never new, it existed!
Does it work with any PrimVar?

···

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
- - - - - - - - - - - - - - - - - -

Dipl.-Hydrol. Marc Walther
Wissenschaftlicher Mitarbeiter

Technische Universität Dresden
Fakultät Umweltwissenschaften
Fachrichtung Hydrowissenschaften
Institut für Grundwasserwirtschaft
01062 Dresden

Besucheradresse
Bergstr. 66, 01069 Dresden
Neubau Chemie 2. Bauabschnitt, CHE-BA2, Raum E59

Tel +49 351 463 425 56
Fax +49 351 463 425 52
Mobil +49 178 334 18 19

- - - - - - - - - - - - - - - - - -

i am not sure. please ask Wenqing as he has implemented it some time ago

···

On 01/16/2014 03:59 PM, Marc Walther wrote:

This is cool; never new, it existed!
Does it work with any PrimVar?

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
Dr.-Ing. Norihiro Watanabe
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

Yes, you can do like this but keep in mind that only + and * are operators, and x, y, z denote coordinates. For BC, the syntax is that in Nori's email. For IC: it takes the form:

$DIS_TYPE
   FUNCTION
  $GEO_TYPE
   SUBDOMAIN
   2
    0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the constant to be input.

Wenqing

···

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

It should work for any primary variable.

···

On 16.01.2014 16:30, Norihiro Watanabe wrote:

i am not sure. please ask Wenqing as he has implemented it some time ago

On 01/16/2014 03:59 PM, Marc Walther wrote:

This is cool; never new, it existed!
Does it work with any PrimVar?

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

typo in 'For IC':
$DIS_TYPE
   FUNCTION
  $GEO_TYPE
   SUBDOMAIN
   2
    0 1.693567e5+4000*x+-11559.45*y
    2 1.693567e5+4000*x+-11559.45*y

The integer before the expression is the material index of a domain.

···

On 16.01.2014 16:52, Wenqing Wang wrote:

Yes, you can do like this but keep in mind that only + and * are operators, and x, y, z denote coordinates. For BC, the syntax is that in Nori's email. For IC: it takes the form:

$DIS_TYPE
  FUNCTION
$GEO_TYPE
  SUBDOMAIN
  2
   0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the constant to be input.

Wenqing

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for increase
of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below]
[point_ID1] [PrimVar_value] [name of point with ID1]
[point_ID2] [PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are only
defined for z-direction. I remember I looked this up in the code,
because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might get
"wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least direction
looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still some

Dear Hua,

Could you upload the setup to your branch and give a link, that we
can take a look?

I don't understand what you mean by 1. "it doesn't work" - can you
give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on some
Benchmark-example, e.g. Foryth, we got a very quick start. However,
we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system
2. Definition of coordinate system is quite confuse, i.e. 2D
problem in x-z section, output of velocity is ux, uy
3. The calculation of velocity is with great uncertainty: in
some test cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com <mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains, please
patch your source code with my today's commits , which contain the
bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

My BC definition using GRADIENT is the same as Georg, as said that it is not correct if I use it in xy section.

Both Version 5.4 and 5.5 doesn't work if I used the function described by Marc. I attached the .BC in this mail.

We have not understood the function described by Nori and Wenqing. What are a, b, c ....?

If I need only a distribution
  at y=400 300000 Pa
  at y=700 100000 Pa

How can I define them?

Best
Hua

bedrichov.bc (0 Bytes)

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Donnerstag, 16. Januar 2014 16:57
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: [ogs-users] Bugs fixed --> there are still some

typo in 'For IC':
$DIS_TYPE
   FUNCTION
  $GEO_TYPE
   SUBDOMAIN
   2
    0 1.693567e5+4000*x+-11559.45*y
    2 1.693567e5+4000*x+-11559.45*y

The integer before the expression is the material index of a domain.

On 16.01.2014 16:52, Wenqing Wang wrote:

Yes, you can do like this but keep in mind that only + and * are
operators, and x, y, z denote coordinates. For BC, the syntax is that
in Nori's email. For IC: it takes the form:

$DIS_TYPE
  FUNCTION
$GEO_TYPE
  SUBDOMAIN
  2
   0 1.693567e5+4000*x+-11559.45*y
   0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the
constant to be input.

Wenqing

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

$DIS_TYPE
  FUNCTION
  0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for
increase of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
  $PCS_TYPE
   RICHARDS_FLOW
  $PRIMARY_VARIABLE
   PRESSURE1
  $GEO_TYPE
   POLYLINE LEFT_BOUNDARY
  $DIS_TYPE
   GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below] [point_ID1] [PrimVar_value]
[name of point with ID1] [point_ID2] [PrimVar_value] [name of point
with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are
only defined for z-direction. I remember I looked this up in the
code, because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might
get "wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least
direction looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still
some

Dear Hua,

Could you upload the setup to your branch and give a link, that
we can take a look?

I don't understand what you mean by 1. "it doesn't work" - can
you give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on
some Benchmark-example, e.g. Foryth, we got a very quick start.
However, we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system 2. Definition of coordinate system is quite confuse, i.e.
2D problem in x-z section, output of velocity is ux, uy 3. The
calculation of velocity is with great uncertainty: in some test
cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains,
please patch your source code with my today's commits , which
contain the bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

You can realize it by defining two polyline paralleling to x-z plane at y= 400, ad y = 700, respectively, and assign BC at these lines.

···

On 17.01.2014 09:47, Shao, Hua wrote:

My BC definition using GRADIENT is the same as Georg, as said that it is not correct if I use it in xy section.

Both Version 5.4 and 5.5 doesn't work if I used the function described by Marc. I attached the .BC in this mail.

We have not understood the function described by Nori and Wenqing. What are a, b, c ....?

If I need only a distribution
  at y=400 300000 Pa
  at y=700 100000 Pa

How can I define them?

Best
Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Donnerstag, 16. Januar 2014 16:57
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: [ogs-users] Bugs fixed --> there are still some

typo in 'For IC':
$DIS_TYPE
    FUNCTION
   $GEO_TYPE
    SUBDOMAIN
    2
     0 1.693567e5+4000*x+-11559.45*y
     2 1.693567e5+4000*x+-11559.45*y

The integer before the expression is the material index of a domain.

On 16.01.2014 16:52, Wenqing Wang wrote:

Yes, you can do like this but keep in mind that only + and * are
operators, and x, y, z denote coordinates. For BC, the syntax is that
in Nori's email. For IC: it takes the form:

$DIS_TYPE
   FUNCTION
  $GEO_TYPE
   SUBDOMAIN
   2
    0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the
constant to be input.

Wenqing

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

  $DIS_TYPE
   FUNCTION
   0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for
increase of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
   $PCS_TYPE
    RICHARDS_FLOW
   $PRIMARY_VARIABLE
    PRESSURE1
   $GEO_TYPE
    POLYLINE LEFT_BOUNDARY
   $DIS_TYPE
    GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below] [point_ID1] [PrimVar_value]
[name of point with ID1] [point_ID2] [PrimVar_value] [name of point
with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are
only defined for z-direction. I remember I looked this up in the
code, because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field might
get "wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least
direction looks properly.

The xz coordinate system runs well, however, the velocity changes
direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still
some

Dear Hua,

Could you upload the setup to your branch and give a link, that
we can take a look?

I don't understand what you mean by 1. "it doesn't work" - can
you give some details?

2. For 2D, I think OGS will always write results in the first two
components of the vectors. You can easily transform the vector in
ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on
some Benchmark-example, e.g. Foryth, we got a very quick start.
However, we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system 2. Definition of coordinate system is quite confuse, i.e.
2D problem in x-z section, output of velocity is ux, uy 3. The
calculation of velocity is with great uncertainty: in some test
cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains,
please patch your source code with my today's commits , which
contain the bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

Hallo Wenqing
We don't understand what you mean. My model is in xy-plane. If I define two parallel polylines to xz-plane, what will happen for the nodes in between?
Hua

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Freitag, 17. Januar 2014 10:00
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: AW: [ogs-users] Bugs fixed --> there are still some

You can realize it by defining two polyline paralleling to x-z plane at y= 400, ad y = 700, respectively, and assign BC at these lines.

On 17.01.2014 09:47, Shao, Hua wrote:

My BC definition using GRADIENT is the same as Georg, as said that it is not correct if I use it in xy section.

Both Version 5.4 and 5.5 doesn't work if I used the function described by Marc. I attached the .BC in this mail.

We have not understood the function described by Nori and Wenqing. What are a, b, c ....?

If I need only a distribution
  at y=400 300000 Pa
  at y=700 100000 Pa

How can I define them?

Best
Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im
Auftrag von Wenqing Wang
Gesendet: Donnerstag, 16. Januar 2014 16:57
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: [ogs-users] Bugs fixed --> there are still some

typo in 'For IC':
$DIS_TYPE
    FUNCTION
   $GEO_TYPE
    SUBDOMAIN
    2
     0 1.693567e5+4000*x+-11559.45*y
     2 1.693567e5+4000*x+-11559.45*y

The integer before the expression is the material index of a domain.

On 16.01.2014 16:52, Wenqing Wang wrote:

Yes, you can do like this but keep in mind that only + and * are
operators, and x, y, z denote coordinates. For BC, the syntax is that
in Nori's email. For IC: it takes the form:

$DIS_TYPE
   FUNCTION
  $GEO_TYPE
   SUBDOMAIN
   2
    0 1.693567e5+4000*x+-11559.45*y
    0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the
constant to be input.

Wenqing

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

  $DIS_TYPE
   FUNCTION
   0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for
increase of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
   $PCS_TYPE
    RICHARDS_FLOW
   $PRIMARY_VARIABLE
    PRESSURE1
   $GEO_TYPE
    POLYLINE LEFT_BOUNDARY
   $DIS_TYPE
    GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below] [point_ID1]
[PrimVar_value] [name of point with ID1] [point_ID2]
[PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are
only defined for z-direction. I remember I looked this up in the
code, because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field
might get "wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least
direction looks properly.

The xz coordinate system runs well, however, the velocity
changes direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still
some

Dear Hua,

Could you upload the setup to your branch and give a link, that
we can take a look?

I don't understand what you mean by 1. "it doesn't work" - can
you give some details?

2. For 2D, I think OGS will always write results in the first
two components of the vectors. You can easily transform the
vector in ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on
some Benchmark-example, e.g. Foryth, we got a very quick start.
However, we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system 2. Definition of coordinate system is quite confuse, i.e.
2D problem in x-z section, output of velocity is ux, uy 3. The
calculation of velocity is with great uncertainty: in some test
cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains,
please patch your source code with my today's commits , which
contain the bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

Then, just lines paralleling to x -axis.

···

On 17.01.2014 10:13, Shao, Hua wrote:

Hallo Wenqing
We don't understand what you mean. My model is in xy-plane. If I define two parallel polylines to xz-plane, what will happen for the nodes in between?
Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Freitag, 17. Januar 2014 10:00
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: AW: [ogs-users] Bugs fixed --> there are still some

You can realize it by defining two polyline paralleling to x-z plane at y= 400, ad y = 700, respectively, and assign BC at these lines.

On 17.01.2014 09:47, Shao, Hua wrote:

My BC definition using GRADIENT is the same as Georg, as said that it is not correct if I use it in xy section.

Both Version 5.4 and 5.5 doesn't work if I used the function described by Marc. I attached the .BC in this mail.

We have not understood the function described by Nori and Wenqing. What are a, b, c ....?

If I need only a distribution
  at y=400 300000 Pa
  at y=700 100000 Pa

How can I define them?

Best
Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im
Auftrag von Wenqing Wang
Gesendet: Donnerstag, 16. Januar 2014 16:57
An: ogs-users@googlegroups.com
Betreff: Re: AW: AW: [ogs-users] Bugs fixed --> there are still some

typo in 'For IC':
$DIS_TYPE
     FUNCTION
    $GEO_TYPE
     SUBDOMAIN
     2
      0 1.693567e5+4000*x+-11559.45*y
      2 1.693567e5+4000*x+-11559.45*y

The integer before the expression is the material index of a domain.

On 16.01.2014 16:52, Wenqing Wang wrote:

Yes, you can do like this but keep in mind that only + and * are
operators, and x, y, z denote coordinates. For BC, the syntax is that
in Nori's email. For IC: it takes the form:

$DIS_TYPE
    FUNCTION
   $GEO_TYPE
    SUBDOMAIN
    2
     0 1.693567e5+4000*x+-11559.45*y

The expression should be a+b*x+c*y+d*z, where a, b, c, and d are the
constant to be input.

Wenqing

On 16.01.2014 15:11, Norihiro Watanabe wrote:

you can also try a distribution type FUNCTION, which
f(x,y,z)=ax+by+cz. Example is

   $DIS_TYPE
    FUNCTION
    0+-11453*y

best,
nori

On 01/16/2014 03:07 PM, Georg Kosakowski wrote:

If you only need a gradient in Z and would like to account for
increase of pressure with depth ...you can try this:

#BOUNDARY_CONDITION
    $PCS_TYPE
     RICHARDS_FLOW
    $PRIMARY_VARIABLE
     PRESSURE1
    $GEO_TYPE
     POLYLINE LEFT_BOUNDARY
    $DIS_TYPE
     GRADIENT 0.0 0.0 9.81e+3 ; at z=0 m we have Pressure =0 -> fully
saturated ...saturation is decreasing with increasing z

GRADIENT [reference depth] [reference depth value] [Gradient]

Greetings, Georg

On 01/16/2014 02:50 PM, Marc Walther wrote:

Hi Hua,

yes, I remember Georg having some problems with high gradients of
permeability some time ago.

Concerning the BC, you might try to use this (.bc):

$DIS_TYPE
LINEAR [number of points defined below] [point_ID1]
[PrimVar_value] [name of point with ID1] [point_ID2]
[PrimVar_value] [name of point with ID2]

for this, you will have to give names to the points in the .gli like:

[point_ID1] x y z *$NAME *[name of point with ID1]

Maybe this helps you with applying a gradient.

Best,
Marc

On Do 16 Jan 2014 14:24:30 CET, Georg Kosakowski wrote:

Hoi Hua,
you are correct regarding the boundary conditions. Gradients are
only defined for z-direction. I remember I looked this up in the
code, because I wondered how to define the direction of gradients!
I also found, that it is necessary to change a bit the numerical
settings, as otherwise the solution for the pressure field is not
accurate enough in low-permeability media and the flow field
might get "wrong".
Try setting the solver tolerance to 1e-20 or even smaller ...
Greetings, Georg

On 01/16/2014 02:18 PM, Shao, Hua wrote:

Thank you all for lots tips. I don’t know where is my branch, I
attached my two simplified cases in the zip-form.

The xy coordinate system doesn’t allow me to used gradient as
boundary conditions. The calculation of velocity, at least
direction looks properly.

The xz coordinate system runs well, however, the velocity
changes direction, especially around tunnel in the late time.

Hua

*Von:*ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com]
*Im Auftrag von *Marc Walther
*Gesendet:* Donnerstag, 16. Januar 2014 12:26
*An:* ogs-users@googlegroups.com
*Betreff:* Re: AW: [ogs-users] Bugs fixed --> there are still
some

Dear Hua,

Could you upload the setup to your branch and give a link, that
we can take a look?

I don't understand what you mean by 1. "it doesn't work" - can
you give some details?

2. For 2D, I think OGS will always write results in the first
two components of the vectors. You can easily transform the
vector in ParaView or TecPlot.

3. The setup would help to understand why this happens :slight_smile:

Best,
Marc

Am 16.01.2014 11:56 schrieb "Shao, Hua" <Hua.Shao@bgr.de
<mailto:Hua.Shao@bgr.de>>:

Dear all

We are currently using Richards-model for one project. Based on
some Benchmark-example, e.g. Foryth, we got a very quick start.
However, we have also problems as:
1. 2D problem doesn't work if the problem in y-z coordinate
system 2. Definition of coordinate system is quite confuse, i.e.
2D problem in x-z section, output of velocity is ux, uy 3. The
calculation of velocity is with great uncertainty: in some test
cases, we got even vx=vy (???) in the output file.

Does anyone have an idea or experience on the 'RICHARD's flow'?

Hua

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
[mailto:ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 14. Januar 2014 12:55
An: ogs-users@googlegroups.com
<mailto:ogs-users@googlegroups.com>
Betreff: [ogs-users] Bugs fixed

Dear all,

If you solve THM coupling problem with deactivated domains,
please patch your source code with my today's commits , which
contain the bug fixing for the boundary condition assignment.

Wenqing

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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%2Bunsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out&gt;\.

--
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/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management Paul
Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Fax: +41-56-3102821
Email:georg.kosakowski@psi.ch h323:+41563104743
homepage:Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.

--
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.
For more options, visit https://groups.google.com/groups/opt_out\.