A High-Order Radial Basis Function (RBF) Leray Projection Method for the Solution of the Incompressible Unsteady Stokes Equations

A new projection method based on radial basis functions (RBFs) is presented for discretizing the incompressible unsteady Stokes equations in irregular geometries. The novelty of the method comes from the application of a new technique for computing the Leray-Helmholtz projection of a vector field using generalized interpolation with divergence-free and curl-free RBFs. Unlike traditional projection methods, this new method enables matching both tangential and normal components of divergence-free vector fields on the domain boundary. This allows incompressibility of the velocity field to be enforced without any time-splitting or pressure boundary conditions. Spatial derivatives are approximated using collocation with global RBFs so that the method only requires samples of the field at (possibly scattered) nodes over the domain. Numerical results are presented demonstrating high-order convergence in both space (between 5th and 6th order) and time (up to 4th order) for some model problems in two dimensional irregular geometries.Comment: 34 pages, 8 figure

A High-Order Kernel Method for Diffusion and Reaction-Diffusion Equations on Surfaces

In this paper we present a high-order kernel method for numerically solving diffusion and reaction-diffusion partial differential equations (PDEs) on smooth, closed surfaces embedded in $\mathbb{R}^d$. For two-dimensional surfaces embedded in $\mathbb{R}^3$, these types of problems have received growing interest in biology, chemistry, and computer graphics to model such things as diffusion of chemicals on biological cells or membranes, pattern formations in biology, nonlinear chemical oscillators in excitable media, and texture mappings. Our kernel method is based on radial basis functions (RBFs) and uses a semi-discrete approach (or the method-of-lines) in which the surface derivative operators that appear in the PDEs are approximated using collocation. The method only requires nodes at "scattered" locations on the surface and the corresponding normal vectors to the surface. Additionally, it does not rely on any surface-based metrics and avoids any intrinsic coordinate systems, and thus does not suffer from any coordinate distortions or singularities. We provide error estimates for the kernel-based approximate surface derivative operators and numerically study the accuracy and stability of the method. Applications to different non-linear systems of PDEs that arise in biology and chemistry are also presented

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

Adaptive meshless refinement schemes for RBF-PUM collocation

In this paper we present an adaptive discretization technique for solving elliptic partial differential equations via a collocation radial basis function partition of unity method. In particular, we propose a new adaptive scheme based on the construction of an error indicator and a refinement algorithm, which used together turn out to be ad-hoc strategies within this framework. The performance of the adaptive meshless refinement scheme is assessed by numerical tests

Numerical solutions of a boundary value problem on the sphere using radial basis functions

Boundary value problems on the unit sphere arise naturally in geophysics and oceanography when scientists model a physical quantity on large scales. Robust numerical methods play an important role in solving these problems. In this article, we construct numerical solutions to a boundary value problem defined on a spherical sub-domain (with a sufficiently smooth boundary) using radial basis functions (RBF). The error analysis between the exact solution and the approximation is provided. Numerical experiments are presented to confirm theoretical estimates

Kernel-based stochastic collocation for the random two-phase Navier-Stokes equations

In this work, we apply stochastic collocation methods with radial kernel basis functions for an uncertainty quantification of the random incompressible two-phase Navier-Stokes equations. Our approach is non-intrusive and we use the existing fluid dynamics solver NaSt3DGPF to solve the incompressible two-phase Navier-Stokes equation for each given realization. We are able to empirically show that the resulting kernel-based stochastic collocation is highly competitive in this setting and even outperforms some other standard methods

Combining the radial basis function Eulerian and Lagrangian schemes with geostatistics for modeling of radionuclide migration through the geosphere

To assess the long-term safety of a radioactive waste disposal system, mathematical models are used to describe groundwater flow, chemistry, and potential radionuclide migration through geological formations. A number of processes need to be considered, when predicting the movement of radionuclides through the geosphere. The most important input data are obtained from field measurements, which are not available for all regions of interest. For example, the hydraulic conductivity as an input parameter varies from place to place. In such cases, geostatistical science offers a variety of spatial estimation procedures. Methods for solving the solute transport equation can also be classified as Eulerian, Lagrangian and mixed. The numerical solution of partial differential equations (PDE) has usually been obtained by finite-difference methods (FDM), finite-element methods (FEM), or finite-volume methods (FVM). Kansa introduced the concept of solving partial differential equations using radial basis functions (RBF) for hyperbolic, parabolic, and elliptic PDEs. The aim of this study was to present a relatively new approach to the modeling of radionuclide migration through the geosphere using radial basis function methods in Eulerian and Lagrangian coordinates. In this study, we determine the average and standard deviation of radionuclide concentration with regard to variable hydraulic conductivity, which was modelled by a geostatistical approach. Radionuclide concentrations will also be calculated in heterogeneous and partly heterogeneous 2D porous media. (C) 2004 Elsevier Ltd. All rights reserved

RBF multiscale collocation for second order elliptic boundary value problems

In this paper, we discuss multiscale radial basis function collocation methods for solving elliptic partial differential equations on bounded domains. The approximate solution is constructed in a multi-level fashion, each level using compactly supported radial basis functions of smaller scale on an increasingly fine mesh. On each level, standard symmetric collocation is employed. A convergence theory is given, which builds on recent theoretical advances for multiscale approximation using compactly supported radial basis functions. We are able to show that the convergence is linear in the number of levels. We also discuss the condition numbers of the arising systems and the effect of simple, diagonal preconditioners, now proving rigorously previous numerical observations