# American Institute of Mathematical Sciences

December  2017, 10(4): 1127-1161. doi: 10.3934/krm.2017044

## First-order quarter-and mixed-moment realizability theory and Kershaw closures for a Fokker-Planck equation in two space dimensions

 Technische Universität Kaiserslautern, 67663 Kaiserslautern, Germany

* Corresponding author: Florian Schneider

Received  August 2016 Revised  November 2016 Published  March 2017

Mixed-moment models, introduced in [8,44] for one space dimension, are a modification of the method of moments applied to a (linear) kinetic equation, by choosing mixtures of different partial moments. They are well-suited to handle equations where collisions of particles are modelled with a Laplace-Beltrami operator. We generalize the concept of mixed moments to two dimensions. In the context of minimum-entropy models, the resulting hyperbolic system of equations has desirable properties (entropy-diminishing, bounded eigenvalues), removing some drawbacks of the well-known M1 model. We furthermore provide a realizability theory for a first-order system of mixed moments by linking it to the corresponding quarter-moment theory. Additionally, we derive a type of Kershaw closures for mixed-and quarter-moment models, giving an efficient closure (compared to minimum-entropy models). The derived closures are investigated for different benchmark problems.

Citation: 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
##### References:

show all references

##### References:
Interpolation between two realizability boundaries in the projected three-space $(\phi _{\mathcal{S}_x^ + }^{\left( {1,0} \right)},\phi _{\mathcal{S}_x^ - }^{\left( {1,0} \right)},\phi _{\mathcal{S}_y^ + }^{\left( {0,1} \right)})$ along an isoline of $\mathcal{N}\left( {{\phi ^{\left| 1 \right|}}} \right)$. The realizable set with respect to the quadrant ${\mathcal{S}_{ - + }}$ is plotted in grey
The eigenvector ${\bf{v}}_1({\phi _{{\mathcal{S}_{ + + }}}})$ (blue arrows) of the ${\rm{Q}}{{\rm{M}}_1}$ second moment $\phi _{{\mathcal{S}_{ + + }}}^{\left| 2 \right|}$ (Figure (a)) and its minimal deviation from $\phi _{{\mathcal{S}_{+ +}}}^{\left| 1 \right|}$ in degrees (Figure (b))
Components of the $Q{K_1}$ second moment on ${\mathcal{S}_{ + + }}$
Deviation of $\phi _{{\mathcal{S}_{ + + }}}^{\left( {2,0} \right)}$ and $\phi _{{\mathcal{S}_{ + + }}}^{\left( {2,0} \right)}$ for ${\rm{Q}}{{\rm{M}}_1}$ and $Q{K_1}$, respectively
Eigenvalues of the flux Jacobian in $x$-direction ${\bf{F}}_1^\prime$
Minimal and maximal distance of adjacent eigenvalues of the flux Jacobian ${\bf{F}}_1^\prime$
Visualization of the coefficients $\rho _{\iota \pm}$ for the example $\phi _{{\mathcal{S}_{+ +}}}^{\left| 1 \right|}=\left(0.4,0.6\right)^T$. The length of the arrows represent the maximal value for the respective coefficient
Second moment $\phi _{\mathcal{S}_y^ + }^{\left( {0,2} \right)}$ for the ${\rm{M}}{{\rm{M}}_1}$ and Kershaw model for $\phi _{\mathcal{S}_y^ + }^{\left( {0,1} \right)} = -\phi _{\mathcal{S}_y^ - }^{\left( {0,1} \right)} = \frac14$
Local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 1.2$ in the two-beams test case
Local particle density ${u^{\left( {0,0} \right)}}$ at $t = 1.2$ in the opposing two-beams test case
Local particle density ${u^{\left( {0,0} \right)}}$ at $t = 1.2$ in the rotated two-beams test case
Local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 1$ in the source-beam test case
Local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 2.5$ in the source-beam test case
Horizontal cut ($y = 1.5$) of the local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 1$ and $t = 2.5$, respectively, in the source-beam test case
Local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 0.45$ in the line-source test case
Local particle density ${u^{\left( {0,0} \right)}}$ of some models at $t = 0.45$ in the line-source test case, horizontal and diagonal cuts
Absolute and relative running times (wrt. ${\rm{M}}{{\rm{M}}_1}$) of some models for the source-beam test
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 1.305e+05 3.885e-01 ${\rm{M}}{{\rm{M}}_2}$ 13 3.359e+05 1.0 ${\rm{M}}{{\rm{K}}_1}$ 5 3.408e+01 1.015e-04 ${{\rm{M}}_1}$ 3 9.418e+04 2.804e-01 ${{\rm{M}}_2}$ 6 2.269e+05 6.757e-01 $Q{K_1}$ 12 5.309e+01 1.581e-04
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 1.305e+05 3.885e-01 ${\rm{M}}{{\rm{M}}_2}$ 13 3.359e+05 1.0 ${\rm{M}}{{\rm{K}}_1}$ 5 3.408e+01 1.015e-04 ${{\rm{M}}_1}$ 3 9.418e+04 2.804e-01 ${{\rm{M}}_2}$ 6 2.269e+05 6.757e-01 $Q{K_1}$ 12 5.309e+01 1.581e-04
Absolute and relative running times (wrt. ${\rm{M}}{{\rm{M}}_1}$) of some models for the source-beam test.
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 2.489e+05 1.0 ${\rm{M}}{{\rm{M}}_2}$ 13 9.526e+04 3.828e-01 ${\rm{M}}{{\rm{K}}_1}$ poly 5 1.535e+02 6.169e-04 ${{\rm{M}}_1}$ 3 6.639e+04 2.668e-01 ${{\rm{M}}_2}$ 6 1.189e+05 4.777e-01 ${{\rm{M}}_3}$ 12 2.232e+05 8.967e-01
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 2.489e+05 1.0 ${\rm{M}}{{\rm{M}}_2}$ 13 9.526e+04 3.828e-01 ${\rm{M}}{{\rm{K}}_1}$ poly 5 1.535e+02 6.169e-04 ${{\rm{M}}_1}$ 3 6.639e+04 2.668e-01 ${{\rm{M}}_2}$ 6 1.189e+05 4.777e-01 ${{\rm{M}}_3}$ 12 2.232e+05 8.967e-01
Absolute and relative running times (wrt. ${\rm{M}}{{\rm{M}}_3}$) of some models for the line-source test
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 1.638e+05 2.091e-01 ${\rm{M}}{{\rm{M}}_2}$ 13 3.491e+05 4.456e-01 ${\rm{M}}{{\rm{M}}_3}$ 25 7.834e+05 1.0 ${\rm{M}}{{\rm{K}}_1}$ poly 5 1.888e+02 2.409e-04 ${\rm{M}}{{\rm{K}}_1}$ 5 2.898e+04 3.700e-02 ${{\rm{M}}_1}$ 3 5.458e+04 6.968e-02 ${{\rm{M}}_2}$ 6 9.696e+04 1.238e-01 ${{\rm{M}}_3}$ 10 1.823e+05 2.327e-01
 Model $n$ Abs. runtime [s] Rel. runtime ${\rm{M}}{{\rm{M}}_1}$ 5 1.638e+05 2.091e-01 ${\rm{M}}{{\rm{M}}_2}$ 13 3.491e+05 4.456e-01 ${\rm{M}}{{\rm{M}}_3}$ 25 7.834e+05 1.0 ${\rm{M}}{{\rm{K}}_1}$ poly 5 1.888e+02 2.409e-04 ${\rm{M}}{{\rm{K}}_1}$ 5 2.898e+04 3.700e-02 ${{\rm{M}}_1}$ 3 5.458e+04 6.968e-02 ${{\rm{M}}_2}$ 6 9.696e+04 1.238e-01 ${{\rm{M}}_3}$ 10 1.823e+05 2.327e-01
 [1] A. M. Micheletti, Angela Pistoia. Multiple eigenvalues of the Laplace-Beltrami operator and deformation of the Riemannian metric. Discrete & Continuous Dynamical Systems - A, 1998, 4 (4) : 709-720. doi: 10.3934/dcds.1998.4.709 [2] Micol Amar, Roberto Gianni. Laplace-Beltrami operator for the heat conduction in polymer coating of electronic devices. Discrete & Continuous Dynamical Systems - B, 2018, 23 (4) : 1739-1756. doi: 10.3934/dcdsb.2018078 [3] Andrea Bonito, Roland Glowinski. On the nodal set of the eigenfunctions of the Laplace-Beltrami operator for bounded surfaces in $R^3$: A computational approach. Communications on Pure & Applied Analysis, 2014, 13 (5) : 2115-2126. doi: 10.3934/cpaa.2014.13.2115 [4] 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 [5] Josephus Hulshof, Pascal Noble. Travelling waves for a combustion model coupled with hyperbolic radiation moment models. Discrete & Continuous Dynamical Systems - B, 2008, 10 (1) : 73-90. doi: 10.3934/dcdsb.2008.10.73 [6] 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 [7] Filippo Gazzola. On the moments of solutions to linear parabolic equations involving the biharmonic operator. Discrete & Continuous Dynamical Systems - A, 2013, 33 (8) : 3583-3597. doi: 10.3934/dcds.2013.33.3583 [8] Yong-Jung Kim. A generalization of the moment problem to a complex measure space and an approximation technique using backward moments. Discrete & Continuous Dynamical Systems - A, 2011, 30 (1) : 187-207. doi: 10.3934/dcds.2011.30.187 [9] T. Hillen. On the $L^2$-moment closure of transport equations: The general case. Discrete & Continuous Dynamical Systems - B, 2005, 5 (2) : 299-318. doi: 10.3934/dcdsb.2005.5.299 [10] T. Hillen. On the $L^2$-moment closure of transport equations: The Cattaneo approximation. Discrete & Continuous Dynamical Systems - B, 2004, 4 (4) : 961-982. doi: 10.3934/dcdsb.2004.4.961 [11] John C. Schotland, Vadim A. Markel. Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation. Inverse Problems & Imaging, 2007, 1 (1) : 181-188. doi: 10.3934/ipi.2007.1.181 [12] Shingo Takeuchi. Partial flat core properties associated to the $p$-laplace operator. Conference Publications, 2007, 2007 (Special) : 965-973. doi: 10.3934/proc.2007.2007.965 [13] Matthias Geissert, Horst Heck, Christof Trunk. $H^{\infty}$-calculus for a system of Laplace operators with mixed order boundary conditions. Discrete & Continuous Dynamical Systems - S, 2013, 6 (5) : 1259-1275. doi: 10.3934/dcdss.2013.6.1259 [14] 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 [15] 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 [16] Guoqing Zhang, Jia-yu Shao, Sanyang Liu. Linking solutions for N-laplace elliptic equations with Hardy-Sobolev operator and indefinite weights. Communications on Pure & Applied Analysis, 2011, 10 (2) : 571-581. doi: 10.3934/cpaa.2011.10.571 [17] Peter I. Kogut, Olha P. Kupenko. On optimal control problem for an ill-posed strongly nonlinear elliptic equation with $p$-Laplace operator and $L^1$-type of nonlinearity. Discrete & Continuous Dynamical Systems - B, 2019, 24 (3) : 1273-1295. doi: 10.3934/dcdsb.2019016 [18] Julián López-Gómez. On the structure of the permanence region for competing species models with general diffusivities and transport effects. Discrete & Continuous Dynamical Systems - A, 1996, 2 (4) : 525-542. doi: 10.3934/dcds.1996.2.525 [19] Pengwen Chen, Changfeng Gui. Alpha divergences based mass transport models for image matching problems. Inverse Problems & Imaging, 2011, 5 (3) : 551-590. doi: 10.3934/ipi.2011.5.551 [20] Martin Burger, Jan-Frederik Pietschmann, Marie-Therese Wolfram. Identification of nonlinearities in transport-diffusion models of crowded motion. Inverse Problems & Imaging, 2013, 7 (4) : 1157-1182. doi: 10.3934/ipi.2013.7.1157

2018 Impact Factor: 1.38