Investigation of 2-Dimensional Fluid Flow using Finite Difference Flow Method of Navier Stokes Equation

Investigation of 2-Dimensional Fluid Flow using Finite Difference Flow Method of Navier Stokes Equation

Rakshit Kaushik, Dr. Mohit Kumar Srivastava

Under the guidance of

Shri Vishwakarma Skill University

Abstract:- The presented research paper illustrates the development of a new methodology to solve 2-dimentional (2D) Navier-Stoke equations, which Pukhnachev proposed through introducing unknown functions in the stream and pressure functions of fluid flow. The proposed novel method is distinguishable from the common vorticity-stream given in Navier-Stokes expression because it has a stream function that corresponds to the unknown function in the elliptic expression. The equation represents a couple scheme in algorithmic considerations because it enables two situations to be solved using one function of the subject stream without putting new conditions on the innovative function. Here, the concept of numerical algorithm is applied in a flow under cavity to represent a benchmark task to be solved. The benchmark task gives enough representation of the subject flow.

Key Words: Incompressible viscous flow, Navier-Stoke equations, finite-different scheme

INTRODUCTION

The Navier-Stoke equations can be solved using various numerical schemes when dealing with incompressible viscous flows. Some the numerous schemes are those that use variables commonly known as primitive or velocity-pressure and stream function. Others in the same category uses vorticity-steam and other formulations. Usually, finding numerical solutions when using the primitive variable poses a great challenge because the evolution expression will be missing in most cases, which can help solve the variable related to pressure. Therefore, the arising problem with using approach of primitive variables can be avoided by using vorticity-velocity and stream function vorticity formulations of the Navier-Stoke equation. However, the right values for vorticity boundary are not easy to determine. This study uses a new form of Navier-Stokes equation.

Pukhnachev proposed a new form of Navier-Stoke equation to solve 2D and axisymmetric flows. The novel Navier- Stoke expression in the 2D stream considering the unknown function has one expression for transportation that represent an elliptic mathematical expression and the functionality of the stream for the unknown. The stream function corresponds to both the stream and vorticity functions. However, the coupling function has a different physical meaning from the vorticity. A finite difference scheme is structed for the Navier-Stoke equation in the new methodology development. The algorithm considers the mathematical expression as a coupled organization and allows satisfaction of double situation for the function of the stream without putting condition on the undetermined function. It is possible to extend the proposed scheme to solve axisymmetric Navier-Stoke equations. The performance of this proposed method can be determined using a popular benchmark problem. One case that is applicable in the test is the simulation of a 2D lid-driven cavity flow under Reynolds number Re 1000, under strict steady and laminar motion [3]. So far, research on numerical has revealed the properties of this cavity flow. Additionally, the flow of viscous stream in a forced cavity is one of the common test concepts in determining and verifying numerical problems.[3]

In the paper, several sections are included to illustrate various parts. The first part deals with the formulation of the proposed novel Navier-Stokes expression with no-slip boundary limits.

After that, a brief account of a task applied in the trial case is given with arithmetical algorithm. The next section gives the authentication findings of finite-different design with detailed comparison to tentative and numerical statistics.

FORMULATION OF NOVEL NAVIER-STOKES EQUATION

It is essential to illustrate the steps involved in the formulation of novel equation for Navier-Stokes regarding incompressible viscosity in 2D, which makes the paper comprehensive. The Navier-Stokes expression governs the incompressible viscosity of the flow in a Cartesian coordinate presentation (x, y),

………………………………………………………………..Equation (1)

…………………………………………………………………Equation (2)

…………………………………………………………………………………………………Equation (3)

In the above equations, the parts of velocity are represented by u and in x-and y- routes in that order. = pressure, = fluid density, and = kinematic viscosity. Notably, the fluid is under potential external forces[4]. Considering the flow of stream under 2D, the compressibility limit can be satisfied through expressing vector of speed regarding the purpose of stream according to:

……………………………………………………………………………………….Equation (4)

From here, a new design of the popular Navier-Stokes equation can be pegged to the following scrutiny. When equation (4) is substituted in equation (1), it gives:

………………………………………………………………Equation (5)

Where

Therefore, a function or ( ) satisfies the equation.

and

………………………………………………………………………………………………..Equation (6)

…………………………………………………………………………………….Equation (7)

When equations (6) and (7) are differentiated with respect to y and x, and the result substituted into equation (2), where u and are expressed in terms of , it yields

…………………………………………………………………………………………………Equation (8)

Here, the case that correspond the no-slip state is considered as the limit of the flow sphere. Considering only as a function of the stream, the conditions of the limits are

…………………………………………………………………………………………..Equation (9)

Where represent a vector with normal derivative towards the direction of the boundary. The formulation of the task can be completed through identifying the original forms

…………………………………………………………………………………………Equation (10)

The key objective here is to develop and confirm a finite-difference design for solving the system (7) to (10) [1][2].

Numerical Technique

The lid-driven cavity flow is the standard problem for benchmarking verifying 2D Navier-Stokes equations shown below (Figure 1). The upper wall of the cavity pushes the fluid and results in a steady motion towards the right from left [4]. Let L represent the characteristic scale related to the moving boundary.

Figure 1: Configuration of lid-driven square cavity flow

The non-dimensional parameter of the problem is

……………………………………………………………………………………………………Equation (11)

The system of equations (1) to (3) is rendered dimensionless as follows

………………………………………………………Equation (12)

The domain is covered using a uniform grid

with spacing

