Batch file steady state flow

It is possible to execute the steady state flow directly in the command line without using the dialog box, if a batch-file *.bge (for GEONEU) or *.bsi (for SITRA) exists in the directory. For this you have to change to the directory in which the calculation is to be executed. The input of geoneu (or sitra) and confirming with the Enter key starts the flow calculation. The default file name is geoneu.bge or sitra.bsi. This file is searched in the directory.

To use another batch file, type sitra (or geoneu) file name (the extension ".bge" or “.bsi” is automatically attached). Batch files can be created, viewed and modified with the help of any text editor.

Example of a batch file for a steady-state flow computation with the module GEONEU

out.g # Output file

1 1.0e-8 1.0e-6 # Equation solver (1-iterative,0-direct), residual, iter. difference

1 0 0 # Control lines? Velocity (0 average, 1 right, -1 left), vertical flow rate? (1-yes,0-no)

0 # Calculation of velocities? (1-yes, 0-no)

0 # Import K-values from EICHEN? (1-yes, 0-no)

1 5 1 0.5 # sat/unsat? (1/-1/0)(yes/contin/no),numb.,init.potentials EICH, atten.

0 # Calculation of pathlines ? (1-yes, 0-no)

Example of a batch file for a steady-state flow computation with the module SITRA

out.s # Output file

1 # detailed protocol (1-yes/0-no)

1 0 # Equation solver P/U (1-iterative, 0-direct, -1=multigrid)

1 0 # Calculation of control lines (1-yes, 0-no), IMITTEL

0 0 # Peclet numbers, hydr. cond. from EICHEN (1-yes, 0-no)

STROEMUNG #

1 # satur./unsat. computation (1-yes,0-no)

0 0 0 # equal conc., density-related calc., p or h (1-yes/0-no)

5 0.5 # Number of iteration steps, factor of attenuation

0 0 # Flow/transport steady st./0 trans./1 trans.+conf./unconf./2

0 0 # Initial values POTE/KONZ from ASAT(3) null(2) EICH/AKON(1) =0(0)

