Next Article in Journal
Dynamical Analysis of Hyper-ILSR Rumor Propagation Model with Saturation Incidence Rate
Next Article in Special Issue
Correspondence between the Energy Equipartition Theorem in Classical Mechanics and Its Phase-Space Formulation in Quantum Mechanics
Previous Article in Journal
Process and Time
Previous Article in Special Issue
Uniform Error Estimates of the Finite Element Method for the Navier–Stokes Equations in R2 with L2 Initial Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radial Basis Function Finite Difference Method Based on Oseen Iteration for Solving Two-Dimensional Navier–Stokes Equations

College of Mathematics and System Sciences, Xinjiang University, Urumqi 830017, China
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(5), 804; https://doi.org/10.3390/e25050804
Submission received: 20 March 2023 / Revised: 16 April 2023 / Accepted: 18 April 2023 / Published: 16 May 2023
(This article belongs to the Collection Foundations of Statistical Mechanics)

Abstract

:
In this paper, the radial basis function finite difference method is used to solve two-dimensional steady incompressible Navier–Stokes equations. First, the radial basis function finite difference method with polynomial is used to discretize the spatial operator. Then, the Oseen iterative scheme is used to deal with the nonlinear term, constructing the discrete scheme for Navier–Stokes equation based on the finite difference method of the radial basis function. This method does not require complete matrix reorganization in each nonlinear iteration, which simplifies the calculation process and obtains high-precision numerical solutions. Finally, several numerical examples are obtained to verify the convergence and effectiveness of the radial basis function finite difference method based on Oseen Iteration.

1. Introduction

The Navier–Stokes equations are a set of equations used to describe fluid substance such as liquid and air, and can be used to simulate weather, ocean currents, water flows in pipes, the motion of stars in galaxies, etc. Therefore, they are of significant research value; this paper discusses the numerical solution method of steady incompressible Navier–Stokes equations.
Scholars have studied many numerical methods for Navier–Stokes equations [1]. For unsteady Navier–Stokes problems [2] and steady Navier–Stokes problems [3], a variety of traditional numerical methods have been proposed, including the finite element method, finite difference method, and finite volume method. In response to incompressible conditions, nonlinearity, long-time integration, and other difficulties encountered in numerical solutions of three-dimensional unsteady Navier–Stokes equations, this paper discusses the research status and latest research results of highly efficient and fully discrete finite element methods, which can be used to overcome these difficulties. In addition, this paper illustrates stability and error estimation of finite element-space discrete solutions and optimal error estimation of efficient fully-discrete finite element solutions for solving three-dimensional unsteady Navier–Stokes equations [4]. These methods are mesh-based, and the generation of meshes increases computational cost; moreover, the computational domain and the technology used for mesh quality and stabilization affect the computational accuracy. Especially for high-dimensional problems, it is difficult and costly to generate good meshes. Therefore, in recent years many scholars have introduced meshless methods [5,6] for solving partial differential equations, among which the radial basis function method [7] is increasingly popular.
The radial basis function method has been widely used for solving various partial differential equations, with the radial basis function finite difference method (RBF-FD) proposed considering ill-conditioned problems with a high occurrence rate in the matrix generated by the global radial basis function method. In this way, it is easy to realize RBF-FD discretization without needing to use a mesh at all. The RBF-FD method can be regarded as an extension of the standard finite difference and local refinement of global radial basis function method. This method uses the radial basis function to create weights for the RBF-FD formula, and the differential matrix obtained in this way is as sparse as that obtained by the standard difference method, which combines many advantages of radial basis function and traditional finite difference approximation. Thus, RBF-FD approximation [8] is an attractive substitute for the global radial basis function method.
In order to control the stagnation error, Natasha Flyer et al. introduced higher-order polynomials under the background of Gaussian (GA) and Polyharmonic Spline (PHS) interpolation [9,10]. PHS eliminates the stagnation error in combination with the polynomial to capture the basic physical properties of the problem, and achieves high-order accuracy without the need to adjust the shape parameters. The present paper chooses the PHS basis function and applies the radial basis function finite difference method with the polynomial to discretize the spatial operator of the equation. It is well known that the iterative method [11] is effective for Navier–Stokes equations with strong nonlinear terms, e.g., the Newton iterative method, Stokes iterative method, and Oseen iterative method. This paper focuses on the Oseen iterative method to deal with the nonlinear terms of the equation. The radial basis function finite difference method based on Oseen iteration proposed in this paper does not need a complete set of matrices in each nonlinear iteration, which is more suitable for nonlinear problems, and obtains high-precision numerical solutions.

2. Problem Setup

This paper studies the following two-dimensional steady incompressible Navier–Stokes equations [12].
ν Δ u + ( u · ) u + p = f on Ω , · u = 0 on Ω , u = 0 on Ω .
where  Ω  is the bounded region on  R 2 u = ( u , v ) : Ω R 2  is the velocity vector,  p : Ω R  is the pressure,  ν > 0  is the viscosity coefficient, and  f = ( f 1 , f 2 )  is the external force acting on a unit volume of fluid, assuming that the uniqueness of p can be achieved by imposing the condition,  Ω p d x = 0 .
We define the 2-norm and -norm of the vector  x  as
x 2 = ( i = 1 n x i 2 ) 1 2 , x = m a x 1 i n | x i | .

3. Numerical Method

In this part, we focus on the radial basis function finite difference method with polynomial.

3.1. Radial Basis Function Finite Difference Method with Polynomial

First, this part provides a brief introduction to the radial basis function (RBF). The RBF, expressed as  Φ ( r ) : R d R , is a d-dimensional radially symmetric function, and is only dependent on  r = x x k 2 , where  · 2  represents the Euclidean distance between two points,  x  is the point to be solved, and  x k  is the center position of the RBF. Common radial basis functions include Gauss (GA) ( e ( ε r ) 2 ), multiquadratic function ( 1 + ( ε r ) 2 ), and Polyharmonic Spline (PHS) ( r m , m = 1 , 3 , 5 , ), where  ε  is the shape parameter that determines the radial basis function. Next, we introduce the radial basis function finite difference method with polynomial. First, L is defined as a linear operator, which can be  Δ x y , etc. We can approximate the value of the operator L at node  x c  by linear combination of the function value  { u k } k = 1 n . For this, we need to select n nodes near the central node  x c  to form a node template, which is recorded as  X = { x k } k = 1 n , including  x c , as follows:
L u x = x c = k = 1 n w k u k ,
where n is the size of the node template and  w k  is the differential weight. In order to calculate the weight  w k , we can specify that the linear combination of function values be accurate for interpolation  s ( x ) :
s x = k = 1 n λ k Φ x x k 2 + k = 1 μ k p k x ,
The constraint condition is
k = 1 n λ k p j x k = 0 j = 1 , 2 ,
where  p l ( x )  is degree l of the binary polynomial and  Φ x x k 2  is the radial basis function.
When  d = 2  and the degree of the polynomial is  l = 1 , Equation (3) is
s x = k = 1 n λ k Φ x x k 2 + μ 1 + μ 2 x + μ 3 y ,
with the constraint condition (4)
k = 1 n λ k = k = 1 n λ k x k = k = 1 n λ k y k = 0 .
Suppose that the function value at  x k  is  f k , k = 1 , , n ; then, the matrix form of Equations (5) and (6) is
Φ ( x 1 x 1 2 ) Φ ( x 1 x n 2 ) 1 x 1 y 1 Φ ( x 2 x 1 2 ) Φ ( x 2 x n 2 ) 1 x 2 y 2 Φ ( x n x 1 2 ) Φ ( x n x n 2 ) 1 x n y n 1 1 0 0 0 x 1 x n 0 0 0 y 1 y n 0 0 0 λ 1 λ 2 λ n μ 1 μ 2 μ 3 = f 1 f 2 f n 0 0 0 ,
In this formula,  A ^  is used to represent the matrix at the left end  ( n + 3 ) × ( n + 3 ) ; thus,
s ( x ) = [ Φ ( x x 1 2 ) Φ x x n 2 1 x y ] [ λ 1 λ n μ 1 μ 2 μ 3 ] T = [ Φ ( x x 1 2 ) Φ x x n 2 1 x y ] A ^ 1 [ f 1 f n 0 0 0 ] T ,
The linear operator L is used to calculate the result at  x = x c :
L s ( x ) | x = x c = [ L Φ ( x x 1 2 ) | x = x c L Φ x x n 2 | x = x c L 1 | x = x c L x | x = x c L y | x = x c ] A ^ 1 [ f 1 f n 0 0 0 ] T .
If  f i = 1  and  f j = 0 , then  j i ; thus,  L s ( x ) | x = x c = w i  and  i = 1 , , n , according to Formula (2). Therefore,
[ w 1 , , w n ] = [ L Φ ( x x 1 2 ) | x = x c L Φ x x n 2 | x = x c
L 1 | x = x c L x | x = x c L y | x = x c ] A ^ 1 D O ( n + 3 ) × n ,
where D is the identity matrix of  ( n × n )  and O is the null matrix of  ( 3 × n ) . An identity matrix is formed by adding three columns to the right matrix, which can be ignored. Thus, the left-hand matrix is added with  [ w n + 1 , w n + 2 , w n + 3 ] . The right-hand matrix is multiplied by  A ^  and then transposed, obtaining the following linear equations:
Φ ( x 1 x 1 2 ) Φ ( x 1 x n 2 ) 1 x 1 y 1 Φ ( x 2 x 1 2 ) Φ ( x 2 x n 2 ) 1 x 2 y 2 Φ ( x n x 1 2 ) Φ ( x n x n 2 ) 1 x n y n 1 1 0 0 0 x 1 x n 0 0 0 y 1 y n 0 0 0 w 1 w 2 w n w n + 1 w n + 2 w n + 3 = L Φ ( x x 1 2 ) x = x c L Φ ( x x 2 2 ) x = x c L Φ ( x x n 2 ) x = x c L 1 x = x c L x x = x c L y x = x c .
Therefore, the weights  w 1 w n  are obtained from the above equations. In the following numerical examples, the resulting linear system changes with the varying degree of the polynomial. In this paper, the PHS radial basis function is mainly considered;  Φ ( r ) = r m , m = 1 , 3 , 5 , , PHS is a piecewise smooth function, where  r = x x k 2 , and it is unnecessary to select a shape parameter  ε .

