Related papers: Numerical evolutions of fields on the 2-sphere usi…
In this paper, we present a stochastic method for the simulation of Laplace's equation with a mixed boundary condition in planar domains that are polygonal or bounded by circular arcs. We call this method the Reflected Walk-on-Spheres…
Within the well-established optical response function formalism, a new strategy with the central idea of employing the forward-backward stochastic Schr\"{o}dinger equations in a segmented way to accurately obtain the two-dimensional (2D)…
3D image processing constitutes nowadays a challenging topic in many scientific fields such as medicine, computational physics and informatics. Therefore, development of suitable tools that guaranty a best treatment is a necessity.…
In a previous article [J. Comp. Phys. $\mathbf{357}$ (2018) 282-304], the mixed mimetic spectral element method was used to solve the rotating shallow water equations in an idealized geometry. Here the method is extended to a smoothly…
We present a robust algorithm for the computation of electromagnetic fields radiated by point sources (Hertzian dipoles) in cylindrically stratified media where each layer may exhibit material properties (permittivity, permeability, and…
We present a new formulation of Einstein's equations for an axisymmetric spacetime with vanishing twist in vacuum. We propose a fully constrained scheme and use spherical polar coordinates. A general problem for this choice is the…
In this paper we introduce a reproducing kernel Hilbert space defined on $\mathbb{R}^{d+1}$ as the tensor product of a reproducing kernel defined on the unit sphere $\mathbb{S}^{d}$ in $\mathbb{R}^{d+1}$ and a reproducing kernel defined on…
Spectral measures arise in numerous applications such as quantum mechanics, signal processing, resonances, and fluid stability. Similarly, spectral decompositions (pure point, absolutely continuous and singular continuous) often…
This paper introduces filtered finite difference methods for numerically solving a dispersive evolution equation with solutions that are highly oscillatory in both space and time. We consider a semiclassically scaled nonlinear Schr\"odinger…
We have presented some practical consequences on the molecular-dynamics simulations arising from the numerical algorithm published recently in paper Int. J. Mod. Phys. C 16, 413 (2005). The algorithm is not a finite-difference method and…
We present a numerical study of the Einstein equations, according to the Arnowitt-Deser-Misner (ADM) formalism, in order to simulate the dynamics of gravitational fields. We took in consideration the original $3+1$ decomposition of the ADM…
Odd-type spin 2 perturbations of Einstein's equation can be reduced to the scalar Regge-Wheeler equation. We show that the weighted norms of solutions are in L^2 of time and space. This result uses commutator methods and applies uniformly…
We develop an efficient numerical scheme for the 3D mean-field spherical dynamo equation. The scheme is based on a semi-implicit discretization in time and a spectral method in space based on the divergence-free spherical harmonic…
In the present manuscript, we consider the problem of dispersive wave simulation on a rotating globally spherical geometry. In this Part IV, we focus on numerical aspects while the model derivation was described in Part III. The algorithm…
Understanding the complex dynamics and structure of the upper solar atmosphere benefits strongly from the use of a combination of several diagnostics. Frequently, such diverse diagnostics can only be obtained from telescopes and/or…
An unsteady problem is considered for a space-fractional equation in a bounded domain. A first-order evolutionary equation involves the square root of an elliptic operator of second order. Finite element approximation in space is employed.…
Stochastic differential equations are widely used in various fields; in particular, the usefulness of duality relations has been demonstrated in some models such as population models and Brownian momentum processes. In this study, a…
We develop numerical homotopy algorithms for solving systems of polynomial equations arising from the classical Schubert calculus. These homotopies are optimal in that generically no paths diverge. For problems defined by hypersurface…
A new method for solving inverse spectral problems on quantum star graphs is proposed. The method is based on Neumann series of Bessel functions representations for solutions of Sturm-Liouville equations. The representations admit estimates…
Exact diagonalization and other numerical studies of quantum spin systems are notoriously limited by the exponential growth of the Hilbert space dimension with system size. A common and well-known practice to reduce this increasing…