To numerically test the turbulent dynamo theory in a weakly ionized medium, it is necessary to carry out
a two-fluid simulation that treats the dynamics of neutrals and ions separately.
[1]}
performed a 3D high-resolution two-fluid dynamo simulation to
solve the following equations
[2]}
by using the two-fluid RIEMANN code
[3]}, [4]},
\(\begin{aligned}&\frac{\partial \rho _i}{\partial t} + \nabla \cdot (\rho _i v_i) = 0, \\&\frac{\partial v_i}{\partial t} + (v_i \cdot \nabla )v_i =-c_s^2 \nabla \ln \rho _i - \frac{1}{4\pi }B\times (\nabla \times B) \\& ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - \gamma _d \rho _n (v_i - v_n), \\&\frac{\partial B}{\partial t} = \nabla \times ( v_i \times B), \\&\frac{\partial \rho _n}{\partial t} + \nabla \cdot (\rho _n v_n) = 0, \\&\frac{\partial v_n}{\partial t} + (v_n \cdot \nabla )v_v =-c_s^2 \nabla \ln \rho _n - \gamma _d \rho _i (v_n - v_i),\end{aligned}\) 