3.2. Radial Basis Function Finite Difference Method Based on Oseen Iteration

In this part, we introduce discretization of the Navier–Stokes equations using the radial basis function finite difference method based on Oseen iteration. First, we use the radial basis function finite difference method with polynomial to discretize the component equations. Then, we linearize the equations based on Oseen iteration starting from this scheme, and finally obtain the fully discrete scheme of the equations.

3.2.1. Discrete Scheme of Equation

First, the radial basis function finite difference method with polynomial is used to discretize the spatial operator of the equation. Equation (1) is expressed as the following components:
ν ( 2 u x 2 + 2 u y 2 ) + u u x + v u y + p x = f 1 ,
ν ( 2 v x 2 + 2 v y 2 ) + u v x + v v y + p y = f 2 ,
u x + v y = 0 .
For Equation (9), the differential value at a point  x c  is
L u x c = j = 1 n w j u x c j ,
where  x c j j n  is the n points around  x c  and  w j j n  is the weight coefficient corresponding to these points obtained by applying the interpolation condition  f x c j = f j .
First, we assume that there are N nodes in a zone. At a certain point  x c , Formula (9) is discretized according to RBF-FD interpolation and the following form is obtained:
ν j = 1 n 1 w j Δ u ( x c j ) + u ( x c ) j = 1 n 1 w j x u ( x c j ) + v ( x c ) j = 1 n 1 w j y u ( x c j ) + j = 1 n 2 w j x p ( x c j ) = f 1 ( x c ) ,
where  n 1  represents the number of nearest nodes at the central node  x c , used for approximating the weight coefficients of  Δ u  and  u , and  n 2  represents the number of nodes at the central node  x c , used to form a template of nearest nodes for approximating the weight coefficients of  p ; moreover,  w j Δ  is a weight coefficient for approximating  Δ u w j x  is a weight coefficient for approximating operator  x , and  w j y  is a weight coefficient for approximating operator  y . Next, we use the Oseen iterative scheme to deal with the nonlinear terms of the steady Navier–Stokes equation. Iterating according to Discrete Equation (13), we obtain the following form:
ν j = 1 n 1 w j Δ u ( x c j k + 1 ) + u ( x c k ) j = 1 n 1 w j x u ( x c j k + 1 ) + v ( x c k ) j = 1 n 1 w j y u ( x c j k + 1 ) + j = 1 n 2 w j x p ( x c j k + 1 ) = f 1 ( x c k + 1 ) ,
then
ν j = 1 n 1 w j Δ + u ( x c k ) j = 1 n 1 w j x + v ( x c k ) j = 1 n 1 w j y u ( x c j k + 1 ) + j = 1 n 2 w j x p ( x c j k + 1 ) = f 1 ( x c k + 1 ) ,
where k denotes the number of iterations,  u ( x c k )  denotes the value of u at the kth iteration,  u ( x c k + 1 )  denotes the value of u at the (k + 1)th iteration, and the initial value of the velocity  u 0  is obtained from the corresponding Stokes equation. In Equation (15), the value of u at the kth level is known; thus, the nonlinear term in Equation (9) is linearized. Likewise, there are similar discrete forms for Equations (10) and (11):
ν j = 1 n 1 w j Δ + u ( x c k ) j = 1 n 1 w j x + v ( x c k ) j = 1 n 1 w j y v ( x c j k + 1 ) + j = 1 n 2 w j y p ( x c j k + 1 ) = f 2 ( x c k + 1 ) ,
j = 1 n 1 w j x u ( x c j k + 1 ) + j = 1 n 1 w j y v ( x c j k + 1 ) = 0 ,
There are N nodes in a zone including  N B  boundary nodes and  N I  internal nodes, where  N = N I + N B . The central node  x c  in Formulas (15) and (16) is considered on  N I  internal nodes, and the central node  x c  in Formula (17) is considered on N total nodes. Considering the boundary condition  u = 0 , Equations (15)–(17) and the boundary conditions form the following equations:
W N I × N I Δ + W N I × N B Δ + O N I × N I O N I × N B W N I × N x O N B × N I E N B × N B O N B × N I O N B × N B O N B × N O N I × N I O N I × N B W N I × N I Δ + W N I × N B Δ + W N I × N y O N B × N I O N B × N B O N B × N I E N B × N B O N B × N W N × N I x W N × N B x W N × N I y W N × N B y O N × N U I U B V I V B P = F 1 0 F 2 0 0 ,
We can now write (18) in the form of a block matrix:
W u Δ + O W p x O W v Δ + W p y W u x W v y O U V P = F 1 F 2 0 .
The stabilization parameter  ϵ  is added to Equation (19) to obtain a system of equations in the form of BU = F:
W u Δ + O W p x O W v Δ + W p y W u x W v y ϵ D U V P = F 1 F 2 0 .
where B is the left-end weight matrix,  W u Δ +  and  W v Δ +  are  ( N × N )  matrices composed of weight coefficients and boundary condition coefficients used to approximate u and v in Equations (15) and (16), respectively,  W u x  and  W v y  are a matrix of  ( N × N )  composed of weighted coefficients for approximating  u x  and  v y W p x  and  W p y  are a matrix of  ( N × N )  composed of weighted coefficients for approximating  p x  and  p y F 1  and  F 2  are known right-end terms, U is the solution of velocity u on point set X, V is the solution of velocity v on point set X, P is the solution of pressure p on point set X, D is the unit matrix of  ( N × N ) , and  ϵ  is the stabilization parameter,  ϵ  = 1.0 × 10 6 .

