An economical method for direct numerical simulation studies of transitional round jets

A pseudo-spectral method based on Fourier expansions in a Cartesian coordinate system is shown to be an economical method for direct numerical simulation studies of transitional round jets, Several characteristics of the solutions are presented to establish the validity of the solutions in spite of the unnatural choices. We show that neither periodicity, nor the use of a Cartesian system have adversely affected the simulations, Instead, there are benefits in terms of ease of computing and lack of the usual restrictions due to grid structure near the jet axis. By computing the simultaneous evolution of passive scalers, the process of reaction in round jet burners, between a fuel-laden jet and an ambient oxidizer, was also simulated. Some typical solutions are shown and then the results of analysis of these data are summarized. (C) 2001 Elsevier Science Ltd, All rights reserved.

Toward a Mean Field Formulation of the Babcock-Leighton Type Solar Dynamo. I. α-Coefficient versus Durney's Double-Ring Approach

We develop a model of the solar dynamo in which, on the one hand, we follow the Babcock-Leighton approach to include surface processes, such as the production of poloidal field from the decay of active regions, and, on the other hand, we attempt to develop a mean field theory that can be studied in quantitative detail. One of the main challenges in developing such models is to treat the buoyant rise of the toroidal field and the production of poloidal field from it near the surface. A previous paper by Choudhuri, Schüssler, & Dikpati in 1995 did not incorporate buoyancy. We extend this model by two contrasting methods. In one method, we incorporate the generation of the poloidal field near the solar surface by Durney's procedure of double-ring eruption. In the second method, the poloidal field generation is treated by a positive α-effect concentrated near the solar surface coupled with an algorithm for handling buoyancy. The two methods are found to give qualitatively similar results.

Very large array detection of radio recombination lines from the radio nucleus of NGC 253: Ionization by a weak active galactic nucleus, an obscured super star cluster, or a compact supernova remnant?

We have imaged the H92alpha and H75alpha radio recombination line (RRL) emissions from the starburst galaxy NGC 253 with a resolution of similar to4 pc. The peak of the RRL emission at both frequencies coincides with the unresolved radio nucleus. Both lines observed toward the nucleus are extremely wide, with FWHMs of similar to200 km s(-1). Modeling the RRL and radio continuum data for the radio nucleus shows that the lines arise in gas whose density is similar to10(4) cm(-3) and mass is a few thousand M., which requires an ionizing flux of (6-20) x 10(51) photons s(-1). We consider a supernova remnant (SNR) expanding in a dense medium, a star cluster, and also an active galactic nucleus (AGN) as potential ionizing sources. Based on dynamical arguments, we rule out an SNR as a viable ionizing source. A star cluster model is considered, and the dynamics of the ionized gas in a stellar-wind driven structure are investigated. Such a model is only consistent with the properties of the ionized gas for a cluster younger than similar to10(5) yr. The existence of such a young cluster at the nucleus seems improbable. The third model assumes the ionizing source to be an AGN at the nucleus. In this model, it is shown that the observed X-ray flux is too weak to account for the required ionizing photon flux. However, the ionization requirement can be explained if the accretion disk is assumed to have a big blue bump in its spectrum. Hence, we favor an AGN at the nucleus as the source responsible for ionizing the observed RRLs. A hybrid model consisting of an inner advection-dominated accretion flow disk and an outer thin disk is suggested, which could explain the radio, UV, and X-ray luminosities of the nucleus.

Progressively More Prolate Dark Matter Halo In the Outer Galaxy As Traced by Flaring Hi gas

A galactic disk in a spiral galaxy is generally believed to be embedded in an extended dark matter halo, which dominates its dynamics in the outer parts. However, the shape of the halo is not clearly understood. Here we show that the dark matter halo in the Milky Way is prolate in shape. Further, it is increasingly more prolate at larger radii, with the vertical-to-planar axis ratio monotonically increasing to 2.0 at 24 kpc. This is obtained by modeling the observed steeply flaring atomic hydrogen gas layer in the outer Galactic disk, where the gas is supported by pressure against the net gravitational field of the disk and the halo. The resulting prolate-shaped halo can explain several long-standing puzzles in galactic dynamics, for example, it permits long-lived warps thus explaining their ubiquitous nature.

Regimes of spray formation in gas-centered swirl coaxial atomizers

