Related papers: Computing Nearly Singular Solutions Using Pseudo-S…
The high-order accuracy of Fourier method makes it the method of choice in many large scale simulations. We discuss here the stability of Fourier method for nonlinear evolution problems, focusing on the two prototypical cases of the…
In this paper, we perform a careful numerical study of nearly singular solutions of the 3D incompressible Euler equations with smooth initial data. We consider the interaction of two perturbed antiparallel vortex tubes which was previously…
Pseudo-spectral methods are widely used for direct numerical simulations of turbulence, but the standard 2/3 truncation rule for dealiasing is computationally expensive -- accounting for up to 80% of the total cost in three dimensions.…
Efficient and accurate spectral solvers for nonlocal models in any spatial dimension are presented. The approach we pursue is based on the Fourier multipliers of nonlocal Laplace operators introduced in a previous work. It is demonstrated…
Despite their ubiquity throughout science and engineering, only a handful of partial differential equations (PDEs) have analytical, or closed-form solutions. This motivates a vast amount of classical work on numerical simulation of PDEs and…
We present a pseudo-spectral method for solving the three-dimensional Boussinesq equations in unbounded cylindrical domains, specifically tailored for rotating, stably stratified flows subject to strong azimuthal shear. To effectively…
Pseudospectral numerical schemes for solving the Dirac equation in general static curved space are derived using a pseudodifferential representation of the Dirac equation along with a simple Fourier-basis technique. Owing to the presence of…
Three pseudospectral algorithms are described (Euler, leapfrog and trapez) for solving numerically the time dependent nonlinear Schroedinger equation in one, two or three dimensions. Numerical stability regions in the parameter space are…
We discuss the rigorous justification of the spatial discretization by means of Fourier spectral methods of quasilinear first-order hyperbolic systems. We provide uniform stability estimates that grant spectral convergence of the…
Whether the 3D incompressible Euler equations can develop a singularity in finite time from smooth initial data is one of the most challenging problems in mathematical fluid dynamics. This work attempts to provide an affirmative answer to…
In this paper we are interested in the numerical solution of stochastic differential equations with non negative solutions. Our goal is to construct explicit numerical schemes that preserve positivity, even for super linear stochastic…
A quasi-second order scheme is developed to obtain approximate solutions of the shallow water equationswith bathymetry. The scheme is based on a staggered finite volume scheme for the space discretization:the scalar unknowns are located in…
We present a new numerical scheme which combines the Spectral Difference (SD) method up to arbitrary high order with \emph{a-posteriori} limiting using the classical MUSCL-Hancock scheme as fallback scheme. It delivers very accurate…
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…
We present an open-source Python implementation of an idealized high-order pseudo-spectral solver for the one-dimensional nonlinear Schr\"odinger equation (NLSE). The solver combines Fourier spectral spatial discretization with an adaptive…
A high-performance shooting algorithm is developed to compute time-periodic solutions of the free-surface Euler equations with spectral accuracy in double and quadruple precision. The method is used to study resonance and its effect on…
In this paper, we propose a linearized Fourier pseudo-spectral method, which preserves the total mass and energy conservation laws, for the damped nonlinear Schr\"{o}dinger equation in three dimensions. With the aid of the semi-norm…
In this paper, we propose Fourier pseudospectral methods to solve the variable-order space fractional wave equation and develop an accelerated matrix-free approach for its effective implementation. In constant-order cases, our methods can…
The diffuse-domain, or smoothed boundary, method is an attractive approach for solving partial differential equations in complex geometries because of its simplicity and flexibility. In this method the complex geometry is embedded into a…
Many physical problems such as Allen-Cahn flows have natural maximum principles which yield strong point-wise control of the physical solutions in terms of the boundary data, the initial conditions and the operator coefficients.…