3.2.2. Selection of Degrees of Added Polynomial and Number of Template Nodes

According to [13,14], without depending on the dimension, when using radial basis function PHS:  Φ ( r ) = r p , when p is an odd number and the l-degree polynomials are supplemented to approximate the k-order derivative the convergence rate is not determined by the order of PHS, and is instead determined by the highest degree of the polynomial used. The solution accuracy is  O ( h l k + 1 ) . For the PDE with only the first-order spatial derivative, the convergence speed can reach  O ( h l ) , while for the PDE with the second-order spatial derivative, the convergence speed can reach  O ( h l 1 ) . The number of terms of the l-degree polynomial is  ( l   +   1 ) ( l   +   2 ) 2 . When approximating the value of an operator acting on the center point of the template, symmetry of the node templates is beneficial, as symmetrical node templates provide uniform information. In this paper, the numerical examples are considered in two layouts, namely, a right-angle node layout and a hexagonal node layout, with the node layouts shown in Figure 1. The right-angled node layout has completely symmetrical templates on numbers 5, 9, 13, 21, 25, 29, 37, 45, etc., while the hexagonal node layout has completely symmetrical templates on numbers 7, 13, 19, 31, 37, 43, 55, 61, etc. In this paper, the minimum distance between two points under the right-angled node layout is  h = 1 N     1 , while N is the total number of nodes.

