# American Institute of Mathematical Sciences

October 2018, 11(5): 1255-1276. doi: 10.3934/krm.2018049

## Second-order mixed-moment model with differentiable ansatz function in slab geometry

 Fachbereich Mathematik, TU Kaiserslautern, Erwin-Schrödinger-Str., 67663 Kaiserslautern, Germany

Received  January 2017 Revised  October 2017 Published  May 2018

Mixed-moment minimum-entropy models (${\rm{M}}{{\rm{M}}_N}$ models) are known to overcome the zero net-flux problem of full-moment minimum entropy ${{\rm{M}}_N}$ models but lack regularity. We study differentiable mixed-moment models (full zeroth and first moment, half higher moments, called ${\rm{DM}}{{\rm{M}}_N}$ models) for a Fokker-Planck equation in one space dimension. Realizability theory for these modification of mixed moments is derived for second order. Numerical tests are performed with a kinetic first-order finite volume scheme and compared with ${{\rm{M}}_N}$, classical ${\rm{M}}{{\rm{M}}_N}$ and a ${{\rm{P}}_N}$ reference scheme.

Citation: Florian Schneider. Second-order mixed-moment model with differentiable ansatz function in slab geometry. Kinetic & Related Models, 2018, 11 (5) : 1255-1276. doi: 10.3934/krm.2018049
##### References:

show all references