Spray formation in ambient atmosphere from gas-centered swirl coaxial atomizers is described by carrying out experiments in a spray test facility. The atomizer discharges a circular air jet and an axisymmetric swirling water sheet from its coaxially arranged inner and outer orifices. A high-speed digital imaging system along with a backlight illumination arrangement is employed to record the details of liquid sheet breakup and spray development. Spray regimes exhibiting different sheet breakup mechanisms are identified and their characteristic features presented. The identified spray regimes are wave-assisted sheet breakup, perforated sheet breakup, segmented sheet breakup, and pulsation spray regime. In the regime of wave-assisted sheet breakup, the sheet breakup shows features similar to the breakup of two-dimensional planar air-blasted liquid sheets. At high air-to-liquid momentum ratios, the interaction process between the axisymmetric swirling liquid sheet and the circular air jet develops spray processes which are more specific to the atomizer studied here. The spray exhibits a periodic ejection of liquid masses whose features are dominantly controlled by the central air jet.

Numerical Study Of Heat Transfer From Pin Fin Heat Sink Using Steady And Pulsated Impinging Jets

The work reported in this thesis is an attempt to enhance heat transfer in electronic devices with the use of impinging air jets on pin-finned heat sinks. The cooling per-formance of electronic devices has attracted increased attention owing to the demand of compact size, higher power densities and demands on system performance and re-liability. Although the technology of cooling has greatly advanced, the main cause of malfunction of the electronic devices remains overheating. The problem arises due to restriction of space and also due to high heat dissipation rates, which have increased from a fraction of a W/cm2to 100s of W /cm2. Although several researchers have at-tempted to address this at the design stage, unfortunately the speed of invention of cooling mechanism has not kept pace with the ever-increasing requirement of heat re- moval from electronic chips. As a result, efficient cooling of electronic chip remains a challenge in thermal engineering. Heat transfer can be enhanced by several ways like air cooling, liquid cooling, phase change cooling etc. However, in certain applications due to limitations on cost and weight, eg. air borne application, air cooling is imperative. The heat transfer can be increased by two ways. First, increasing the heat transfer coefficient (forced convec- tion), and second, increasing the surface area of heat transfer (finned heat sinks). From previous literature it was established that for a given volumetric air flow rate, jet im-pingement is the best option for enhancing heat transfer coefficient and for a given volume of heat sink material pin-finned heat sinks are the best option because of their high surface area to volume ratio. There are certain applications where very high jet velocities cannot be used because of limitations of noise and presence of delicate components. This process can further be improved by pulsating the jet. A steady jet often stabilizes the boundary layer on the surface to be cooled. Enhancement in the convective heat transfer can be achieved if the boundary layer is broken. Disruptions in the boundary layer can be caused by pulsating the impinging jet, i.e., making the jet unsteady. Besides, the pulsations lead to chaotic mixing, i.e., the fluid particles no more follow well defined streamlines but move unpredictably through the stagnation region. Thus the flow mimics turbulence at low Reynolds number. The pulsation should be done in such a way that the boundary layer can be disturbed periodically and yet adequate coolant is made available. So, that there is not much variation in temperature during one pulse cycle. From previous literature it was found that square waveform is most effective in enhancing heat transfer. In the present study the combined effect of pin-finned heat sink and impinging slot jet, both steady and unsteady, has been investigated for both laminar and turbulent flows. The effect of fin height and height of impingement has been studied. The jets have been pulsated in square waveform to study the effect of frequency and duty cycle. This thesis attempts to increase our understanding of the slot jet impingement on pin-finned heat sinks through numerical investigations. A systematic study is carried out using the finite-volume code FLUENT (Version 6.2) to solve the thermal and flow fields. The standard k-ε model for turbulence equations and two layer zonal model in wall function are used in the problem Pressure-velocity coupling is handled using the SIMPLE algorithm with a staggered grid. The parameters that affect the heat transfer coefficient are: height of the fins, total height of impingement, jet exit Reynolds number, frequency of the jet and duty cycle (percentage time the jet is flowing during one complete cycle of the pulse). From the studies carried out it was found that: a) beyond a certain height of the fin the rate of enhancement of heat transfer becomes very low with further increase in height, b) the heat transfer enhancement is much more sensitive to any changes at low Reynolds number than compared to high Reynolds number, c) for a given total height of impingement the use of fins and pulsated jet, increases the effective heat transfer coefficient by almost 200% for the same average Reynolds number, d) for all the cases it was observed that the optimum frequency of impingement is around 50 − 100 Hz and optimum duty cycle around 25-33.33%, e) in the case of turbulent jets the enhancement in heat transfer due to pulsations is very less compared to the enhancement in case of laminar jets.

