Changes
Page history
version 4.2
authored
Oct 30, 2025
by
Udo Ziegler
Hide 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$.
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>
:

{width=500px}

{width=7:w
00px}
At
**x**
=(0.3,0,0) a spherical clump with
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:
BC for the gravitational potential
*Φ*
cannot be explicitely set
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
*Φ*
|
|------------------|:----|:---------------------------------------|
| P | | P (periodicity) |
...
...
@@ -1304,7 +1304,7 @@ function `phiUser()` expecting a user-defined gravitational potential.
possible. The solver only supports triple-periodic BC when periodic.
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
apply and any specified BC are ignored.
...
...
@@ -1319,22 +1319,21 @@ respectively.
Fluid viscosity enters the momentum equation and energy equation via the
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
-2/3(
\n
abla
\c
dot
\m
athbf{v})
\m
athrm{I}
\r
ight]$$
*τ*
=
*ν*
\[
∇
**v**
+(∇
**v**
)
<sup>
⊤
</sup>
−2/3(∇⋅
**v**
)
*I*
\]
where $
\n
u$ in units
\[
$
\m
athtt{kg}
\c
dot
\m
athtt{m}^{-1}
\c
dot
\m
athtt{s}^{-1}$
\]
is the
*dynamic*
viscosity coefficient, $
\m
athbf{v}$ the fluid velocity and
$
\m
athrm{I}$ the identity operator.
where
*ν*
in units
\[
kg⋅m
<sup>
−1
</sup>
⋅s
<sup>
−1
</sup>
\]
is the
*dynamic*
viscosity coefficient,
**v**
the fluid velocity and
*I*
the
identity operator.
**Important:**
The viscosity coefficient
$
\n
u$ 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
related to the dynamic coefficient by
$
\n
u=
\v
arrho
\n
u_{
\m
athrm{kinetic}}$
.
**Important:**
The viscosity coefficient
*ν*
here should not be
confused with the kinetic viscosity coefficient,
*ν*
<sub>
*kinetic*
</sub>
with unit
\[
m
<sup>
2
</sup>
⋅s
<sup>
−1
</sup>
\]
,
which is
related to the dynamic coefficient by
*ν*
= 𝜚
*ν*
<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
viscosityCoeffUser(g,vis);
...
...
...
...