Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Micro-shock waves generated inside a fluid jet impinging on plane wall

An Nd:glass laser pulse (18 ns, 1.38 J) is focused in a tiny area of about 100-mum diam under ambient conditions to produce micro-shock waves. The laser is focused above a planar surface with a typical standoff distance of about 4 mm, The laser energy is focused inside a supersonic circular jet of carbon dioxide gas produced by a nozzle with internal diameter of 2.9 mm and external diameter of 8 mm, Nominal value of the Mach number of the jet is around 2 with the corresponding pressure ratio of 7.5 (stagnation pressure/static pressure at the exit of the nozzle), The interaction process of the micro-shock wave generated inside the supersonic jet with the plane wall is investigated using double-pulse holographic interferometry. A strong surface vortex field with subsequent generation of a side jet propagating outward along the plane wail is observed. The interaction of the micro-shock wave with the cellular structure of the supersonic jet does not seem to influence the near surface features of the flowfield. The development of the coherent structures near the nozzle exit due to the upstream propagation of pressure waves seems to be affected by the outward propagating micro-shock wave. Mach reflection is observed when the micro-shock wave interacts with the plane wall at a standoff distance of 4 mm, The Mach stem is slightly deflected, indicating strong boundary-layer and viscous effects near the wall. The interaction process is also simulated numerically using an axisymmetric transient laminar Navier-Stokes solver. Qualitative agreement between experimental and numerical results is good.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Modeling Study on the Characteristics of Laminar and Turbulent Argon Plasma Jets Impinging Normally upon a Flat Plate in Ambient Air

Modeling study is performed to compare the flow and heat transfer characteristics of laminar and turbulent argon thermal-plasma jets impinging normally upon a flat plate in ambient air. The combined-diffusion-coefficient method and the turbulence-enhanced combined-diffusion-coefficient method are employed to treat the diffusion of argon in the argon-air mixture for the laminar and the turbulent cases, respectively. Modeling results presented include the flow, temperature and argon concentration fields, the air mass flow-rates entrained into the impinging plasma jets, and the distributions of the heat flux density on the plate surface. It is found that the formation of a radial wall jet on the plate surface appreciably enhances the mass flow rate of the ambient air entrained into the laminar or turbulent plasma impinging-jet. When the plate standoff distance is comparatively small, there exists a significant difference between the laminar and turbulent plasma impinging-jets in their flow fields due to the occurrence of a large closed recirculation vortex in the turbulent plasma impinging-jet, and no appreciable difference is found between the two types of jets in their maximum values and distributions of the heat flux density at the plate surface. At larger plate standoff distances, the effect of the plate on the jet flow fields only appears in the region near the plate, and the axial decaying-rates of the plasma temperature, axial velocity and argon mass fraction along the axis of the laminar plasma impinging-jet become appreciably less than their turbulent counterparts.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

超声速大攻角旋成体迎风激波数值模拟

旋成弹体在大攻角超声速流动中其迎风面上有时会出现一种迎风激波.采用非定常雷诺平均的Navier-Stokes方程以及两方程kω-sst和代数B-L湍流模型成功地模拟了该现象,验证了该迎风激波是由于流场主旋涡导致超声速流动偏移造成的涡诱导激波的一部分的结论.

Evaluation of the Cumulative Formation of High-Speed Liquid Jets

This paper describes the generation of pulsed, high-speed liquid jets using the cumulation method. This work mainly includes (1) the design of the nozzle assembly, (2) the measurement of the jet velocity and (3) flow visualization of the injection sequences. The cumulation method can be briefly described as the liquid being accelerated first by the impact of a moving projectile and then further after it enters a converging section. The experimental results show that the cumulation method is useful in obtaining a liquid jet with high velocity. The flow visualization shows the roles of the Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the breakup of the liquid depend on the jet diameter and the downstream distance. When the liquid jet front is far downstream from the nozzle exit, the jet is decelerated by air drag. Meanwhile, large coherent vortex structures are formed surrounding the jet. The liquid will break up totally by the action of these vortices. Experimental results showing the effect of the liquid volume on the jet velocity are also included in this paper. Finally, a method for measuring the jet velocity by cutting two carbon rods is examined.

Three-Dimensional Propagation of the Transmitted Shock Wave in a Square Cross-Sectional Chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.

Numerical Simulation of Particle Separation in An Oil-Sand Separator

The gathering systems of crude oil are greatly endangered by the fine sand and soil in oil. Up to now , how to separate sand from the viscid oil is still a technical problem for oil production home or abroad. Recently , Institute of Mechanics in Chinese Academy of Sciences has developed a new type of oil-sand separator , which has been applied successfully in oil field in situ. In this paper, the numerical method of vortex-stream function is used to predict the liquid-solid separating course and the efficiency for this oil-sand separator. Results show that the viscosity and particle diameter have much influence on the particle motion. The calculating separating efficiency is compared with that of experiment and indicates that this method can be used to model the complex two-phase flow in the separator.

Steady Current Induced Seabed Scour around a Vibrating Pipeline

Most of the existing researches either focus on vortex-induced vibrations (VIVs) of a pipeline near a rigid boundary, or on seabed scour around a fixed pipeline. In this study, the coupling effects between pipeline vibration and sand scour are investigated experimentally. Experimental results indicate that there often exist two phases in the process of sand scouring around the pipeline with an initial embedment, i.e. Phase I: scour beneath pipe without VIV, and Phase II: scour with VIV of pipe. During Phase II, the amplitude of pipe vibration gets larger and its frequency gets smaller while the sand beneath the pipe is being scoured, and finally the pipe vibration and sand scour get into an equilibrium state. This indicates that sand scouring has an influence upon not only the amplitude of pipe vibration but also its frequency. Moreover, the equilibrium scour depth decreases with increasing initial gap-to-diameter ratio for both the fixed pipes and vibrating pipes. For a given value of initial gapto- diameter ratio (e0/D), the vibrating pipe may induce a deeper scour hole than the fixed pipe in the examined range of initial gap-to-diameter ratios (−0.25 < e0/D < 0.75).

DNS of Compressible Turbulent Boundary Layer Over a Blunt Wedge

Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.

Statistics of dissipation and enstrophy induced by localized vortices

The ratios of enstrophy and dissipation moments induced by localized vorticity are inferred to be finite. It follows that the scaling exponents for locally averaged dissipation and enstrophy are equal. However, enstrophy and dissipation exponents measured over finite ranges of scales may be different. The cylindrical vortex profile that yields maximal moment ratios is determined. The moment ratios for cylindrical vortices are used to interpret differences in scale dependence of enstrophy and dissipation previously found in numerical simulations.

Control of Mean and Fluctuating Forces on a Circular Cylinder at High Reynolds Numbers

A narrow strip is used to control mean and fluctuating forces on a circular cylinder at Reynolds numbers from 2.0 x 10(4) to 1.0 x 10(5). The axes of the strip and cylinder are parallel. The control parameters are strip width ratio and strip position characterized by angle of attack and distance from the cylinder. Wind tunnel tests show that the vortex shedding from both sides of the cylinder can be suppressed, and mean drag and fluctuating lift on the cylinder can be reduced if the strip is installed in an effective zone downstream of the cylinder. A phenomenon of mono-side vortex shedding is found. The strip-induced local changes of velocity profiles in the near wake of the cylinder are measured, and the relation between base suction and peak value in the power spectrum of fluctuating lift is studied. The control mechanism is then discussed from different points of view.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.