A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

On separated shear layer of blunt circular cylinder

Separated shear layer of blunt circular cylinder has been experimentally investigated for the Reynolds numbers (based on the diameter) ranging from 2.8×10³ to 1.0×10⁵, with emphasis on evolution of separated shear layer, its structure and distribution of Reynolds shear stress and turbulence kinetic energy. The results demonstrate that laminar separated shear layer experiences 2–3 times vortex merging before it reattaches, and turbulence separated shear layer takes 5–6 times vortex merging. In addition, relationship between dimensionless initial frequencies of K-H instability and Reynolds numbers is identified, and reasons for the decay of turbulence kinetic energy and Reynolds shear stress in reattachment region are discussed. The project supported by the National Natural Science Foundation of China and the Key Laboratory for Hydrodynamics of NDCST.     

The numerical solution of the flow in a general bifurcating channel at moderately high Reynolds number using boundary-fitted co-ordinates,primitive variables and Newton iteration

A numerical code has been implemented for the numerical solution of the steady, incompressible Navier–Stokes equations using primitive variables in a bifurcating channel. A boundary-fitted, numerically generated grid is placed onto the domain of the channel which is transformed into either a rectilinear C- or T-shaped region. The differenced equations are solved using Newton's iteration which makes upwinding at high Reynolds number unnecessary. Practical implications of inverting the huge Jacobian matrix of Newton's method are discussed. The results have relative error of 2–3 × 10^?3 at Reynolds number 100, with T-geometry being marginally but significantly more accurate than C-geometry. Results have been obtained for Reynolds numbers up to 1000 for three bifurcations one of which models the carotid arterial bifurcation in the human head. For this latter bifurcation the wall shear stress is calculated in connection with the onset of atherosclerosis. Finally, the results of flows having different daughter tube end pressures are presented.     

Experimental investigation of mixed convection in a rectangular duct with a heated plate in the middle of cross section

The flow and heat transfer in an inclined and horizontal rectangular duct with a heated plate longitudinally mounted in the middle of cross section was experimentally investigated. The heated plate and rectangular duct were both made of highly conductive materials, and the heated plate was subjected to a uniform heat flux. The heat transfer processes through the test section were under various operating conditions: Pr ≈ 0.7, inclination angle ϕ = −60° to +60°, Reynolds number Re = 334–1,911, Grashof number Gr = 5.26 × 10²–5.78 × 10⁶. The experimental results showed that the average Nusselt number in the entrance region was 1.6–2 times as large as that in the fully developed region. The average Nusselt numbers and pressure drops increased with the Reynolds number. The average Nusselt numbers and pressure drops decreased with an increase in the inclination angle from −60° to +60° when the Reynolds number was less than 1,500. But when the Reynolds number increased to over about 1,800, the heat transfer coefficients and pressure drops were independent of inclination angles.     

Numerical investigation of turbulent flow,heat transfer and entropy generation in a helical coiled tube with larger curvature ratio

Three-dimensional turbulent forced convective heat transfer and flow characteristics, and the non-dimensional entropy generation number in a helical coiled tube subjected to uniform wall temperature are simulated using the k–ε standard turbulence model. A finite volume method is employed to solve the governing equations. The effects of Reynolds number, curvature ratio, and coil pitch on the average friction factor and Nusselt number are discussed. The results presented in this paper cover a Reynolds number range of 2 × 10⁴ to 6 × 10⁴, a pitch range of 0.1–0.2 and a curvature ratio range of 0.1–0.3. The results show that the coil pitch, curvature ratio and Reynolds number have different effects on the average friction factor and Nusselt number at different cross-sections. In addition, the flow and heat transfer characteristics in a helical coiled tube with a larger curvature ratio for turbulent flow are different from that of smaller curvature ratio for laminar and turbulent flow in certain ways. Some new features that are not obtained in previous researches are revealed. Moreover, the effects of Reynolds number, curvature ratio, and coil pitch on the non-dimensional entropy generation number of turbulent forced convection in a helical coiled tube are also discussed.     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 2. Validation and application at Reynolds numbers up to <Emphasis Type="Italic">Re</Emphasis> ≤ 10<Superscript>4</Superscript>

Hybrid unsteady-flow simulation combining particle tracking velocimetry (PTV) and direct numerical simulation (DNS) is introduced in the series of two papers. Particle velocities on a laser-light sheet acquired with time-resolved PTV in a water tunnel are supplied to two-dimensional DNS with time intervals corresponding to the frame rate of the PTV. Hybrid velocity fields then approach those representing the PTV data in the course of time, and the reconstructed velocity fields satisfy the governing equations with the resolution comparable to numerical simulation. In part 2, by extending the capabilities of the hybrid simulation to higher Reynolds numbers, we simulate flows past the NACA0012 airfoil over ranges of Reynolds numbers (Re ≤ 10⁴) and angles of attack (−5° ≤ α ≤ 20°) and validate the proposed technique by comparing with experimental results in terms of the lift and drag coefficients. We also compare the results with unsteady Reynolds-averaged Navier–Stokes (URANS) simulation in two-dimensions and show the advantages of the hybrid simulation against two-dimensional URANS.     

