Related papers: A structure-preserving implicit exponential time d…
Machine learning simulations of open quantum dynamics often rely on recursive predictors that accumulate error. We develop a non-recursive convolutional neural networks (CNNs) that maps system parameters and a redundant time encoding…
Lorentz invariant structure-preserving algorithms possess reference-independent secular stability, which is vital for simulating relativistic multi-scale dynamical processes. The splitting method has been widely used to construct…
An efficient numerical scheme for solving transport equations for tokamak plasmas within an integrated modelling framework is presented. The plasma transport equations are formulated as diffusion-advection equations in two coordinates (a…
This work presents a novel formulation and numerical strategy for the simulation of geometrically nonlinear structures. First, a non-canonical Hamiltonian (Poisson) formulation is introduced by including the dynamics of the stress tensor.…
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.…
Numerical solution of equations governing time domain simulations in computational electromagnetics, is usually based on grid methods in space and on explicit schemes for the time evolution. A predefined grid in the problem domain and a…
In this paper we propose a computational framework for the investigation of the correlated motion between positive and negative ions exposed to the attraction of a bubble surface that mimics the (oscillating) cell membrane. The correlated…
Stability is an important aspect of numerical methods for hyperbolic conservation laws and has received much interest. However, continuity in time is often assumed and only semidiscrete stability is studied. Thus, it is interesting to…
We develop a Macroscopic Auxiliary Asymptotic-Preserving Neural Network (MA-APNN) method to solve the time-dependent linear radiative transfer equations (LRTEs), which have a multi-scale nature and high dimensionality. To achieve this, we…
The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions and have applications in a wide variety of fields. Using an adaptive time-stepper based…
We explore the dynamics of active elements performing persistent random motion with fluctuating active speed and in the presence of translational noise in a $d$-dimensional harmonic trap, modeling active speed generation through an…
For the Maxwell's equations in a Havriliak-Negami (H-N) dispersive medium, the associated energy dissipation law has not been settled at both continuous level and discrete level. In this paper, we rigorously show that the energy of the H-N…
An energy preserving reduced order model is developed for two dimensional nonlinear Schr\"odinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty…
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…
In this paper, we introduce and analyse numerical schemes for the homogeneous and the kinetic L\'evy-Fokker-Planck equation. The discretizations are designed to preserve the main features of the continuous model such as conservation of…
The stability of classical semi-implicit scheme, and some more advanced iterative schemes recently proposed for Numerical Weather Prediction (NWP) purpose is examined. In all these schemes, the solution of the centred-implicit non-linear…
The paper focuses on the development of numerical methods for the compressible Euler equations. It is well-known that if the Mach number is small, the system becomes stiff and hence explicit schemes suffer from severe time-step…
In this paper, we utilize the maximum-principle-preserving flux limiting technique, originally designed for high order weighted essentially non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to develop a class of high…
This paper is concerned with the numerical approximation of the isothermal Euler equations for charged particles subject to the Lorentz force. When the magnetic field is large, the so-called drift-fluid approximation is obtained. In this…
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces…