Related papers: pizza: An open-source pseudo-spectral code for sph…
We investigate the computation of stable fracture paths in brittle thin films using one-dimensional damage models with an elastic foundation. The underlying variational formulation is non-convex, making the evolution path sensitive to…
This paper presents a novel spatial discretisation method for the reliable and efficient simulation of Bose-Einstein condensates modelled by the Gross-Pitaevskii equation and the corresponding nonlinear eigenvector problem. The method…
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an…
We have developed a cosmic ray (CR) shock code in one dimensional spherical geometry with which the particle distribution, the gas flow and their nonlinear interaction can be followed numerically in a frame comoving with an expanding shock.…
Large regions of giant planets are thought to possess unstable thermal gradients stabilised by gradients in heavy-element composition. The fluid can then develop semi-convection, a double-diffusive instability driven by the unequal…
The proliferation of quantum fluctuations and long-range entanglement presents an outstanding challenge for the numerical simulation of interacting spin systems with exotic ground states. Here, we present a toolset of Chebyshev…
We present a general formalism for performing a time-dependent Bogoliubov analysis of a dynamically unstable Bose-Einstein condensate, which extends the quasiparticle projection method of Morgan et al. [Phys. Rev. A 57, 3818 (1998)] to…
Let $S$ be a compact, orientable surface of hyperbolic type. Let $(k_+,k_-)$ be a pair of negative numbers and let $(g_+, g_-)$ be a pair of marked metrics over $S$ of constant curvature equal to $k_+$ and $k_-$ respectively. Using a…
In this paper, we consider stochastic second-order methods for minimizing a finite summation of nonconvex functions. One important key is to find an ingenious but cheap scheme to incorporate local curvature information. Since the true…
We present a new implementation of the geometric method of Cullen & Purser (1984) for solving the semi-geostrophic Eady slice equations which model large scale atmospheric flows and frontogenesis. The geometric method is a Lagrangian…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
In most current debris disc models, the dynamical and the collisional evolutions are studied separately, with N-body and statistical codes, respectively, because of stringent computational constraints. We present here LIDT-DD, the first…
In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization (SP) algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent…
We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three…
We study equal weight numerical integration, or Quasi Monte Carlo (QMC) rules, for functions in a Sobolev space $H^s(S^d)$ with smoothness parameter $s>d/2$ defined over the unit sphere $S^d$ in $R^{d+1}$. Focusing on $N$-point sets that…
Spherical confinement in 3D harmonic, quartic and other higher oscillators of even order is studied. The generalized pseudospectral method is employed for accurate solution of relevant Schr\"odinger equation in an \emph{optimum,…
A novel second order family of explicit stabilized Runge-Kutta-Chebyshev methods for advection-diffusion-reaction equations is introduced. The new methods outperform existing schemes for relatively high Peclet number due to their favorable…
We study the ubiquitous super-resolution problem, in which one aims at localizing positive point sources in an image, blurred by the point spread function of the imaging device. To recover the point sources, we propose to solve a convex…
Typical multispecies compressible Navier-Stokes computations employ conservative equations for mass fraction transport. Upwind discretisations of these governing equations produce spurious pressure oscillations at diffuse contact surfaces…
Spectral methods are now common in the solution of ordinary differential eigenvalue problems in a wide variety of fields, such as in the computation of black hole quasinormal modes. Most of these spectral codes are based on standard…