Related papers: Positivity-preserving methods for population model…
High-order entropy-stable discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations require the positivity of thermodynamic quantities in order to guarantee their well-posedness. In this work, we introduce a…
Since almost twenty years, modified Patankar--Runge--Kutta (MPRK) methods have proven to be efficient and robust numerical schemes that preserve positivity and conservativity of the production-destruction system irrespectively of the time…
We develop a new method for equality constrained optimization problems based on a sequential cubic programming framework. Each iteration utilizes a step decomposition based on the Jacobian of the constraints into a normal and a tangential…
Optimal Strong Stability Preserving (SSP) Runge--Kutta methods has been widely investegated in the last decade and many open conjectures have been formulated. The iterated implicit midpoint rule has been observed numerically optimal in…
In this paper, we are concerned with arbitrarily high-order momentum-preserving and energy-preserving schemes for solving the generalized Rosenau-type equation, respectively. The derivation of the momentum-preserving schemes is made within…
Ordinary differential equations (ODEs) and ordinary difference systems (O$\Delta$Ss) invariant under the actions of the Lie groups $\mathrm{SL}_x(2)$, $\mathrm{SL}_y(2)$ and $\mathrm{SL}_x(2)\times\mathrm{SL}_y(2)$ of projective…
In this paper, we present two stochastic positive-preserving symplectic methods for the stochastic Lotka-Volterra predator-prey model driven by a multiplicative noise. To inherit the intrinsic characteristic of the original system, the…
In this paper, we are interested in constructing a scheme solving compressible Navier--Stokes equations, with desired properties including high order spatial accuracy, conservation, and positivity-preserving of density and internal energy…
This survey reviews variational and iterative methods for reconstructing non-negative solutions of ill-posed problems in infinite-dimensional spaces. We focus on two classes of methods: variational methods based on entropy-minimization or…
We further develop a simple modification of Runge--Kutta methods that guarantees conservation or stability with respect to any inner-product norm. The modified methods can be explicit and retain the accuracy and stability properties of the…
Stochastic Hamiltonian partial differential equations, which possess the multi-symplectic conservation law, are an important and fairly large class of systems. The multi-symplectic methods inheriting the geometric features of stochastic…
Efficient high order numerical methods for evolving the solution of an ordinary differential equation are widely used. The popular Runge--Kutta methods, linear multi-step methods, and more broadly general linear methods, all have a global…
This work considers multirate generalized-structure additively partitioned Runge-Kutta (MrGARK) methods for solving stiff systems of ordinary differential equations (ODEs) with multiple time scales. These methods treat different partitions…
Inverse problems arise in situations where data is available, but the underlying model is not. It can therefore be necessary to infer the parameters of the latter starting from the former. Statistical mechanics offers a toolbox of…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
We design high order accurate methods that exploit low rank structure in the density matrix while respecting the essential structure of the Lindblad equation. Our methods preserves complete positivity and are trace preserving.
Classical convergence theory of Runge-Kutta methods assumes that the time step is small relative to the Lipschitz constant of the ordinary differential equation (ODE). For stiff problems, that assumption is often violated, and a problematic…
Differentially private mechanisms are increasingly used to publish tables of counts, where each entry represents the number of individuals belonging to a particular category. A distribution of counts summarizes the information in the count…
We propose a novel non-compact, positivity-preserving scheme for linear non-divergence form elliptic equations. Based on the Feynman--Kac formula, the solution is represented as a conditional expectation associated with a diffusion…
We derive a second-order realizability-preserving scheme for moment models for linear kinetic equations. We apply this scheme to the first-order continuous and discontinuous models in slab and three-dimensional geometry derived in a…