Abstract
In this article, we propose high order, semi-discrete, entropy stable, finite difference schemes for Ten-Moment Gaussian Closure equations. The crucial components of these schemes are an entropy conservative flux and suitable high order entropy dissipative operators to ensure entropy stability. We design two numerical fluxes, one is approximately entropy conservative, and another is entropy conservative flux. For the construction of appropriate entropy dissipative operators, we also derive entropy scaled right eigenvectors. This is used for sign preserving reconstruction of scaled entropy variables, which results in second and third order entropy stable schemes. We also extend these schemes to a plasma flow model with source term. Several numerical results are presented for homogeneous and non-homogeneous cases to demonstrate stability and performance of these schemes.








Similar content being viewed by others
References
Barth, T.J.: Numerical methods for gas-dynamics systems on unstructured meshes. In: Kroner, D., Ohlberger, M., Rohde, C. (eds.) An Introduction to Recent Developments in Theory and Numerics of Conservation Laws. Lecture Notes in Computational Science and Engineering, vol. 5, pp. 195–285. Springer, Berlin (1999)
Berthon, C.: Numerical approximations of the 10-moment Gaussian-closure. Math. Comput. 75(256), 1809–1831 (2006)
Berthon, C., Desveaux, V.: An entropy preserving MOOD scheme for the Euler equations. Int. J. Finite 11, 1–39 (2014)
Berthon, C., Dubroca, B., Sangam, A.: An entropy preserving relaxation scheme for the Ten-Moments equations with source terms. Commun. Math. Sci. 13(8), 2119–2154 (2015)
Chandrashekar, P.: Kinetic energy preserving and entropy stable finite volume schemes for compressible euler and Navier–Stokes equations. Commun. Comput. Phys. 14, 1252–1286 (2013)
Brown, S.L.,Roe, P.L., Groth, C.P.: Numerical solution of a 10-moment model for nonequilibrium gasdynamics. In: 12th AIAA Computational Fluid Dynamics Conference (1995)
Dubroca, B., Tchong, M., Charrier, P., Tikhonchuk, V.T., Morreeuw, J.P.: Magnetic field generation in plasmas due to anisotropic laser heating. Phys. Plasm. 11, 3830–3839 (2004)
Farzad Ismail, R., Philip, L.: Affordable, entropy-consistent euler flux functions II: entropy production at shocks. J. Comput. Phys. 228(15), 5410–5436 (2009)
Fjordholm, U.S., Mishra, S., Tadmor, E.: Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM J. Numer. Anal. 50(2), 544–573 (2012)
Fjordholm, U.S., Mishra, S., Tadmor, E.: ENO reconstruction and ENO interpolation are stable. Found. Comput. Math. 13(2), 139–159 (2012)
Fjordholm, U.S.: High-order accurate entropy stable numerical schemes for hyperbolic conservation laws. Diss. No. 21025, Ph.D. Thesis, ETH Zurich (2013)
Gottlieb, S., Shu, C.W., Tadmor, E.: Strong stability-preserving high-order time discretization methods. SIAM Rev. 43, 89–112 (2001)
Hakim, A.: Extended MHD modelling with ten-moment equations. J. Fusion Energy 27, 36–43 (2008)
Johnson, E.A.: Gaussian-Moment relaxation closures for verifiable numerical simulation of fast magnetic reconnection in plasma. Ph.D. Thesis, UW-Madison (2011)
LeFloch, P.G., Mercier, J.M., Rohde, C.: Fully discrete entropy conservative schemes of arbitrary order. SIAM J. Numer. Anal. 40(5), 1968–1992 (2002)
Levermore, C.D.: Moment closure hierarchies for kinetic theories. J. Stat. Phys. 83, 1021–1065 (1996)
Levermore, C.D., Morokoff, W.J.: The Gaussian moment closure for gas dynamics. SIAM J. Appl. Math. 59(1), 72–96 (1996)
Meena, A.K., Kumar, H.: A robust MUSCL scheme for ten-moment Gaussian closure equations with source terms. Submitted
Meena, A.K., Kumar, H., Chnadrashekar, P.: Positivity-preserving high-order discontinuous Galerkin schemes for ten-moment Gaussian closure equations. J. Comput. Phys. 339, 370–395 (2017)
Ray, D., Chandrashekar, P., Fjordholm, U.S., Mishra, S.: Entropy stable scheme on two-dimensional unstructured grids for Euler equations. Commun. Comput. Phys. 19, 1111–1140 (2016)
Sangam, A., Morreeuw, J.P., Tikhonchuk, V.T.: Anisotropic instability in a laser heated plasma. Phys. Plasm. 14, 053111 (2007)
Sangam, A.: An HLLC scheme for ten-moments approximation coupled with magnetic field. Int. J. Comput. Sci. Math. 2, 73–109 (2008)
Shu, C.W.: TVD time discretizations. SIAM J. Math. Anal. 14, 1073–1084 (1988)
Struchtrup, H., Torrilhon, M.: Regularization of grad’s 13-moment-equations: derivation and linear analysis. Phys. Fluids 15(9), 2668–2680 (2003)
Struchtrup, H.: Macroscopic Transport Equations for Rarefied Gas Flows. Springer, Berlin (2005)
Tadmor, E.: The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Math. Comp. 49, 91–103 (1987)
Tadmor, E.: Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems. Act. Numer. 12, 451–512 (2004)
Author information
Authors and Affiliations
Corresponding author
Additional information
Ms. Chhanda Sen is supported by CSIR Fellowship at Dept. of Mathematics, IIT Delhi. Harish Kumar has been funded in part by SERB, DST Grant with file No. YSS/2015/001663.
Appendix: Symmetrization and Entropy Scaled Eigenvectors for Ten-Moment Equations
Appendix: Symmetrization and Entropy Scaled Eigenvectors for Ten-Moment Equations
In this section we will derive the entropy scaled right eigenvectors. For simplicity, let us consider the x-directional homogenous part of the system (1), i.e.
The entropy-entropy flux pair for this system is same as for the complete system. So, the results presented here can be easily extended to the complete system. Let us recall that,
and
Theorem A.1
(Barth [1]) (Eigenvector scaling)
Let \(A\in \mathbb {R}^{n\times n}\) be an arbitrary diagonalizable matrix and S the set of all right symmetrizers:
Further, let \(\mathbf{R}\in \mathbb {R}^{n\times n}\) denote the right eigenvector matrix which diagonalizes A, \(A=\mathbf{R}\Lambda \mathbf{R}^{-1}\) with r distinct eigenvalues, \(\Lambda =Diag(\lambda _1I_{m_1\times m_1}, \lambda _2I_{m_2\times m_2}, \ldots , \lambda _rI_{m_r\times m_r}).\) Then for each \(B\in S\) there exists a symmetric block diagonal matrix \(T=Diag(T_{m_1\times m_1}, T_{m_2\times m_2},\ldots T_{m_r\times m_r} )\), such that the block scales columns of \(\mathbf{R}\), \(\tilde{\mathbf{R}}=\mathbf{R}T\), with
which imply
By change of variable from \(\mathbf{u}\) to \(\mathbf{v}\) we can rewrite (26) as,
As we have \((e,q_{x})\) as entropy pair for the system (26), the above change of variable should symmetrize the system (see [1]). Furthermore, the matrix \(\frac{\partial \mathbf{u }}{\partial \mathbf{v }}\) is symmetric positive definite. We now want to find scaled entropy right eigenvectors \(\tilde{\mathbf{R}}\) such that,
To simplify the calculation of the entropy scaling, let us introduce primitive variables \(\mathbf{w}=\{\rho , \vec {v},\mathbf{p}\}^{\top }.\) Rewriting the system (26) in primitive variables, we get eigenvalues as,
and corresponding right eigenvector matrix as,
Once the entropy scaled right eigenvector for primitive variable system, \(\tilde{\mathbf{R}}_{\mathbf{w}}^x\) is found, entropy scaled right eigenvectors for the conservative system can be easily calculated using,
Here, the jacobian matrix for the change of variable is,
Using Eqs.(28) and (29), we note that,
which implies,
To calculate the matrix \(\frac{\partial \mathbf{w }}{\partial \mathbf{v }}\) we need to write primitive variables in terms of entropy variables \(\mathbf{v}\). Let us define the notations,

