Numerical study of turbulent fluid flow and heat transfer in lateral perforated extended surfaces

Numerical study has been performed in this study to investigate the turbulent convection heat transfer on a rectangular plate mounted over a flat surface. Thermal and fluid dynamic performances of extended surfaces having various types of lateral perforations with square, circular, triangular and hexagonal cross sections are investigated. RANS (Reynolds averaged Navier–Stokes) based modified k–ω turbulence model is used to calculate the fluid flow and heat transfer parameters. Numerical results are compared with the results of previously published experimental data and obtained results are in reasonable agreement. Flow and heat transfer parameters are presented for Reynolds numbers from 2000 to 5000 based on the fin thickness.

Simultaneous measurement of aerodynamic and heat transfer data for large angle blunt cones in hypersonic shock tunnel

Aerodynamic forces and fore-body convective surface heat transfer rates over a 60 degrees apex-angle blunt cone have been simultaneously measured at a nominal Mach number of 5.75 in the hypersonic shock tunnel HST2. An aluminum model incorporating a three-component accelerometer-based balance system for measuring the aerodynamic forces and an array of platinum thin-film gauges deposited on thermally insulating backing material flush mounted on the model surface is used for convective surface heat transfer measurement in the investigations. The measured value of the drag coefficient varies by about +/-6% from the theoretically estimated value based on the modified Newtonian theory, while the axi-symmetric Navier-Stokes computations overpredict the drag coefficient by about 9%. The normalized values of measured heat transfer rates at 0 degrees angle of attack are about 11% higher than the theoretically estimated values. The aerodynamic and the heat transfer data presented here are very valuable for the validation of CFD codes used for the numerical computation of How fields around hypersonic vehicles.

Effect of uncertainty on helicopter performance predictions

The effect of uncertainties on performance predictions of a helicopter is studied in this article. The aeroelastic parameters such as the air density, blade profile drag coefficient, main rotor angular velocity, main rotor radius, and blade chord are considered as uncertain variables. The propagation of these uncertainties in the performance parameters such as thrust coefficient, figure of merit, induced velocity, and power required are studied using Monte Carlo simulation and the first-order reliability method. The Rankine-Froude momentum theory is used for performance prediction in hover, axial climb, and forward flight. The propagation of uncertainty causes large deviations from the baseline deterministic predictions, which undoubtedly affect both the achievable performance and the safety of the helicopter. The numerical results in this article provide useful bounds on helicopter power requirements.

Shock tunnel study of spiked aerodynamic bodies flying at hypersonic Mach numbers

A three-component accelerometer balance system is used to study the drag reduction effect of an aerodisc on large angle blunt cones flying at hypersonic Mach numbers. Measurements in a hypersonic shock tunnel at a freestream Mach number of 5.75 indicate more than 50% reduction in the drag coefficient for a 120degrees apex angle blunt cone with a forward facing aerospike having a flat faced aerodisc at moderate angles of attack. Enhancement of drag has been observed for higher angles of attack due to the impingement of the flow separation shock on the windward side of the cone. The flowfields around the large angle blunt cone with aerospike assembly flying at hypersonic Mach numbers are also simulated numerically using a commercial CFD code. The pressure and density levels on the model surface, which is under the aerodynamic shadow of the flat disc tipped spike, are found very low and a drag reduction of 64.34% has been deduced numerically.

Uncertainty Quantification in Helicopter Performance Using Monte Carlo Simulations

The effect of structural and aerodynamic uncertainties on the performance predictions of a helicopter is investigated. An aerodynamic model based on blade element and momentum theory is used to predict the helicopter performance. The aeroelastic parameters, such as blade chord, rotor radius, two-dimensional lift-curve slope, blade profile drag coefficient, rotor angular velocity, blade pitch angle, and blade twist rate per radius of the rotor, are considered as random variables. The propagation of these uncertainties to the performance parameters, such as thrust coefficient and power coefficient, are studied using Monte Carlo Simulations. The simulations are performed with 100,000 samples of structural and aerodynamic uncertain variables with a coefficient of variation ranging from 1 to 5%. The scatter in power predictions in hover, axial climb, and forward flight for the untwisted and linearly twisted blades is studied. It is found that about 20-25% excess power can be required by the helicopter relative to the determination predictions due to uncertainties.

