Coupling Eulerian‐Lagrangian method of air‐particle two‐phase flow with population balance equations to simulate the evolution of vehicle exhaust plume

In this paper, we present a new numerical scheme to describe the dynamic evolution of multiphase polydisperse systems in terms of time, space, and properties by coupling the Eulerian‐Lagrangian method for air‐particle two‐phase flow and population balance equations to describe particle property evolution due to microbehaviors (eg, aggregation, breakage, and growth). This coupling scheme was used to comprehensively simulate the two‐phase flow structure, particle size spectrum, particle number, and volume concentrations. These were characterized by a high‐resolution particle tracking using the Lagrangian approach and the high precision of moments of the particle size spectrum by solving the population balance equation with the quadrature method of moments. The algorithm of the coupling scheme was incorporated into the open source computational fluid dynamics software OpenFOAM to simulate the dynamic evolution of vehicle exhaust plume. The impacts of vehicle velocity, exhaust temperature, and aggregation efficiency on the distribution of auto exhaust particles in space and changes in their properties were analyzed. The results indicate that the particle number concentration, volume concentration, and average diameter of particles in the vehicle exhaust plume could be strongly affected by the plume structure and flow properties.     

Large eddy simulations of a turbulent thermal plume

Large eddy simulations of a three-dimensional turbulent thermal plume in an open environment have been carried out using a self-developed parallel computational fluid dynamics code SMAFS (smoke movement and flame spread) to study the thermal plume’s dynamics including its puffing, self-preserving and air entrainment. In the simulation, the sub-grid stress was modeled using both the standard Smagorinsky and the buoyancy modified Smagorinsky models, which were compared. The sub-grid scale (SGS) scalar flux in the filtered enthalpy transport equation was modeled based on a simple gradient transport hypothesis with constant SGS Prandtl number. The effect of the Smagorinsky model constant and the SGS Prandtl number were examined. The computation results were compared with experimental measurements, thermal plume theory and empirical correlations, showing good agreement. It is found that both the buoyancy modification and the SGS turbulent Prandtl number have little influence on simulation. However, the SGS model constant C _s has a significant effect on the prediction of plume spreading, although it does not affect much the prediction of puffing.     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Rayleigh scattering in supersonic high-temperature exhaust plumes

A supersonic exhaust plume test rig and a Rayleigh scattering system were developed. Molecular number densities in the supersonic high-temperature exhaust plume with and without an annular base flow were investigated. The physical meaning of the inferred mean temperature from the number density measurement in turbulent flows is clarified. For both flows, the potential core extends up to about six nozzle diameters, and self-similarity of the radial density distributions is observed at downstream sections Z/d=10–50. The recovery of the flow density deficit (or the decay of temperature) with the annular flow is faster than that without the annular flow at upstream sections Z/d ≤ 10. Received: 16 August 2000 / Accepted: 20 November 2001     

Temperature spectra from a turbulent thermal plume by ultrasound scattering

This paper reports the results of a study on temperature inhomogeneities conducted on a thermal plume by using ultrasound scattering as a non-intrusive measurement technique. The plume rises from a metallic disk which can be heated up to 800 °C. The working fluid is air at atmospheric pressure. In the measurement technique, an incoming ultrasound wave is emitted towards the thermal plume. The incident wave is scattered because of non-linear couplings with the flow instabilities present in the measurement region. The scattered wave carries information about those flow instabilities. The technique allows for the retrieving of this information. The shape of the obtained spectrum of temperature fluctuations as a function of wave vector modulus is consistent with previous theoretical analysis. Three qualitatively different regions were identified: first, a production region characterized by a q² law; secondly, a region with behavior as per q⁻³ associated with a buoyancy region and; finally, a dissipation region associated with a q⁻⁷ law. These spectral regions characterize the energy transfers mechanisms among the length scales of flow investigated here. A coefficient of anisotropy γ was defined to analyze anisotropic features of the flow.     

Direct Numerical Simulation of the Puffing Phenomenon of an Axisymmetric Thermal Plume

