Related papers: An efficient semi-Lagrangian algorithm for simulat…

Safety is essential for reinforcement learning (RL) applied in the real world. Adding chance constraints (or probabilistic constraints) is a suitable way to enhance RL safety under uncertainty. Existing chance-constrained RL methods like…

Numerical solution of partial differential equations on parallel computers using domain decomposition usually requires synchronization and communication among the processors. These operations often have a significant overhead in terms of…

Particle swam optimization (PSO) is a popular stochastic optimization method that has found wide applications in diverse fields. However, PSO suffers from high computational complexity and slow convergence speed. High computational…

Model Order Reduction (MOR) based on Proper Orthogonal Decomposition (POD) and Smooth Particle Hydrodynamics (SPH) has proven effective in various applications. Most MOR methods utilizing POD are implemented within a pure Eulerian…

Spatial reaction-diffusion models have been employed to describe many emergent phenomena in biological systems. The modelling technique most commonly adopted in the literature implements systems of partial differential equations (PDEs),…

We introduce improved Reduced Order Models (ROM) for convection-dominated flows. These non-linear closure models are inspired from successful numerical stabilization techniques used in Large Eddy Simulations (LES), such as Local Projection…

Suspensions with fiber-like particles in the low Reynolds number regime are modeled by two different approaches that both use a Lagrangian representation of individual particles. The first method is the well-established formulation based on…

The Lorentz equations describe the motion of electrically charged particles in electric and magnetic fields and are used widely in plasma physics. The most popular numerical algorithm for solving them is the Boris method, a variant of the…

In this paper, we propose an efficient high order semi-Lagrangian (SL) discontinuous Galerkin (DG) method for solving linear convection-diffusion equations. The method generalizes our previous work on developing the SLDG method for…

We propose a kinetic Ising model to study phase separation driven by surface diffusion. This model is referred to as "Model S", and consists of the usual Kawasaki spin-exchange kinetics ("Model B") in conjunction with a kinetic constraint.…

We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the…

Several stochastic simulation algorithms (SSAs) have been recently proposed for modelling reaction-diffusion processes in cellular and molecular biology. In this paper, two commonly used SSAs are studied. The first SSA is an on-lattice…

This study proposes an extension of the high-order compact gas-kinetic scheme (CGKS) to compressible flow simulation in an arbitrary Lagrangian-Eulerian (ALE) formulation in unstructured mesh. The ALE method is achieved by subdividing…

We propose a class of semi-Lagrangian methods of high approximation order in space and time, based on spectral element space discretizations and exponential integrators of Runge-Kutta type. We discuss the extension of these methods to the…

The inclusion of convection in stellar evolution models lacks realism, especially near convective-radiative interfaces. Furthermore, the interaction of convection with oscillations prevent us from accurately predicting seismic frequencies,…

This paper presents a stabilization scheme which addresses the rank-deficiency problem in meshless collocation methods for solid mechanics. Specifically, Smooth-Particle Hydrodynamics (SPH) in the Total Lagrangian formalism is considered.…

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the Discrete Ion Stochastic Continuum Overdamped Solvent (DISCOS) algorithm was derived for triply…

Existing approaches to solving the Vlasov equation treat the system as a partial differential equation on a phase space grid, and track in either an Eulerian, Lagrangian, or semi-Lagrangian picture. We present an alternative approach, which…

In this paper we propose a new inexact dual decomposition algorithm for solving separable convex optimization problems. This algorithm is a combination of three techniques: dual Lagrangian decomposition, smoothing and excessive gap. The…

We describe an exact and highly efficient numerical algorithm for solving a special but important class of convection-diffusion equations. These equations occur in many problems in physics, chemistry, or biology, and they are usually hard…