Related papers: pizza: An open-source pseudo-spectral code for sph…
We present investigations of rapidly-rotating convection in a thick spherical shell geometry relevant to planetary cores, comparing results from Quasi-Geostrophic, 3D and hybrid QG-3D models. The 170 reported calculations span Ekman…
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…
A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere.…
We examine the dynamics associated with weakly compressible convection in a spherical shell by running 3D direct numerical simulations using the Boussinesq formalism [1]. Motivated by problems in astrophysics, we assume the existence of a…
We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are…
We present a new high-order accurate computational fluid dynamics model based on the incompressible Navier-Stokes equations with a free surface for the accurate simulation of nonlinear and dispersive water waves in the time domain. The…
Buoyancy-induced (Rayleigh-Benard) convection of a fluid between two horizontal plates is a central paradigm for studying the transition to complex spatiotemporal dynamics in sustained nonequilibrium systems. To improve the analysis of…
We present a new hybrid discrete exterior calculus (DEC) and finite difference (FD) method to simulate fully three-dimensional Boussinesq convection in spherical shells subject to internal heating and basal heating, relevant in the…
A three-dimensional finite-difference solver has been developed and implemented for Boussinesq convection in a spherical shell. The solver transforms any complex curvilinear domain into an equivalent Cartesian domain using Jacobi…
Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar…
Rapidly rotating spherical kinematic dynamos are computed using the combination of a quasi geostrophic (QG) model for the velocity field and a classical spectral 3D code for the magnetic field. On one hand, the QG flow is computed in the…
Models of inhomogeneous universes constructed with exact solutions of Einstein's General Relativity have been proposed in the literature with the aim of reproducing the cosmological data without any need for a dark energy component. Besides…
We study rapidly-rotating Boussinesq convection driven by internal heating in a full sphere. We use a numerical model based on the quasi-geostrophic approximation for the velocity field, whereas the temperature field is three-dimensional.…
Cohesive fracture is among the few techniques able to model complex fracture nucleation and propagation with a sharp (nonsmeared) representation of the crack. Implicit time-stepping schemes are often favored in mechanics due to their…
We generalize ideas in the recent literature and develop new ones in order to propose a general class of contour integral methods for linear convection-diffusion PDEs and in particular for those arising in finance. These methods aim to…
We present a fully pseudo-spectral scheme to solve axisymmetric hyperbolic equations of second order. With the Chebyshev polynomials as basis functions, the numerical grid is based on the Lobbato (for two spatial directions) and Radau (for…
We present a major update to QuSpin, SciPostPhys.2.1.003 -- an open-source Python package for exact diagonalization and quantum dynamics of arbitrary boson, fermion and spin many-body systems, supporting the use of various (user-defined)…
We present a novel numerical solver for the systems of coupled non-linear elliptical differential equations. The solver partitions the computational domain into a set of rectangular pseudo-spectral collocation subdomains and is especially…
Numerical approximation of a general class of nonlinear unidirectional wave equations with a convolution-type nonlocality in space is considered. A semi-discrete numerical method based on both a uniform space discretization and the discrete…
We apply the ultraspherical spectral method to solving time-dependent PDEs by proposing two approaches to discretization based on the method of lines and show that these approaches produce approximately same results. We analyze the…