Adaptive change of basis in entropy-based moment closures for linear kinetic equations
From MaRDI portal
Publication:348610
DOI10.1016/J.JCP.2013.10.049zbMATH Open1349.82066arXiv1306.2881OpenAlexW1971599530MaRDI QIDQ348610FDOQ348610
Authors: Graham Alldredge, André L. Tits, Dianne P. O'Leary, Cory D. Hauck 
Publication date: 5 December 2016
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: Entropy-based (M_N) moment closures for kinetic equations are defined by a constrained optimization problem that must be solved at every point in a space-time mesh, making it important to solve these optimization problems accurately and efficiently. We present a complete and practical numerical algorithm for solving the dual problem in one-dimensional, slab geometries. The closure is only well-defined on the set of moments that are realizable from a positive underlying distribution, and as the boundary of the realizable set is approached, the dual problem becomes increasingly difficult to solve due to ill-conditioning of the Hessian matrix. To improve the condition number of the Hessian, we advocate the use of a change of polynomial basis, defined using a Cholesky factorization of the Hessian, that permits solution of problems nearer to the boundary of the realizable set. We also advocate a fixed quadrature scheme, rather than adaptive quadrature, since the latter introduces unnecessary expense and changes the computationally realizable set as the quadrature changes. For very ill-conditioned problems, we use regularization to make the optimization algorithm robust. We design a manufactured solution and demonstrate that the adaptive-basis optimization algorithm reduces the need for regularization. This is important since we also show that regularization slows, and even stalls, convergence of the numerical simulation when refining the space-time mesh. We also simulate two well-known benchmark problems. There we find that our adaptive-basis, fixed-quadrature algorithm uses less regularization than alternatives, although differences in the resulting numerical simulations are more sensitive to the regularization strategy than to the choice of basis.
Full work available at URL: https://arxiv.org/abs/1306.2881
Recommendations
- High-Order Entropy-Based Closures for Linear Transport in Slab Geometry II: A Computational Study of the Optimization Problem
- A regularized entropy-based moment method for kinetic equations
- Entropic quadrature for moment approximations of the Boltzmann-BGK equation
- Positive filtered \(\mathrm P_N\) moment closures for linear kinetic equations
- On the convergence of the regularized entropy-based moment method for kinetic equations
Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Transport processes in time-dependent statistical mechanics (82C70)
Cites Work
- Title not available (Why is that?)
- Title not available (Why is that?)
- Title not available (Why is that?)
- Title not available (Why is that?)
- Entropy optimization and mathematical programming
- Moment closure hierarchies for kinetic theories.
- Minimization of functions having Lipschitz continuous first partial derivatives
- High-Order Entropy-Based Closures for Linear Transport in Slab Geometry II: A Computational Study of the Optimization Problem
- A realizability-preserving discontinuous Galerkin method for the \(M_1\) model of radiative transfer
- Title not available (Why is that?)
- Etude théorique et numérique d'une hiérarchie de modèles aux moments pour le transfert radiatif
- Strong stability-preserving high-order time discretization methods
- Duality Relationships for Entropy-Like Minimization Problems
- High-order entropy-based closures for linear transport in slab geometry
- Dual Methods in Entropy Maximization. Application to Some Problems in Crystallography
- Domain of definition of Levermore's five-moment system
- MAXIMUM ENTROPY FOR REDUCED MOMENT PROBLEMS
- Entropic approximation in kinetic theory
- Convex Duality and Entropy-Based Moment Closures: Characterizing Degenerate Densities
- The multidimensional moment-constrained maximum entropy problem: A BFGS algorithm with constraint scaling
- Modified moments and Gaussian quadratures
- Heuristic solutions to polynomial moment problems with some convex entropic objectives
- When Not to Use an Automatic Quadrature Routine
- A Fast Heuristic Method for Polynomial Moment Problems with Boltzmann–Shannon Entropy
- The multidimensional maximum entropy moment problem: a review on numerical methods
Cited In (33)
- A realizability-preserving high-order kinetic scheme using WENO reconstruction for entropy-based moment closures of linear kinetic equations in slab geometry
- A consistent BGK model with velocity-dependent collision frequency for gas mixtures
- Approximating the \(M_2\) method by the extended quadrature method of moments for radiative transfer in slab geometry
- Wasserstein-penalized entropy closure: a use case for stochastic particle methods
- Machine Learning Moment Closure Models for the Radiative Transfer Equation II: Enforcing Global Hyperbolicity in Gradient-Based Closures
- Entropy–Based Methods for Uncertainty Quantification of Hyperbolic Conservation Laws
- A numerical comparison of the method of moments for the population balance equation
- On the Renormalization Maps for the φ -Divergence Moment Closures Applied in Radiative Transfer
- Gaussian process regression for maximum entropy distribution
- Maximum-principle-satisfying second-order Intrusive Polynomial Moment scheme
- Entropy-based moment closures for kinetic equations
- A positive and stable L2-minimization based moment method for the Boltzmann equation of gas dynamics
- Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency
- First-order quarter- and mixed-moment realizability theory and Kershaw closures for a Fokker-Planck equation in two space dimensions
- First-order continuous- and discontinuous-Galerkin moment models for a linear kinetic equation: realizability-preserving splitting scheme and numerical analysis
- Kershaw closures for linear transport equations in slab geometry. II: High-order realizability-preserving discontinuous-Galerkin schemes
- Discretizing distributions with exact moments: error estimate and convergence analysis
- Machine learning moment closure models for the radiative transfer equation. III: enforcing hyperbolicity and physical characteristic speeds
- TheM2Model for Dose Simulation in Radiation Therapy
- Optimization and large scale computation of an entropy-based moment closure
- Moment methods for the radiative transfer equation based on \(\varphi\)-divergences
- A comparison of moment closures for linear kinetic transport equations: the line source benchmark
- The 35-moment system with the maximum-entropy closure for rarefied gas flows
- Title not available (Why is that?)
- A realizability-preserving discontinuous Galerkin scheme for entropy-based moment closures for linear kinetic equations in one space dimension
- Efficient algorithms and implementations of entropy-based moment closures for rarefied gases
- First-order continuous- and discontinuous-Galerkin moment models for a linear kinetic equation: model derivation and realizability theory
- A regularized entropy-based moment method for kinetic equations
- An approximation of the \(M_2\) closure: application to radiotherapy dose simulation
- A new entropy-variable-based discretization method for minimum entropy moment approximations of linear kinetic equations
- Machine learning moment closure models for the radiative transfer equation. I: Directly learning a gradient based closure
- Some recent advances on the method of moments in kinetic theory
- Second-order mixed-moment model with differentiable ansatz function in slab geometry
This page was built for publication: Adaptive change of basis in entropy-based moment closures for linear kinetic equations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q348610)
