Related papers: Divergence preservation in the ADI algorithms for …
Fourth-order accurate compact schemes for variable coefficient convection diffusion equations are considered. A sufficient condition for the stability of the fully discrete problem is derived using a difference equation based approach. The…
Many modern discontinuous Galerkin (DG) methods for conservation laws make use of summation by parts operators and flux differencing to achieve kinetic energy preservation or entropy stability. While these techniques increase the robustness…
We emphasize that the stabilizing symmetry for dark matter (DM) particles does not have to be the commonly used parity (Z_2) symmetry. We therefore examine the potential of the colliders to distinguish models with parity stabilized DM from…
We highlight a striking difference in behavior between two widely used variants of coordinate ascent variational inference: the sequential and parallel algorithms. While such differences were known in the numerical analysis literature in…
Atomistic/continuum coupling methods aim to achieve optimal balance between accuracy and efficiency. Adaptivity is the key for the efficient implementation of such methods. In this paper, we carry out a rigorous a posteriori analysis of the…
We consider a class of multiscale parabolic problems with diffusion coefficients oscillating in space at a possibly small scale $\varepsilon$. Numerical homogenization methods are popular for such problems, because they capture efficiently…
Second-order Volume-preserving algorithms (VPAs) for simulating charged particle motion in electromagnetic fields have been generalized to a rotating angle formulation by using the matrix decomposition methods. Based on this method, the…
An initial-boundary value problem for the 1D self-adjoint parabolic equation on the half-axis is solved. We study a broad family of two-level finite-difference schemes with two parameters related to averagings both in time and space.…
The aim of this paper is to develop a stability theory and a numerical model for the three density-stratified electrically conductive liquid layers. Using regular perturbation methods to reduce the full 3d problem to the shallow layer…
We note that in extensions of the Standard Model that allow for a varying fine structure constant, alpha, all matter species, apart from right-handed neutrinos, will gain an intrinsic electric dipole moment (EDM). In a large subset of…
We present and analyze two stabilized finite element methods for solving numerically the Poisson--Nernst--Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for…
A key characteristic of the anomalous sub-solution equation is that the solution exhibits algebraic decay rate over long time intervals, which is often refered to the Mittag-Leffler type stability. For a class of power nonlinear…
A novel Douglas alternating direction implicit (ADI) method is proposed in this work to solve a two-dimensional (2D) heat equation with interfaces. The ADI scheme is a powerful finite difference method for solving parabolic equations, due…
Structural properties of a finite number ($N = 2 - 20$) of point charges (classical electrons) confined laterally in a two-dimensional two-minima potential are calculated as a function of the distance ($d$) between the minima. The particles…
The present article is devoted to well-posedness by noise for the continuity equation. Namely, we consider the continuity equation with non-linear and partially degenerate stochastic perturbations in divergence form. We prove the existence…
We propose two algorithms for discrete-time parameter estimation, one for time-varying parameters under persistent excitation (PE) condition, another for constant parameters under no PE condition. For the first algorithm, we show that in…
Perturbation theory is an indispensable tool in quantum mechanics and electrodynamics that handles weak effects on particle motion or fields. However, its extension to plasmons involving complex motion of {\it both} particles and fields…
There is a qualitative difference between one-dimensional and multi-dimensional solutions to the Euler equations: new features that arise are vorticity and a nontrivial incompressible (low Mach number) limit. They present challenges to…
This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A…
We introduce a novel method to investigate the stability of wave packet dynamics under perturbations of the Hamiltonian. Our approach relies on semiclassical approximations, but is non-perturbative. Two separate contributions to the quantum…