Related papers: Divergence preservation in the ADI algorithms for …
Beyond the choice of wavefront control systems or coronographs, advanced data processing methods play a crucial role in disentangling potential planetary signals from bright quasi-static speckles. Among these methods, angular differential…
Divergence and vorticity damping, which operate upon horizontal divergence and relative vorticity, are explicit diffusion mechanisms used in dynamical cores to ensure stability. To avoid numerical blow-up from excessively strong diffusion,…
In this work, we present a positivity-preserving entropy-based adaptive filtering method for shock capturing in discontinuous spectral element methods. By adapting the filter strength to enforce positivity and a local discrete minimum…
The presence of global conserved quantities in interacting systems generically leads to diffusive transport at late times. Here, we show that systems conserving the dipole moment of an associated global charge, or even higher moment…
In this paper, by using Strang's second-order splitting method, the numerical procedure for the three-dimensional (3D) space fractional Allen-Cahn equation can be divided into three steps. The first and third steps involve an ordinary…
Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design.…
A thermodynamically consistent particle-based model for fluid dynamics with continuous velocities and a non-ideal equation of state is presented. Excluded volume interactions are modeled by means of biased stochastic multiparticle…
Three dimensional computations of self consistent three species gyrofluid turbulence are carried out for tokamak edge conditions. Profiles as well as disturbances in dependent variables are followed, running the dynamical system to…
We prove that adapted entropy solutions of scalar conservation laws with discontinuous flux are stable with respect to changes in the flux under the assumption that the flux is strictly monotone in u and the spatial dependency is piecewise…
Charge-transfer excited states are highly relevant for applications in molecular electronics. However, the accurate calculation of these states in large systems is challenging since wave function methods are prohibitively expensive,…
A direct numerical simulation of the three-dimensional elektrokinetic instability near a charge selective surface (electric membrane, electrode, or system of micro-/nanochannels) is carried out and analyzed. A special finite-difference…
In this article we describe a stable partitioned algorithm that overcomes the added mass instability arising in fluid-structure interactions of light rigid bodies and inviscid compressible flow. The new algorithm is stable even for bodies…
This article presents an entropy stability analysis of smoothed dissipative particle dynamics (SDPD) to review the validity of particle discretization of entropy equations. First, we consider the simplest SDPD system: a simulation of…
In this manuscript, we present the development of implicit and implicit-explicit ADER and DeC methodologies within the DeC framework using the two-operators formulation, with a focus on their stability analysis both as solvers for ordinary…
We present an accurate and efficient algorithm to calculate the electrostatic interaction of charged point particles with partially periodic boundary conditions that are confined along the nonperiodic direction by two metallic parallel…
In the present work, a novel particle merging scheme is proposed for PIC-DSMC simulations, based on the solution of a Non-negative Least Squares problem. The merging algorithm conserves arbitrary moments of the velocity distribution…
The paper studies numerical methods that preserve a Lyapunov function of a dynamical system, i.e. numerical approximations whose energy decreases, just like in the original differential equation. With this aim, a discrete gradient method is…
This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal…
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…
A consistency/latency tradeoff arises as soon as a distributed storage system replicates data. For low latency, modern storage systems often settle for weak consistency conditions, which provide little, or even worse, no guarantee for data…