CFD Based Design Optimizations of the Diffuser of a Gas Dynamically Driver Laser Cavity

Based on the an earlier CFD analysis of the performance of the gas-dynamically controlled laser cavity [1]it was found that there is possibility of optimizing the geometry of the diffuser that can bring about reductions in both size and cost of the system by examining the critical dimensional requirements of the diffuser. Consequently,an extensive CFD analysis has been carried out for a range of diffuser configurations by simulating the supersonic flow through the arrangement including the laser cavity driven by a bank of converging – diverging nozzles and the diffuser. The numerical investigations with 3D-RANS code are carried out to capture the flow patterns through diffusers past the cavity that has multiple supersonic jet interactions with shocks leading to complex flow pattern. Varying length of the diffuser plates is made to be the basic parameter of the study. The analysis reveals that the pressure recovery pattern during the flow through the diffuser from the simulation, being critical for the performance of the laser device shows its dependence on the diffuser length is weaker beyond a critical lower limit and this evaluation of this limit would provide a design guideline for a more efficient system configuration.The observation based on the parametric study shows that the pressure recovery transients in the near vicinity of the cavity is not affected for the reduction in the length of the diffuser plates up to its 10% of the initial size, indicating the design in the first configuration that was tested experimentally has a large factor of margin. The flow stability in the laser cavity is found to be unaffected since a strong and stable shock is located at the leading edge of the diffuser plates while the downstream shock and flow patterns are changed, as one would expect. Results of the study for the different lengths of diffusers in the range of 10% to its full length are presented, keeping the experimentally tested configuration used in the earlier study [1] as the reference length. The conclusions drawn from the analysis is found to be of significance since it provides new design considerations based on the understanding of the intricacies of the flow, allowing for a hardware optimization that can lead to substantial size reduction of the device with no loss of performance.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Direct Infrared Absorption Spectroscopy of Benzene Dimer

The direct infrared (IR) absorption spectrum of benzene dimer formed in a free-jet expansion was recorded in the 3.3 mu m region for the first time. This has led to the observation of the C H stretching fundamental mode nu(13) (B(1u)), which is both IR and Raman forbidden in the monomer. Moreover, the IR forbidden and Raman allowed nu(7) (E(2g)) mode has been observed as well. These two modes were found to be red-shifted along with the IR allowed nu(20) (E(1u)) mode, as previously reported by Erlekam et al. [Erlekam; Frankowski; Meijer; Gert von Helden J. Chem. Phys. 2006, 124, 171101], using ion-dip spectroscopy, contrary to the blue-shift predicted earlier by theoretical studies. The observation of the nu(13) band indicates that the symmetry is reduced in the dimer, confirming the T-shaped structure observed by Erlekam et al. Our experimental results have not provided any direct evidence for the presence of the parallel displaced geometry, the main objective of the present work, as predicted by theoretical calculations.

Computational Analysis of Flow through a Multiple Nozzle Driven Laser Cavity and Diffuser

Computations have been carried out for simulating supersonic flow through a set of converging-diverging nozzles with their expanding jets forming a laser cavity and flow patterns through diffusers, past the cavity. A thorough numerical investigation with 3-D RANS code is carried out to capture the flow distribution which comprises of shock patterns and multiple supersonic jet interactions. The analysis of pressure recovery characteristics during the flow through the diffusers is an important parameter of the simulation and is critical for the performance of the laser device. The results of the computation have shown a close agreement with the experimentally measured parameters as well as other established results indicating that the flow analysis done is found to be satisfactory.

Dispersion and vaporization of biofuels and conventional fuels in a crossflow pre-mixer

