A RBF partition of unity collocation method based on finite difference for initial-boundary value problems

Meshfree radial basis function (RBF) methods are popular tools used to numerically solve partial differential equations (PDEs). They take advantage of being flexible with respect to geometry, easy to implement in higher dimensions, and can also provide high order convergence. Since one of the main disadvantages of global RBF-based methods is generally the computational cost associated with the solution of large linear systems, in this paper we focus on a localizing RBF partition of unity method (RBF-PUM) based on a finite difference (FD) scheme. Specifically, we propose a new RBF-PUM-FD collocation method, which can successfully be applied to solve time-dependent PDEs. This approach allows to significantly decrease ill-conditioning of traditional RBF-based methods. Moreover, the RBF-PUM-FD scheme results in a sparse matrix system, reducing the computational effort but maintaining at the same time a high level of accuracy. Numerical experiments show performances of our collocation scheme on two benchmark problems, involving unsteady convection-diffusion and pseudo-parabolic equations

A Hybrid Radial Basis Function - Pseudospectral Method for Thermal Convection in a 3-D Spherical Shell

A novel hybrid spectral method that combines radial basis function (RBF) and Chebyshev pseudospectral (PS) methods in a “2+1” approach is presented for numerically simulating thermal convection in a 3-D spherical shell. This is the first study to apply RBFs to a full 3D physical model in spherical geometry. In addition to being spectrally accurate, RBFs are not defined in terms of any surface based coordinate system such as spherical coordinates. As a result, when used in the lateral directions, as in this study, they completely circumvent the pole issue with the further advantage that nodes can be “scattered” over the surface of a sphere. In the radial direction, Chebyshev polynomials are used, which are also spectrally accurate and provide the necessary clustering near the boundaries to resolve boundary layers. Applications of this new hybrid methodology are given to the problem of convection in the Earth’s mantle,which is modeled by a Boussinesq fluid at infinite Prandtl number. To see whether this numerical technique warrants further investigation, the study limits itself to an isoviscous mantle.Benchmark comparisons are presented with other currently used mantle convection codes for Rayleigh number 7 · 103 and 105. The algorithmic simplicity of the code (mostly due to RBFs)allows it to be written in less than 400 lines of Matlab and run on a single workstation. We find that our method is very competitive with those currently used in the literature

Wave model for longitudinal dispersion: Application to the laminar-flow tubular reactor

The wave model for longitudinal dispersion, published elsewhere as an alternative to the commonly used dispersed plug-flow model, is applied to the classic case of the laminar-flow tubular reactor. The results are compared in a wide range of situations to predictions by the dispersed plug-flow model as well as to exact numerical calculations with the 2-D model of the reactor and to other available methods. In many practical cases, the solutions of the wave model agree closely with the exact data. The wave model has a much wider region of validity than the dispersed plug-flow model, has a distinct physical background, and is easier to use for reactor calculations. This provides additional support to the theory developed elsewhere. The properties and the applicability of the wave model to situations with rapidly changing concentration fields are discussed. Constraints to be satisfied are established to use the new theory with confidence for arbitrary initial and boundary conditions

DNS of bifurcations in an air-filled rotating baroclinic annulus

Three-dimensional Direct Numerical Simulation (DNS) on the nonlinear dynamics and a route to chaos in a rotating fluid subjected to lateral heating is presented here and discussed in the context of laboratory experiments in the baroclinic annulus. Following two previous preliminary studies by Maubert and Randriamampianina, the fluid used is air rather than a liquid as used in all other previous work. This study investigated a bifurcation sequence from the axisymmetric flow to a number of complex flows. The transition sequence, on increase of the rotation rate, from the axisymmetric solution via a steady, fully-developed baroclinic wave to chaotic flow followed a variant of the classical quasi-periodic bifurcation route, starting with a subcritical Hopf and associated saddle-node bifurcation. This was followed by a sequence of two supercritical Hopf-type bifurcations, first to an amplitude vacillation, then to a three-frequency quasi-periodic modulated amplitude vacillation (MAV), and finally to a chaotic MAV\@. In the context of the baroclinic annulus this sequence is unusual as the vacillation is usually found on decrease of the rotation rate from the steady wave flow. Further transitions of a steady wave with a higher wave number pointed to the possibility that a barotropic instability of the side wall boundary layers and the subsequent breakdown of these barotropic vortices may play a role in the transition to structural vacillation and, ultimately, geostrophic turbulence.Comment: 31 page

Linear and Nonlinear Convection in Porous Media between Coaxial Cylinders

In this thesis we develop a mathematical model for describing three-dimensional natural convection in porous media filling a vertical annular cylinder. We apply a linear stability analysis to determine the onset of convection and the preferred convective mode when the annular cylinder is subject to two different types of boundary conditions: heat insulated sidewalls and heat conducting sidewalls. The case of an annular cylinder with insulated sidewalls has been investigated earlier, but our results reveal more details than previously found. We also investigate the case of the radius of the inner cylinder approaching zero and the results are compared with previous work for non-annular cylinders. Using pseudospectral methods we have built a high-order numerical simulator to uncover the nonlinear regime of the convection cells. We study onset and geometry of convection modes, and look at the stability of the modes with respect to different types of perturbations. Also, we examine how variations in the Rayleigh number affects the convection modes and their stability regimes. We uncover an increased complexity regarding which modes that are obtained and we are able to identify stable secondary and mixed modes. We find the different convective modes to have overlapping stability regions depending on the Rayleigh number. The motivation for studying natural convection in porous media is related to geothermal energy extraction and we attempt to determine the effect of convection cells in a geothermal heat reservoir. However, limitations in the simulator do not allow us to make any conclusions on this matter.Master i Anvendt og beregningsorientert matematikkMAMN-MABMAB39

Poloidal-toroidal decomposition in a finite cylinder. II. Discretization, regularization and validation

The Navier-Stokes equations in a finite cylinder are written in terms of poloidal and toroidal potentials in order to impose incompressibility. Regularity of the solutions is ensured in several ways: First, the potentials are represented using a spectral basis which is analytic at the cylindrical axis. Second, the non-physical discontinuous boundary conditions at the cylindrical corners are smoothed using a polynomial approximation to a steep exponential profile. Third, the nonlinear term is evaluated in such a way as to eliminate singularities. The resulting pseudo-spectral code is tested using exact polynomial solutions and the spectral convergence of the coefficients is demonstrated. Our solutions are shown to agree with exact polynomial solutions and with previous axisymmetric calculations of vortex breakdown and of nonaxisymmetric calculations of onset of helical spirals. Parallelization by azimuthal wavenumber is shown to be highly effective