Related papers: A structure-preserving implicit exponential time d…
The density matrix is a widely used tool in quantum mechanics. In order to determine its evolution with respect to time, the Liouville-von Neumann equation must be solved. However, analytic solutions of this differential equation exist only…
We propose a fully decoupled, structure-preserving relaxation Crank--Nicolson finite element method (FEM) for the coupled Gross--Pitaevskii--Poisson (GPP) system modeling ultracold plasmas. By introducing suitable auxiliary variables to…
In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it…
We investigate a two-state conformational conversion system and introduce a novel structure-preserving numerical scheme that couples a local discontinuous Galerkin space discretization with the backward Euler time-integration method. The…
The paper is devoted to a constitutive solution, limit load analysis and Newton-like methods in elastoplastic problems containing the Mohr-Coulomb yield criterion. Within the constitutive problem, we introduce a self-contained derivation of…
In this paper we consider unconditionally energy stable numerical schemes for the nonstationary 3D magneto-micropolar equations that describes the microstructure of rigid microelements in electrically conducting fluid flow under some…
We propose a new pressure-based structure-preserving (SP) and quasi asymptotic preserving (AP) staggered semi-implicit finite volume scheme for the unified first order hyperbolic formulation of continuum mechanics. The unified model is…
In this work, we carry out a study of the conserved quantities and dynamical structure of the four-dimensional modified axion electrodynamics theory described by the axion-photon coupling. In the first part of the analysis, we employ the…
We consider coupled models for particulate flows, where the disperse phase is made of particles with distinct sizes. We are thus led to a system coupling the incompressible Navier-Stokes equations to the multi-component Vlasov-Fokker-Planck…
Techniques are developed for decoupling dissipative differential equations. The approach considered is based upon obtaining a sufficient gap in the time dependent linear portion of the equation that corresponds to the linear variational…
We develop structure preserving schemes for a class of nonlinear mobility continuity equation. When the mobility is a concave function, this equation admits a form of gradient flow with respect to a Wasserstein-like transport metric. Our…
We develop a structure-preserving computational framework for optimal mixing control in incompressible flows. Our approach exactly conserves the continuous system's key invariants (mass and $L^2$-energy), while also maintaining discrete…
Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure…
We present a finite volume scheme for modeling the diffusion of charged particles, specifically ions, in constrained geometries using a degenerate Poisson-Nernst-Planck system with size exclusion yielding cross-diffusion. Our method…
In this paper, we investigate the use of a mass lumped fully explicit time stepping scheme for the discretisation of the wave equation with underlying material parameters that vary at arbitrarily fine scales. We combine the leapfrog scheme…
This research demonstrates analytical time-dependent non-equilibrium green function (TD-NEGF) algorithms to investigate dynamical functionalities of quantum devices, especially for photon-assisted transports. Together with the lumped…
We propose an arbitrarily higher (even) order implicit leapfrog scheme for time discretization of a three-field formulation of Maxwell's equations. We use this in conjunction with an arbitrarily higher-order and compatible discretization…
In this paper we develop a neural network for the numerical simulation of time-dependent linear transport equations with diffusive scaling and uncertainties. The goal of the network is to resolve the computational challenges of…
Two fast L1 time-stepping methods, including the backward Euler and stabilized semi-implicit schemes, are suggested for the time-fractional Allen-Cahn equation with Caputo's derivative. The time mesh is refined near the initial time to…
A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions.…