In lean premixed pre-vaporized (LPP) combustion, controlled atomization, dispersion and vaporization of different types of liquid fuel in the premixer are the key factors required to stabilize the combustion process and improve the efficiency. The dispersion and vaporization process for biofuels and conventional fuels sprayed into a crossflow pre-mixer have been simulated and analyzed with respect to vaporization rate, degree of mixedness and homogeneity. Two major biofuels under investigation are Ethanol and Rapeseed Methyl Esters (RME), while conventional fuels are gasoline and jet-A. First, the numerical code is validated by comparing with the experimental data of single n-heptane and decane droplet evaporating under both moderate and high temperature convective air now. Next, the spray simulations were conducted with monodispersed droplets with an initial diameter of 80 mu m injected into a turbulent crossflow of air with a typical velocity of 10 m/s and temperature of around 800K. Vaporization time scales of different fuels are found to be very different. The droplet diameter reduction and surface temperature rise were found to be strongly dependent on the fuel properties. Gasoline droplet exhibited a much faster vaporization due a combination of higher vapor pressure and smaller latent heat of vaporization compared to other fuels. Mono-dispersed spray was adopted with the expectation of achieving more homogeneous fuel droplet size than poly-dispersed spray. However, the diameter histogram in the zone near the pre-mixer exit shows a large range of droplet diameter distributions for all the fuels. In order to improve the vaporization performance, fuels were pre-heated before injection. Results show that the Sauter mean diameter of ethanol improved from 52.8% of the initial injection size to 48.2%, while jet-A improved from 48.4% to 18.6% and RME improved from 63.5% to 31.3%. The diameter histogram showed improved vaporization performance of jet-A. (C) 2011 Elsevier Ltd. All rights reserved.

Flow-acoustic Characterisation of a Cavity-based Combustor Configuration

This study concerns the flow-acoustic characterisation of a cavity-based combustor configuration. A well-validated numerical tool has been used to simulate the unsteady, two-dimensional reacting flow. Initially, a conventional flow over a cavity with dimensions and conditions corresponding to a compact cavity combustor was studied. Cavity mass injections in the form of fuel and air injections required for trapped vortex formation were then employed and the resonance features of this configuration were studied. The results indicate that the cavity depth mode resonance mechanism is dominant at the conditions studied in this work and that the oscillation frequencies do not change with cavity air injection. This observation is important since it implies that the only important variable which can alter resonant frequencies is the cavity depth. With combustion, the pressure oscillation amplitude was observed to increases significantly due to periodic entrainment of the cavity air jet and fluctuation of fuel-air mixture composition to produce highly fluctuating heat-release rates. The underlying mechanisms of the unsteady flow in the cavity combustor identified in this study indicate the strong dependence of the acoustics on the cavity injection strategies.

Fluid Mechanics of Aquatic Locomotion at Large Reynolds Numbers

Abstract | There exist a huge range of fish species besides other aquatic organisms like squids and salps that locomote in water at large Reynolds numbers, a regime of flow where inertial forces dominate viscous forces. In the present review, we discuss the fluid mechanics governing the locomotion of such organisms. Most fishes propel themselves by periodic undulatory motions of the body and tail, and the typical classification of their swimming modes is based on the fraction of their body that undergoes such undulatory motions. In the angulliform mode, or the eel type, the entire body undergoes undulatory motions in the form of a travelling wave that goes from head to tail, while in the other extreme case, the thunniform mode, only the rear tail (caudal fin) undergoes lateral oscillations. The thunniform mode of swimming is essentially based on the lift force generated by the airfoil like crosssection of the fish tail as it moves laterally through the water, while the anguilliform mode may be understood using the “reactive theory” of Lighthill. In pulsed jet propulsion, adopted by squids and salps, there are two components to the thrust; the first due to the familiar ejection of momentum and the other due to an over-pressure at the exit plane caused by the unsteadiness of the jet. The flow immediately downstream of the body in all three modes consists of vortex rings; the differentiating point being the vastly different orientations of the vortex rings. However, since all the bodies are self-propelling, the thrust force must be equal to the drag force (at steady speed), implying no net force on the body, and hence the wake or flow downstream must be momentumless. For such bodies, where there is no net force, it is difficult to directly define a propulsion efficiency, although it is possible to use some other very different measures like “cost of transportation” to broadly judge performance.

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Real source-sink pair

This work is concerned with the interaction of a source-sink pair. The main parameters of the problem are source and sink flow rates, the axial and lateral separations of the source and sink, and the angle between the axes of source and sink. Of concern is the percentage of source fluid that enters the sink as a function of these parameters. The experiments have been carried using the source nozzle diameter of 6 mm and the sink pipe diameter of two sizes: 10 mm and 20 mm. The Reynolds numbers of the source jet is about 3200. The main diagnostics are flow visualization using dye, laser induced fluorescence (LIF), particle streak photographs and particle image velocimetry (Ply). To obtain the removal effectiveness (that is percentage of source fluid that is going through the sink pipe), titration method is used. The sink diameter and the angle between source and the sink axes do not influence efficiencies as do the sink flow rate and the lateral separation. Data from experiments have been consolidated so that these results can be used for designing sinks for removal of heat and pollutants. (C) 2011 Elsevier Ltd. All rights reserved.