Related papers: Efficient numerical solution to vacuum decay with …
A numerical study of the Einstein field equations, based on the 3+1 foliation of the spacetime, is presented. A pseudo-spectral technique has been employed for simulations in vacuum, within two different formalisms, namely the…
Many astrophysical simulations involve extreme dynamic range of timescales around 'special points' in the domain (e.g. black holes, stars, planets, disks, galaxies, shocks, mixing interfaces), where processes on small scales couple strongly…
We propose a numerical method to solve the three-dimensional static Maxwell equations in a singular axisymmetric domain, generated by the rotation of a singular polygon around one of its sides. The mathematical tools and an in-depth study…
We present a numerical method to deal efficiently with large numbers of particles in incompressible fluids. The interactions between particles and fluid are taken into account by a physically motivated ansatz based on locally defined drag…
Explicit simulations of fluid mixtures of highly size-dispersed particles are constrained by numerical challenges associated with identifying pair-interaction neighbors. Recent algorithmic developments have ameliorated these difficulties to…
We develop a numerical method for solving a system of nonlinear integral equations involving two integral terms: at the current time t, one integral is taken from 0 to t, and a different integral is taken from t to infinity. We prove the…
Modelling real world systems frequently requires the solution of systems of nonlinear equations. A number of approaches have been suggested and developed for this computational problem. However, it is also possible to attempt solutions…
We develop two numerical methods to solve the differential equations with deviating arguments for the motion of two charges in the action-at-a-distance electrodynamics. Our first method uses St\"urmer's extrapolation formula and assumes…
This paper presents a numerical function optimization framework designed for constrained optimization problems in robotics. The tool is designed with real-time considerations and is suitable for online trajectory and control input…
We consider the problem of numerically integrating functions with hyperplane discontinuities over the entire Euclidean space in many dimensions. We describe a simple process through which the Euclidean space is partitioned into simplices on…
We propose a nonlinear Discrete Duality Finite Volume scheme to approximate the solutions of drift diffusion equations. The scheme is built to preserve at the discrete level even on severely distorted meshes the energy / energy dissipation…
The scaled boundary finite element method is known for its capability in reproducing highly-detailed solution fields. This, however, is only attainable in those cases where analytical solutions exist. Many others invoke the use of numerical…
This work focuses on the development of a self adjusting multirate strategy based on an implicit time discretization for the numerical solution of hyperbolic equations, that could benefit from different time steps in different areas of the…
Simulations of biophysical systems inevitably include steps that correspond to time integrations of ordinary differential equations. These equations are often related to enzyme action in the synthesis and destruction of molecular species,…
An important class of problems exhibits smooth behaviour on macroscopic space and time scales, while only a microscopic evolution law is known. For such time-dependent multi-scale problems, an "equation-free" framework has been proposed, of…
"Bubble Ball" is a game built on a 2D physics engine, where a finite set of objects can modify the motion of a bubble-like ball. The objective is to choose the set and the initial configuration of the objects, in order to get the ball to…
Modeling of microlensing events poses computational challenges for the resolution of the lens equation and the high dimensionality of the parameter space. In particular, numerical noise represents a severe limitation to fast and efficient…
Finite difference method was extended to unstructured meshes to solve Euler equations. The spatial discretization is made of two steps. First, numerical fluxes are computed at the middle point of each edge with high order accuracy. In this…
We present a novel algorithm to compute multi-scale curvature fields on triangle meshes. Our algorithm is based on finding robust mean curvatures using the ball neighborhood, where the radius of a ball corresponds to the scale of the…
In this paper, we develop a class of robust numerical methods for solving dynamical systems with multiple time scales. We first represent the solution of a multiscale dynamical system as a transformation of a slowly varying solution. Then,…