##### References:
Two ansatz functions and their derivatives for the ${\rm{M}}{{\rm{M}}_2}$ and ${\rm{DM}}{{\rm{M}}_2}$ model, respectively.
Left: ${\hat \psi }(\mu ) = \exp\left(\left(\mu -\mu ^2\right)1_{[0, 1]}-2\mu 1_{[-1, 0]}\right)$,
Right: ${\hat \psi }(\mu ) = \exp\left(-2\mu -\mu ^21_{[0, 1]}\right)$
The normalized realizable set for the differentiable mixed-moment basis of order $N = 2$.
Online version: Press to activate 3D view ($x$-axis (red): ${\phi _{2 + }}$, $y$-axis (green): ${\phi _{2 - }}$, $z$-axis (blue): ${{\phi }_{1}}$)
Eigenvalues of the ${\rm{DM}}{{\rm{M}}_2}$ flux Jacobian $\frac{\partial \mathbf{F}(\mathbf{u} )}{\partial \mathbf{u}}$ along ${{\phi }_{1}} = \frac{1}{2}\left({\phi _{2 + }} - \sqrt{{\phi _{2 - }}\, \left(1-{\phi _{2 + }}\right)} +\sqrt{{\phi _{2 + }}\, \left(1-{\phi _{2 - }}\right)} - {\phi _{2 - }}\right)$
Minimal and maximal distance between adjacent eigenvalues of the ${\rm{DM}}{{\rm{M}}_2}$ flux Jacobian $\frac{\partial \mathbf{F}(\mathbf{u})}{\partial \mathbf{u}}$ along the cut ${{\phi }_{1}} = \frac{1}{2}\left({\phi _{2 + }} - \sqrt{{\phi _{2 - }}\, \left(1-{\phi _{2 + }}\right)} +\sqrt{{\phi _{2 + }}\, \left(1-{\phi _{2 - }}\right)} - {\phi _{2 - }}\right)$
Eigenvalues at $5\%$ regularized boundary moments.
Online version: Press to activate 3D view ($x$-axis (red): ${\phi _{2 + }}$, $y$-axis (green): ${\phi _{2 - }}$, $z$-axis (blue): ${{\phi }_{1}}$)
Minimal eigenvalue distance for $5\%$ regularized boundary moments.
Online version: Press to activate 3D view ($x$-axis (red): ${\phi _{2 + }}$, $y$-axis (green): ${\phi _{2 - }}$, $z$-axis (blue): ${{\phi }_{1}}$)
Results for the plane-source test at the final time ${t_f} = 1$
Results for the source-beam test at the final time ${t_f} = 2.5$
DMM2 approximation of the non-differentiable MM2 ansatz (27)
Quantitative investigation of the minimal and maximal gaps of the eigenvalues at selected parts $S$ of the realizable set for fixed ${{\phi }_{1}} = \frac{1}{2}\left({\phi _{2 + }} - \sqrt{{\phi _{2 - }}\, \left(1-{\phi _{2 + }}\right)} +\sqrt{{\phi _{2 + }}\, \left(1-{\phi _{2 - }}\right)} - {\phi _{2 - }}\right)$
 $S$ $\min\limits_{\mathbf{u} \in S} \lambda _{\min}$ $\max\limits_{\mathbf{u} \in S} \lambda _{\max}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{min}}}}}\limits_{\mathbf{u} \in S}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{max}}}}}\limits_{\mathbf{u} \in S}$ $\min\limits_{\mathbf{u} \in S} \lambda _{1}$ $\max\limits_{\mathbf{u} \in S} \lambda _{1}$ $\min\limits_{\mathbf{u} \in S} \lambda _{4}$ $\max\limits_{\mathbf{u} \in S} \lambda _{4}$ $\mathcal{R}_\mathbf{b}{}$ 0.0260 1.9058 0.3287 1.0283 -0.9948 -0.1326 0.1326 0.9948 ${\phi _{2 + }}+{\phi _{2 - }}\geq \frac{19}{20}$ 0.0260 1.9058 0.0438 1.7811 -0.9948 -0.3778 0.3778 0.9948 ${\phi _{2 + }}\leq \frac{1}{120}$ 0.0314 1.2696 0.2066 0.5884 -0.9943 -0.1326 0.1326 0.4376 ${\phi _{2 - }}\leq \frac{1}{120}$ 0.0314 1.2696 0.2066 0.5884 -0.4376 -0.1326 0.1326 0.9943 ${\phi _{2 + }}={\phi _{2 - }}$ 0.0392 1.9058 0.3949 0.8791 -0.9921 -0.1326 0.1326 0.9921
 $S$ $\min\limits_{\mathbf{u} \in S} \lambda _{\min}$ $\max\limits_{\mathbf{u} \in S} \lambda _{\max}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{min}}}}}\limits_{\mathbf{u} \in S}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{max}}}}}\limits_{\mathbf{u} \in S}$ $\min\limits_{\mathbf{u} \in S} \lambda _{1}$ $\max\limits_{\mathbf{u} \in S} \lambda _{1}$ $\min\limits_{\mathbf{u} \in S} \lambda _{4}$ $\max\limits_{\mathbf{u} \in S} \lambda _{4}$ $\mathcal{R}_\mathbf{b}{}$ 0.0260 1.9058 0.3287 1.0283 -0.9948 -0.1326 0.1326 0.9948 ${\phi _{2 + }}+{\phi _{2 - }}\geq \frac{19}{20}$ 0.0260 1.9058 0.0438 1.7811 -0.9948 -0.3778 0.3778 0.9948 ${\phi _{2 + }}\leq \frac{1}{120}$ 0.0314 1.2696 0.2066 0.5884 -0.9943 -0.1326 0.1326 0.4376 ${\phi _{2 - }}\leq \frac{1}{120}$ 0.0314 1.2696 0.2066 0.5884 -0.4376 -0.1326 0.1326 0.9943 ${\phi _{2 + }}={\phi _{2 - }}$ 0.0392 1.9058 0.3949 0.8791 -0.9921 -0.1326 0.1326 0.9921
