Defining z-dependend pressure boundary condition

bastian.graupner · March 19, 2013, 3:00pm

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left -> the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Mauro_Cacace · March 19, 2013, 3:18pm

Are you using a plane for the interpolation? If yes, try to set a
relevant espilon for the interpolation. If not, try to export a tin
of the boundary.

···

Am 19/03/2013 16:00, schrieb Graupner
Bastian:

Hi everyone,

        I am using OGS 64bit on Windows, ver 5.3.7 JOD .

        I would like to define a
pressure boundary condition on the sites of a cube depending
on the z-coordinate. Therefore I am going to use the
distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

        SURFACE

SURF_Left

$DIS_TYPE

      LINEAR  4                          ; Pa
Den Knoten der Fläche werden linear die Druckdifferenzen
zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is
used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

        Unfortunately the result
of the pressure distribution for the first time step looks a
little strange. I got heterogeneous pressure definitions at
the side and obviously no pressure definition at the corners
of the surface. What should I do to get a proper pressure
distribution as boundary condition?

Best regards,

Bastian

** Bastian
Graupner**

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–
  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 .
For more options, visit .


-- Dr. Mauro Cacace _________________________________________
Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum GFZ (German Research Centre for Geosciences)
- Sect. 4.4: Basin Analysis -
Telegrafenberg Haus C4 Raum 30/1.08
D-14473 Potsdam Germany
_________________________________________
Tel. : +49 (0)331 288 1783
Fax : +49 (0)331 288 1349
e-mail: _________________________________________

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

cacace@gfz-potsdam.de

Marc_Walther · March 19, 2013, 4:06pm

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

···

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

  Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

        I am using OGS 64bit on Windows, ver 5.3.7 JOD .

        I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

        SURFACE SURF_Left

$DIS_TYPE

      LINEAR  4                          ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

          PLY_Left                           -> the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

        Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

** Bastian Graupner**

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

  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](https://groups.google.com/groups/opt_out).

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.


-- Dr. Mauro Cacace _________________________________________ Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum GFZ (German Research Centre for Geosciences) - Sect. 4.4: Basin Analysis - Telegrafenberg Haus C4 Raum 30/1.08 D-14473 Potsdam Germany _________________________________________ Tel. : +49 (0)331 288 1783
Fax : +49 (0)331 288 1349
e-mail: cacace@gfz-potsdam.de
_________________________________________

–

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.

bastian.graupner · March 19, 2013, 4:12pm

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

···

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax : +49 (0)331 288 1349

e-mail:  cacace@gfz-potsdam.de

–

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.

Marc_Walther · March 19, 2013, 4:20pm

Did you try it? I remember having another issue where I had to switch to a surface bc from a polyline and I think I used gradient there… But let me check later, I’m mobile currently…

···

On 19 Mar 2013 17:12, “Graupner Bastian” Bastian.Graupner@ensi.ch wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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.

Marc_Walther2 · March 19, 2013, 5:04pm

Hi Sebastian,

we were both right... in the trunk, there is no gradient option.

But, Jens-Olaf once included that in a version and this is why it
worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD ;)

rf_bc_new.h (7.31 KB)

rf_bc_new.cpp (65.6 KB)

···

On 19.03.2013 17:12, Graupner Bastian
wrote:

Hi

Marc,

        unfortunately GRADIENT is only implemented for

Polylines but not for surfaces.

Thanks,

Bastian

Von:
        []
Marc Walther
Dienstag, 19. März 2013 17:06
Re: [ogs-users] Defining z-dependend
pressure boundary condition

Hej Sebastian,

      I'm not sure, but maybe you can use gradient instead of

linear.

      Cheers,

      Marc

On 19 Mar 2013 16:18, “Mauro Cacace” <cacace@gfz-potsdam.de >
wrote:

            Am 19/03/2013 16:00, schrieb Graupner

Bastian:

Hi

everyone,

                I am using OGS 64bit on Windows, ver 5.3.7 JOD .

                I would like to define a pressure
