Large eddy simulation of a premixed flame with approximate deconvolution modeling

Approximate deconvolution modeling is a very recent approach to large eddy simulation of turbulent flows. It has been applied to compressible flows with success. Here, a premixed flame which forms in the wake of a flameholder has been selected to examine the subgrid-scale modeling of reaction rate by this new method because a previous plane two-dimensional simulation of this wake flame, using a wrinkling function and artificial flame thickening, had revealed discrepancies when compared with experiment. The present simulation is of the temporal evolution of a round wakelike flow at two Reynolds numbers, Re = 2000 and 10,000, based on wake defect velocity and wake diameter. A Fourier-spectral code has been used. The reaction is single-step and irreversible, and the rate follows an Arrhenius law. The reference simulation at the lower Reynolds number is fully resolved. At Re = 10,000, subgrid-scale contributions are significant. It was found that subgrid-scale modeling in the present simulation agrees more closely with unresolved subgrid-scale effects observed in experiment. Specifically, the highest contributions appeared in thin folded regions created by vortex convection. The wrinkling function approach had not selected subgrid-scale effects in these regions.

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Equilibrium and relaxation in turbulent wakes

In order to study the memory of the larger eddies in turbulent shear flow, experiments have been conducted on plane turbulent wakes undergoing transition from an initial (carefully prepared) equilibrium state to a different final one, as a result of a nearly impulsive pressure gradient. It is shown that under the conditions of the experiments the equations of motion possess self-preserving solutions in the sense of Townsend (1956), but the observed behaviour of the wake is appreciably different when the pressure gradient is not very small, as the flow goes through a slow relaxation process before reaching final equilibrium. Measurements of the Reynolds stresse show that the approach to a new equilibrium state is exponential, with a relaxation length of the order of 103 momentum thicknesses. It is suggested that a flow satisfying the conditions required by a self-preservation analysis will exhibit equilibrium only if the relaxation length is small compared with a characteristic streamwise length scale of the flow.

Turbulent non-equlibrium wakes

We consider here the detailed application of a model Reynolds stress equation (Narasimha 1969) to plane turbulent wakes subjected to pressure gradients. The model, which is a transport equation for the stress exhibiting relaxation and diffusion, is found to be consistent with the observed response of a wake to a nearly impulsive pressure gradient (Narasimha & Prabhu 1971). It implies in particular that a wake can be in equilibrium only if the longitudinal strain rate is appreciably less than the wake shear. We then describe a further series of experiments, undertaken to investigate the range of validity of the model. It is found that, with an appropriate convergence correction when necessary, the model provides excellent predictions of wake development under favourable, adverse and mixed pressure gradients. Furthermore, the behaviour of constant-pressure distorted wakes, as reported by Keffer (1965, 1967), is also explained very well by the model when account is taken of the effective flow convergence produced by the distortion. In all these calculations, only a simple version of the model is used, involving two non-dimensional constants both of which have been estimated from a single relaxation experiment.

Influence of matrix characteristics on fracture toughness of high volume fraction Al2O3/Al–AlN composites

The role of matrix microstructure on the fracture of Al-alloy composites with 60 vol% alumina particulates was studied. The matrix composition and microstructure were systematically varied by changing the infiltration temperature and heat treatment. Characterization was carried out by a combination of metallography, hardness measurements, and fracture studies conducted on compact tension specimens to study the fracture toughness and crack growth in the composites. The composites showed a rise in crack resistance with crack extension (R curves) due to bridges of intact matrix ligaments formed in the crack wake. The steady-state or plateau toughness reached upon stable crack growth was observed to be more sensitive to the process temperature rather than to the heat treatment. Fracture in the composites was predominantly by particle fracture, extensive deformation, and void nucleation in the matrix. Void nucleation occurred in the matrix in the as-solutionized and peak-aged conditions and preferentially near the interface in the underaged and overaged conditions. Micromechanical models based on crack bridging by intact ductile ligaments were modified by a plastic constraint factor from estimates of the plastic zone formed under indentations, and are shown to be adequate in predicting the steady-state toughness of the composite.

On an oscillatory point force in a rotating stratified fluid

The forced oscillations due to a point forcing effect in an infinite or contained, inviscid, incompressible, rotating, stratified fluid are investigated taking into account the density variation in the inertia terms in the linearized equations of motion. The solutions are obtained in closed form using generalized Fourier transforms. Solutions are presented for a medium bounded by a finite cylinder when the oscillatory forcing effect is acting at a point on the axis of the cylinder. In both the unbounded and bounded case, there exist characteristic cones emanating from the point of application of the force on which either the pressure or its derivatives are discontinuous. The perfect resonance existing at certain frequencies in an unbounded or bounded homogeneous fluid is avoided in the case of a confined stratified fluid.

