Changes
Page history
version 4.2
authored
Oct 30, 2025
by
Udo Ziegler
Show whitespace changes
Inline
Side-by-side
3-NIRVANA-user-guide/3.2-User-interfaces.md
View page @
231a5906
...
@@ -1097,7 +1097,8 @@ algebraic expression $a^2+b^2+c^2$.
...
@@ -1097,7 +1097,8 @@ algebraic expression $a^2+b^2+c^2$.
IC for the shock-cloud interaction problem simulated in a 3D Cartesian
IC for the shock-cloud interaction problem simulated in a 3D Cartesian
box given by (x,y,z)=[-1/2,1/2]
<sup>
3
</sup>
:
box given by (x,y,z)=[-1/2,1/2]
<sup>
3
</sup>
:

{width=500px}

{width=7:w
00px}
At
**x**
=(0.3,0,0) a spherical clump with
At
**x**
=(0.3,0,0) a spherical clump with
radius 0.15 and density of 10 is embedded and co-moving with its
radius 0.15 and density of 10 is embedded and co-moving with its
...
@@ -1279,9 +1280,8 @@ depends on the type of variable as listed in the following table:
...
@@ -1279,9 +1280,8 @@ depends on the type of variable as listed in the following table:
BC for the gravitational potential
*Φ*
cannot be explicitely set
BC for the gravitational potential
*Φ*
cannot be explicitely set
except when specifying U in
`nirvana.par`
. Otherwise, the following rule
except when specifying U in
`nirvana.par`
. Otherwise, the following rule
is adopted to transform MHD BC types into BC types for
$
\P
hi$
:
is adopted to transform MHD BC types into BC types for
*Φ*
:
*Φ*
| selected BC type | ⇒ | BC type for
*Φ*
|
| selected BC type | ⇒ | BC type for
*Φ*
|
|------------------|:----|:---------------------------------------|
|------------------|:----|:---------------------------------------|
| P | | P (periodicity) |
| P | | P (periodicity) |
...
@@ -1304,7 +1304,7 @@ function `phiUser()` expecting a user-defined gravitational potential.
...
@@ -1304,7 +1304,7 @@ function `phiUser()` expecting a user-defined gravitational potential.
possible. The solver only supports triple-periodic BC when periodic.
possible. The solver only supports triple-periodic BC when periodic.
When using the screening mass approach for cylindrical/spherical
When using the screening mass approach for cylindrical/spherical
problems with
$2
\p
i$
2
*π*
-periodicity BC are computed to discretization error
problems with 2
*π*
-periodicity BC are computed to discretization error
accuracy by the gravity solver itself. In this case above rules do not
accuracy by the gravity solver itself. In this case above rules do not
apply and any specified BC are ignored.
apply and any specified BC are ignored.
...
@@ -1319,22 +1319,21 @@ respectively.
...
@@ -1319,22 +1319,21 @@ respectively.
Fluid viscosity enters the momentum equation and energy equation via the
Fluid viscosity enters the momentum equation and energy equation via the
viscous stress tensor $
\t
au$ given by
viscous stress tensor $
\t
au$ given by
viscous stress tensor
*τ*
given by
$$
\t
au=
\n
u
\l
eft[
\n
abla
\m
athbf{v}+(
\n
abla
\m
athbf{v})^
\t
op
*τ*
=
*ν*
\[
∇
**v**
+(∇
**v**
)
<sup>
⊤
</sup>
−2/3(∇⋅
**v**
)
*I*
\]
-2/3(
\n
abla
\c
dot
\m
athbf{v})
\m
athrm{I}
\r
ight]$$
where $
\n
u$ in units
where
*ν*
in units
\[
kg⋅m
<sup>
−1
</sup>
⋅s
<sup>
−1
</sup>
\]
is the
\[
$
\m
athtt{kg}
\c
dot
\m
athtt{m}^{-1}
\c
dot
\m
athtt{s}^{-1}$
\]
is the
*dynamic*
viscosity coefficient,
**v**
the fluid velocity and
*I*
the
*dynamic*
viscosity coefficient, $
\m
athbf{v}$ the fluid velocity and
identity operator.
$
\m
athrm{I}$ the identity operator.
**Important:**
The viscosity coefficient
$
\n
u$ here should not be
**Important:**
The viscosity coefficient
*ν*
confused with the kinetic viscosity coefficient,
here should not be
confused with the kinetic viscosity coefficient,
$
\n
u_{
\m
athrm{kinetic}}$ with unit
\[
m$^2
\c
dot$s$^{-1}$
\]
, which is
*ν*
<sub>
*kinetic*
</sub>
with unit
\[
m
<sup>
2
</sup>
⋅s
<sup>
−1
</sup>
\]
,
related to the dynamic coefficient by
which is
related to the dynamic coefficient by
$
\n
u=
\v
arrho
\n
u_{
\m
athrm{kinetic}}$
.
*ν*
= 𝜚
*ν*
<sub>
*kinetic*
</sub>
.
A user-defined coefficient,
$
\n
u$
, has to be assigned in module
A user-defined coefficient,
*ν*
, has to be assigned in module
`viscosityUser.c`
in the function
`viscosityUser.c`
in the function
viscosityCoeffUser(g,vis);
viscosityCoeffUser(g,vis);
...
...
...
...