The program ignores the characters behind the sign `#'. This space can be used for individual remarks.

The presented batch file for the steady state flow can be downloaded on our Homepage.

Here you find a detailed

description of the batch file for a steady state flow computation with GEONEU

1st line: name of the output file (starting at the first position!)

2nd line: flag for the solver: 0=direct, 1=iterative, -1=multigrid

3rd line: three flags: 1st flag for calculation of control lines (yes=1), 2nd flag for velocities for calculation of control lines (0=mean velocity, 1=velocity on the right hand side of the control line, -1=velocity on the left hand side of the control line), 3rd flag for calculation of vertical mass flow rates (3D), (only for 3D models without ceased elements! )

4th line: flag for calculation of velocities (yes=1), (automatically in the case of calculating control lines, mass flow rates or path lines)

5th line: flag for using calibrated hydraulic conductivity (yes=1), (only if calibration was running before!)

6th line: flag and parameter for iteration of thickness / free surface:

Three-dimensional and vertical models:

flag for iteration of the free surface(1 or -1 = yes/0=no), number of iteration (if flaf = 1 or -1), flag for initial values from EICH (flag=1), damping factor for iterations

Horizontal models with UNTE:

Flag for iteration of thickness (1 or -1=yes, 0=no), number of iterations, flag for protocol of iter. Thickness (1=yes/0=no), flag for protocol of iter. Potentials (1=yes/0=no), attenuation factor for iteration

With flag=-1 the iteration is performed, too. If possible, a former iteration is continued! For this case the use of calibrated hydraulic conductivities from EICHEN is not possible anymore!)

7th line: flag for pathline calculation (1=yes/0=no)

from the 8th line: input parameter for the path line calculation:

8th line: maximum length (m) of a path line

9th line: maximum flowing time (unit of time) of a path line

10th line: path lines equal (1) or against (0) the flow direction

11th line: porosity

12th line: snapping radius near wells (between 0.0 and 1.0)

13th line: 3 flags for number of pathlines with starting points in: in the centre of an element, at nodes, at any location:

From 14^th line:

Start in the centres of the elements: element numbers (maximum 10 per line)

Start in nodes: first line: number of starting points per node, radius (m), next lines: node numbers (maximum 10 per line)

Start at any location: line with x, y (and z) coordinates (3D model) for each starting point

Here you find a detailed

description of the batch file for a flow computation with SITRA

1st line: name of the output file (starting at the 1st column!)

2nd line: flag for complete record (flag=1)

3rd line: flag for the solver: 0=direct, 1=iterative

4th line: flag for computing control lines (1=yes/0=no), if flag = 1:

Flag for computing the velocities for the control lines:

0 = average between left and right velocity

1 = velocity from the right of the control line

-1 = velocity from the left of the control line

5th line: flag for the computing of Peclet numbers and (transient calculation only) Courant numbers (1=yes/0=no), flag for using the calibrated hydraulic conductivity from calibration (1=yes/0=no)

6th line: 'STROEMUNG' (only flow modelling) or 'TRANSPORT' (flow and transport modelling)

7th line: 3D model or vertical model: flag for a saturated/unsaturated calculation (1=yes/0=no)

Horizontal model using UNTE: flag for iteration of thickness (1=yes/0=no), else irrelevant line

8th line: 1. flag for transport: considering equal concentration at equal potential head (1=yes/0=no), 2. flag only important in 3D model or vertical model: density dependent calculation (1=yes/0=no), 3. flag: considering pressure or potential equation (1=pressure/0=potential)

9th line: Saturated/unsaturated calculation or iteration of thickness and/or density dependent calculation: number of iterations, attenuation factor, else irrelevant

Remark: For transient calculations with an iteration of thickness the number of iterations per time step is always 1 (no iteration!). Then, the attenuation factor should be equal 1.0 (thickness is only updated).

10th line: flag for steady state or transient flow and flag for steady state or transient transport: 3D model or vertical model: 0= steady state, 1= transient,

horizontal model with groundwater filled thickness:

0=steady state, 1=transient (without consideration of confined/unconfined conditions), 2=transient (with consideration of confined/unconfined conditions)

11th line: 2 flags for initial values of potential head and concentration/temperature (potential head for sat./unsat. calculation or transient flow) and concentration/temperature for density dep. calculation or transient transport: 2= data from file null (if file null not found then initial value = EICH/AKON), 1= initial potential head EICH or concentration/temperature AKON (if no values for EICH or AKON found then initial value = 0), 0= start value = 0.0

from the 12th line: input parameter for transient calculation

12th line: name of the transient input file (starting in 1st column!)

13th line: flag for warm start using the data in the file null (1=yes/0=no), (if file null not found then input is ignored!)

14th line: flag for definition of transient time steps:

0= flag for using the time steps from the transient input file, further input: factor for scaling down the time steps (integer), add. factor for scaling down time steps in transport eq. (optional)

1= constant time step length, further input (separated by blank or tab): time step length (double), unit of time (integer), (0=sec. 1=min. 2=hour. 3=day 4=month 5=year), add. factor for scaling down time steps in transport eq. (optional)

15th line: interpolation flags for transient boundary conditions for POTE, KNOT, FLAE, KONZ, 1KON,VORF and EFLA

For POTE, KNOT, FLAE, KONZ, 1KON, VORF and EFLA: 0= implicit use without interpolation, 1= implicit use with interpolation

For FLAE, EFLA and KNOT: -1= explicit use without interpolation, -2= explicit use with interpolation, flag for calculation with transient terrain settings as a result of subsurface mining (1=yes/0=no)

16th line: number of time steps or input for flooding: GRUB POTE or MENGE

17th line: flag for a complete record of the final results (file out66) to continue the transient calculation (flag=1)

18th line: flag for saving intermediate results to binary files: 0= flag for saving the results at every n-th time step, further input: n, 1= flag for saving the results at n defined time steps, further input: n, 2= flag for saving complete time-series plots at single nodes, further input: n

if Flag=1 in the 18th line:

from the 19th line: time steps in increasing order, maximum 10 per line, without the starting (0) and the last time step!

if Flag=2 in the 18th line:

from the 19th line: node numbers, maximum 10 per line