Abstract:Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. Combined with explicit time integration schemes, the finite cell method introduces a lower bound for the critical time step size. Explicit transient analyses commonly use the spectral element method due to its natural way of obtaining diagonal mass matrices through nodal lumping. Its combination with the finite cell method is called the spectral cell method. Unfortunately, a direct application of nodal lumping in the spectral cell method is impossible due to the special quadrature necessary to treat the discontinuous integrand inside the cut cells. We analyze an implicit-explicit (IMEX) time integration method to exploit the advantages of the nodal lumping scheme for uncut cells on one side and the unconditional stability of implicit time integration schemes for cut cells on the other. In this hybrid, immersed Newmark IMEX approach, we use explicit second-order central differences to integrate the uncut degrees of freedom that lead to a diagonal block in the mass matrix and an implicit trapezoidal Newmark method to integrate the remaining degrees of freedom (those supported by at least one cut cell). The immersed Newmark IMEX approach preserves the high-order convergence rates and the geometric flexibility of the finite cell method. We analyze a simple system of spring-coupled masses to highlight some of the essential characteristics of Newmark IMEX time integration. We then solve the scalar wave equation on two- and three-dimensional examples with significant geometric complexity to show that our approach is more efficient than state-of-the-art time integration schemes when comparing accuracy and runtime.

An efficient implicit boundary integral solver for the Vlasov-Maxwell system

Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems

A novel implicit solver for the Vlasov-Maxwell system

A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations

Conservative semi-Lagrangian kinetic scheme coupled with implicit finite element field solver for multidimensional Vlasov Maxwell system

Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

Canonical Symplectic Particle-In-Cell Method For Long-Term Large-Scale Simulations Of The Vlasov-Maxwell Equations

Implicit and Coupled Multi-Fluid Solver for Collisional Low-Temperature Plasma

A high-order integral equation-based solver for the time-dependent Schrodinger equation

Discontinuous Galerkin Methods for the Vlasov-Maxwell Equations

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

Kinetic Solvers with Adaptive Mesh in Phase Space for Low-Temperature Plasmas

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

A discontinuous Galerkin method for the Vlasov-Poisson system

An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations

Implicit-Explicit Time Integration for the Immersed Wave Equation

A 3D Nonlinear Maxwell's Equations Solver Based On A Hybrid Numerical Method