A spatial direct numerical simulation of an axisymmetric buoyant thermal plume is presented. The governing flow field equations at the centerline are put into a special form to circumvent the axis singularity associated with the cylindrical coordinates and the high order accuracy of the numerical scheme is preserved at the centerline. Boundary conditions associated with the spatial DNS of open-boundary buoyant flows and compatible with the modern nondissipative high-order finite difference schemes have been developed. The fluid exhibits a periodic oscillatory motion known as the puffing phenomenon, which is the formation and convection of vortex at the near field of the plume. Budgets of the vorticity transport are determined to examine the mechanisms leading to the puffing phenomenon. The analysis on vorticity transport shows that vorticity is created mainly by the gravitational term which is due to the interaction between the radial density gradients and gravity at the initial stage of the establishment of the puffing structure, while the baroclinic torque dominates the vorticity transport when the flow is established. Density stratification in the radial direction close to the plume base is found to be essential to the development of the buoyant flow instability. Simulations with different initial temperature ratios reveal that entrainment close to the plume base is enhanced at a higher temperature ratio despite the fact that the puffing structures and the plume pulsation frequency only vary very weakly with the initial temperature ratio. The predicted puffing frequencies are in agreement with the values from experimental correlations for fire and isothermal helium/air plumes. Received 12 May 1999 and accepted 9 December 1999     

Design of an exhaust manifold to improve transient performance of a high-speed turbocharged diesel engine

In this paper a study to analyse the influence of exhaust gases thermal energy saving and reduction of pressure pulses interference on engine dynamic performance during the load transient of high speed direct injection turbocharged diesel engines is described. The analysis has been performed thanks to the use of a dual walled air gap exhaust manifold and a 4-2-1 pulse exhaust manifold, which have not only been modeled but also manufactured and tested on an engine test bench able to carry out constant speed load transient tests. Finally, in order to estimate maximum improvement in engine performance, an optimum design of a synthesis exhaust manifold has been proposed and modeled showing that exist a potential not only to improve engine transient performance at low engine speeds but also the gas temperature at the catalyst inlet and the steady effective torque.     

Prediction of entrance length and mass suction rate for a cylindrical sucking funnel

Conservation equations for mass, momentum and energy have been solved numerically for a cylindrical funnel with louvers (lateral openings on the side wall of the cylindrical funnel through which air can come into it) to compute the suction rate of air into the funnel. The nozzle placed centrally at the bottom of the cylinder ejects high‐velocity hot gaseous products so that atmospheric air gets sucked into the funnel. The objective of the work is to compute the ratio of the rate of mass suction to that of the mass ejected by the nozzle for different operating conditions and geometrical size of the funnel. From the computation it has been found that there exists optimum funnel diameter and optimum funnel height for which the mass suction is the highest. The protruding length of the nozzle into the funnel has almost no effect on the mass suction rate after a certain funnel height. The louvers opening area has a very high impact on the mass suction rate. The entrance length for such a sucking funnel is strikingly much lower compared with a simple cylindrical pipe having uniform flow at the inlet at same Reynolds number. A new correlation has been developed to propose the entrance length for a sucking pipe, the rate of mass suction into it and the exhaust plume temperature over a wide range of operating parameters that are normally encountered in a general funnel operations of naval or merchant ship. Copyright © 2009 John Wiley & Sons, Ltd.     

Semi‐coupled air/water immersed boundary approach for curvilinear dynamic overset grids with application to ship hydrodynamics

For many problems in ship hydrodynamics, the effects of air flow on the water flow are negligible (the frequently called free surface conditions), but the air flow around the ship is still of interest. A method is presented where the water flow is decoupled from the air solution, but the air flow uses the unsteady water flow as a boundary condition. The authors call this a semi‐coupled air/water flow approach. The method can be divided into two steps. At each time step the free surface water flow is computed first with a single‐phase method assuming constant pressure and zero stress on the interface. The second step is to compute the air flow assuming the free surface as a moving immersed boundary (IB). The IB method developed for Cartesian grids (Annu. Rev. Fluid Mech. 2005; 37 :239–261) is extended to curvilinear grids, where no‐slip and continuity conditions are used to enforce velocity and pressure boundary conditions for the air flow. The forcing points close to the IB can be computed and corrected under a sharp interface condition, which makes the computation very stable. The overset implementation is similar to that of the single‐phase solver (Comput. Fluids 2007; 36 :1415–1433), with the difference that points in water are set as IB points even if they are fringe points. Pressure–velocity coupling through pressure implicit with splitting of operators or projection methods is used for water computations, and a projection method is used for the air. The method on each fluid is a single‐phase method, thus avoiding ill‐conditioned numerical systems caused by large differences of fluid properties between air and water. The computation is only slightly slower than the single‐phase version, with complete absence of spurious velocity oscillations near the free surface, frequently present in fully coupled approaches. Validations are performed for laminar Couette flow over a wavy boundary by comparing with the analytical solution, and for the surface combatant model David Taylor Model Basin (DTMB) 5512 by comparing with Experimental Fluid Dynamics (EFD) and the results of two‐phase level set computations. Complex flow computations are demonstrated for the ONR Tumblehome DTMB 5613 with superstructure subject to waves and wind, including 6DOF motions and broaching in SS7 irregular waves and wind. Copyright © 2008 John Wiley & Sons, Ltd.     