4. Numerical Method

In this section, several numerical examples are used to demonstrate the effectiveness of the proposed methods; the relative error  L 2  and error  L  are applied for comparison of the numerical solution and true solution. The relative error  L 2  is defined as  e r r o r L 2 = u u h 2 u 2 , while the error  L  is defined as  e r r o r L = u u h L ( Ω ) = m a x x X | u ( x ) u h ( x ) | , where  u h  is the obtained numerical solution and  u  is the true solution  x X .

4.1. Convergence Test

This example studies the two-dimensional Navier–Strokes equations under the layout of right-angled nodes. We assume that Equation (1) has the true solutions
u = x 2 2 x 3 + x 4 2 y 6 y 2 + 4 y 3 on Ω , v = y 2 2 y 3 + y 4 2 x 6 x 2 + 4 x 3 on Ω , p = cos π x sin π y on Ω .
where  Ω = 0 , 1 × 0 , 1 , the viscosity coefficient  ν = 1 , and the nodes are uniformly distributed on  Ω  as shown in Figure 1 left. For purpose of the convergence test, the basis function  Φ ( r ) = r 7  is applied to add third-order and second-order polynomials for approximation of the velocity  u = ( u , v ) , while the basis function  Φ ( r ) = r 5  is applied to add second-order and first-order polynomials for approximation of the pressure  p . For the number of node templates approximated at a central node of the velocity and pressure in Equations (15) and (16), we select  n 1 = 21  and  n 2 = 9 , respectively.
Figure 2 left shows the addition of the third-order polynomial for approximation of the velocity and the error results for pressure approximation when adding the second-order polynomial. The figure shows the relative error  L 2  when the number of nodes N is 25, 81, 289, 1089, and 4225. The error decreases with increasing number of nodes, and the convergence order of velocity and pressure reaches third-order and second-order. When the number of nodes N is 4225, the velocity error is  4.4583 × 10 5  and the pressure error is  2.2254 × 10 2 . The right0hand side of Figure 2 shows the error results of the approximation when adding the second-order polynomial to the velocity and adding the first-order polynomial to the pressure. When the number of nodes N is 4225, the relative errors  L 2  of the velocity and pressure are  4.2393 × 10 5  and  2.8000 × 10 2 , respectively.
By comparing the two figures, the error between the numerical solution and the true solution obtained from approximation by adding third-order and second-order polynomials to the velocity and pressure is far less than that obtained from approximation by adding second-order and first-order polynomials to the velocity and pressure. When the number of nodes is less than 1089, the former reaches the corresponding convergence order. As the number of nodes approaches 4225, the order of velocity is less than 2 due to the reduction in point density. This example illustrates that under the right-angled node layout, the radial basis function difference method based on Oseen iteration is convergent and effective for solving incompressible Navier–Stokes equations.

4.2. Small Viscosity Problem

In order to test the suitability of the method proposed in this paper for small viscosity problems, we begin by assuming that Equation (1) has a true solution (21). By setting the viscosity coefficient  ν = 0.00001 , the nodes are distributed in a right-angled node layout. In this example, the basis function  Φ ( r ) = r 7  is applied to add the third-order polynomial for approximation of the velocity  u = ( u , v ) , and the basis function  Φ ( r ) = r 5  is applied to add the second-order polynomial for approximation of the pressure p. Table 1 shows the error results when  ν  is 0.00001 and number of nodes N is 25, 81, 289, 1089, and 4225. When the number of nodes N is 4225, the relative error  L 2  of the velocity is  6.7100 × 10 5  and the relative error  L 2  of the pressure is  2.0271 × 10 2 . Figure 3 shows the corresponding error order.