Defining the ‘modified Griffin plot’ in vortex-induced vibration: revealing the effect of Reynolds number using controlled damping

In the present work, we study the transverse vortex-induced vibrations of an elastically mounted rigid cylinder in a fluid flow. We employ a technique to accurately control the structural damping, enabling the system to take on both negative and positive damping. This permits a systematic study of the effects of system mass and damping on the peak vibration response. Previous experiments over the last 30 years indicate a large scatter in peak-amplitude data ($A^*$) versus the product of mass–damping ($\alpha$), in the so-called ‘Griffin plot’. A principal result in the present work is the discovery that the data collapse very well if one takes into account the effect of Reynolds number ($\mbox{\textit{Re}}$), as an extra parameter in a modified Griffin plot. Peak amplitudes corresponding to zero damping ($A^*_{{\alpha}{=}0}$), for a compilation of experiments over a wide range of $\mbox{\textit{Re}}\,{=}\,500-33000$, are very well represented by the functional form $A^*_{\alpha{=}0} \,{=}\, f(\mbox{\textit{Re}}) \,{=}\, \log(0.41\,\mbox{\textit{Re}}^{0.36}$). For a given $\mbox{\textit{Re}}$, the amplitude $A^*$ appears to be proportional to a function of mass–damping, $A^*\propto g(\alpha)$, which is a similar function over all $\mbox{\textit{Re}}$. A good best-fit for a wide range of mass–damping and Reynolds number is thus given by the following simple expression, where $A^*\,{=}\, g(\alpha)\,f(\mbox{\textit{Re}})$: \[ A^* \,{=}\,(1 - 1.12\,\alpha + 0.30\,\alpha^2)\,\log (0.41\,\mbox{\textit{Re}}^{0.36}). \] In essence, by using a renormalized parameter, which we define as the ‘modified amplitude’, $A^*_M\,{=}\,A^*/A^*_{\alpha{=}0}$, the previously scattered data collapse very well onto a single curve, $g(\alpha)$, on what we refer to as the ‘modified Griffin plot’. There has also been much debate over the last three decades concerning the validity of using the product of mass and damping (such as $\alpha$) in these problems. Our results indicate that the combined mass–damping parameter ($\alpha$) does indeed collapse peak-amplitude data well, at a given $\mbox{\textit{Re}}$, independent of the precise mass and damping values, for mass ratios down to $m^*\,{=}\,1$.

Tuning the structure of surfactant complexes with DNA and other polyelectrolytes

We have carried out small-angle X-ray diffraction studies on complexes formed by the anionic polyelectrolytes, namely, sodium salts of double and single stranded (ds and ss) DNA, poly( glutamic acid) ( PGA), poly( acrylic acid) (PAA), and poly( styrene sulfonate) (PSS) with a cationic surfactant system consisting of cetyltrimethylammonium bromide ( CTAB) and sodium 3-hydroxy-2-naphthoate (SHN). All complexes have a two-dimensional (2D) hexagonal structure at low SHN concentrations. DNA-CTAB-SHN complexes exhibit a hexagonal to lamellar transition near the SHN concentration at which CTAB-SHN micelles show a cylinder to bilayer transformation. On the other hand, PGA and PAA complexes form a 2D centered rectangular phase at higher SHN concentrations, and PSS complexes show a primitive rectangular structure. These results provide a striking example of polyion specificity in polyelectrolytesurfactant interactions.

Continuous Recording of Retardation and Intensity of Echoes from the Ionosphere

Pulse retardation method of Breit and Tuve has been modified to record continuously the equivalent height as well as the intensity of reflections from the ionosphere. Synchronized pulses are transmitted, and the received ground pulse and the reflected pulses, after amplification and suitable distortion, are applied to the focusing cylinder of a cathode ray tube the horizontal deflecting plates of which are connected to a synchronized linear time base circuit. The pattern on the screen is composed of a bright straight line corresponding to the time base with dark gaps corresponding to the received pulses. The distance between the initial points of the gaps represents retardation while the widths of the gaps correspond to the intensity of the pulses. The pattern is photographed on a vertically moving film. One of the first few records taken at Bangalore on 4 megacycles is reproduced. It shows, among other things, that the less retarded component of magneto-ionic splitting from the F layer is present most of the time. Whenever the longer retardation component does occur, it has stronger intensity than the former. Towards the late evening hours, just before disappearing, when the F layer rises and exhibits magnetoionic splitting, the intensity of the less retarded component is extremely low compared with the other component.