boundary condition on the sites of a cube depending
on the z-coordinate. Therefore I am going to use the
distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

              LINEAR  4                          ; Pa Den Knoten der
Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

                  PLY_Left

→ the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

                Unfortunately the result of the
pressure distribution for the first time step looks
a little strange. I got heterogeneous pressure
definitions at the side and obviously no pressure
definition at the corners of the surface. What
should I do to get a proper pressure distribution as
boundary condition?

              Best

regards,

Bastian

** Bastian
Graupner**

Dr.

Safety Analyst

                Swiss Federal
Nuclear Safety Inspectorate ENSI

                Secton

Disposal & Analyses

                Industriestrasse

19

                CH-5200

Brugg

                Tel.:

+41 56 460 85 52

                Fax:

+41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

            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](https://groups.google.com/groups/opt_out).

          Are you using a plane for the
interpolation? If yes, try to set a relevant espilon for
the interpolation. If not, try to export a tin of the
boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

        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](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](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 .
For more options, visit .


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

ogs-users@googlegroups.commailto:ogs-users@googlegroups.comIm Auftrag von
Gesendet:
**An:**ogs-users@googlegroups.com
Betreff:

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

http://www.tu-dresden.de/fghhigw

Cmcdermo · March 20, 2013, 8:32am

Hi
Just one tip I found in the past. Make sure if you are defining along a polyline that all the points mentioned as belonging to the polyline in the gli file are specifically addressed under the LINEAR description.
Regards
Chris McDermott

···

On 19/03/2013 15:18, Mauro Cacace wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a
cube depending on the z-coordinate. Therefore I am going to use the
distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Flï¿½che werden linear die Druckdifferenzen
zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left -> the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first
time step looks a little strange. I got heterogeneous pressure
definitions at the side and obviously no pressure definition at the
corners of the surface. What should I do to get a proper pressure
distribution as boundary condition?

Best regards,

Bastian

imap://cacace@imap.gfz-potsdam.de:143/fetch%3EUID%3E/INBOX%3E34860?header=quotebody&part=1.2&filename=image001.png

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

*Bastian Graupner*

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch <mailto:Bastian.Graupner@ensi.ch>__

Swiss Federal Nuclear Safety Inspectorate ENSI » ENSI

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

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

Are you using a plane for the interpolation? If yes, try to set a
relevant espilon for the interpolation. If not, try to export a tin of
the boundary.

--
Dr. Mauro Cacace
_________________________________________

Helmholtz-Zentrum Potsdam
Deutsches GeoForschungsZentrum GFZ
(German Research Centre for Geosciences)
- Sect. 4.4: Basin Analysis -
Telegrafenberg
Haus C4 Raum 30/1.08
D-14473 Potsdam
Germany
_________________________________________

Tel. : +49 (0)331 288 1783
Fax : +49 (0)331 288 1349
e-mail:cacace@gfz-potsdam.de
_________________________________________

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

--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.

jens-olaf_Delfs · March 20, 2013, 8:47am

Hi,
thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE
DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3
nodenr3 value3
. .
. .
. .
#STOP

···

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right... in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison...



Marc



PS: All credits go to JOD ;)





  On 19.03.2013 17:12, Graupner Bastian wrote:

        Hi Marc,

        unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com ] Im Auftrag von Marc Walther
        **Gesendet:** Dienstag, 19. März 2013 17:06

        **An:** ogs-users@googlegroups.com

        **Betreff:** Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

      I'm not sure, but maybe you can use gradient instead of linear.

      Cheers,

      Marc

On 19 Mar 2013 16:18, “Mauro Cacace” <cacace@gfz-potsdam.de > wrote:

            Am 19/03/2013 16:00, schrieb Graupner Bastian:

              Hi everyone,

                I am using OGS 64bit on Windows, ver 5.3.7 JOD .

                I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

                SURFACE SURF_Left

$DIS_TYPE

              LINEAR  4                          ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

                  PLY_Left                           -> the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

                Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

              Best regards,

Bastian

** Bastian Graupner**

Dr.

