version 4.2 authored by Udo Ziegler's avatar Udo Ziegler
...@@ -1021,20 +1021,8 @@ square domain of length L (given by `_C.up[0]-_C.lo[0]`): ...@@ -1021,20 +1021,8 @@ square domain of length L (given by `_C.up[0]-_C.lo[0]`):
**B** = *B*<sub>0</sub>(−sin(2*π* *y*/*L*),sin(4*π* *x*/*L*),0) **B** = *B*<sub>0</sub>(−sin(2*π* *y*/*L*),sin(4*π* *x*/*L*),0)
*e* = *p*<sub>0</sub>/(*γ*−1) + **m**<sup>2</sup>/(2𝜚) + **B**<sup>2</sup>/(2*μ*) *e* = *p*<sub>0</sub>/(*γ*−1) + **m**<sup>2</sup>/(2𝜚) + **B**<sup>2</sup>/(2*μ*)
where 𝜚<sub>0</sub>=2.77778, *p*<sub>0</sub>=5/3 and where 𝜚<sub>0</sub>=2.77778, *p*<sub>0</sub>=5/3 and *B*<sub>0</sub>=1.
*B*<sub>0</sub>=1. The adiabatic index *γ*=5/3 and the magnetic The adiabatic index *γ*=5/3 and the magnetic permeability *μ*=1.
permeability *μ*=1.
$\varrho = \varrho_0$
$\mathbf{m}=\varrho_0 \left(-\sin (2\pi\,\!y/L),\sin (2\pi\,\!x/L),0\right)$
$\mathbf{B}=B_0 \left(-\sin (2\pi\,\!y/L),\sin (4\pi\,\!x/L),0\right)$
$e=p_0/(\gamma -1)+\mathbf{m}^2/(2\varrho)+\mathbf{B}^2/(2\mu)$
where $\varrho_0=2.77778$, $p_0=5/3$ and $B_0=1$. The adiabatic index
$\gamma =5/3$ and the magnetic permeability $\mu=1$.
void configUser(GRD **gm) void configUser(GRD **gm)
/** /**
...@@ -1107,15 +1095,19 @@ algebraic expression $a^2+b^2+c^2$. ...@@ -1107,15 +1095,19 @@ algebraic expression $a^2+b^2+c^2$.
[Z05](#references)) [Z05](#references))
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)\in [-1/2,1/2]^3$: box given by (x,y,z)=[-1/2,1/2]<sup>3</sup>:
![ic_example2](uploads/e31ebefeb64b91e0ef74d839ebac4de8/ic_example2.png)
$$(\varrho,p,v_x,v_y,v_z,B_x,B_y,B_z)=\left\{\begin{array}{ll} $$(\varrho,p,v_x,v_y,v_z,B_x,B_y,B_z)=\left\{\begin{array}{ll}
(3.86859,167.345,0,0,0,0,2.1826182,-2.1826182) & x<0.1\\ (3.86859,167.345,0,0,0,0,2.1826182,-2.1826182) & x<0.1\\
(1,1,-11.2536,0,0,0,0.56418958,0.56418958) & x\ge 0.1 (1,1,-11.2536,0,0,0,0.56418958,0.56418958) & x\ge 0.1
\end{array}\right.$$ At ${\mathbf x}=(0.3,0,0)$ a spherical clump with \end{array}\right.$$
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
surrounding flow under the assumption of pressure equilibrium. The surrounding flow under the assumption of pressure equilibrium.
adiabatic index $\gamma =5/3$ and magnetic permeability $\mu=1$. The The adiabatic index *γ*=5/3 and the magnetic permeability *μ*=1.
grid is initially refined with 3 refinement levels, `_C.level`=3, in the The grid is initially refined with 3 refinement levels, `_C.level`=3, in the
vicinity of the shock and cloud (cf. [User-defined initial mesh vicinity of the shock and cloud (cf. [User-defined initial mesh
refinement and refinement refinement and refinement
control](#user-defined-initial-mesh-refinement-and-refinement-control) control](#user-defined-initial-mesh-refinement-and-refinement-control)
... ...
......