in the directions of x- and the y- axes accordingly. The given grid allows the use of median difference in approximating the conditions of the boundary with second-order on double-point patterns.

The sugested algorithm has one essential element that is considered as a pair of systems given in mathematic expressions

(7) and (8) for and The new formulation used here is based on the concept of the two border situations for as real situations for the – system. Notably, and are assessed using permanent procedure. The approximation for second- order central difference is used for the operators in equation (7) and (8). The difference equation for the system becomes

………………………………………………………………Equation (13)

……………………………Equation (14)

The conditions for the boundary are written as given below:

…………………………………Equation (15)

It is possible to add (13) to (15) and obtain one linear scheme that is banded with matrix. The matrix enable an introduction of a new scheme of indices.

Each of the nodules of grid Qh, especially (i, j) is associates with indices including k(i,j) and even m(i,j). Note that the index k(i,j)

represents an odd number while m(i,j) represents an even number. Therefore,

…………………………………………………………Equation (16a)

………………………………………………………Equation (16b)

A new grid function can be introduced as defined in the composite grid. Taking the components of the grid function to be with even indices representing i,j and components with odd indices

to represent can be substituted instead of i,j and substituted instead of in equation (13) and (14), which

allows recasting of algebraic system as follows

……………………..Equation (17a)

…………………….Equation (17b)

Where

A simple implementation of the algorithm results in a problem with a single numerical matrix. The resultant problem

can be standardized in different ways. Notably, addition of a small element to the boundary gives the best results [5]. Using this idea, equation (15) can be rearranged for function

………………………………Equation (18a)

…………………………….Equation (18b)

…………………………….Equation (18c)

…………………………..Equation (18d)

Where is a small number. For steady flow. The algorithm is considered as an interactive step and the interaction is terminated at n-N satisfying the following criteria.

Notably, the coupling of the system results in the formation of a linear system , which can be rearranged based on multi- diagonal arrangement for the function of grid composite

…………………………..Equation (19)

Where Equation 19 gives a linear system with matrix that is bound with on both upper side and lower side. The problem given in equation 19 can be solved using routings of LAPACK.

RESULTS AND DISCUSSION

A confirmation test involve flow in 2D cavity at Reynolds number of up to 1000, where they give steady and laminar flow. Calculations were done for the lid-driven cavity task on grid from 32 * 32 to 102 * 102. , which represents a small number, has a sensible impact on the findings of the test, which were evaluated using the numerical analysis. The analysis showed that . The results prove that the solution was compatible with the results of the test cases. The validation of the subject scheme was done through comparing values of external points and those of locations in space functions of the target stream. The target stream was having [2,3,4,5, and 6]. Consider table 1. The top row shows values form simulation while the bottom row shows values from other sources. For the primary vortex, the results showed agreement within 5% as compared to values from other sources. For the secondary vortex, the 52 * 52 grid, portrayed results agree within tolerable 5% [5]. However, the values were different as compared to (Sellountos et al., 2019)[6]. The data from 102 * 102 grid revealed perfect match with those data from [2,3,4,6]. The structures for geometry of the flow are shown in figures 2 and 3 using the 102 * 102 grid. Notably, the values of velocity (u) along x- direction, velocity (v) along y- axis, and w were calculated from the after the convergence of the interaction. The values of u, v, and w in the domain were approximated using central difference and the same values on the boundary were approximated using difference of one side first-order. Figure 4 reveals a comparison between centerline u- and v- profiles for velocities.

Table 1: Strength and location of vortices at Re 100 to 1000

It is notable that understanding the behavior of function for various Re. Figure 5 illustrates the contour of the function for several Re and the set of figures given from 52 * 52 grid with parameters .

Figure 2: Functions of stream , vorticity w, and velocity components u and v contours for the lid-driven flow at Re-100.

CONCLUSION

Remarkably, the development of finite difference sachem was successful and it confirmed the validity of the proposed novel formulation using the 2D Navier-Stokes equation. The development revealed that the stream function was reliable enough to be used in evaluating the equation, which is a distinguishing property of the novel development in the application of Navier- Stokes equation. It is notable that the developed algorithm and its counterpart scheme considered the subject equation as a coupled system, and it allows satisfaction of two conditions when dealing with functions when dealing with function streams without putting new conditions on the derived function. The formulation portrayed high accuracy and reliable efficiency for the benchmark problem in lid-driven cavity flow.

REFERENCES

1. Anon, Vorticitystream function formulation of navierstokes equation. Krylov Subspace Methods with Application in Incompressible Fluid Flow Solvers, 2020, pp.177217.

2. Anon, Navierstokes flow. Introductory Incompressible Fluid Mechanics, 2021, pp.143181.

3. Michael Quinn, Quinns Law of Fluid Dynamics: Supplement #3 a unique solution to the Navier-stokes equation for fluid flow in closed conduits. Fluid Mechanics, 2020, p.30.

4. Peter, E.E. & Bala, A., Application of navier stokes equation to solve fluid flow problems. Journal of Advances in Mathematics and Computer Science, 2020, pp.101114.

5. Rouzbar, R. & Eyi, S., Scramjet combustor flow analysis using coupled Navier-Stokes and finite rate chemical equations. 46th AIAA Fluid Dynamics Conference. 2016.

6. Sellountos, E.J., Tiago, J. & Sequeira, A., Meshless velocity vorticity local boundary integral equation (LBIE) method for two dimensional incompressible navier-stokes equations. International Journal of Numerical Methods for Heat & Fluid Flow, 2019, 29(11), pp.40344073.