Other Mozart Inputs
Input which is not related to mesh/node/grids is collected under the "stupid inputs" categories (Philippe's terminology). These inputs include parameters related to linear systems, physical properties of materials, boundary conditions, and so on.
The input files for this category have the prefix "allstupids." in the file name. The middle part of the filename indicates what the input is for:
"aleg" for ALE/regridding"Xlsys" for linear systems solutions
"Xbcs" for boundary conditions
"Xers" for errors
"Xins" for initializations
"Xmedia" for "media" (physical properties, e.g. rheology)
where X can be t (thermal) for advection-diffusion production problems, v for velocity problems, or m for mesh.
|
Name of file |
Input related to... |
|---|---|
allstupids.aleg.i |
ALE/regridding |
allstupids.ters.i |
thermal errors |
allstupids.tmedia.i |
definition of physical materials for thermal system |
allstupids.vins.i |
velocity initializations |
allstupids.mlsys.i |
mesh linear systems solution |
allstupids.tins.i |
thermal initializations |
allstupids.vbcs.i |
velocity boundary conditions |
allstupids.vlsys.i |
velocity linear systems |
allstupids.tbcs.i |
thermal boundary conditions |
allstupids.tlsys.i |
thermal linear systems solution |
allstupids.vers.i |
velocity errors |
allstupids.vmedia.i |
definition of physical materials for velocity system |
The file Syimpt.01.i defines which of the above inputs will be read by Mozart. The Syimpt.01.i file is read when the "ALL" button in the "IOS" menu is pushed.
[.... Format for Syimpt.01.i here ...]
|
itype in Syimpt.01.i |
Meaning |
|---|---|
|
1 |
initialization of ADP equations (temperature initialization) - initializes fields on grids, etc. [input file allstupids.tins.i] |
|
2 |
initialization of VF equations (viscous flow, velocity) [input file allstupids.vins.i] |
|
3 |
import thermal media (physics/dynamics of materials for thermal) [input file allstupids.tmedia.i] |
|
4 |
import VF media (materials, material properties, rheologies, etc.) [input file allstupids.vmedia.i] |
|
5 |
linear systems for temperature [input file allstupids.tlsys.i] |
|
6 |
linear systems for velocities [input file allstupids.vlsys.i] |
|
7 |
linear systems for meshes [input file allstupids.mlsys.i] |
|
8 |
boundary conditions for thermal problems [input file allstupids.tbcs.i] |
|
9 |
boundar conditions for velocity problems [input file allstupids.vbcs.i] |
|
10 |
errors for thermal problems [input file allstupids.ters.i] |
|
11 |
errors for velocity problems [input file allstupids.vers.i] |
|
12 |
importing ALE graphs (defines which part of Eulerian grid is related to which part of Lagrangian grid or cloud) [input file allstupids.aleg.i] |
FILE : mlsys.i
1 ! template number
2 ! total number of parameters in this template
2 1 ! number of parameters in the following section; number of lines of documentation
ndof = number of degrees of freedom; mlinks = maximum number of nonzeros in the columns of a matrix
2 18 ! ndof mlinks
0 0 ! number of parameters in the following section; number of lines of documentation
FILE : tbcs.i
5 ! Template number [NOTE: template 5 is just a silly example]
88 ! Total number of parameters in this template
11 0 ! Number of parameters in the next section; number of lines of documentation
0 0 0 0 0 0 0 0 0 0 1 ! These are the 11 parameters for this section
11 1 ! Number of parameters in the next section; number of lines of documentation
mozart ! one line of documentation
0 0 0 0 0 0 0 0 0 0 1 ! These are the 11 parameters for this section
11 2 ! Number of parameters in the next section; number of lines of documentation
is ! Two lines
nice ! of documentation
0 0 0 0 0 0 0 0 0 0 1 ! These are the 11 parameters for this section
11 3 ! Number of parameters in the next section; number of lines of documentation
mozart ! Three
is ! lines
nice ! of documentation
0 0 0 0 0 0 0 0 0 0 1 ! These are the 11 parameters for this section
11 0 ! Number of parameters in the next section; number of lines of documentation
0 0 0 0 0 0 0 0 0 0 0 ! These are the 11 parameters for this section
11 0 ! Number of parameters in the next section; number of lines of documentation
0 0 0 0 0 0 0 0 0 0 0 ! These are the 11 parameters for this section
11 0 ! Number of parameters in the next section; number of lines of documentation
0 0 0 0 0 0 0 0 0 0 0 ! These are the 11 parameters for this section
11 0 ! Number of parameters in the next section; number of lines of documentation
0 0 0 0 0 0 0 0 0 0 0 ! These are the 11 parameters for this section
0 0
1 ! voulu Template number
88 ! Total number of parameters for this template
11 2 ! # of parameters in next section; # of lines of doc
Lt1(y) p1<y<p2 t1=p3.y+p4+p5.exp(-|y-p7|/p6)+p8.degrecosinus((y-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 1 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Lf1(y) p1<y<p2 t1=p3.y+p4+p5.exp(-|y-p7|/p6)+p8.degrecosinus((y-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 1 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Rt1(y) p1<y<p2 t1=p3.y+p4+p5.exp(-|y-p7|/p6)+p8.degrecosinus((y-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 1 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Rf11(y) p1<y<p2 t1=p3.y+p4+p5.exp(-|y-p7|/p6)+p8.degrecosinus((y-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 1 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Tt1(x) p1<x<p2 t1=p3.x+p4+p5.exp(-|x-p7|/p6)+p8.degrecosinus((x-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 0 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Tf1(x) p1<x<p2 t1=p3.x+p4+p5.exp(-|x-p7|/p6)+p8.degrecosinus((x-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 0 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Bt1(x) p1<x<p2 t1=p3.x+p4+p5.exp(-|x-p7|/p6)+p8.degrecosinus((x-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 0 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
11 2 ! # of parameters in next section; # of lines of doc
Bf1(x) p1<x<p2 t1=p3.x+p4+p5.exp(-|x-p7|/p6)+p8.degrecosinus((x-p10)/p9)
p11 :1 ignore(=free) 0 use
0 0 0 0 0 0 0 0 0 0 0 ! Parameters p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11
0 0 ! # of parameters in next section; # of lines of doc
0 ! termination template
FILE : ters.i
1 ! Template number
2 ! Total number of parameters for this template
2 1 ! # of parameters in next section; # of lines of doc
maxiter = maximum number of iterations, maxerr = maximum error on t or f
10 1.d-2 ! maxiter ,maximum error on t or f
0 0 ! # of parameters in next section; # of lines of doc
FILE : tins.i
1 ! Template number
8 ! Total number of parameters for this template
4 1 ! # of parameters in next section; # of lines of doc
Initial temperature t0(y)= a1+a2.y+a3.y.y+a4.y.y.y
0. 1. 0. 0. ! a1 a2 a3 a4
4 1 ! # of parameters in next section; # of lines of doc
f0(y)= a1+a2.exp(-(y-a3)/a4)
1 1 1 1 ! a1 a2 a3 a4
0 0 ! # of parameters in next section; # of lines of doc
FILE : vins.i
1 ! Template number
21 ! Total number of parameters for this template
1 1 ! # of parameters in next section; # of lines of doc
Initial dummy uniform viscosity
1.d0
7 2 ! # of parameters in next section; # of lines of doc
cell material numbers pass 1 grid-aligned layers (1d boxes aligned with x)
<#layers> foreach < iy1 iy2 boxvalue> last hit wins
2 1 15 1 15 51 2 !
13 5 ! # of parameters in next section; # of lines of doc
cell material numbers : rectangular boxes
<#Rboxes> ,<type>
foreach of them < IFTYPE=0 x1,y1,x2,y2,boxvalue rectangle >
< IFTYPE=1 xc lx yc ly,boxvalue ellipsis >
. first hit wins
2 0 1.5 .5 1.8 .8 3. 1 3. .8 .5 .3 4.
0 0 ! # of parameters in next section; # of lines of doc
FILE : vmedia.i
2 ! <MEDIUM type=VF name=ZOZO>
2
9
3 2
ZOZO:VF:MODEL:
ZOZO:VF:DATA :SUBMEDIUM 1
1 2 3
1 1
us flag
4
3 0
5 6 7
2 0
8 9
0 0
9
8 2
ZOZO:VF:MODEL:
ZOZO:VF:DATA :SUBMEDIUM 2
1 1 1 1 1 1 1 1
1 1
us flag
1
0 0
1 ! <MEDIUM type=VF name=SIMPLE>
4
13
13 2
SIMPLE:VF:MODEL:submedium 1
fx fy vist sc0 edot0 n ts tc vmin vmax coh friction minp
0 0 2.d7 1.d-1 1.d0 4 0 0 1.d0 1.d8 1.d-2 30.d0 1.d-2
0 0
13
13 2
SIMPLE:VF:MODEL:submedium 2
fx fy vist sc0 edot0 n ts tc vmin vmax coh friction minp
0 0 2.d7 1.d0 1.d0 4 0 0 1.d0 1.d8 1.d-2 30.d0 1.d-2
0 0
13
13 2
SIMPLE:VF:MODEL:submedium 3
fx fy vist sc0 edot0 n ts tc vmin vmax coh friction minp
0 0 2.d7 1.d0 1.d0 4 0 0 1.d0 1.d8 1.d-2 30.d0 1.d-2
0 0
13
13 2
SIMPLE:VF:MODEL:submedium 4
fx fy vist sc0 edot0 n ts tc vmin vmax coh friction minp
0 5 2.d7 1.d-2 1.d0 4 0 0 1.d0 1.d8 1.d-2 30.d0 1.d-2
0 0
10 ! <MEDIUM type=VF name=BOEUF-BOURGUIGNON>
2 ! --> #submedia in medium boeuf-bourguignon
116 ! #parms in submedium -------> boeuf-bourguignon-1
1 1
bulk viscosity
2.d7
1 2
UCRUST:VF-NLCREEP:MODEL:
NLCREEP-1 : stress scale
1.84374d-6
8 1
NLCREEP-2 : strain control on stress scale
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
NLCREEP-3 : pressure control on stress scale
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 3
NLCREEP-4 : temperature control on stress scale
this is additional to the exponential boltzmann factor
which will be seen later
1. -1. 1.d10 2 -1. 1. 1.d10 1.
2 2
NLCREEP-5 : strain rate control on stress scale
edot and power in power law
1. 4.
1 1
NLCREEP-6 : shift to absolute temperatures: Ts
0
1 2
NLCREEP-7 : boltzmann energy (in shifted temperature units) T part :
E=(Ts+Tc)
.989872
2 3
NLCREEP-8 : boltzmann energy (in shifted temperature units) P part
expressed as E*a*(p/pc).
so Etotal=E(1+a.(p/pc)) . enter a (e.g) 0 or 1 and pc (pressure scale)
0. 1.
1 1
DRUCKER-PRAGER-1 : friction angle scale
1.
8 1
DRUCKER-PRAGER-2 : strain control on friction
1. -1. 1.d10 2 -1. 30. 1.d10 30.
8 1
DRUCKER-PRAGER-3 : pressure control on friction
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-4 : temperature control on friction
1. -1. 1.d10 2 -1. 1. 1.d10 1.
1 1
DRUCKER-PRAGER-5 : cohesion scale
1.d-2
8 1
DRUCKER-PRAGER-6 : strain control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-7 : pressure control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-8 : temperature control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
1 2
BODY-FORCES-1
#force ndofs (or axis)
1
6 3
BODY-FORCES-2
on axis naxis volumic source and taylor t p expansion
naxis r0 t0 p0 drdt drdp
2 -1 0 0 0 0
9 1
BODY-FORCES-3 : strain control on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
9 1
BODY-FORCES-4 : pressure control on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
9 1
BODY-FORCES-5 : temperature control on on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
0 0
116 ! #parms in submedium -------> boeuf-bourguignon-2
1 1
bulk viscosity
2.d7
1 2
UCRUST:VF-NLCREEP:MODEL:
NLCREEP-1 : stress scale
1.84374d-6
8 1
NLCREEP-2 : strain control on stress scale
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
NLCREEP-3 : pressure control on stress scale
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 3
NLCREEP-4 : temperature control on stress scale
this is additional to the exponential boltzmann factor
which will be seen later
1. -1. 1.d10 2 -1. 1. 1.d10 1.
2 2
NLCREEP-5 : strain rate control on stress scale
edot and power in power law
1. 4.
1 1
NLCREEP-6 : shift to absolute temperatures: Ts
0
1 2
NLCREEP-7 : boltzmann energy (in shifted temperature units) T part :
E=(Ts+Tc)
.989872
2 3
NLCREEP-8 : boltzmann energy (in shifted temperature units) P part
expressed as E*a*(p/pc).
so Etotal=E(1+a.(p/pc)) . enter a (e.g) 0 or 1 and pc (pressure scale)
0. 1.
1 1
DRUCKER-PRAGER-1 : friction angle scale
1.
8 1
DRUCKER-PRAGER-2 : strain control on friction
1. -1. 1.d10 2 -1. 30. 1.d10 30.
8 1
DRUCKER-PRAGER-3 : pressure control on friction
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-4 : temperature control on friction
1. -1. 1.d10 2 -1. 1. 1.d10 1.
1 1
DRUCKER-PRAGER-5 : cohesion scale
1.d-2
8 1
DRUCKER-PRAGER-6 : strain control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-7 : pressure control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
8 1
DRUCKER-PRAGER-8 : temperature control on cohesion
1. -1. 1.d10 2 -1. 1. 1.d10 1.
1 2
BODY-FORCES-1
#force ndofs (or axis)
1
6 3
BODY-FORCES-2
on axis naxis volumic source and taylor t p expansion
naxis r0 t0 p0 drdt drdp
2 -1 0 0 0 0
9 1
BODY-FORCES-3 : strain control on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
9 1
BODY-FORCES-4 : pressure control on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
9 1
BODY-FORCES-5 : temperature control on on bodyforces
2 1. -1. 1.d10 2 -1. 1. 1.d10 1.
0 0
5 ! <MEDIUM type=VF name=RIRI>
2
9
8 2
ZOZO:T-ADP:MODEL:
ZOZO:T-ADP:DATA :SUBMEDIUM 1
1 1 1 1 1 1 1 1
1 1
us flag
1
0 0
9
8 2
ZOZO:T-ADP:MODEL:
ZOZO:T-ADP:DATA :SUBMEDIUM 2
1 1 1 1 1 1 1 1
1 1
us flag
1
0 0