Filling Characteristics for an Activated Carbon Based Adsorbed Natural Gas Storage System

The storage capacity of an activated carbon bed is studied using a 2D transport model with constant inlet flow conditions. The predicted filling times and variation in bed pressure and temperature are in good agreement with experimental observations obtained using a 1.82 L prototype ANG storage cylinder. Storage efficiencies based on the maximum achievable V/V (volume of gas/volume of container) and filling times are used to quantify the performance of the charging process. For the high permeability beds used in the experiments, storage efficiencies are controlled by the rate of heat removal. Filling times, defined as the time at which the bed pressure reaches 3.5 MPa, range from 120 to 3.4 min for inlet flow rates of 1.0 L min(-1) and 30.0 L min(-1), respectively. The corresponding storage efficiencies, eta(s), vary from 90% to 76%, respectively. Simulations with L/D ratios ranging from 0.35 to 7.8 indicate that the storage efficiencies can be improved with an increase in the LID ratios and/or with water cooled convection. Thus for an inlet flow rate of 30.0 L min(-1), an eta(s) value of 90% can be obtained with water cooling for an L/D ratio of 7.8 and a filling time of a few minutes. In the absence of water cooling the eta(s) value reduces to 83% at the same L/D ratio. Our study suggests that with an appropriate choice of cylinder dimensions, solutions based on convective cooling during adsorptive storage are possible with some compromise in the storage capacity.

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.

Viscous effects in the inception of cavitation

The inception of cavitation in the steady flow of liquids around bodies is seen to depend upon the real fluid flow around the bodies as well as the supply of nucleating cavitation sources—or nuclei—within the fluid. A primary distinction is made between bodies having a laminar separation or not having a laminar separation. The former group is relatively insensitive to the nuclei concentration whereas the latter is much more sensitive. Except for the case of fully separated wake flows and for gaseous cavitation by diffusion the cavitation inception index tends always to be less than the magnitude of the minimum pressure coefficient and only approaches that value for high Reynolds numbers in flows well supplied with nuclei.

Nonthermal plasma approach in direct methanol synthesis from CH4

Direct methanol synthesis from CH4 and O2 has been experimentally studied using pulsed discharge plasma in concentric-cylinder-type reactors. The methanol production becomes efficient with an increase in the average electric field strength of the reactor. A combination of the pulsed discharge and catalysts was tested and was proved to be effective in increasing both the production and selectivity of methanol. In the present stage, about 2% of CH4 can be converted into other hydrocarbons, and a methanol yield of around 0.5% and selectivity of 38% can be obtained when a catalyst of V2O5+SiO2 is combined with the pulsed discharge plasma

Treatment of NOx from Diesel Engine Exhaust by Dielectric Barrier Discharge Method

This paper reports improved performance of advantages when compared to its counterpart as it is cost discharge plasma in filtered engine exhaust treatment. Our effective, low capital and operation costs, salable by- paper deals about the removal of NOX emissions from the diesel products, and integration with the existing systems. In this exhaust by electric discharge plasma. For the treatment of diesel paper we describe an alternate reactor geometry referred to exhaust a new type of reactor referred to as cross-flow dielectric as cross-flow DBD reactor, where the exhaust gas flow barrier discharge reactor has been used, where the gas flow is perpendicular to the wire-cylinder reaction chamber. This perpendicular to the corona electrode. Experiments were reactor is used to treat the actual exhaust of a 3.75 kW diesel- conducted at different flow rates ranging from 2 l/min to 10 l/ generator set. The main emphasis is laid on the NOX treatment min. The discharge plasma assisted barrier discharge reactor of diesel engine exhaust. Experiments were carried out at has shown promising results in NOX removal at high flow rates.

A core study on NOX removal in diesel exhaust by pulsed/ac/dc electric discharge plasma

In this paper a study on effect of different energization on removal of NOX in diesel engine exhaust has been presented. Here we made a detailed qualitative study of effect of pulsed/ac/dc voltage energizations on the NOX treatment of using conventional wire-cylinder reactor configuration. It was observed that amongst different energizations, pulse energization exhibits maximum NOX removal efficiency when compared to ac and dc energizations. For a given specific energy density, wire-cylinder reactor filled with BaTiO3 pellet gives higher NOX removal efficiency when compared to reactor without pellets under both pulse and ac energization. The dc energization does not have much impact on the removal processes. The paper further discusses the individual energization cases in detail.