Related papers: A Blob Method for the Aggregation Equation
Entropy stabilization of the compressible Euler system is achieved by adapting the averages that are applied to the density and internal energy variables. The approach achieves non-linear robustness despite the use of simplified symmetric…
We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general…
Most fluid flow problems that are vital in engineering applications involve at least one of the following features: turbulence, shocks, and/or material interfaces. While seemingly different phenomena, these flows all share continuous…
We propose a new method to accelerate the convergence of optimization algorithms. This method simply adds a power coefficient $\gamma\in[0,1)$ to the gradient during optimization. We call this the Powerball method and analyze the…
In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations.…
We present a subdivision method to solve systems of congruence equations. This method is inspired in a subdivision method, based on Bernstein forms, to solve systems of polynomial inequalities in several variables and arbitrary degrees. The…
The evaluation of the electrostatic potential is fundamental to the study of condensed phase systems. We discuss the calculation of the relevant lattice summations by Ewald-type techniques. A model charge density is introduced, that cancels…
We study the strong approximation of a rough volatility model, in which the log-volatility is given by a fractional Ornstein-Uhlenbeck process with Hurst parameter $H<1/2$. Our methods are based on an equidistant discretization of the…
Although it is relatively easy to apply, the gradient method often displays a disappointingly slow rate of convergence. Its convergence is specially based on the structure of the matrix of the algebraic linear system, and on the choice of…
Given a large ensemble of interacting particles, driven by nonlocal interactions and localized repulsion, the mean-field limit leads to a class of nonlocal, nonlinear partial differential equations known as aggregation-diffusion equations.…
We introduce vertex block descent, a block coordinate descent solution for the variational form of implicit Euler through vertex-level Gauss-Seidel iterations. It operates with local vertex position updates that achieve reductions in global…
Grouped data are commonly encountered in applications. The Bernstein polynomial model is proposed as an approximate model in this paper for estimating a univariate density function based on grouped data. The coefficients of the Bernstein…
We consider the classical point vortex model in the mean-field scaling regime, in which the velocity field experienced by a single point vortex is proportional to the average of the velocity fields generated by the remaining point vortices.…
Finding the stationary states of a free energy functional is an important problem in phase field crystal (PFC) models. Many efforts have been devoted for designing numerical schemes with energy dissipation and mass conservation properties.…
The paper studies the convergence of some parallel multisplitting block iterative methods for the solution of linear systems arising in the numerical solution of Euler equations. Some sufficient conditions for convergence are proposed. As…
We consider in this paper random batch interacting particle methods for solving the Poisson-Nernst-Planck (PNP) equations, and thus the Poisson-Boltzmann (PB) equation as the equilibrium, in the external unbounded domain. To justify the…
In \cite{NRxx}, we proposed a numerical regularized moment method of arbitrary order (abbreviated as NRxx method) for Boltzmann-BGK equation, which makes numerical simulation using very large number of moments possible. In this paper, we…
Context: There is increasing need for good algorithms for modeling the aggregation and fragmentation of solid particles (dust grains, dust aggregates, boulders) in various astrophysical settings, including protoplanetary disks, planetary-…
The Wasserstein barycenter extends the Euclidean mean to the space of probability measures by minimizing the weighted sum of squared 2-Wasserstein distances. We develop a free-support algorithm for computing Wasserstein barycenters that…
We present the method for computation of fluid flows that are characterized by the large degree of expansion/contraction and in which the fluid velocity is dominated by the bulk component associated with the expansion/contraction and/or…