input files for liquid

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

CemBench-conservative_liquid.ogs (184 KB)

CemBench-conservative_groundwater.ogs (184 KB)

路路路

--
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: http://people.web.psi.ch/kosakowski/

Hi Georg,

try using this num:

$OVERALL_COUPLING
2 5
#NUMERICS
  $PCS_TYPE
   LIQUID_FLOW
  $ELE_UPWINDING
   0.5
  $ELE_MASS_LUMPING
   1
  $LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
   2 1 1e-12 500 1 100 4
  $NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
  PICARD ERNORM 100 0.0 1e-2
  $ELE_GAUSS_POINTS
   2
  $COUPLING_CONTROL
  ERNORM 1e-2
#NUMERICS
  $PCS_TYPE
   MASS_TRANSPORT
  $ELE_GAUSS_POINTS
   2
  $LINEAR_SOLVER
   2 1 1e-12 500 1 100 4
  $NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
  PICARD ERNORM 100 0.0 1e-2
  $COUPLING_CONTROL
  ERNORM 1e-2
#STOP

路路路

On 04.02.2013 13:25, Georg Kosakowski wrote:

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

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

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

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

Hoi Marc,
thanks and good try, but the wine is not yet yours :wink:
The coupling is not really necessary and (with or without coupling) the Darcy fluxes are still much to high. They will introduce some dispersion and too much advective transport.

Greetings, Georg

路路路

On 02/04/2013 01:47 PM, Marc Walther wrote:

Hi Georg,

try using this num:

$OVERALL_COUPLING
2 5
#NUMERICS
$PCS_TYPE
聽聽LIQUID_FLOW
$ELE_UPWINDING
聽聽0.5
$ELE_MASS_LUMPING
聽聽1
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$ELE_GAUSS_POINTS
聽聽2
$COUPLING_CONTROL
ERNORM 1e-2
#NUMERICS
$PCS_TYPE
聽聽MASS_TRANSPORT
$ELE_GAUSS_POINTS
聽聽2
$LINEAR_SOLVER
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$COUPLING_CONTROL
ERNORM 1e-2
#STOP

On 04.02.2013 13:25, Georg Kosakowski wrote:

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

--
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: http://people.web.psi.ch/kosakowski/

Hoi,
the num file below seems to work. But can someone explain me why the former numerical settings do not work with LIQUID flow?
In any case thanks to Marc and I hope Pinot Noir (Blauburgunder) is ok.
Greetings, Georg

GeoSys-NUM: Numerical Parameter ----------------------------------------
#NUMERICS
聽聽$PCS_TYPE
聽聽聽LIQUID_FLOW
聽聽$ELE_GAUSS_POINTS
聽聽聽3
聽聽$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
; 2 6 1.e-12 200000 1.0 1 2
聽聽聽2 1 1.e-12 1000 1.0 100 4
#NUMERICS
聽聽$PCS_TYPE
聽聽聽MASS_TRANSPORT
聽聽$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
聽聽聽2 6 1.e-014 2000 1.0 1 2
聽聽$ELE_GAUSS_POINTS
聽聽聽3
#STOP

路路路

On 02/04/2013 02:11 PM, Georg Kosakowski wrote:

Hoi Marc,
thanks and good try, but the wine is not yet yours :wink:
The coupling is not really necessary and (with or without coupling) the Darcy fluxes are still much to high. They will introduce some dispersion and too much advective transport.

Greetings, Georg

On 02/04/2013 01:47 PM, Marc Walther wrote:

Hi Georg,

try using this num:

$OVERALL_COUPLING
2 5
#NUMERICS
$PCS_TYPE
聽聽LIQUID_FLOW
$ELE_UPWINDING
聽聽0.5
$ELE_MASS_LUMPING
聽聽1
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$ELE_GAUSS_POINTS
聽聽2
$COUPLING_CONTROL
ERNORM 1e-2
#NUMERICS
$PCS_TYPE
聽聽MASS_TRANSPORT
$ELE_GAUSS_POINTS
聽聽2
$LINEAR_SOLVER
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$COUPLING_CONTROL
ERNORM 1e-2
#STOP

On 04.02.2013 13:25, Georg Kosakowski wrote:

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

--
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: http://people.web.psi.ch/kosakowski/

Hoi,
last message for today....the settings below do not seem to work too. As soon as I set the permeability contrast between fracture and rest the solution does not converge to the correct pressure profile.....sorry.
Greetings, Georg

路路路

On 02/04/2013 03:52 PM, Georg Kosakowski wrote:

Hoi,
the num file below seems to work. But can someone explain me why the former numerical settings do not work with LIQUID flow?
In any case thanks to Marc and I hope Pinot Noir (Blauburgunder) is ok.
Greetings, Georg

