We present a Fourier Continuation-based parallel pseudospectral method for incompressible fluids in cuboid non-periodic domains. The method produces dispersionless and dissipationless derivatives with fast spectral convergence inside the domain, and with very high order convergence at the boundaries. Incompressibility is imposed by solving a Poisson equation for the pressure. Being Fourier-based, the method allows for fast computation of spectral transforms. It is compatible with uniform grids (although refined or nested meshes can also be implemented), which in turn allows for explicit time integration at sufficiently high Reynolds numbers. Using a new parallel code named SPECTER we illustrate the method with two problems: channel flow, and plane Rayleigh-Bénard convection under the Boussinesq approximation. In both cases the method yields results compatible with previous studies using other high-order numerical methods, with mild requirements on the time step for stability.

No abstract

A parametric study of the magnetic dipole behavior in resistive incompressible MHD inside a rotating sphere is performed, using direct numerical simulations and considering Reynolds and Ekman numbers as controlling parameters. The tendency is to obtain geodynamo-like magnetic dipole reversal regimes for sufficiently small Ekman and large Reynolds numbers. The typical dipole latitude obtained in the reversal regime is around 40 degrees (with respect to the rotation axis of the sphere). A statistical analysis of waiting times between dipole reversals is also performed, obtaining a non-Poissonian distribution of waiting times, indicating long-term memory effects. We also report the presence of a 1/f frequency power spectrum in the magnetic dipole time-series, which also shows a tendency to grow toward lower frequencies as the Ekman number is decreased. arXiv:1808.01370v2 [physics.flu-dyn]