Finite element analysis of vortex shedding using equal order interpolations

An operator splitting and element‐by‐element conjugated gradient solver, and equal order interpolations are applied for solving time dependent Navier–Stokes (NS) equations to simulate flow induced vortex shedding in the present study. In addition, the convection term is corrected by balanced tensor diffusivity, which can stabilize the numerical simulation and overcome the numerical oscillations. The evolution of the interested flowing properties with time is analyzed by using spectral analysis. The developed code has been validated by the application of two examples: a driven cavity flow and a flow induced vortex vibration. Results from the first example for Reynolds number Re=10³ and Re=10⁴ are compared with other numerical simulations. Results from the second example, uniform flow past a square rod over a wide range of high Reynolds numbers from Re=10³～10⁵, are compared with experimental data and other numerical studies. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical study of mixed convection and conduction in a 2-D square ventilated cavity with an inlet at the vertical glazing wall and outlet at the top surface

A steady state numerical study of combined laminar mixed convection and conduction heat transfer in a ventilated square cavity is presented. The air inlet gap is located at the bottom of a vertical glazing wall and air exits the cavity via a gap located at the top surface. Three locations for the opening at the top surface: left (case a), center (case b) and right side (case c) are considered. All the remaining surfaces are considered adiabatic. The mass, momentum and energy conservation equations were solved using the finite volume method for different Rayleigh numbers in the interval of 10⁴ < Ra < 10⁶ and Reynolds number in the interval of 100 < Re < 700. Temperature, flow field, and heat transfer rates are analyzed. The effect of the interaction between ambient conditions outside the glazing and the air inlet gap at the bottom for different air outlet gap positions at the top surface modifies the flow structure and temperature distribution of the air inside the cavity. The Nusselt number as a function of the Reynolds number was determined for the three cases. It was found that configuration for case (a) removes a higher amount of heat entering the cavity compared to cases (b) and (c). This is due to the short distance between the main stream and the glass wall surface. Thus, the forced airflow entering the cavity is assisted by the buoyancy forces, and most of the cavity remains at the inlet flow temperature, which should be appropriate for warm climates. These results may provide useful information about the heat transfer and fluid flow for future studies.     

Developing a physical influence upwind scheme for pressure-based cell-centered finite volume methods

In the framework of a cell-centered finite volume method (FVM), the advection scheme plays the most important role in developing FVMs to solve complicated fluid flow problems for a wide range of Reynolds numbers. Advection schemes have been widely developed for FVMs employing pressure-velocity coupling methodology in the incompressible flow limit. In this regard, the physical influence upwind scheme (PIS) is developed for a cell-centered finite volume coupled solver (FVCS) using a pressure-weighted interpolation method for linking the pressure and velocity fields. The well-known exponential differencing scheme and skew upwind differencing scheme are also deployed in the current FVCS and their numerical results are presented. The accuracy and convergence of the present PIS are evaluated solving flow in a lid-driven square cavity, a lid-driven skewed cavity, and over a backward-facing step (BFS). The flow within the lid-driven square cavity is numerically solved at Reynolds numbers from 400 to 10 000 on a relatively coarse mesh with respect to other reported solutions. The lid-driven skewed cavity test case at Reynolds number of 1000 demonstrates the numerical performance of the present PIS on nonorthogonal grids. The flow over a BFS at Reynolds number of 800 is numerically solved to examine capabilities of current FVCS employing the current PIS in inlet-outlet flow computations. The numerical results obtained by the current PIS are in excellent agreement with those of benchmark solutions of corresponding test cases. Incorporating implicit role of pressure terms in a pressure-weighted interpolation method and development of PIS provides satisfactory solution convergence alongside the numerical accuracy for the current FVCS. A particular numerical verification is performed for the V velocity calculation within the BFS flow field, which confirms the reliability of present PIS.     

Mixed convection cooling of a heated circular cylinder by laminar upward-directed slot jet impingement

An experimental and numerical study has been carried out to investigate the heat transfer characteristics of a horizontal circular cylinder exposed to a slot jet impingement of air. A square-edged nozzle is mounted parallel with the cylinder axis and jet flow impinges on the bottom of the cylinder. The study is focused on low Reynolds numbers ranging from 120 to 1,210, Grashof numbers up to Gr = 10Re ² and slot-to-cylinder spacing from 2 to 8 of the slot width. The flow field is greatly influenced by the slot exit velocity and the buoyancy force due to density change. A Mach–Zehnder Interferometer is used for measurement of local Nusselt number around the cylinder at 10° interval. It is observed that the average Nusselt number decreases with increasing the jet spacing and increases with rising the Reynolds number. A finite volume method utilizing a curvilinear coordinate transformation is used for numerical modeling. The numerical results show good agreement with the experimental results. The flow and thermal field are seen to be stable and symmetric around the cylinder over the range of parameters studied.     

Robust semidirect finite difference methods for solving the Navier–Stokes and energy equations

Semidirect solution techniques can be an effective alternative to the more conventional iterative approaches used in many finite difference methods. This paper summarizes several semidirect techniques which generally have not been applied to the Navier–Stokes and energy equations in finite difference form. The methods presented use both successive substitution and Jacobian-based updates as well as two variations of Broyden's full matrix update. A hybrid method is also presented, as is a norm-reducing search technique that can be used to enhance the convergence characteristics of any semidirect approach. These methods have been compared with the well known iterative methods SIMPLE and SIMPLER. The comparison was performed on the natural convection and driven cavity problems. The semidirect methods proved to be reliably convergent without the need for a priori specification of variable under-relaxation factors, which was necessary with the iterative methods. Natural convection and driven cavity solutions have been readily obtained with the proposed methods for Rayleigh and Reynolds numbers up to 10⁹ and 10⁶ respectively. Of the semidirect techniques, the hybrid approach was the most robust. From an arbitrary zero initial guess this method was able to obtain a solution to the natural convection problem for Rayleigh numbers three orders of magnitude larger than was possible with the Newton-Raphson update. The computational effort required by the semidirect methods is comparable to that required by the iterative methods; however, the memory requirements can be significantly greater.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

Flow around a NACA0018 airfoil with a cavity and its dynamical response to acoustic forcing

Trapping of vortices in a cavity has been explored in recent years as a drag reduction measure for thick airfoils. If, however, trapping fails, then oscillation of the cavity flow may couple with elastic vibration modes of the airfoil. To examine this scenario, the effect of small amplitude vertical motion on the oscillation of the shear layer above the cavity is studied by acoustic forcing simulating a vertical translation of a modified NACA0018 profile. At low Reynolds numbers based on the chord (O(10⁴)), natural instability modes of this shear layer are observed for Strouhal numbers based on the cavity width of order unity. Acoustic forcing sufficiently close to the natural instability frequency induces a strong non-linear response due to lock-in of the shear layer. At higher Reynolds numbers (above 10⁵) for Strouhal number 0.6 or lower, no natural instabilities of the shear layer and only a linear response to forcing were observed. The dynamical pressure difference across the airfoil is then dominated by added mass effects, as was confirmed by numerical simulations.     

Finite element analysis of the axially symmetric motion of an incompressible viscous fluid in a spherical annulus

This paper presents results obtained by employing a modified Galerkin finite element method to analyse the steady state flow of a fluid contained between two concentric, rotating spheres. The spheres are assumed to be rigid and the cavity region between the spheres is filled with an incompressible, viscous, Newtonian fluid. The inner sphere is constrained to rotate about a vertical axis with a prescribed angular velocity, while the outer sphere is fixed. Results for the circumferential function Ω, streamfunction ψ, vorticity function ζ and inner boundary torque T₁ are presented for Reynolds numbers Re ? 2000 and radius ratios 0.1 ? α ? 0.9. The method proved effective for obtaining results for a wide range of radius ratios (0.1 ? α ? 0.9) and Reynolds numbers (0 ? Re ? 2000). Previous investigators who employed the finite difference method experienced difficulties in obtaining results for cases with radius ratios α ? 0.2, except for small Reynolds numbers (Re ? 100). Results for Ω, Ψ, ζ and T₁ obtained in this study for radius ratios 0.8 ≤ α ≤ 0.9 verified the development of Taylor vortices reported by other investigators. The research indicates that the method may be useful for analysing other non-linear fluid flow problems.     

A LES study of the flow interference between tandem square cylinder pairs

The results of an investigation on the interference effects of the tandem square cylinders exposed to a uniform flow are presented in this paper. Time-dependent and three-dimensional flow simulations are carried out using large eddy simulation with a one-equation subgrid model. An incompressible three-dimensional finite volume code with a collocated grid arrangement is used for solving filtered Navier–Stokes equations. These equations are solved with an implicit fractional two-step method. Simulations are conducted with different Reynolds numbers between 10³ and 10⁵. The longitudinal spacing between the cylinders is selected 4D for the chosen Reynolds numbers, where D is the side of the cylinders. Also the effect of the spacing between cylinders, ranging from 1D to 12D, is studied for the selected Reynolds numbers. The instantaneous flow field is studied by analyzing the vortices, pressure, streamlines and Q-criterion to assist understanding of the various flow patterns, vortical structures and Kelvin–Helmholtz vortices in the separating shear layers. The hysteresis is observed in a certain range of the gap spacing, which this range depends on the selected Reynolds number. The global results are also computed and compared with available experimental results. The results indicate that there is a satisfactory agreement between the predictions and available experimental data considering the fine grid adopted.     

Aerodynamic Characteristics of Transonic and Supersonic Flow over Rectangular Cavities

An experimental and computational investigation has been performed to investigate flow characteristics and flow-field structures for three types of rectangular cavities. The data presented herein was obtained with cavity length to depth ratio of 6, 10 and 15 at 0° of attack, yawing and rolling angles of 0° over free-stream Mach numbers of 0.6, 0.8, 1.2 and 1.5 at Reynolds numbers of 1.23 ×10⁷, 1.55 ×10⁷, 2.01 ×10⁷ and 2.26 ×10⁷ per meter. The results indicate that the shear-layer expands over the cavity leading edge and impinges on the cavity floor for closed cavity flow, whereas it bridges the open cavity. The static pressure distributions are relatively uniform with the exception of a small adverse gradient occurring ahead of the rear face inside open cavity. Cavity length to depth ratio is a key geometrical parameter to define cavity flow types and influent pressure distributions inside cavities, and its decrease induces a decrease in pressure gradient. Increase in free-stream Mach numbers results in the trend that cavity flow types transform from closed to transitional cavity flow and from transitional to open cavity flow.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

Aspect ratio effects on three-dimensional incompressible flow in a two-sided non-facing lid-driven parallelepiped cavity

A numerical study of the three-dimensional fluid flow has been carried out to determine the effects of the transverse aspect ratio, Ay, on the flow structure in two-sided non-facing lid-driven cavities. The flow is complex, unstable and can undergo bifurcation. The numerical method is based on the finite volume method and multigrid acceleration. Computations have been investigated for several Reynolds numbers and various aspect ratio values. At a fixed Reynolds number, Re=500, the three-dimensional flow characteristics are analyzed considering four transverse aspect ratios, Ay=1,0.75,0.5 and 0.25. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. An analysis of the flow evolution shows that, at Ay=0.75, the flow bifurcates to a periodic regime at (Re=600) with a frequency f=0.093 less than the predicted value in the cubical cavity. A correlation is established when Ay=0.5 and gives the critical Reynolds number value. At Ay=0.25, the periodic regime occurs at high Re value beyond 3500, after which the flow becomes chaotic. It is shown that, when increasing Ay over the unit, the flow in the cavity exhibits a complex behavior. The kinetic energy transmission from the driven walls to the cavity center is reduced at low Ay values.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Numerical solutions of 2‐D steady incompressible driven cavity flow at high Reynolds numbers

Numerical calculations of the 2‐D steady incompressible driven cavity flow are presented. The Navier–Stokes equations in streamfunction and vorticity formulation are solved numerically using a fine uniform grid mesh of 601 × 601. The steady driven cavity flow solutions are computed for Re ? 21 000 with a maximum absolute residuals of the governing equations that were less than 10^?10. A new quaternary vortex at the bottom left corner and a new tertiary vortex at the top left corner of the cavity are observed in the flow field as the Reynolds number increases. Detailed results are presented and comparisons are made with benchmark solutions found in the literature. Copyright © 2005 John Wiley & Sons, Ltd.     

Time-dependent natural convection in a square cavity: Application of a new finite volume method

A new finite volume (FV) approach with adaptive upwind convection is used to predict the two-dimensional unsteady flow in a square cavity. The fluid is air and natural convection is induced by differentially heated vertical walls. The formulation is made in terms of the vorticity and the integral velocity (induction) law. Biquadratic interpolation formulae are used to approximate the temperature and vorticity fields over the finite volumes, to which the conservation laws are applied in integral form. Image vorticity is used to enforce the zero-penetration condition at the cavity walls. Unsteady predictions are carried sufficiently forward in time to reach a steady state. Results are presented for a Prandtl number (Pr) of 0-71 and Rayleigh numbers equal to 10³, 10⁴ and 10⁵. Both 11 × 11 and 21 × 21 meshes are used. The steady state predictions are compared with published results obtained using a finite difference (FD) scheme for the same values of Pr and Ra and the same meshes, as well as a numerical bench-mark solution. For the most part the FV predictions are closer to the bench-mark solution than are the FD predictions.     

Multiplicity of Steady Two-Dimensional Flows in Two-Sided Lid-Driven Cavities

The two-dimensional steady incompressible flow in rectangular cavities is calculated numerically by a finite volume method. The flow is driven by two opposing cavity side walls which move with constant velocities tangentially to themselves. Depending on the cavity aspect ratio and the two side-wall Reynolds numbers different flow states exist. Their range of existence and the bifurcations between different states are investigated by a continuation method accurately locating the bifurcation points. When both side walls move in opposite directions up to seven solutions are found to exist for the same set of parameters. Three of these are point-symmetric and four are asymmetric with respect to the center of the cavity, if the side-wall Reynolds numbers have the same magnitude. When the walls move in the same direction, up to five different flow states are found. In this case only a single mirror symmetric solution exists for equal Reynolds numbers. Received 9 February 2000 and accepted 9 October 2000     

Flow structure and heat transfer characteristics of an unconfined impinging air jet at high jet Reynolds numbers

The flow and heat transfer characteristics of an unconfined air jet that is impinged normally onto a heated flat plate have been experimentally investigated for high Reynolds numbers ranging from 30,000 to 70,000 and a nozzle-to-plate spacing range of 1–10. The mean and turbulence velocities by using hot-wire anemometry and impingement surface pressures with pressure transducer are measured. Surface temperature measurements are made by means of an infrared thermal imaging technique. The effects of Reynolds number and nozzle-to-plate spacing on the flow structure and heat transfer characteristics are described and compared with similar experiments. It was seen that the locations of the second peaks in Nusselt number distributions slightly vary with Reynolds number and nozzle-to-plate spacing. The peaks in distributions of Nusselt numbers and radial turbulence intensity are compatible for spacings up to 3. The stagnation Nusselt number was correlated for the jet Reynolds number and the nozzle-to-plate spacing as Nu _st∝Re ^0.69(H/D)^0.019.     

A least-squares finite element method for incompressible Navier-Stokes problems

A least-squares finite element method based on the velocity–pressure–vorticity formulation was proposed for solving steady incompressible Navier-Stokes problems. This method leads to a minimization problem rather than to the saddle point problem of the classic mixed method and can thus accommodate equal-order interpolations. The method has no parameter to tune. The associated algebraic system is symmetric and positive definite. In order to show the validity of the method for high-Reynolds-number problems, this paper provides numerical results for cavity flow at Reynolds number up to 10 000 and backward-facing step flow at Reynolds number up to 900.     

Direct numerical simulation of turbulent flow and heat transfer in a square duct with natural convection

In this paper, a direct numerical simulation of a fully developed turbulent flow and heat transfer are studied in a square duct with an imposed temperature difference between the vertical walls and the perfectly insulated horizontal walls. The natural convection is considered on the cross section in the duct. The numerical scheme employs a time-splitting method to integrate the three dimensional incompressible Navier-Stokes equation. The unsteady flow field was simulated at a Reynolds number of 400 based on the Mean friction velocity and the hydraulic diameter (Re _m = 6200), while the Prandtl number (Pr) is assumed 0.71. Four different Grashof numbers (Gr = 10⁴, 10⁵, 10⁶ and 10⁷) are considered. The results show that the secondary flow and turbulent characteristics are not affected obviously at lower Grashof number (Gr ≤ 10⁵) cases, while for the higher Grashof number cases, natural convection has an important effect, but the mean flow and mean temperature at the cross section are also affected strongly by Reynolds stresses. Compared with the laminar heat transfer at the same Grashof number, the intensity of the combined heat transfer is somewhat decreased.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.