Numerical solution of Navier–Stokes equations using multiquadric radial basis function networks

A numerical method based on radial basis function networks (RBFNs) for solving steady incompressible viscous flow problems (including Boussinesq materials) is presented in this paper. The method uses a 'universal approximator' based on neural network methodology to represent the solutions. The method is easy to implement and does not require any kind of 'finite element-type' discretization of the domain and its boundary. Instead, two sets of random points distributed throughout the domain and on the boundary are required. The first set defines the centres of the RBFNs and the second defines the collocation points. The two sets of points can be different; however, experience shows that if the two sets are the same better results are obtained. In this work the two sets are identical and hence commonly referred to as the set of centres. Planar Poiseuille, driven cavity and natural convection flows are simulated to verify the method. The numerical solutions obtained using only relatively low densities of centres are in good agreement with analytical and benchmark solutions available in the literature. With uniformly distributed centres, the method achieves Reynolds number Re = 100 000 for the Poiseuille flow (assuming that laminar flow can be maintained) using the density of 11x11, Re = 400 for the driven cavity flow with a density of 33x33 and Rayleigh number Ra = 1 000 000 for the natural convection flow with a density of 27x27.

[1]  E. Kansa.  MULTIQUADRICS--A SCATTERED DATA APPROXIMATION SCHEME WITH APPLICATIONS TO COMPUTATIONAL FLUID-DYNAMICS-- II SOLUTIONS TO PARABOLIC, HYPERBOLIC AND ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS , 1990 .

[2]  Leopold Škerget,et al.  Mixed boundary elements for laminar flows , 1999 .

[3]  D. Broomhead,et al.  Radial Basis Functions, Multi-Variable Functional Interpolation and Adaptive Networks , 1988 .

[4]  J. Reddy,et al.  The Finite Element Method in Heat Transfer and Fluid Dynamics , 1994 .

[5]  Simon Haykin,et al.  Neural Networks: A Comprehensive Foundation , 1998 .

[6]  M. J. D. Powell,et al.  Radial basis functions for multivariable interpolation: a review , 1987 .

[7]  C. Fletcher.  Computational Galerkin Methods , 1983 .

[8]  G. D. Davis.  Natural convection of air in a square cavity: A bench mark numerical solution , 1983 .

[9]  J. Z. Zhu,et al.  The finite element method , 1977 .

[10]  E. J. Kansa,et al.  Application of the Multiquadric Method for Numerical Solution of Elliptic Partial Differential Equations , 2022 .

[11]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[12]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[13]  F. Girosi,et al.  Networks for approximation and learning , 1990, Proc. IEEE.

[14]  E. Kansa.  Multiquadrics—A scattered data approximation scheme with applications to computational fluid-dynamics—I surface approximations and partial derivative estimates , 1990 .

[15]  B. P. Leonard,et al.  Why You Should Not Use 'Hybrid', 'Power-Law' or Related Exponential Schemes for Convective Modelling—There Are Much Better Alternatives , 1995 .

[16]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[17]  John Moody,et al.  Fast Learning in Networks of Locally-Tuned Processing Units , 1989, Neural Computation.

[18]  R. Cook,et al.  Concepts and Applications of Finite Element Analysis , 1974 .

[19]  G. J. Moridis,et al.  The Laplace Transform Multiquadric Method: A Highly Accurate Scheme for the Numerical Solution of Partial Differential Equations , 1993 .

[20]  Nam Mai-Duy,et al.  Numerical solution of differential equations using multiquadric radial basis function networks , 2001, Neural Networks.

[21]  Patrick J. Roache,et al.  Verification and Validation in Computational Science and Engineering , 1998 .