Effect of Channel Confinement on Mixed Convective Flow Past an Equilateral Triangular Cylinder

The present work investigates the mixed convective flow and heat transfer characteristics past a triangular cylinder placed symmetrically in a vertical channel. At a representative Reynolds number, Re = 100, simulations are carried out for the blockage ratios beta = 1/3; 1/4; and 1/6. Effect of aiding and opposing buoyancy is brought about by varying the Richardson number in the range -1.0 <= Ri <= 1.0. At a blockage ratio of 1/3, suppression of vortex shedding is found at Ri = 1, whereas von Karman vortex street is seen both at beta = 1/4 and 1/6, respectively. This is the first time that such behavior of blockage ratio past a triangular cylinder in the present flow configuration is reported. Drag coefficient increases progressively with increasing Ri and a slightly higher value is noticed at beta = 1/3. For all b, heat transfer increases with increasing Ri. Flattening of Nu(avg)-Ri curve beyond Ri > 0: 75 is observed at beta = 1/3.

The Effect of Conical Flowfields on the Performance of Waveriders at Mach 6

The performance of 23 kinds of waveriders, derived from different conical flowfields, is analyzed by the numerical computation under the conditions of fight speed of Mach 6, attack angle of 0° and flight altitude of 30 km. These results indicate that the performance is influenced by the shapes and the width to height ratios (W/H ) of generating cones. The geometrical parameter and the lift coefficient are proportional to W/H, while the drag coefficient and the lift to drag ratio (L/D ) have extreme values. Considering the base drag and the computation errors, the waverider with the highest L/D is cut from the elliptical cone’s flowfield (W/H = 1.5―1.618), and the configuration with the lowest drag can also be obtained at W/H = 1:1.5. Accordingly, good suggestions are proposed for practical design based on these computational results.

Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

Behavior of spherical particles at low Reynolds numbers in a fluctuating translational flow

The behavior of spheres in non-steady translational flow has been studied experimentally for values of Reynolds number from 0.2 to 3000. The aim of the work was to improve our qualitative understanding of particle transport in turbulent gaseous media, a process of extreme importance in power plants and energy transfer mechanisms.

Particles, subjected to sinusoidal oscillations parallel to the direction of steady translation, were found to have changes in average drag coefficient depending upon their translational Reynolds number, the density ratio, and the dimensionless frequency and amplitude of the oscillations. When the Reynolds number based on sphere diameter was less than 200, the oscillation had negligible effect on the average particle drag.

For Reynolds numbers exceeding 300, the coefficient of the mean drag was increased significantly in a particular frequency range. For example, at a Reynolds number of 3000, a 25 per cent increase in drag coefficient can be produced with an amplitude of oscillation of only 2 per cent of the sphere diameter, providing the frequency is near the frequency at which vortices would be shed in a steady flow at the mean speed. Flow visualization shows that over a wide range of frequencies, the vortex shedding frequency locks in to the oscillation frequency. Maximum effect at the natural frequency and lock-in show that a non-linear interaction between wake vortex shedding and the oscillation is responsible for the increase in drag.

Turbulent buoyant jets into stratified or flowing ambient fluids

Theoretical and experimental studies were made on two classes of buoyant jet problems, namely:

1) an inclined, round buoyant yet in a stagnant environment with linear density-stratification;

2) a round buoyant jet in a uniform cross stream of homogenous density.

Using the integral technique of analysis, assuming similarity, predictions can be made for jet trajectory, widths, and dilution ratios, in a density-stratified or flowing environment. Such information is of great importance in the design of disposal systems for sewage effluent into the ocean or waste gases into the atmosphere.

The present study of a buoyant jet in a stagnant environment has extended the Morton type of analysis to cover the effect of the initial angle of discharge. Numerical solutions have been presented for a range of initial conditions. Laboratory experiments were conducted for photographic observations of the trajectories of dyed jets. In general the observed jet forms agreed well with the calculated trajectories and nominal half widths when the value of the entrainment coefficient was taken to be α = 0.082, as previously suggested by Morton.

