Related papers: A second-order numerical method for the aggregatio…
A quasi-second order scheme is developed to obtain approximate solutions of the shallow water equationswith bathymetry. The scheme is based on a staggered finite volume scheme for the space discretization:the scalar unknowns are located in…
The sum-of-squares method can give rigorous lower bounds on the energy of quantum Hamiltonians. Unfortunately, typically using this method requires solving a semidefinite program, which can be computationally expensive. Further, the…
We propose a second order finite volume scheme for nonlinear degenerate parabolic equations. For some of these models (porous media equation, drift-diffusion system for semiconductors, ...) it has been proved that the transient solution…
In this paper, we establish the optimal convergence result of a second order exponential-type integrator from (136, Numer. Math., 2017) for solving the KdV equation under rough initial data. The scheme is explicit and efficient to…
Formally second-order correct, mathematical descriptions of long-crested water waves propagating mainly in one direction are derived. These equations are analogous to the first-order approximations of KdV- or BBM-type. The advantage of…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
A general procedure to construct a class of simple and efficient high resolution Total Variation Diminishing (TVD) schemes for non-linear hyperbolic conservation laws by introducing anti-diffusive terms with the flux limiters is presented.…
Conservative symmetric second-order one-step schemes are derived for dynamical systems describing various many-body systems using the Discrete Multiplier Method. This includes conservative schemes for the $n$-species Lotka-Volterra system,…
The scaling of the exact solution of a hyperbolic balance law generates a family of scaled problems in which the source term does not depend on the current solution. These problems are used to construct a sequence of solutions whose…
A new second-order method for approximating the compressible Euler equations is introduced. The method preserves all the known invariant domains of the Euler system: positivity of the density, positivity of the internal energy and the local…
An "exact" method for scalar one-dimensional hyperbolic conservation laws is presented. The approach is based on the evolution of shock particles, separated by local similarity solutions. The numerical solution is defined everywhere, and is…
This study presents an efficient, accurate, effective and unconditionally stable time stepping scheme for the Darcy-Brinkman equations in double-diffusive convection. The stabilization within the proposed method uses the idea of stabilizing…
We focus in this work on the numerical discretization of the one dimensional aggregation equation $\pa_t\rho + \pa_x (v\rho)=0$, $v=a(W'*\rho)$, in the attractive case. Finite time blow up of smooth initial data occurs for potential $W$…
Successive quadratic approximations, or second-order proximal methods, are useful for minimizing functions that are a sum of a smooth part and a convex, possibly nonsmooth part that promotes regularization. Most analyses of iteration…
In this paper, we develop reliable a posteriori error estimates for numerical approximations of scalar hyperbolic conservation laws in one space dimension. Our methods have no inherent small-data limitations and are a step towards error…
A singularly perturbed linear system of second order ordinary differential equations of reaction-diffusion type with given boundary conditions is considered. The leading term of each equation is multiplied by a small positive parameter.…
We introduce a numerical method, based on finite elements and lattice gauge theory, to compute approximate solutions to Schr\"odinger and Pauli equations. The crucial geometric property of the method is discrete gauge invariance. The main…
The convergence analysis of a third-order scheme for the highly nonlinear Landau-Lifshitz-Gilbert equation with a non-convex constraint is considered. In this paper, we first present a fully discrete semi-implicit method for solving the…
Langevin simulation provides an effective way to study collisional effects in beams by reducing the six-dimensional Fokker-Planck equation to a group of stochastic ordinary differential equations. These resulting equations usually have…
We introduce new adaptive schemes for the one- and two-dimensional hyperbolic systems of conservation laws. Our schemes are based on an adaption strategy recently introduced in [{\sc S. Chu, A. Kurganov, and I. Menshov}, Appl. Numer. Math.,…