4.3. Hexagonal Node Layout

In this section, in order to test the suitability of the method proposed in this paper for different node layouts, we consider the two-dimensional equation under the hexagonal node layout, as shown in Figure 1 right, assuming that Equation (1) has a true solution (21) where  Ω = 0 , 1 × 0 , 1  and the viscosity coefficient  ν = 1 . In this example, the basis function  Φ ( r ) = r 7  is applied to add the third-order polynomial for approximation of the velocity  u  and the basis function  Φ ( r ) = r 5  is applied to add the second-order polynomial for approximation of the pressure p. For the velocity and pressure in a discrete equation, the number of node templates approximated at a central node ( n 1  and  n 2 ) are separately selected as  n 1 = 31  and  n 2 = 13 . When number of nodes is 4699, the relative error  L 2  of the velocity is  2.8638 × 10 5  and the relative error  L 2  of the pressure is  1.3900 × 10 2 . Figure 4 shows the relative error  L 2  between the numerical solution and the true solution when the number of nodes N is 68, 279, 1166, and 4699.
In summation, these numerical experiments show that the radial basis function finite difference method based on Oseen iteration has good performance for solving Navier–Strokes equations under the hexagonal node layout.

5. Conclusions

This paper proposes a radial basis function finite difference method based on Oseen iteration for solving two-dimensional steady Navier–Stokes equations with discretization of the spatial operators of the Navier–Stokes equations using the radial basis function finite difference method with polynomial, then linearizing the equation based on Oseen iteration. In this paper, we provide numerical solutions for the Navier–Stokes equations under right-angled and hexagonal node layouts, analyze and compare the influence of the number of polynomial additions on the accuracy and convergence of the solution, and verify the effectiveness of the proposed method. Under the two node layouts studied in this paper, the proposed method offers high-precision numerical solutions for solving Navier–Stokes equations while demonstrating good performance. Its extension to unsteady problems will be the focus of our future research work.

Author Contributions

Conceptualization, X.F.; methodology, L.M.; software, L.M.; validation, X.F.; formal analysis, L.M.; investigation, L.M.; resources, X.F.; data curation, L.M.; writing—original draft preparation, L.M.; writing—review and editing, X.F. and L.M.; visualization, L.M.; supervision, X.F.; project administration, X.F.; funding acquisition, X.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Teman, R. Navier-Stokes Equations: Theory and Numerical Analysis; Springer: Amsterdam, The Netherland; New York, NY, USA; Oxford, UK, 1984. [Google Scholar]
  2. He, Y.; Miao, H.; Mattheij, R. Numerical analysis of a modified finite element nonlinear Galerkin method. Numer. Math. 2004, 97, 725–756. [Google Scholar] [CrossRef]
  3. He, Y.; Wang, A.; Mei, L. Stabilized finite-element method for the stationary Navier-Stokes equations. J. Eng. Math. 2005, 51, 367–380. [Google Scholar] [CrossRef]
  4. He, Y.; Feng, X. Research progress of efficient fully discrete finite element method for 3D unsteady Navier-Stokes equations. J. Xinjiang Univ. 2022, 39, 257–265. [Google Scholar]
  5. Raissi, M.; Perdikaris, P.; Karniadakis, G. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. J. Comput. Phys. 2019, 378, 686–707. [Google Scholar] [CrossRef]
  6. Bayona, V.; Moscoso, M.; Carretero, M. RBF-FD formulas and convergence properties. J. Comput. Phys. 2010, 229, 8281–8295. [Google Scholar] [CrossRef]
  7. Bayona, V.; Moscoso, M.; Kindelan, M. Optimal constant shape parameter for multiquadric based RBF-FD method. J. Comput. Phys. 2011, 230, 7384–7399. [Google Scholar] [CrossRef]
  8. Bayona, V.; Kindelan, M. Propagation of premixed laminar flames in 3D narrow open ducts using RBF-generated finite differences. Combust. Theory Model. 2013, 17, 789–803. [Google Scholar] [CrossRef]
  9. Bayona, V.; Flyer, N.; Fornberg, B. On the role of polynomials in RBF-FD approximations, II: Numerical solution of elliptic PDEs. J. Comput. Phys. 2017, 332, 257–273. [Google Scholar] [CrossRef]
  10. Bayona, V.; Flyer, N.; Fornberg, B. On the role of polynomials in RBF-FD approximations, III: Behavior near domain boundaries. J. Comput. Phys. 2019, 38, 378–399. [Google Scholar] [CrossRef]
  11. He, Y. Stability and convergence of iterative methods related to viscosities for the 2D/3D steady Navier-Stokes equations. J. Math. Anal. Appl. 2015, 423, 1129–1149. [Google Scholar] [CrossRef]
  12. He, Y.; Li, J. Convergence of three iterative methods based on the finite element discretization for the stationary Navier-Stokes equations. Comput. Methods Appl. Mech. Eng. 2009, 198, 1351–1359. [Google Scholar] [CrossRef]
  13. Flyer, N.; Fornberg, B.; Bayona, V. On the role of polynomials in RBF-FD approximations, I: Interpolation and accuracy. J. Comput. Phys. 2016, 321, 21–38. [Google Scholar] [CrossRef]
  14. Flyer, N.; Barnett, G.; Wicker, L. Enhancing finite differences with radial basis functions: Experiments on the Navier-Stokes equations. J. Comput. Phys. 2016, 316, 39–62. [Google Scholar] [CrossRef]