Safety Analyst

                Swiss Federal Nuclear Safety Inspectorate ENSI

                Secton Disposal & Analyses

                Industriestrasse 19

                CH-5200 Brugg

                Tel.: +41 56 460 85 52

                Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

            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](https://groups.google.com/groups/opt_out).

          Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

        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](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](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](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 [http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

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

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

–

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.

bastian.graupner · March 20, 2013, 11:38am

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

···

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax +49 351 463 425 52

Mobil     +49 178 334 18 19

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

--

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](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](https://groups.google.com/groups/opt_out).

 

</details>

jens-olaf_Delfs · March 20, 2013, 12:14pm

Please go with the debuger and your example through the function

void CBoundaryConditionsGroup::Set(CRFProcess* pcs, int ShiftInNodeVector,
const std::string& this_pv_name)

Search for FiniteElement::GRADIENT and the comment // 6/2012 JOD

I am not sure if current implementation works or not, because there are many brackets.

Thanks

Jens

···

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 09:48

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE

DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3

nodenr3 value3

. .

. .

. .

#STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax       +49 351 463 425 52

Mobil     +49 178 334 18 19

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

–

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.

Marc_Walther2 · March 20, 2013, 12:22pm

the implementation worked at least with JOD’s 5.1.0 version… this is were the files are from. I used it for PRESSURE1 BC in RICHARDS_FLOW at several surfaces.

···

On 20.03.2013 13:14, jens-olaf Delfs wrote:

  Please go with the debuger and your example through the function



  void CBoundaryConditionsGroup::Set(CRFProcess* pcs, int ShiftInNodeVector,

                                     const std::string& this_pv_name)



  Search for FiniteElement::GRADIENT and the comment // 6/2012 JOD

   

  I am not sure if current implementation works or not, because there are many brackets.



  Thanks



  Jens

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

              Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

              Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von: ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com ] Im Auftrag von jens-o laf Delfs
              **Gesendet:** Mittwoch, 20. März 2013 09:48

              **An:** ogs-users@googlegroups.com

              **Betreff:** Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

            thanks for the credits                 . GRADIENT is the right way. Can you switch to the trunk version? I saw that it's already implemented for surfaces, but there is a bug                 . We haven't tested it for the surface yet.

              If you prefer to stick to your version, you can use WWs well-tested



              $DIS_TYPE

                DIRECT   filename.tex



              Bad think is, that you have to know the node numbers and calculate the values for a file:



              nodenr1  value1

              nodenr2  value3

              nodenr3  value3

               .             .

               .             . 

               .             .

              #STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

                    we were both right... in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison...



                    Marc



                    PS: All credits go to JOD ;)

                          On 19.03.2013 17:12, Graupner Bastian wrote:

                              Hi Marc,

                              unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com
[mailto:ogs-users@googlegroups.com ] Im Auftrag von Marc Walther

                              **Gesendet:** Dienstag, 19. März 2013 17:06

                              **An:** ogs-users@googlegroups.com

                              **Betreff:** Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

                            I'm not sure, but maybe you can use gradient instead of linear.

                            Cheers,

                            Marc

                              On 19 Mar 2013 16:18, "Mauro Cacace" <cacace@gfz-potsdam.de                                  > wrote:

                                  Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

                                      I am using OGS 64bit on Windows, ver 5.3.7 JOD .

                                      I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

                                      SURFACE SURF_Left

$DIS_TYPE

                                    LINEAR  4                          ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

                                      The gli-file contains:

#POLYLINE

$NAME

                                        PLY_Left                           -> the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

                                      Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

** Bastian Graupner**

Dr.

Safety Analyst

                                      Swiss Federal Nuclear Safety Inspectorate ENSI

                                      Secton Disposal & Analyses

                                      Industriestrasse 19