Numerical investigation of ethylene flame bubble instability induced by shock waves

In this paper, the ethylene/oxygen/nitrogen premixed flame instabilities induced by incident and reflected shock wave were investigated numerically. The effects of grid resolutions and chemical mechanisms on the flame bubble deformation process are evaluated. In the computational frame, the 2D multi-component Navier–Stokes equations with second-order flux-difference splitting scheme were used; the stiff chemical source term was integrated using an implicit ordinary differential equations (ODEs) solver. The two ethylene/oxygen/nitrogen chemical mechanisms, namely 3-step reduced mechanism and 35-step elementary skeletal mechanism, were used to examine the reliability of chemistry. On the other hand, the different grid sizes, Δx × Δy = 0.25 × 0.5mm and Δx × Δy = 0.15 × 0.2mm, were implemented to examine the accuracy of the grid resolution. The computational results were qualitatively validated with experimental results of Thomas et al. (Combust Theory Model 5:573–594, 2001). Two chemical mechanisms and two grid resolutions used in present study can qualitatively reproduce the ethylene spherical flame instability process generated by an incident shock wave of Mach number 1.7. For the case of interaction between the flame and reflected shock waves, the 35-steps mechanism qualitatively predicts the physical process and is somewhat independent on the grid resolutions, while the 3-steps mechanism fails to reproduce the instability of ethylene flame for the two selected grid resolutions. It is concluded that the detailed chemical mechanism, which includes the chain elementary reactions of fuel combustion, describes the flame instability induced by shock wave, in spite of the fact that the flame thickness (reaction zone) is represented by 1–2 grids only.      

Numerical investigation of corner angle and wing number effects on fluid flow characteristics of a turbulent stellar jet

In this research the fluid dynamics characteristics of a stellar turbulent jet flow is studied numerically and the results of three dimensional jet issued from a stellar nozzle are presented. A numerical method based on control volume approach with collocated grid arrangement is employed. The turbulent stresses are approximated using k–ε and k–ω models with four different inlet conditions. The velocity field is presented and the rate of decay at jet centerline is noted. Special attention is drawn on the influence of corner angle and number of wings on mixing in stellar cross section jets. Stellar jets with three; four and five wings and 15–65° corner angles are studied. Also the effect of Reynolds number (based on hydraulic diameter) as well as the inflow conditions on the evolution of the stellar jet is studied. The Numerical results show that the jet entrains more with corner angle 65° and five wings number. The jet is close to a converged state for high Reynolds numbers. Also the influence of the inflow conditions on the jet characteristics is so strong.     

Flow past a plate with a moving surface

On the basis of a numerical solution of the unsteady Navier-Stokes equations, the flow past a finite plate with an upstream-moving surface is investigated. For the Reynolds numbers Re =10²−10⁴, the flow past the plate is analyzed as a function of the relative plate surface velocity. On the basis of this analysis a limiting mathematical model of the flow as Re → ∞ is proposed.     

Experimental study on heat transfer of thermally developing and developed, turbulent, horizontal pipe flow of dilute air-solids suspensions

Heat transfer characteristics of a turbulent, dilute air-solids suspension flow in thermally developing/developed regions were experimentally studied, using a uniformly heated, horizontal 54.5 mm-ID pipe and 43-μm-diameter glass beads. The local heat transfer was measured at 27 locations from the inlet to 120-dia downstream of the heated section over a range of Reynolds numbers 3×10⁴−1.2×10⁵ and solids loading ratio 0–3, and the fully developed profiles of air velocity/temperature and particle mass flux were measured at a location 140-dia downstream of the heated section using specially designed probes, inserted into the suspension flow. The effects of the Reynolds number, solids loading ratio, and azimuthal/longitudinal locations on the heat transfer characteristics and their interactions are discussed through comparison of the present results with the data obtained by other investigators. Received on 14 October 1996     

Dynamic NMR microscopy measurement of the dynamics and flow partitioning of colloidal particles in a bifurcation

The flow and distribution of Newtonian, polymeric and colloid suspension fluids at low Reynolds numbers in bifurcations has importance in a wide range of disciplines, including microvascular physiology and microfluidic devices. A bifurcation consisting of circular capillaries laser etched into a hard polymer with inlet diameter 2.50 ± 0.01 mm, bifurcating to a small diameter outlet of 0.76 ± 0.01 mm and a large diameter outlet of 1.25 ± 0.01 mm is examined. Four distinct fluids (water, 0.25%wt xanthan gum, 8 and 22%vol hard-sphere colloidal suspensions) are flowed at flow rates from 10 to 30 ml/h corresponding to Reynolds numbers based on the entry flow from 0.001 to 8. PGSE NMR techniques are applied to obtain dynamic images of the fluids inside the bifurcation with spatial resolution of 59 × 59 μm/pixel in plane over a 200-μm-thick slice. Velocity in all three spatial directions is examined to determine the impact of secondary flows and characterize the transport in the bifurcation. The velocity data provide direct measurement of the volumetric distribution of the flow between the two channels as a function of flow rate. Water and the 8% colloidal suspension show a constant distribution with increasing flow rate, the xanthan gum shows an increase in fluid going into the larger outlet with higher flow rate, and the 22% colloidal suspension shows a decrease in fluid entering the larger channel with higher flow rate. For the colloidal particle flow, the distribution of colloid particles down the capillary is determined by examining the spectrally resolved propagator for the oil inside the core–shell particles in a direction perpendicular to the axial flow. Using dynamic magnetic resonance microscopy, the potential for using magnetic resonance for “particle counting” in a microscale bifurcation is thus demonstrated.     

A numerical study of the steady forced convection heat transfer from an unconfined circular cylinder

Forced convection heat transfer from an unconfined circular cylinder in the steady cross-flow regime has been studied using a finite volume method (FVM) implemented on a Cartesian grid system in the range as 10 ≤ Re ≤ 45 and 0.7 ≤ Pr ≤ 400. The numerical results are used to develop simple correlations for Nusselt number as a function of the pertinent dimensionless variables. In addition to average Nusselt number, the effects of Re, Pr and thermal boundary conditions on the temperature field near the cylinder and on the local Nusselt number distributions have also been presented to provide further physical insights into the nature of the flow. The rate of heat transfer increases with an increase in the Reynolds and/or Prandtl numbers. The uniform heat flux condition always shows higher value of heat transfer coefficient than the constant wall temperature at the surface of the cylinder for the same Reynolds and Prandtl numbers. The maximum difference between the two values is around 15–20%.     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 1. Algorithm of hybrid simulation and hybrid velocity field at <Emphasis Type="Italic">Re</Emphasis> ≈ 10<Superscript>3</Superscript>

We develop a hybrid unsteady-flow simulation technique combining direct numerical simulation (DNS) and particle tracking velocimetry (PTV) and demonstrate its capabilities by investigating flows past an airfoil. We rectify instantaneous PTV velocity fields in a least-squares sense so that they satisfy the equation of continuity, and feed them to the DNS by equating the computational time step with the frame rate of the time-resolved PTV system. As a result, we can reconstruct unsteady velocity fields that satisfy the governing equations based on experimental data, with the resolution comparable to numerical simulation. In addition, unsteady pressure distribution can be solved simultaneously. In this study, particle velocities are acquired on a laser-light sheet in a water tunnel, and unsteady flow fields are reconstructed with the hybrid algorithm solving the incompressible Navier–Stokes equations in two dimensions. By performing the hybrid simulation, we investigate nominally two-dimensional flows past the NACA0012 airfoil at low Reynolds numbers. In part 1, we introduce the algorithm of the proposed technique and discuss the characteristics of hybrid velocity fields. In particular, we focus on a vortex shedding phenomenon under a deep stall condition (α = 15°) at Reynolds numbers of Re = 1000 and 1300, and compare the hybrid velocity fields with those computed with two-dimensional DNS. In part 2, the extension to higher Reynolds numbers is considered. The accuracy of the hybrid simulation is evaluated by comparing with independent experimental results at various angles of attack and Reynolds numbers up to Re = 10⁴. The capabilities of the hybrid simulation are also compared with two-dimensional unsteady Reynolds-Averaged Navier–Stokes (URANS) solutions in part 2. In the first part of these twin papers, we demonstrate that the hybrid velocity field approaches the PTV velocity field over time. We find that intensive alternate vortex shedding past the airfoil, which is predicted by the two-dimensional DNS, is substantially suppressed in the hybrid simulation and the resultant flow field is similar to the PTV velocity field, which is projection of the three-dimensional velocity field on the streamwise plane. We attempt to identify the motion that originates three-dimensional flow patterns by highlighting the deviation of the PTV velocity field from the two-dimensional governing equations at each snapshot. The results indicate that the intensive spots of the deviation appear in the regions in which three-dimensional instabilities are induced in the shear layer separated from the pressure side.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

Numerical simulation of viscous flow interaction with an elastic membrane

A numerical fluid–structure interaction model is developed for the analysis of viscous flow over elastic membrane structures. The Navier–Stokes equations are discretized on a moving body‐fitted unstructured triangular grid using the finite volume method, taking into account grid non‐orthogonality, and implementing the SIMPLE algorithm for pressure solution, power law implicit differencing and Rhie–Chow explicit mass flux interpolations. The membrane is discretized as a set of links that coincide with a subset of the fluid mesh edges. A new model is introduced to distribute local and global elastic effects to aid stability of the structure model and damping effects are also included. A pseudo‐structural approach using a balance of mesh edge spring tensions and cell internal pressures controls the motion of fluid mesh nodes based on the displacements of the membrane. Following initial validation, the model is applied to the case of a two‐dimensional membrane pinned at both ends at an angle of attack of 4° to the oncoming flow, at a Reynolds number based on the chord length of 4 × 10³. A series of tests on membranes of different elastic stiffness investigates their unsteady movements over time. The membranes of higher elastic stiffness adopt a stable equilibrium shape, while the membrane of lowest elastic stiffness demonstrates unstable interactions between its inflated shape and the resulting unsteady wake. These unstable effects are shown to be significantly magnified by the flexible nature of the membrane compared with a rigid surface of the same average shape. Copyright © 2007 John Wiley & Sons, Ltd.     

Non-symmetric bi-stable flow around the Ahmed body

The flow around the Ahmed body at varying Reynolds numbers under yawing conditions is investigated experimentally. The body geometry belongs to a regime subject to spanwise flow instability identified in symmetric flow by Cadot and co-workers (Grandemange et al., 2013b). Our experiments cover the two slant angles 25° and 35° and Reynolds numbers up to 2.784 × 10⁶. Special emphasis lies on the aerodynamics under side wind influence. For the 35° slant angle, forces and moments change significantly with the yawing angle in the range 10° ≤ |β| ≤ 15°. The lift and the pitching moment exhibit strong fluctuations due to bi-stable flow around a critical angle β of ±12.5°, where the pitching moment changes sign. Time series of the forces and moments are studied and explained by PIV measurements in the flow field near the rear of the body.     

Simulation of Richtmyer–Meshkov instability by sixth-order filter methods

Simulation of a 2-D Rightmyer–Meshkov instability (RMI), including inviscid, viscous and magnetic field effects was conducted comparing recently developed sixth-order filter schemes with various standard shock-capturing methods. The suppression of the inviscid gas dynamics RMI in the presence of a magnetic field was investigated by Samtaney and Wheatley et al. Numerical results illustrated here exhibit behavior similar to the work of Samtaney. Due to the different amounts and different types of numerical dissipation contained in each scheme, the structures and the growth of eddies for the chaotic-like inviscid gas dynamics RMI case are highly grid size and scheme dependent, even with many levels of refinement. The failure of grid refinement for all studied numerical methods extends to the viscous gas dynamics case for high Reynolds number. For lower Reynolds number, grid convergence has been achieved by all studied methods. To achieve similar resolution, standard shock-capturing methods require more grid points than filter schemes and yet the CPU times using the same grid for all studied methods are comparable. This paper is based on work that was presented at the 17th International Shock Interaction Symposium (ISIS17), Rome, Italy, 4–8 September 2006.     

Finite-Amplitude Equilibrium Solutions for Plane Poiseuille–Couette Flow

Two-dimensional nonlinear equilibrium solutions for the plane Poiseuille–Couette flow are computed by directly solving the full Navier–Stokes equations as a nonlinear eigenvalue problem. The equations are solved using the two-point fourth-order compact scheme and the Newton–Raphson iteration technique. The linear eigenvalue computations show that the combined Poiseuille–Couette flow is stable at all Reynolds numbers when the Couette velocity component σ₂ exceeds 0.34552. Starting with the neutral solution for the plane Poiseuille flow, the nonlinear neutral surfaces for the combined Poiseuille–Couette flow were mapped out by gradually increasing the velocity component σ₂. It is found that, for small σ₂, the neutral surfaces stay in the same family as that for the plane Poiseuille flow, and the nonlinear critical Reynolds number gradually increases with increasing σ₂. When the Couette velocity component is increased further, the neutral curve deviates from that for the Poiseuille flow with an appearance of a new loop at low wave numbers and at very low energy. By gradually increasing the σ₂ values at a constant Reynolds number, the nonlinear critical Reynolds numbers were determined as a function of σ₂. The results show that the nonlinear neutral curve is similar in shape to a linear case. The critical Reynolds number increases slowly up to σ₂∼ 0.2 and remains constant until σ₂∼ 0.58. Beyond σ₂ > 0.59, the critical Reynolds number increases sharply. From the computed results it is concluded that two-dimensional nonlinear equilibrium solutions do not exist beyond a critical σ₂ value of about 0.59. Received: 26 November 1996 and accepted 12 May 1997     

Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.