Figure 1. Right-angled node layout (left) and hexagonal node layout (right).
Figure 1. Right-angled node layout (left) and hexagonal node layout (right).
Entropy 25 00804 g001
Figure 2. Relative error  L 2  of the velocity u and pressure p (right-angled node layout).
Figure 2. Relative error  L 2  of the velocity u and pressure p (right-angled node layout).
Entropy 25 00804 g002
Figure 3. Relative error  L 2  of the velocity  u  (left) and the pressure p (right) ( ν = 0.00001 ).
Figure 3. Relative error  L 2  of the velocity  u  (left) and the pressure p (right) ( ν = 0.00001 ).
Entropy 25 00804 g003
Figure 4. Relative error  L 2  of the velocity  u  and pressure p (hexagonal node layout).
Figure 4. Relative error  L 2  of the velocity  u  and pressure p (hexagonal node layout).
Entropy 25 00804 g004
Table 1. Relative error  L 2  and error  L  between numerical solution and true solution of the velocity u and pressure p ( ν = 0.00001 ).
Table 1. Relative error  L 2  and error  L  between numerical solution and true solution of the velocity u and pressure p ( ν = 0.00001 ).
NRelative Error  L 2  of uRelative Error  L 2  of p  L  Error of u L  Error of p
25   2.9942 × 10 2   2.4962 × 10 1   9.9977 × 10 1   1.0191 × 10 1
81   1.1883 × 10 2   1.5148 × 10 1   1.9150 × 10 1   3.3872 × 10 2
289   1.5680 × 10 3   7.9769 × 10 2   1.2605 × 10 2   9.4008 × 10 3
1089   1.5564 × 10 4   4.0410 × 10 2   6.2549 × 10 4   2.4502 × 10 3
4225   6.7100 × 10 5   2.0271 × 10 2   1.3483 × 10 4   6.2380 × 10 4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mu, L.; Feng, X. Radial Basis Function Finite Difference Method Based on Oseen Iteration for Solving Two-Dimensional Navier–Stokes Equations. Entropy 2023, 25, 804. https://doi.org/10.3390/e25050804

AMA Style

Mu L, Feng X. Radial Basis Function Finite Difference Method Based on Oseen Iteration for Solving Two-Dimensional Navier–Stokes Equations. Entropy. 2023; 25(5):804. https://doi.org/10.3390/e25050804

Chicago/Turabian Style

Mu, Liru, and Xinlong Feng. 2023. "Radial Basis Function Finite Difference Method Based on Oseen Iteration for Solving Two-Dimensional Navier–Stokes Equations" Entropy 25, no. 5: 804. https://doi.org/10.3390/e25050804

APA Style

Mu, L., & Feng, X. (2023). Radial Basis Function Finite Difference Method Based on Oseen Iteration for Solving Two-Dimensional Navier–Stokes Equations. Entropy, 25(5), 804. https://doi.org/10.3390/e25050804

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop