Navier-Stokes equations numerical solutions

There are a few analytical solutions to Navier-Stokes equations. This explains why solving Navier-Stokes equations numerically has become a discipUne per se in the field of fluid mechanics. However, notion of a solution to Navier-Stokes equations calls for further discussion. 

Taylor, C. and Hood, P., 1973. A numerical solution of the Navier-Stokes equations using the finite element technique. Comput. Fluids 1, 73-100. 

Gresho, P. M., Lee, R. L. and Sani, R. L., 1980. On the time-dependent solution of the incompressible Navier-Stokes equations in two and three dimensions. In Recent Advances in Numerical Methods in fluids, Ch. 2, Pineridge Press, Swansea, pp. 27-75. 

Because the Navier-Stokes equations are first-order in pressure and second-order in velocity, their solution requires one pressure bound-aiy condition and two velocity boundaiy conditions (for each velocity component) to completely specify the solution. The no sBp condition, whicn requires that the fluid velocity equal the velocity or any bounding solid surface, occurs in most problems. Specification of velocity is a type of boundary condition sometimes called a Dirichlet condition. Often boundary conditions involve stresses, and thus velocity gradients, rather than the velocities themselves. Specification of velocity derivatives is a Neumann boundary condition. For example, at the boundary between a viscous liquid and a gas, it is often assumed that the liquid shear stresses are zero. In numerical solution of the Navier-... 

Computational fluid dynamics (CFD) emerged in the 1980s as a significant tool for fluid dynamics both in research and in practice, enabled by rapid development in computer hardware and software. Commercial CFD software is widely available. Computational fluid dynamics is the numerical solution of the equations or continuity and momentum (Navier-Stokes equations for incompressible Newtonian fluids) along with additional conseiwation equations for energy and material species in order to solve problems of nonisothermal flow, mixing, and chemical reaction. 

These values can be compared to predicted values for numerical solutions of the incompressible Navier-Stokes equations. For d = 2, for example, we have the lower bounds Sa=2 (TZ/M) and Wd=2 TZ /M for a LG and the bounds S num, d=2 and d=2 where the bounds for the numerical solutions... 

The strategies discussed in the previous chapter are generally applicable to convection-diffusion equations such as Eq. (32). If the function O is a component of the velocity field, the incompressible Navier-Stokes equation, a non-linear partial differential equation, is obtained. This stands in contrast to O representing a temperature or concentration field. In these cases the velocity field is assumed as given, and only a linear partial differential equation has to be solved. The non-linear nature of the Navier-Stokes equation introduces some additional problems, for which special solution strategies exist. Corresponding numerical techniques are the subject of this section. 

Perhaps even more important is die fact that LEM does not require a numerical solution to die Navier-Stokes equation. Indeed, even a three-dimensional diffusion equation is generally less computationally demanding than the Poisson equation needed to find die pressure field. 

Thus, analytic solutions for flow around a spherical particle have little value for Re > 1. For Re somewhat greater than unity, the most accurate representation of the flow field is given by numerical solution of the full Navier-Stokes equation, while empirical forms should be used for C. These results are discussed in Chapter 5. 

At larger Re and for more marked deformation, theoretical approaches have had limited success. There have been no numerical solutions to the full Navier-Stokes equation for steady flow problems in which the shape, as well as the flow, has been an unknown. Savic (S3) suggested a procedure whereby the shape of a drop is determined by a balance of normal stresses at the interface. This approach has been extended by Pruppacher and Pitter (P6) for water drops falling through air and by Wairegi (Wl) for drops and bubbles in liquids. The drop or bubble adopts a shape where surface tension pressure increments, hydrostatic pressures, and hydrodynamic pressures are in balance at every point. Thus... 

One important use of the stream function is for the visualization of flow fields that have been determined from the solution of Navier-Stokes equations, usually by numerical methods. Plotting stream function contours (i.e., streamlines) provides an easily interpreted visual picture of the flow field. Once the velocity and density fields are known, the stream function field can be determined by solving a stream-function-vorticity equation, which is an elliptic partial differential equation. The formulation of this equation is discussed subsequently in Section 3.13.1. Solution of this equation requires boundary values for l around the entire domain. These can be evaluated by integration of the stream-function definitions, Eqs. 3.14, around the boundaries using known velocities on the boundaries. For example, for a boundary of constant z with a specified inlet velocity u(r),... 

In all the above derivations in this section, the influence of viscosity is neglected so that analytical solutions for velocity and pressure profiles can be obtained. When the viscosity of fluid is taken into account, it is difficult to obtain any analytical solution. Kuts and Dolgushev [35] solved numerically the flow field in the impingement of two axial round jets of a viscous impressible liquid ejected at the same velocity from conduits with the same diameter and located very close to each other. The mathematical formulation incorporated the complete Navier-Stokes equations transformed into stream and velocity functions in cylindrical coordinates r and z, with the assumption that the velocity profiles at the entrance and the exit of the conduit were parabolic. The continuity equation is given by Eq. (1.22) and the equations for motion in dimensionless form are ... 

A numerical analysis using FlumeCAD was made, solving the incompressible Navier-Stokes equation for the velocity and pressure fields [70], The steady-state velocity field was then used in the coupled solution of three species transport equations (two reagents and one product). Further details are given in [70],... 

The numerical results are based on the solution of the incompressible Navier-Stokes equation ... 

P 61] The numerical simulations were based on the solution of the incompressible Navier-Stokes equation and a convection-diffusion equation for a concentration field by means of the finite-volume method [152], The Einstein convention of summation over repeated indices was used. For pressure-velocity coupling, the SIMPLEC algorithm and for discretization of the species concentration equation the QUICK differencing scheme were applied. Hybrid and the central differencing schemes referred to velocities and pressure, respectively (commercial flow solvers CFX4 and CFX5). 

We should note that the Navier-Stokes equation holds only for Newtonian fluids and incompressible flows. Yet this equation, together with the equation of continuity and with proper initial and boundary conditions, provides all the equations needed to solve (analytically or numerically) any laminar, isothermal flow problem. Solution of these equations yields the pressure and velocity fields that, in turn, give the stress and rate of strain fields and the flow rate. If the flow is nonisothermal, then simultaneously with the foregoing equations, we must solve the thermal energy equation, which is discussed later in this chapter. In this case, if the temperature differences are significant, we must also account for the temperature dependence of the viscosity, density, and thermal conductivity. 

---