Quantitative investigation of the minimal and maximal gaps of the eigenvalues at selected parts $S$ of the realizable set for the $5\%$-regularized boundary moments
 $S$ $\min\limits_{\mathbf{u} \in S} \lambda _{\min}$ $\max\limits_{\mathbf{u} \in S} \lambda _{\max}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{min}}}}}\limits_{\mathbf{u} \in S}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{max}}}}}\limits_{\mathbf{u} \in S}$ $\min\limits_{\mathbf{u} \in S} \lambda _{1}$ $\max\limits_{\mathbf{u} \in S} \lambda _{1}$ $\min\limits_{\mathbf{u} \in S} \lambda _{4}$ $\max\limits_{\mathbf{u} \in S} \lambda _{4}$ $\mathcal{R}_\mathbf{b}{}$ 0.0334 1.8812 0.1921 1.0945 -0.9936 -0.2315 0.2315 0.9936 ${\phi _{2 + }}+{\phi _{2 - }}\geq \frac{19}{20}$ 0.0334 1.8812 0.0393 1.8611 -0.9936 -0.9709 0.9709 0.9936 ${\phi _{2 + }}\leq \frac{1}{120}$ 0.0603 1.4064 0.2291 0.8726 -0.9917 -0.2806 0.2315 0.9855 ${\phi _{2 - }}\leq \frac{1}{120}$ 0.0603 1.4064 0.2291 0.8726 -0.9855 -0.2315 0.2806 0.9917 ${\phi _{2 + }}={\phi _{2 - }}$ 0.0527 1.8138 0.1792 1.0526 -0.9897 -0.2789 0.2789 0.9897
 $S$ $\min\limits_{\mathbf{u} \in S} \lambda _{\min}$ $\max\limits_{\mathbf{u} \in S} \lambda _{\max}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{min}}}}}\limits_{\mathbf{u} \in S}$ $\mathop {{\rm{mean}}{\lambda _{{\rm{max}}}}}\limits_{\mathbf{u} \in S}$ $\min\limits_{\mathbf{u} \in S} \lambda _{1}$ $\max\limits_{\mathbf{u} \in S} \lambda _{1}$ $\min\limits_{\mathbf{u} \in S} \lambda _{4}$ $\max\limits_{\mathbf{u} \in S} \lambda _{4}$ $\mathcal{R}_\mathbf{b}{}$ 0.0334 1.8812 0.1921 1.0945 -0.9936 -0.2315 0.2315 0.9936 ${\phi _{2 + }}+{\phi _{2 - }}\geq \frac{19}{20}$ 0.0334 1.8812 0.0393 1.8611 -0.9936 -0.9709 0.9709 0.9936 ${\phi _{2 + }}\leq \frac{1}{120}$ 0.0603 1.4064 0.2291 0.8726 -0.9917 -0.2806 0.2315 0.9855 ${\phi _{2 - }}\leq \frac{1}{120}$ 0.0603 1.4064 0.2291 0.8726 -0.9855 -0.2315 0.2806 0.9917 ${\phi _{2 + }}={\phi _{2 - }}$ 0.0527 1.8138 0.1792 1.0526 -0.9897 -0.2789 0.2789 0.9897
 [1] Florian Schneider, Andreas Roth, Jochen Kall. First-order quarter-and mixed-moment realizability theory and Kershaw closures for a Fokker-Planck equation in two space dimensions. Kinetic & Related Models, 2017, 10 (4) : 1127-1161. doi: 10.3934/krm.2017044 [2] Michael Herty, Lorenzo Pareschi. Fokker-Planck asymptotics for traffic flow models. Kinetic & Related Models, 2010, 3 (1) : 165-179. doi: 10.3934/krm.2010.3.165 [3] Michael Herty, Christian Jörres, Albert N. Sandjo. Optimization of a model Fokker-Planck equation. Kinetic & Related Models, 2012, 5 (3) : 485-503. doi: 10.3934/krm.2012.5.485 [4] José Antonio Alcántara, Simone Calogero. On a relativistic Fokker-Planck equation in kinetic theory. Kinetic & Related Models, 2011, 4 (2) : 401-426. doi: 10.3934/krm.2011.4.401 [5] Sylvain De Moor, Luis Miguel Rodrigues, Julien Vovelle. Invariant measures for a stochastic Fokker-Planck equation. Kinetic & Related Models, 2018, 11 (2) : 357-395. doi: 10.3934/krm.2018017 [6] Marco Torregrossa, Giuseppe Toscani. On a Fokker-Planck equation for wealth distribution. Kinetic & Related Models, 2018, 11 (2) : 337-355. doi: 10.3934/krm.2018016 [7] Andreas Denner, Oliver Junge, Daniel Matthes. Computing coherent sets using the Fokker-Planck equation. Journal of Computational Dynamics, 2016, 3 (2) : 163-177. doi: 10.3934/jcd.2016008 [8] Ioannis Markou. Hydrodynamic limit for a Fokker-Planck equation with coefficients in Sobolev spaces. Networks & Heterogeneous Media, 2017, 12 (4) : 683-705. doi: 10.3934/nhm.2017028 [9] Giuseppe Toscani. A Rosenau-type approach to the approximation of the linear Fokker-Planck equation. Kinetic & Related Models, 2018, 11 (4) : 697-714. doi: 10.3934/krm.2018028 [10] John W. Barrett, Endre Süli. Existence of global weak solutions to Fokker-Planck and Navier-Stokes-Fokker-Planck equations in kinetic models of dilute polymers. Discrete & Continuous Dynamical Systems - S, 2010, 3 (3) : 371-408. doi: 10.3934/dcdss.2010.3.371 [11] Ludovic Dan Lemle. $L^1(R^d,dx)$-uniqueness of weak solutions for the Fokker-Planck equation associated with a class of Dirichlet operators. Electronic Research Announcements, 2008, 15: 65-70. doi: 10.3934/era.2008.15.65 [12] Roberta Bosi. Classical limit for linear and nonlinear quantum Fokker-Planck systems. Communications on Pure & Applied Analysis, 2009, 8 (3) : 845-870. doi: 10.3934/cpaa.2009.8.845 [13] Jessy Mallet, Stéphane Brull, Bruno Dubroca. General moment system for plasma physics based on minimum entropy principle. Kinetic & Related Models, 2015, 8 (3) : 533-558. doi: 10.3934/krm.2015.8.533 [14] Hyung Ju Hwang, Juhi Jang. On the Vlasov-Poisson-Fokker-Planck equation near Maxwellian. Discrete & Continuous Dynamical Systems - B, 2013, 18 (3) : 681-691. doi: 10.3934/dcdsb.2013.18.681 [15] Linjie Xiong, Tao Wang, Lusheng Wang. Global existence and decay of solutions to the Fokker-Planck-Boltzmann equation. Kinetic & Related Models, 2014, 7 (1) : 169-194. doi: 10.3934/krm.2014.7.169 [16] Florian Méhats, Olivier Pinaud. A problem of moment realizability in quantum statistical physics. Kinetic & Related Models, 2011, 4 (4) : 1143-1158. doi: 10.3934/krm.2011.4.1143 [17] Axel Klar, Florian Schneider, Oliver Tse. Approximate models for stochastic dynamic systems with velocities on the sphere and associated Fokker--Planck equations. Kinetic & Related Models, 2014, 7 (3) : 509-529. doi: 10.3934/krm.2014.7.509 [18] Florian Schneider, Jochen Kall, Graham Alldredge. A realizability-preserving high-order kinetic scheme using WENO reconstruction for entropy-based moment closures of linear kinetic equations in slab geometry. Kinetic & Related Models, 2016, 9 (1) : 193-215. doi: 10.3934/krm.2016.9.193 [19] Zbigniew Banach, Wieslaw Larecki. Entropy-based mixed three-moment description of fermionic radiation transport in slab and spherical geometries. Kinetic & Related Models, 2017, 10 (4) : 879-900. doi: 10.3934/krm.2017035 [20] Renjun Duan, Shuangqian Liu. Cauchy problem on the Vlasov-Fokker-Planck equation coupled with the compressible Euler equations through the friction force. Kinetic & Related Models, 2013, 6 (4) : 687-700. doi: 10.3934/krm.2013.6.687

2016 Impact Factor: 1.261

## Tools

Article outline

Figures and Tables