The transient flow calculation can also be called up directly without the graphical user interface via the command line, provided that a batch file *.bsi for the SITRA module or a batch file *.bis for the INSTAT module is available in the directory. When entering the command line, the navigate to the directory in which the calculation is to be carried out. Entering sitra and confirming with the Enter key starts a calculation with the SITRA module. The default batch file name for the SITRA module is sitra.bsi. This file is searched for in the directory. However, another batch file can also be specified with the sitra file name call (the extension ".bsi" is automatically appended if required).
Entering instat and confirming with the Enter key starts a calculation with the INSTAT module. The default batch file name for the INSTAT module is instat.bis. This file is searched for in the directory. However, another batch file can also be specified with the instat file name call (the ".bis" extension is automatically appended if required).
The batch files can also be created or modified manually using any editor.
Example of a batch file for a transient flow computation with the module INSTAT
out.s # Output file
1 # Equation solver (1-iterative, 0-direct)
0 0 # Control lines? Velocity (0 average, 1 right, -1 left)
1 # Calculation of velocities? (1-yes, 0-no)
1 15 1.0 # Iter. thickn.? (1/-1/0)(yes/contin/no), numb., atten.
instat.txt # Transient input file
0 # Detailed protocol (1-yes/0-no)
0 # (1/-1/0)-(warm start, cold start-null, cold start-EICH)
0 3 # Time step 0->from file+factor, 1->time step+time step width+time unit
1 -1 -1 1 -1 0 #Transient consideration of POTE,KNOT,FLAE,VORF,EFLA, transient mining subsidence?
10 # Number of time steps
0 # Saving for optional continuation of the transient computation? (1->yes, 0->no)
0 1 # Save intermediary results 0-> every n.th step, 1 -> for n follow. TS, 2-> hydrographs for n nodes
1 # Computation of pathlines ? (1-yes, 0-no)
0.2 # Porosity
0.2 # Capture radius
6 0 0 # Number of starting points in elements, around nodes, indiv. points
1210 1275 1551 1664 1665 1857 # Element numbers
Example of a batch file for a transient flow computation with the module SITRA
out.s # Output file
0 # 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 # Courant 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
1 0 # Flow/transport steady st./0 trans./1 trans.+conf./unconf./2
1 0 # Initial values POTE/KONZ from ASAT(3) null(2) EICH/AKON(1) =0(0)
instat.txt # Transient input file
0 # Warm start with null-file (1-yes, 0-no)
0 1 1 # time step 0->from file+factor(flow)+factor(transport), 1->time step+time step width+time unit
1 -1 -1 0 0 1 -1 1 #transient consideration of POTE,KNOT,FLAE,KONZ,1KON,VORF,EFLA, mining subsidence?
GRUB POTE # Mine with fixed potential
10 # Number of time steps
0 # Saving for optional continuation of the transient computation? (1->yes, 0->no)
0 1 # Save intermediary results 0-> every n.th step, 1 -> for n follow. TS, 2-> hydrographs for n nodes
SPEITERM VG
EXPORT STRING
The program ignores the characters behind the sign `#'. This space can be used for individual remarks.
The presented batch file for the transient flow can be downloaded on our Homepage.
.
Here you find a detailed description of the batch file for a transient flow computation with INSTAT
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 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!)
6th line: name of the transient input file (starting in 1st column!)
7th line: flag for detailed protocol (flag=1)
8th line: flag for warm/cold start and initial potential head: flag = 1 or -1: warm/cold start with data of the null-file, flag =0: cold start with EICH
9th 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),
1= constant time step length, further input (separated by blank or tab): time step length, unit of time, (0=sec. 1=min. 2=hour. 3=day 4=month 5=year),
10th line: interpolation flags for transient boundary conditions for POTE, KNOT, FLAE, VORF and EFLA
For POTE, KNOT, FLAE, 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)
11th line: number of time steps
12th line: flag for a complete record of the final results (file out66) to continue the transient calculation (1=yes/0=no)
13th 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 13th line:
from the 14th line: time steps in increasing order, maximum 10 per line, without the starting (0) and the last time step!
if Flag=2 in the 13th line:
from the 14th line: node numbers, maximum 10 per line
14+x. line: pathline calculation (1=yes/0=no)
From 15+x. line: input parameter for the path line calculation:
15+x. line: porosity
16+x. line: snapping radius near wells (between 0.0 and 1.0)
17+x. line: 3 flags for number of pathlines with starting points in: in the centre of an element, at nodes, at any location:
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 transient 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
(whithout node numbers) in the 19th line: SPEITERM: calculation of the storage coefficient by Van Genuchten (VG)
in 20th line: EXPORT STRING: Saving the velocity fields for each time step
Stand: 11.06.2025 | 