NUMERICAL SIMULATION OF WING-BODY JUNCTION TURBULENCE FLOW

Wing-body junction turbulence flow is simulated by using RANS equation and boundary fitted coordinate technique.Three order differential scheme is used in the computation of convection term and two layers turbulence model are employed in the calculation.     

High velocity flow through fractured and porous media: the role of flow non-periodicity

In the present paper we examine the evolution of the macroscopic flow law in a crenellated channel, representing an element of fractured or porous medium and in function of the Reynolds number Re. A numerical analysis based on the Navier–Stokes equations is applied. We focus on the influence of the flow periodicity or non-periodicity upon the macroscopic law. The physical explanation of the non-linear deviation from Darcy's law is still an issue, as the Ergun–Forchheimer law admitted for high Reynolds numbers comes up against some theoretical problems. In the periodic case, three non-linear flow regimes were revealed: a cubic flow with respect to velocity at low Re, an intermediate non-quadratic law, and a self-similar mode independent of Re at very high Re. The Forchheimer law is not confirmed. The case of a non-periodic flow clearly highlights the link between the flow non-periodicity and the quadratic law. The quadratic deviation becomes all the more important as the non-periodicity degree is high.     

Effect of upstream velocity profile and integral flow straighteners on turbine flowmeters

One of the most important factors determining the shape of the calibration curve for a given turbine meter is the change in the upstream velocity distribution with flowrate. A theoretical model is evolved which can be used to predict the effects of velocity profile, viscosity and swirl on the calibration curve. It has also been used to explain the calibration curve of a commercial meter having a geometry very different from that for which the theory was developed. The effect of different types of integral flow straighteners on turbine meters is also investigated and found to depend on both the number of vanes and their length. A correlation is suggested for radial-vaned flow-straighteners     

Lattice-Boltzmann simulations of flow through Fontainebleau sandstone

We report preliminary results from simulations of single-phase and two-phase flow through three-dimensional tomographic reconstructions of Fontainebleau sandstone. The simulations are performed with the lattice-Boltzmann method, a variant of lattice-gas cellular-automation models of fluid mechanics. Simulations of single-phase flow on a sample of linear size 0.2 cm yield a calculated permeability in the range 1.0–1.5 darcys, depending on direction, which compares qualitatively well with a laboratory measurement of 1.3 darcys on a sample approximately an order of magnitude larger. The sensitivity of permeability calculations to sample size, grid resolution, and choice of model parameters is quantified empirically. We also present a qualitative study of immiscible two-phase flow in a sample of linear size 0.05 cm; simulations of both drainage and imbibition are presented.     

Electrical fluctuations and entry flow instabilities during capillary flow of aqueous poly (ethylene oxide) solutions

Capillary flow of poly (ethylene oxide) solutions generates voltage fluctuations (noise) between electrodes placed on both sides of the capillary. The noise has a 1/f type spectral distribution, the value of increasing with (shear rate) to a limiting level. Within certain ranges, two sets of harmonic peaks appear in the spectra. It is demonstrated that these peaks are related to the frequency components of the instabilities in the entry flow region. The lower frequency set of harmonics corresponds to axial oscillations (pulsations) of the flow, while the higher frequency peaks are associated with the transverse oscillation of the stream lines in the vicinity of the entry. The corresponding frequencies were measured by visual counting and by spectral analysis of laser light transmitted through the entry region during injection of a coloured solution.The noise measurements were carried out using both platinum and reversible Ag/AgCl electrodes. In the latter case the streaming potential was also measured; its variation due to pressure fluctuations in the instable flow region appears to provide a plausible explanation of flow-induced noise phenomena as observed in elastic solutions.Dedicated to Prof. Dr. J. Schurz on the occasion of his sixtieth birthday.     

Downstream pressure and elastic wall reflection of droplet flow in a T-junction microchannel

This paper discusses pressure variation on a wall during the process of liquid flow and droplet formation in a T-junction microchannel. Relevant pressure in the chan-nel, deformation of the elastic wall, and responses of the droplet generation are analyzed using a numerical method. The pressure difference between the continuous and dis-persed phases can indicate the droplet-generation period. The pressure along the channel of the droplet flow is affected by the position of droplets, droplet-generation period, and droplet escape from the outlet. The varying pressures along the channel cause a nonuniform deformation of the wall when they are elastic. The deformation is a vibration and has the same period as the droplet generation arising from the process of droplet formation.     

On the accurate prediction of the wall-normal velocity in compressible boundary-layer flow

We consider a problem which arises in the numerical solution of the compressible two-dimensional or axisymmetric boundary-layer equations. Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x, y) plane to a computational (ξ, η) plane in which the evolution of the flow is ‘slow’ in the time-like ξ direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently non-linear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which satisfies the continuity equation nearly to machine precision. As demonstration of the utility of the method, the boundary layers of three prototypical high-speed flows are investigated and compared: the flat plate, the hollow cylinder, and the cone. An important implication for classical linear stability theory is also briefly discussed.     