GeoSys-NUM: Numerical Parameter ----------------------------------------
#NUMERICS
$PCS_TYPE
聽聽LIQUID_FLOW
$ELE_GAUSS_POINTS
聽聽3
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
; 2 6 1.e-12 200000 1.0 1 2
聽聽2 1 1.e-12 1000 1.0 100 4
#NUMERICS
$PCS_TYPE
聽聽MASS_TRANSPORT
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
聽聽2 6 1.e-014 2000 1.0 1 2
$ELE_GAUSS_POINTS
聽聽3
#STOP

On 02/04/2013 02:11 PM, Georg Kosakowski wrote:

Hoi Marc,
thanks and good try, but the wine is not yet yours :wink:
The coupling is not really necessary and (with or without coupling) the Darcy fluxes are still much to high. They will introduce some dispersion and too much advective transport.

Greetings, Georg

On 02/04/2013 01:47 PM, Marc Walther wrote:

Hi Georg,

try using this num:

$OVERALL_COUPLING
2 5
#NUMERICS
$PCS_TYPE
聽聽LIQUID_FLOW
$ELE_UPWINDING
聽聽0.5
$ELE_MASS_LUMPING
聽聽1
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$ELE_GAUSS_POINTS
聽聽2
$COUPLING_CONTROL
ERNORM 1e-2
#NUMERICS
$PCS_TYPE
聽聽MASS_TRANSPORT
$ELE_GAUSS_POINTS
聽聽2
$LINEAR_SOLVER
聽聽2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$COUPLING_CONTROL
ERNORM 1e-2
#STOP

On 04.02.2013 13:25, Georg Kosakowski wrote:

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

--
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: http://people.web.psi.ch/kosakowski/

what do you set as permeability contrast (ie how does the input file look like)?

I was also experimenting with only one point (POINT0) as the pressure boundary, which seemed to convert much faster...

路路路

On 04.02.2013 16:15, Georg Kosakowski wrote:

Hoi,
last message for today....the settings below do not seem to work too. As soon as I set the permeability contrast between fracture and rest the solution does not converge to the correct pressure profile.....sorry.
Greetings, Georg

On 02/04/2013 03:52 PM, Georg Kosakowski wrote:

Hoi,
the num file below seems to work. But can someone explain me why the former numerical settings do not work with LIQUID flow?
In any case thanks to Marc and I hope Pinot Noir (Blauburgunder) is ok.
Greetings, Georg

GeoSys-NUM: Numerical Parameter ----------------------------------------
#NUMERICS
$PCS_TYPE
  LIQUID_FLOW
$ELE_GAUSS_POINTS
  3
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
; 2 6 1.e-12 200000 1.0 1 2
  2 1 1.e-12 1000 1.0 100 4
#NUMERICS
$PCS_TYPE
  MASS_TRANSPORT
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
  2 6 1.e-014 2000 1.0 1 2
$ELE_GAUSS_POINTS
  3
#STOP

On 02/04/2013 02:11 PM, Georg Kosakowski wrote:

Hoi Marc,
thanks and good try, but the wine is not yet yours :wink:
The coupling is not really necessary and (with or without coupling) the Darcy fluxes are still much to high. They will introduce some dispersion and too much advective transport.

Greetings, Georg

On 02/04/2013 01:47 PM, Marc Walther wrote:

Hi Georg,

try using this num:

$OVERALL_COUPLING
2 5
#NUMERICS
$PCS_TYPE
  LIQUID_FLOW
$ELE_UPWINDING
  0.5
$ELE_MASS_LUMPING
  1
$LINEAR_SOLVER
; method error_tolerance max_iterations theta precond storage
  2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$ELE_GAUSS_POINTS
  2
$COUPLING_CONTROL
ERNORM 1e-2
#NUMERICS
$PCS_TYPE
  MASS_TRANSPORT
$ELE_GAUSS_POINTS
  2
$LINEAR_SOLVER
  2 1 1e-12 500 1 100 4
$NON_LINEAR_ITERATIONS
;type -- error_method -- max_iterations -- relaxation -- tolerance(s)
PICARD ERNORM 100 0.0 1e-2
$COUPLING_CONTROL
ERNORM 1e-2
#STOP

On 04.02.2013 13:25, Georg Kosakowski wrote:

Hoi,
I send it to the group..so everyone has the chance to find my mistake :wink:
I will send a bottle of local wine to the one who finds the problem first.

Please change the file endings to tgz and gunzip and untar the files. Google does not like certain file endings inside zipped files. OGS is a already a virus, no need to hide another one in the input files :wink:

Greetings, Georg

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

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

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