CFD Analysis of a Thermal Plume and the Indoor Air Flow Using k-ε Models with Buoyancy Effects

Turbulence models based on the eddy viscosity concept perform poorly for simulation of non-isothermal flows, which are characterized by strongly non-isotropic features of turbulence. In the present study, a thermal plume is calculated using the WET model and severalk-ε models. Comparison with the experiments reveals that the WET model predicts the flow and temperature fields most accurately. Furthermore, the WET model has been modified to account for the damping effect near the wall and applied to the indoor air flow. The computed flow field agrees well with the measurement, except for higher temperatures appearing near the walls. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical study on heat/mass transfer by shear/thermal/solutal-convection in a cavity under counteracting buoyancies

This study provides physical insights into a lid-driven square cavity filled with a mixture of a solvent vapor and non-condensable gas, subjected to the vertically parallel thermal and solutal gradients. The top lid is maintained at constant speed while bottom lid and the other two walls are kept fixed. Zero heat and mass fluxes are imposed on the vertical side walls. The transport equations are solved numerically through a pressure-correction-based iterative algorithm (SIMPLE) with the QUICK scheme for convective terms. The diffusivities of heat and salt are assumed to be equal throughout this investigation. The essential details of flow, temperature and concentration fields are presented for the opposing buoyancy forces ratio (B < 0) with special attention being given for the values of parameters for which either the flow inside the cavity is operated by the mechanically induced convection; or the flow structure inside the cavity is akin to a single-diffusive thermal or solutal convection. The variations of average rates of heat and mass transfer are uniform with the Reynolds number, while the variations of these quantities against the solutal Richardson number (Ri _C ) and thermal Richardson number (Ri _T ) point to the existence of the local minimum/maximum. Finally, two linear relations between Ri _C and Ri _T at a constant sliding speed are proposed to identify the above points of high and low transport phenomena and justified with the exhibition of the flow structures inside the cavity.     

Stability of a melt/solid interface with reference to a plume channel

The structure of free-convection flow in a plume channel formed as a result of melting above a local heat source placed on the basement of a solid mass is experimentally investigated. The channel shape and the flow pattern in it are functions of the relative power Ka = N/N ₁, where N is the plume source power and N ₁ is the heat removed to the surrounding mass. When the heat is withdrawn from the plume channel by heat conduction, the channel represents a system of convective cells on whose boundaries there are channel constrictions. The temperature fields and the cell flow patterns are investigated. For mantle plumes, such as the Hawaiian, Iceland, and Bouvet plumes and extended igneous provinces, the basement diameter and the values of the criterion Ka are determined.     

Effect of swirl intensity on the flow and combustion of a turbulent non-premixed flat flame

An experiment in a turbulent non-premixed flat flame was carried out in order to investigate the effect of swirl intensity on the flow and combustion characteristics. First, stream lines and velocity distribution in the flow field were obtained using PIV (Particle Image Velocimetry) method in a model burner. In contrast with the axial flow without swirl, highly swirled air induced streamlines going along the burner tile, and its backward flow was generated by recirculation in the center zone of the flow field. In the combustion, the flame shape with swirled air also became flat and stable along the burner tile with increment of the swirl number. Flame structure was examined by measuring OH and CH radicals intensity and by calculating Damkohler number (Da) and turbulence Reynolds number (Re _T ). It appeared that luminescence intensity decreased at higher swirl number due to the recirculated flue gas, and the flat flames were comprised in the wrinkled laminar-flame regime. Backward flow by recirculation of the flue gas widely contacted on the flame front, and decreased the flame temperature and emissions concentration as thermal NO. The homogeneous temperature field due to the widely flat flame was obtained, and the RMS in the high temperature region was rather lower at higher swirl number. Consequently, the stable flat flame with low NO concentration was achieved.     

PIV, 2D-LIF and 1D-Raman measurements of flow field,composition and temperature in premixed gas turbine flames

