Solving hyperbolic-elliptic problems on singular mapped disk-like domains with the method of characteristics and spline finite elements
From MaRDI portal
Publication:2222580
Poisson-like equationpolar splinespseudo-Cartesian coordinatessemi-Lagrangiansingular domainsVlasov-like equation
Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) First-order hyperbolic equations (35L02) Numerical aspects of the method of characteristics for initial value and initial-boundary value problems involving PDEs (65M25) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Ionized gas flow in electromagnetic fields; plasmic flow (76X05)
Abstract: A common strategy in the numerical solution of partial differential equations is to define a uniform discretization of a tensor-product multi-dimensional logical domain, which is mapped to a physical domain through a given coordinate transformation. By extending this concept to a multi-patch setting, simple and efficient numerical algorithms can be employed on relatively complex geometries. The main drawback of such an approach is the inherent difficulty in dealing with singularities of the coordinate transformation. This work suggests a comprehensive numerical strategy for the common situation of disk-like domains with a singularity at a unique pole, where one edge of the rectangular logical domain collapses to one point of the physical domain (for example, a circle). We present robust numerical methods for the solution of Vlasov-like hyperbolic equations coupled to Poisson-like elliptic equations in such geometries. We describe a semi-Lagrangian advection solver that employs a novel set of coordinates, named pseudo-Cartesian coordinates, to integrate the characteristic equations in the whole domain, including the pole, and a finite element elliptic solver based on globally smooth splines (Toshniwal et al., 2017). The two solvers are tested both independently and on a coupled model, namely the 2D guiding-center model for magnetized plasmas, equivalent to a vorticity model for incompressible inviscid Euler fluids. The numerical methods presented show high-order convergence in the space discretization parameters, uniformly across the computational domain, without effects of order reduction due to the singularity. Dedicated tests show that the numerical techniques described can be applied straightforwardly also in the presence of point charges (equivalently, point-like vortices), within the context of particle-in-cell methods.
Recommendations
- A Finite-Element Method for Laplace- and Helmholtz-Type Boundary Value Problems with Singularities
- scientific article; zbMATH DE number 3963131
- Spline solution methods for a class of singular boundary value problems
- On the finite element method for a singular problem
- scientific article; zbMATH DE number 2051947
- scientific article; zbMATH DE number 749500
- Characteristic finite element methods for nonlinear Sobolev equations
- The method of auxiliary mapping for the finite element solutions of elliptic problems containing singularities
- Finite-element approximation for singularly perturbed nonlinear elliptic boundary-value problems
- The finite element method for a boundary value problem with strong singularity
Cites work
- A splitting scheme for the numerical solution of a one-dimensional Vlasov equation
- An Introduction to Numerical Analysis
- Arbitrarily high order convected scheme solution of the Vlasov-Poisson system
- Bézier surfaces and finite elements for MHD simulations
- Computation of MHD equilibrium of tokamak plasma
- Conservative numerical schemes for the Vlasov equation
- Conservative semi-Lagrangian schemes for Vlasov equations
- Geometric Numerical Integration
- High order multi-dimensional characteristics tracing for the incompressible Euler equation and the guiding-center Vlasov equation
- scientific article; zbMATH DE number 1197366 (Why is no real title available?)
- Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws
- Methods of conjugate gradients for solving linear systems
- Multi-degree smooth polar splines: a framework for geometric modeling and isogeometric analysis
- Semi-Lagrangian schemes for the Vlasov equation on an unstructured mesh of phase space
- Targeting realistic geometry in tokamak code Gysela
- The semi-Lagrangian method for the numerical resolution of the Vlasov equation
- Theory and simulations of two-dimensional vortex motion driven by a background vorticity gradient
- Vortex crystals from 2D Euler flow: Experiment and simulation
Cited in
(9)- Implicitly extrapolated geometric multigrid on disk-like domains for the gyrokinetic Poisson equation from fusion plasma applications
- 2D Fourier finite element formulation for magnetostatics in curvilinear coordinates with a symmetry direction
- Energy-minimizing, symmetric discretizations for anisotropic meshes and energy functional extrapolation
- Solver comparison for Poisson-like equations on tokamak geometries
- Local cubic spline interpolation for Vlasov-type equations on a multi-patch geometry.
- Memory -- and compute-optimized geometric multigrid \textit{GMGPolar} for curvilinear coordinate representations -- applications to fusion plasma
- Highly parallel drift-kinetic semi-Lagrangian simulations in Python
- Metriplectic relaxation to equilibria
- Polar spline discretizations with high-order regularity at the origin
This page was built for publication: Solving hyperbolic-elliptic problems on singular mapped disk-like domains with the method of characteristics and spline finite elements
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222580)