CH-5200 Brugg

                                      Tel.:                                             +41 56 460 85 52

                                      Fax:                                             +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

                                  --

                                  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](https://groups.google.com/groups/opt_out).

                                Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

                              --

                              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](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](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](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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax       +49 351 463 425 52

Mobil     +49 178 334 18 19

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

–

                      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](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](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](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 . For more options, visit .


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

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

http://www.tu-dresden.de/fghhigw

bastian.graupner · March 20, 2013, 12:24pm

The code in the trunk version looks very different to the code I got from Marcs rf_bc.cpp. I will check if it works.

Bastian

···

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 09:48

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE

DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3

nodenr3 value3

. .

#STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax +49 351 463 425 52

Mobil     +49 178 334 18 19

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

--

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](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](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](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](https://groups.google.com/groups/opt_out).

 

</details>

jens-olaf_Delfs · March 21, 2013, 4:03pm

Hi,
I tried with a (almost) 5.3.7 version. The GRADIENT BC works well for surfaces. I also found a test case from that time for 3D Richards flow. There must be another issue. Maybe the poinst on the surface are not found, or too many taken

···

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

The code in the trunk version looks very different to the code I got from Marcs rf_bc.cpp. I will check if it works.

Bastian

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 13:15

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Please go with the debuger and your example through the function

void CBoundaryConditionsGroup::Set(CRFProcess* pcs, int ShiftInNodeVector,
                               const std::string& this_pv_name)
Search for FiniteElement::GRADIENT and the comment // 6/2012 JOD

I am not sure if current implementation works or not, because there are many brackets.

Thanks

Jens

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 09:48

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE

DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3

nodenr3 value3

. .

. .

. .

#STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax       +49 351 463 425 52

Mobil     +49 178 334 18 19

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

–

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.

bastian.graupner · March 21, 2013, 4:10pm

Hi Jens,

you are right. It is working if I use it without MPI. I just noticed this. Unfortunately OGS is making strange things if I try it with MPI.

I am working on that.

Thanks a lot for your help and your time.

Best regards

Bastian

···

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

The code in the trunk version looks very different to the code I got from Marcs rf_bc.cpp. I will check if it works.

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 13:15

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Please go with the debuger and your example through the function

void CBoundaryConditionsGroup::Set(CRFProcess* pcs, int ShiftInNodeVector,

                               const std::string& this_pv_name)

Search for FiniteElement::GRADIENT and the comment // 6/2012 JOD

I am not sure if current implementation works or not, because there are many brackets.

Thanks

Jens

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 09:48

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE

DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3

nodenr3 value3

. .

#STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax +49 351 463 425 52

Mobil     +49 178 334 18 19

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

--

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](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](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](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](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](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](https://groups.google.com/groups/opt_out).

 

</details>

jens-olaf_Delfs · March 21, 2013, 4:23pm

Ok, we didn’t consider parallelization that time. Please report the mpi-guys if there is an issue.

···

2013/3/21 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens,

you are right. It is working if I use it without MPI. I just noticed this. Unfortunately OGS is making strange things if I try it with MPI.

I am working on that.

Thanks a lot for your help and your time.

Best regards

Bastian

Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Donnerstag, 21. März 2013 17:03

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

I tried with a (almost) 5.3.7 version. The GRADIENT BC works well for surfaces. I also found a test case from that time for 3D Richards flow. There must be another issue. Maybe the poinst on the surface are not found, or too many taken

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

The code in the trunk version looks very different to the code I got from Marcs rf_bc.cpp. I will check if it works.

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 13:15

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Please go with the debuger and your example through the function

void CBoundaryConditionsGroup::Set(CRFProcess* pcs, int ShiftInNodeVector,
                               const std::string& this_pv_name)
Search for FiniteElement::GRADIENT and the comment // 6/2012 JOD

I am not sure if current implementation works or not, because there are many brackets.

Thanks

Jens

2013/3/20 Graupner Bastian Bastian.Graupner@ensi.ch

Hi Jens-Olaf, hi Marc

Thanks a lot for your help.

Unfortunately I can not use the trunk version because the mpi is currently not working. That’s why I switched to 5.3.07. Marc send me an email with the code for surfaces. Is there also an bug that I should fix? Did you fix the bug in the trunk version already?

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von jens-olaf Delfs

Gesendet: Mittwoch, 20. März 2013 09:48

An: ogs-users@googlegroups.com

Betreff: Re: AW: [ogs-users] Defining z-dependend pressure boundary condition

Hi,

thanks for the credits . GRADIENT is the right way. Can you switch to the trunk version? I saw that it’s already implemented for surfaces, but there is a bug . We haven’t tested it for the surface yet.

If you prefer to stick to your version, you can use WWs well-tested

$DIS_TYPE

DIRECT filename.tex

Bad think is, that you have to know the node numbers and calculate the values for a file:

nodenr1 value1

nodenr2 value3

nodenr3 value3

. .

. .

. .

#STOP

2013/3/19 Marc Walther marc.walther@tu-dresden.de

Hi Sebastian,

we were both right… in the trunk, there is no gradient option. But, Jens-Olaf once included that in a version and this is why it worked for me! Attached is rf_bc_new.* for comparison…

Marc

PS: All credits go to JOD

On 19.03.2013 17:12, Graupner Bastian wrote:

Hi Marc,

unfortunately GRADIENT is only implemented for Polylines but not for surfaces.

Thanks,

Bastian

Von:
ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Marc Walther

Gesendet: Dienstag, 19. März 2013 17:06

An: ogs-users@googlegroups.com

Betreff: Re: [ogs-users] Defining z-dependend pressure boundary condition

Hej Sebastian,

I’m not sure, but maybe you can use gradient instead of linear.

Cheers,

Marc

On 19 Mar 2013 16:18, “Mauro Cacace” cacace@gfz-potsdam.de wrote:

Am 19/03/2013 16:00, schrieb Graupner Bastian:

Hi everyone,

I am using OGS 64bit on Windows, ver 5.3.7 JOD .

I would like to define a pressure boundary condition on the sites of a cube depending on the z-coordinate. Therefore I am going to use the distribution type LINEAR.

Bc-file:

#BOUNDARY_CONDITION

$PCS_TYPE

LIQUID_FLOW

$PRIMARY_VARIABLE

PRESSURE1

$GEO_TYPE

SURFACE SURF_Left

$DIS_TYPE

LINEAR 4 ; Pa Den Knoten der Fläche werden linear die Druckdifferenzen zugeordnet

3 978480.0

4 978480.0

8 821520.0

7 821520.0

#STOP

The gli-file contains:

#POLYLINE

$NAME

PLY_Left → the PLY_Left is used to define the SURF_Left

$POINTS

3

4

8

7

3

#SURFACE

$NAME

SURF_Left

$POLYLINES

PLY_Left

$MAT_GROUP

0

$TYPE

0

Unfortunately the result of the pressure distribution for the first time step looks a little strange. I got heterogeneous pressure definitions at the side and obviously no pressure definition at the corners of the surface. What should I do to get a proper pressure distribution as boundary condition?

Best regards,

Bastian

Bastian Graupner

Dr.

Safety Analyst

Swiss Federal Nuclear Safety Inspectorate ENSI

Secton Disposal & Analyses

Industriestrasse 19

CH-5200 Brugg

Tel.: +41 56 460 85 52

Fax: +41 56 460 84 99

Bastian.Graupner@ensi.ch

http://www.ensi.ch/en

–

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.

Are you using a plane for the interpolation? If yes, try to set a relevant espilon for the interpolation. If not, try to export a tin of the boundary.

--

Dr. Mauro Cacace

_________________________________________

Helmholtz-Zentrum Potsdam

Deutsches GeoForschungsZentrum GFZ

(German Research Centre for Geosciences)

- Sect. 4.4: Basin Analysis -

Telegrafenberg

Haus C4 Raum 30/1.08

D-14473 Potsdam

Germany

_________________________________________

Tel.  :  +49  (0)331  288  1783

Fax   :  +49  (0)331  288  1349

e-mail:  cacace@gfz-potsdam.de

_________________________________________

–

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

[http://www.tu-dresden.de/fghhigw](http://www.tu-dresden.de/fghhigw)

Tel       +49 351 463 425 56

Fax       +49 351 463 425 52

Mobil     +49 178 334 18 19

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

–

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.