Hoi,
I send it to the group..so everyone has the chance to find my mistake 
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
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 
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
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
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
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
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: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)
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
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
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
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: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)
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
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
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
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: Telefonbuch & E-Mail | Paul Scherrer Institut (PSI)
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
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
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
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
- - - - - - - - - - - - - - - - - -