Several laser diagnostic measurement techniques have been applied to study the lean premixed natural gas/air flames of an industrial swirl burner. This was made possible by equipping the burner with an optical combustion chamber that was installed in the high-pressure test rig facility at the DLR Institute of Combustion Technology in Stuttgart. The burner was operated with preheated air at various operating conditions with pressures up to p = 6 bar and a maximum thermal power of P = 1 MW.The instantaneous planar flow field inside the combustor was studied with particle image velocimetry (PIV). Planar laser induced fluorescence (PLIF) of OH radicals on a single-shot basis was used to determine the shape and the location of the flame front as well as the spatial distribution of reaction products. 1D laser Raman spectroscopy was successfully applied for the measurement of the temperature and the concentration of major species under realistic gas turbine conditions.Results of the flow field analysis show the shape and the size of the main flow regimes: the inflow region, the inner and the outer recirculation zone. The highly turbulent flow field of the inner shear layer is found to be dominated by small and medium sized vortices. High RMS fluctuations of the flow velocity in the exhaust gas indicate the existence of a rotating exhaust gas swirl. From the PLIF images it is seen that the primary reactions happened in the shear layers between inflow and the recirculation zones and that the appearance of the reaction zones changed with flame parameters. The results of the multiscalar Raman measurements show a strong variation of the local mixture fraction allowing conclusions to be drawn about the premix quality. Furthermore, mixing effects of unburnt fuel and air with fully reacted combustion products are studied giving insights into the processes of the turbulence–chemistry interaction.     

Coherent motion of turbulent thermal plume in stably stratified fluid

The purpose of this study is to clarify the existence of an ordered and large scale coherent motion in a turbulent plane thermal plume in a thermally-stable stratified fluid inside a comparatively large enclosure. First, the upper part of the thermal plume was carefully observed by a flow visualization. Secondly, a wave form of plume temperature variation was measured. Thirdly, a spectrum analysis was carried out on time series data of the thermal plume. Finally, physical characteristics were investigated on vortices in the thermal plume based on results of the wave form and the spectrum analysis of the plume temperature. As a result, the main conclusions are obtained as follows. (1) An existence of vortices near the upper part of the thermal plume was firstly found by careful flow visualization. (2) From the wave form of temperature variation and the spectrum analysis of the thermal plume, it was clarified that the vortices are generated in the transition state and are transported to the turbulent state. (3) The vortices are ordered and they behave as a large scale coherent motion in the turbulent thermal plume.     

Numerical simulation of the galvanizing process during GA to GI transition

This paper presents the application of a three‐dimensional finite element solution algorithm for the prediction of the velocity, temperature and species concentration fields in an industrial continuous galvanizing bath. Simulations were carried out using a parallel CFD software developed at IMI‐NRC. The turbulent flow, heat and mass transfer has been solved using a high Reynolds number k–ε model. Simulations were carried out for the case when the density of the molten metal depends only on the temperature and also for the case when both temperature and Al concentration affect the density. When considering the buoyancy effect of the Al concentration, differences are especially apparent during the melting of ingots with high Al content. Otherwise, thermal effects are dominant. The continuous monitoring of the temperature and the Al and Fe content in an industrial bath was used to validate the flow, temperature and compositional variations. A period of three hours, corresponding to three different ingot additions, was simulated successfully, resulting in a good agreement of the temperature and compositional variations. Copyright © 2006 Crown in the right of Canada. Published by John Wiley & Sons, Ltd.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Large-eddy simulation of buoyancy-driven natural ventilation in an enclosure with a point heat source

Rising buoyant plumes from a point heat source in a naturally ventilated enclosure have been investigated using large-eddy simulation (LES). The aim of the work is to assess the performance and the accuracy of LES for modelling buoyancy-driven displacement ventilation of an enclosure and to shed more light on the transitional behaviour of the plume and the coherent structures involved. The Smagorinsky sub-grid scale model is used for the unresolved small-scale turbulence. The Rayleigh number, Ra is chosen to be in the range where spatial transition from laminar to turbulent flow takes place (Ra = 1.5 × 10⁹). The plume properties (source strength and rate of spread) as well as the ventilation properties (stratification height and temperature of stratified layer) estimated using the theory of Linden et al. are found to agree reasonably well with the LES results. The variation of the plume width with height indicates a linear variation of the entrainment coefficient rather than a constant value used by Linden et al. for a fully turbulent thermal plume. Flow visualisation revealed the nature of the large-scale coherent structures involved in the transition to turbulence in the plume. The most excited modes observed in the velocity, pressure and temperature fields spectra correspond to Strouhal number in the range 0.3 ≤ St ≤ 0.55 which is in agreement with those observed by Zhou et al. for a turbulent forced plume. Excited modes less than thisvalue (St = 0.2) were observed and may be due to low-frequency motions felt throughout the flow.     

用于低温风洞的新颖制冷方法

描述了用于低温风洞的新颖制冷系统,利用热交换器回收排气冷量预冷压缩空气,然后再用热分离器将其降至深低温作风洞气源．原理性实验结果证实新制冷方法的可行性．讨论了新制冷方法产生的有一定压力的低温空气作引射气源,引射驱动回流型风洞的特性．其制冷方法与现有低温风洞喷雾液氮制冷相比,由于仅需压缩空气而无需液氮,造价更便宜．更由于能量利用合理,效率高,因而运行成本可显著降低．