Effect of the magnetic field on the pulsatile flow through a rigid tube

I.IntroductionThefunctionofanaPbliedmagneticfield\onhumanbodieshasbeenknownandservedasakindofthemedicaltreatmentforalongtime.ThemagneticstoneusedasamedicineisrdeordedinShennongMaterlbMedica,'which.isthefirstChinesemedicinebookandwaswrittenin200A.D..ItwasalsorecordedtyancientGreecethatholdingmagneticstQnesinhandsandfeetcouldrelievethepainandspasml'l.However,beinglimitedbythelevelQfproductionandmedicaltreatmentsatthattime,themagnetictreatmentwasmainlyusedforlocallydephlogisticating,relievin…     

Concentration and velocity measurements in the flow of droplet suspensions through a tube

Two optical methods, light absorption and LDA, are applied to measure the concentration and velocity profiles of droplet suspensions flowing through a tube. The droplet concentration is non-uniform and has two maxima, one near the tube wall and one on the tube axis. The measured velocity profiles are blunted, but a central plug-flow region is not observed. The concentration of droplets on the tube axis and the degree of velocity profile blunting depend on relative viscosity. These results can be qualitatively compared with the theory of Chan and Leal.List of symbols a particle radius,m - a/R, non-dimensional particle radius - c volume concentration of droplets in suspension, m³/m³ - c ₅ stream-average volume concentration of droplets in suspension, - D 2 R, tube diameter, m - L optical path length, m - L _ij path length of laser beam through thej-th concentric layer when the beam crosses the tube diameter at the point on the inner circumference of thei-th layer, m - N exponent in Eqs. (3) and (4) - Q volumetric flowrate of suspension, - R tube radius, m - Re _{S S} D/µ, flow Reynolds number - r radial position (r = 0 on a tube axis), m - r r/R, non-dimensional radial position - v velocity of suspension, m/s - v v/v _S , non-dimensional velocity - v ₀ centre-line velocity of suspension (r = 0), m/s - v _S Q/ R ², stream-average velocity of suspension, m/s - x streamwise position (x = 0 at tube inlet), m - x x/D, non-dimensional streamwise position - _c density of continuous phase, kg/m³ - _d density of dispersed phase, kg/m³ - _s stream-average density of suspension, kg/m³, equals density when homogenized - _d - _c, phase density difference, kg/m³ - µ_c viscosity of continuous phase, Pa · s - µ_d viscosity of dispersed (droplet) phase, Pa · s - µ_d/_c, viscosity ratio - interfacial tension, N/m This work was financially supported by the National Science Foundation (USA) through an agreement no. J-F7F019P, M. Sklodowska-Curie fund     

On the relationship between drag and vertical velocity fluctuations in flow over riblets and liquid infused surfaces

Direct numerical simulations (DNS) of flow over triangular and rectangular riblets in a wide range of size and Reynolds number have been carried out. The flow within the grooves is directly resolved by exploiting the immersed-boundary method. It is found that the drag reduction property is primarily associated with the capability of inhibiting vertical velocity fluctuations at the plane of the crests, as in liquid-infused surfaces (LIS) devices. This is mimicked in DNS through artificial suppression of the vertical velocity component, which yields large drag decrease, proportionate to the riblets size. A parametrization of the drag reduction effect in terms of the vertical velocity variance is found to be quite successful in accounting for variation of the controlling parameters. A Moody-like friction diagram is thus introduced which incorporates the effect of slip velocity and a single, geometry-dependent parameter. Reduced drag-reduction efficiency of LIS-like riblets is found as compared to cases with artificially imposed slip velocity. Last, we find that simple wall models of riblets and LIS-like devices are unlikely to provide accurate prediction of the flow phenomenon, and direct resolution of flow within the grooves in necessary.     

Effect of collision and velocity model of lattice Boltzmann model on three-dimensional turbulent flow simulation

Various collision and velocity models of the lattice Boltzmann model (LBM) were compared to determine their effects on the efficiency of a three-dimensional homogeneous isotropic decaying turbulent flow simulation. We determined that a decrease in the number of velocities, in particular, 13-velocities, which can be used in the quasi-equilibrium lattice Boltzmann and in the multiple-relaxation time models (MRT), could considerably decrease the computational effort. However, decreasing the number of velocities deteriorates the stability and the accuracy of the results. By comparing the collision models, we also determined that the stability of the entropic lattice Boltzmann model (ELBM), and 19- and 27- velocity MRT is much higher than in other models. However, the numerical viscosity introduced by the ELBM underestimates the enstrophy, and the computational effort increases because of the calculation overhead required to solve the additional equations if special care is not given to the calculation.     

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.     