The problem of a buoyant jet in a uniform cross stream was analyzed by assuming an entrainment mechanism based upon the vector difference between the characteristic jet velocity and the ambient velocity. The effect of the unbalanced pressure field on the sides of the jet flow was approximated by a gross drag term. Laboratory flume experiments with sinking jets which are directly analogous to buoyant jets were performed. Salt solutions were injected into fresh water at the free surface in a flume. The jet trajectories, dilution ratios and jet half widths were determined by conductivity measurements. The entrainment coefficient, α, and drag coefficient, C_d, were found from the observed jet trajectories and dilution ratios. In the ten cases studied where jet Froude number ranged from 10 to 80 and velocity ratio (jet: current) K from 4 to 16, α varied from 0.4 to 0.5 and C_d from 1.7 to 0.1. The jet mixing motion for distance within 250D was found to be dominated by the self-generated turbulence, rather than the free-stream turbulence. Similarity of concentration profiles has also been discussed.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

I. Stokes flow past a thin screen. II. Viscous flows past porous bodies of finite size

Part I

The slow, viscous flow past a thin screen is analyzed based on Stokes equations. The problem is reduced to an associated electric potential problem as introduced by Roscoe. Alternatively, the problem is formulated in terms of a Stokeslet distribution, which turns out to be equivalent to the first approach.

Special interest is directed towards the solution of the Stokes flow past a circular annulus. A "Stokeslet" formulation is used in this analysis. The problem is finally reduced to solving a Fredholm integral equation of the second kind. Numerical data for the drag coefficient and the mean velocity through the hole of the annulus are obtained.

Stokes flow past a circular screen with numerous holes is also attempted by assuming a set of approximate boundary conditions. An "electric potential" formulation is used, and the problem is also reduced to solving a Fredholm integral equation of the second kind. Drag coefficient and mean velocity through the screen are computed.

Part II

The purpose of this investigation is to formulate correctly a set of boundary conditions to be prescribed at the interface between a viscous flow region and a porous medium so that the problem of a viscous flow past a porous body can be solved.

General macroscopic equations of motion for flow through porous media are first derived by averaging Stokes equations over a volume element of the medium. These equations, including viscous stresses for the description, are more general than Darcy's law. They reduce to Darcy's law when the Darcy number becomes extremely small.

The interface boundary conditions of the first kind are then formulated with respect to the general macroscopic equations applied within the porous region. An application of such equations and boundary conditions to a Poiseuille shear flow problem demonstrates that there usually exists a thin interface layer immediately inside the porous medium in which the tangential velocity varies exponentially and Darcy's law does not apply.

With Darcy's law assumed within the porous region, interface boundary conditions of the second kind are established which relate the flow variables across the interface layer. The primary feature is a jump condition on the tangential velocity, which is found to be directly proportional to the normal gradient of the tangential velocity immediately outside the porous medium. This is in agreement with the experimental results of Beavers, et al.

The derived boundary conditions are applied in the solutions of two other problems: (1) Viscous flow between a rotating solid cylinder and a stationary porous cylinder, and (2) Stokes flow past a porous sphere.

自由表面上球形浮体的拖曳运动仿真与试验

浮体在自由表面上的拖曳是一种常见的拖曳形式,由于自由表面上物体的运动情况相当复杂,目前对这种拖曳系统的设计和分析一般采用试验的方法。通过理论分析,对球形浮体在自由表面上的匀速直线拖曳运动建立了数学模型,其中对缆绳的分析采用Ablow-Schechter提出的有限差分模型,对球形浮体的分析采用水动力系数法,整个耦合系统用有限差分方法求解。之后,在拖曳水池进行了试验,试验结果表明,这种方法在拖曳速度比较低的情况下能够模拟真实情况,可以用来指导拖曳系统设计。

Dynamic and Thermodynamic Structure of the Monsoon Boundary Layer over India

The objective of the study is to examine the dynamic and thermodynamic structure and the variations that occur in the surface layer during the pre-monsoon, onset and post-monsoon periods over the Indian region. The variations caused during the occurrence of micro and mesoscale systems, structure and variation in the marine boundary layer over the Indian region is also investigated. The drag coefficient computed indirectly also shows variation during various seasons. The thermodynamic structure of the atmosphere shows variation during the various seasons. The onset monsoon causes lowering of the Lifting Condensation Levels. The outcome of the study is expected to provide a better understanding of the structure and variations in the boundary layer over India, which is useful for many applications especially for numerical modeling studies.