Plugging in the expression for \(v_{4},v_{5}\) and \(v_{6}\) into the equation for \(v_{2}\) and \(v_{3}\) we get,
On solving for \(v^{x}\) and \(v^{y}\), we get,
We now rewrite \(v_{1}\) as,
where we have used
This results in,
Now using expression for \(v_{4}\),
Similarly,
The matrix \(\frac{\partial \mathbf{w }}{\partial \mathbf{v }}\) can now be easily calculated and is found to be,
Furthermore, the matrix \(\frac{\partial \mathbf{w }}{\partial \mathbf{v }}\frac{\partial \mathbf{u }}{\partial \mathbf{w }}^{-\top }\) is,
Following the Theorem (A.1) in [1], the scaling matrix \(\mathbf{T}^x\) is square root of matrix \(\mathbf{Y} =\mathbf{R}_{\mathbf{w}}^{-1x} \frac{\partial \mathbf{w }}{\partial \mathbf{v }}\frac{\partial \mathbf{u }}{\partial \mathbf{w }}^{-\top } \mathbf{R}_{\mathbf{w}}^{-\top x}\), which results in,
where \(s^x=\frac{p^{xx}p^{yy}-p^{xy2}}{\sqrt{3}p^{xx}}\) and \(t^{x2}=\frac{\rho }{3}+\frac{4}{3}\frac{(p^{xx}p^{yy}-p^{xy2})^2}{3p^{xx2}\rho }+\frac{2}{\sqrt{3}}\frac{p^{xx}p^{yy}-p^{xy2}}{p^{xx}}\). So, the entropy scaled right eigenvector matrix in primitive variables is,
where \(A^x=\frac{1}{2\sqrt{3}p^{xx}\sqrt{\rho }}\), \(B^x=\frac{\sqrt{p^{xx}p^{yy}-p^{xy2}}}{2p^{xx}\sqrt{\rho }}\), \(\beta _1=\frac{p^{xx}}{(p^{xx}p^{yy}-p^{xy2})}\). Using (29), we can now calculate entropy scaled right eigenvectors in conservative variables which satisfy (28).
Remark A.1
Proceeding similarly in y-direction, with right eigenvector matrix for system in primitive variable as,
we obtained scaling matrix \(\mathbf{T}^y\) as,
where \(s^y=\frac{p^{xx}p^{yy}-p^{xy2}}{\sqrt{3}p^{yy}}\) and \(t^{y2}=\frac{\rho }{3}+\frac{4}{3}\frac{(p^{xx}p^{yy}-p^{xy2})^2}{3p^{yy2}\rho }+\frac{2}{\sqrt{3}}\frac{p^{xx}p^{yy}-p^{xy2}}{p^{yy}}\). Using this, the entropy scaled right eigenvectors in primitive variable in y-direction are,
with \(A^y=\frac{1}{2\sqrt{3}p^{yy}\sqrt{\rho }}\), \(B^y=\frac{\sqrt{p^{xx}p^{yy}-p^{xy2}}}{2p^{yy}\sqrt{\rho }}\), \(\beta _2=\frac{p^{yy}}{(p^{xx}p^{yy}-p^{xy2})}\).
Rights and permissions
About this article
Cite this article
Sen, C., Kumar, H. Entropy Stable Schemes For Ten-Moment Gaussian Closure Equations. J Sci Comput 75, 1128–1155 (2018). https://doi.org/10.1007/s10915-017-0579-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-017-0579-4