The split of vertical annular flow at a large diameter T junction

Measurements and observations have been made of the split of gas—liquid flow at a T junction where the main pipe was vertical and the side arm horizontal. All three pipes connecting to the junction were of 0.125 m dia. The gas and liquid flow rates were chosen to ensure annular flow in the inlet pipe. The resulting data have been compared with existing models.     

Efficient implementation of the proper outlet flow conditions in blood flow simulations through asymmetric arterial bifurcations

We present an efficient implementation of the proper (in vivo) outlet boundary conditions in detailed, three‐dimensional (3D) and time‐periodic simulations of blood flow through arteries. This is achieved through the intermediate use of an approximate ‘simulant’ model of the outlet pressure/flow relationship corresponding to the full 3D and time‐dependent numerical simulation. This model allows us to efficiently couple the 3D outlet pressure/flow conditions to the equivalent relations due to the downstream arterial network, as obtained from a one‐dimensional approximate model in the form of Fourier frequency impedance coefficients. An adjustable time‐periodic function correction term in the simulant model requires input from the full 3D model that has to run iteratively until convergence. The advantage of the proposed numerical scheme is that it decouples the upstream detailed simulation from the downstream approximate network model offering exceptional versatility. This approach is demonstrated here in a series of detailed 3D simulations of blood flow, performed using the commercial software FLUENT?, through an asymmetric arterial bifurcation. Two cases are considered: first a healthy system patterned after the left main coronary arterial bifurcation, and second a diseased case where an occlusion has developed in one of the daughter vessels, resulting in strengthening the asymmetry of the bifurcation. Rapid convergence of the iterative process was achieved in both cases. Subtle changes occur in the shear patterns of the daughter vessels, whereas the flow distribution is quite different. In the presence of a stenosis additional regions of low shear develop due to inertial effects. Copyright © 2010 John Wiley & Sons, Ltd.     

The possibility of measuring flow velocity using a low-intensity periodic pulse discharge

A direct method of measuring the flow velocity in a supersonic wind tunnel by creating a low-intensity periodic pulse discharge is proposed.     

Influence of the velocity gradient on the stagnation point heating in hypersonic flow

In a number of experimental and numerical publications a deviation has been found between the measured or computed stagnation point heat flux and that given by the theory of Fay and Riddell. Since the formula of Fay and Riddell is used in many applications to yield a reference heat flux for experiments performed in wind tunnels, for flight testing and numerical simulations, it is important that this reference heat flux is as accurate as possible. There are some shortcomings in experiments and numerical simulations which are responsible in some part for the deviations observed. But, as will be shown in the present paper, there is also a shortcoming on the theoretical side which plays a major role in the deviation between the theoretical and experimental/numerical stagnation point heat fluxes. This is caused by the method used so far to determine the tangential velocity gradient at the stagnation point. This value is important for the stagnation point heat flux, which so far has been determined by a simple Newtonian flow model. In the present paper a new expression for the tangential velocity gradient is derived, which is based on a more realistic flow model. An integral method is used to solve the conservation equations and, for the stagnation point, yields an explicit solution of the tangential velocity gradient. The solution achieved is also valid for high temperature flows with real gas effects. A comparison of numerical and experimental results shows good agreement with the stagnation point heat flux according to the theory of Fay and Riddell, if the tangential velocity gradient is determined by the new theory presented in this paper.This article was processed by the author using theLAT_EX style filepljour2 from Springer-Verlag.     

Some features of a turbulent wing–body junction vortical flow

Laser-Doppler velocimeter measurements of a wing/body junction flow field made within a plane to the side of the wing/wall junction and perpendicular both to a 3:2 elliptical nose—NACA 0020 tail wing, and a flat wall are presented. Reynolds number of the approach boundary layer was, Re_θ = 5940, and free-stream air velocity was, U_ref = 27.5 m/s. A large vortical structure residing in the outer region redirects the low-turbulence free-stream flow to the vicinity of the wing/wall junction, resulting in thin boundary layers with velocity magnitudes higher than free-stream flow. Lateral pressure gradients result in a three-dimensional separation on the uplifting side of the vortex. Additionally, a high vorticity vortical structure with opposite sense to the outer-layer vortex forms beneath the outer-layer vortex. Normal and shear stresses increase to attain values an order of magnitude larger compared to values measured in a three-dimensional boundary layer just outside the junction vortex. Bimodal histograms of the w fluctuating velocity occur under the outer-layer vortex near the wall due to the time-dependent nature of the horseshoe vortex. In such a flow the shear-stress angle (SSA) highly lags the flow-gradient angle (FGA), and the turbulence diffusion is highly altered due to presence of vortical structures.