Related papers: Metriplectic particle-in-cell integrators for the …
When underlying probability density functions of nonlinear dynamic systems are unknown, the filtering problem is known to be a challenging problem. This paper attempts to make progress on this problem by proposing a new class of filtering…
Over the last decades, various "non-linear" MCMC methods have arisen. While appealing for their convergence speed and efficiency, their practical implementation and theoretical study remain challenging. In this paper, we introduce a…
Probabilistic cellular automata (PCA) are used to model a variety of discrete spatially extended systems undergoing parallel-updating. We propose an embedding of a number of classical nonequilibrium concepts in the PCA-world. We start from…
We address the approximation of functionals depending on a system of particles, described by stochastic differential equations (SDEs), in the mean-field limit when the number of particles approaches infinity. This problem is equivalent to…
In purely non-dissipative systems, Lagrangian and Hamiltonian reduction have proven to be powerful tools for deriving physical models with exact conservation laws. We have discovered a hint that an analogous reduction method exists also for…
We show convergence of a cell-centered finite volume discretization for linear elasticity. The discretization, termed the MPSA method, was recently proposed in the context of geological applications, where cell-centered variables are often…
A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve…
Starting from a particle model describing self-propelled particles interacting through nematic alignment, we derive a macroscopic model for the particle density and mean direction of motion. We first propose a mean-field kinetic model of…
In multibody simulation, the Gear-Gupta-Leimkuhler method for only persistent contacts enforces constraints on position and velocity level at the same time. It yields a robust numerical discretization of differential algebraic equations…
We propose implicit integrators for solving stiff differential equations on unit spheres. Our approach extends the standard backward Euler and Crank-Nicolson methods in Cartesian space by incorporating the geometric constraint inherent to…
We present a loosely coupled approach for the solution of fluid-structure interaction problems between a compressible flow and a deformable structure. The method is based on staggered Dirichlet-Neumann partitioning. The interface motion in…
In recent years, two important techniques for geometric numerical discretization have been developed. In computational electromagnetics, spatial discretization has been improved by the use of mixed finite elements and discrete differential…
Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given…
We explore the possibilities of applying structure-preserving numerical methods to a plasma hybrid model with kinetic ions and mass-less fluid electrons satisfying the quasi-neutrality relation. The numerical schemes are derived by finite…
In this work, we develop an adaptive, multivariate partitioning algorithm for solving mixed-integer nonlinear programs (MINLP) with multi-linear terms to global optimality. This iterative algorithm primarily exploits the advantages of…
We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a…
In this paper, we introduce MPM Lite, a new hybrid Lagrangian/Eulerian method that eliminates the need for particle-based quadrature at solve time. Standard MPM practices suffer from a performance bottleneck where expensive implicit solves…
On the basis of input-output time-domain data collected from a complex simulator, this paper proposes a constructive methodology to infer a reduced-order linear, bilinear or quadratic time invariant dynamical model reproducing the…
We present a geometric Particle-in-Cell (PIC) algorithm on two-dimensional (2D) unstructured meshes for studying electrostatic perturbations in magnetized plasmas. In this method, ions are treated as fully kinetic particles, and electrons…
We present and analyze a structure-preserving method for the approximation of solutions to nonlinear cross-diffusion systems, which combines a Local Discontinuous Galerkin spatial discretization with the backward Euler time-stepping scheme.…