In-situ measurements of concentration and temperature during transient solidification of aqueous solution of ammonium chloride using laser interferometry

The variation in temperature and concentration plays a crucial role in predicting the final microstructure during solidification of a binary alloy. Most of the experimental techniques used to measure concentration and temperature are intrusive in nature and affect the flow field. In this paper, the main focus is laid on in-situ, non-intrusive, transient measurement of concentration and temperature during the solidification of a binary mixture of aqueous ammonium chloride solution (a metal-analog system) in a top cooled cavity using laser based Mach-Zehnder Interferometric technique. It was found from the interferogram, that the angular deviation of fringe pattern and the total number of fringes exhibit significant sensitivity to refractive index and hence are functions of the local temperature and concentration of the NH4Cl solution inside the cavity. Using the fringe characteristics, calibration curves were established for the range of temperature and concentration levels expected during the solidification process. In the actual solidification experiment, two hypoeutectic solutions (5% and 15% NH4Cl) were chosen. The calibration curves were used to determine the temperature and concentration of the solution inside the cavity during solidification of 5% and 15% NH4Cl solution at different instants of time. The measurement was carried out at a fixed point in the cavity, and the concentration variation with time was recorded as the solid-liquid interface approached the measurement point. The measurement exhibited distinct zones of concentration distribution caused by solute rejection and Rayleigh Benard convection. Further studies involving flow visualization with laser scattering confirmed the Rayleigh Benard convection. Computational modeling was also performed, which corroborated the experimental findings. (C) 2011 Elsevier Ltd. All rights reserved.

Particle image velocimetry and infrared thermography in a levitated droplet with nanosilica suspensions

Preferential accumulation and agglomeration kinetics of nanoparticles suspended in an acoustically levitated water droplet under radiative heating has been studied. Particle image velocimetry performed to map the internal flow field shows a single cell recirculation with increasing strength for decreasing viscosities. Infrared thermography and high speed imaging show details of the heating process for various concentrations of nanosilica droplets. Initial stage of heating is marked by fast vaporization of liquid and sharp temperature rise. Following this stage, aggregation of nanoparticles is seen resulting in various structure formations. At low concentrations, a bowl structure of the droplet is dominant, maintained at a constant temperature. At high concentrations, viscosity of the solution increases, leading to rotation about the levitator axis due to the dominance of centrifugal motion. Such complex fluid motion inside the droplet due to acoustic streaming eventually results in the formation of a ring structure. This horizontal ring eventually reorients itself due to an imbalance of acoustic forces on the ring, exposing larger area for laser absorption and subsequent sharp temperature rise.

Investigations on supersonic combustion in a combustion-driven shock tunnel

This article reports experimental results on supersonic combustion in a new facility. The facility is a combustion-driven shock tunnel, which is cheaper to build than the facilities in which such experiments are carried out conventionally. The observation region is a zone between two parallel flat plates with a 33 degrees wedge attached to the upstream end of the bottom plate. Gaseous hydrogen is injected at an angle of 45 degrees into an oncoming supersonic flow of Mach 2 (approximate) from a port on the bottom plate. The resulting flow field is visualized by a high speed camera in a dark background. Three different test gases, namely nitrogen, air, and oxygen-rich air are used, and the results are compared. A distinct luminosity due to combustion for oxygen-containing test gases is observed. Heat-transfer rates on a probe placed at the downstream end of the observation region and midway between the parallel plates are measured and compared for the three cases. Wall static pressure at 28 mm downstream of the injection port on the bottom plate is also presented.

Full Chemical Kinetic Simulation of Biogas Early Phase Combustion in SI Engines

This paper presents computational work on the biogas early phase combustion in spark ignition (SI) engines using detailed chemical kinetics. Specifically, the early phase combustion is studied to assess the effect of various ignition parameters such as spark plug location, spark energy, and number of spark plugs. An integrated version of the KIVA-3V and CHEMKIN codes was developed and used for the simulations utilizing detailed kinetics involving 325 reactions and 53 species The results show that location of the spark plug and local flow field play an important role. A central plug configuration, which is associated with higher local flow velocities in the vicinity of the spark plug, showed faster initial combustion. Although a dual plug configuration shows the highest rate of fuel consumption, it is comparable to the rate exhibited by the central plug case. The radical species important in the initiation of combustion are identified, and their concentrations are monitored during the early phase of combustion. The concentration of these radicals is also observed to correlate very well with the above-mentioned trend.Thus, the role of these radicals in promoting faster combustion has been clearly established. It is also observed that the minimum ignition energy required to initiate a self-sustained flame depends on the flow field condition in the vicinity of the spark plug.Increasing the methane content in the biogas has shown improved combustion.

Mixing enhancement in a compact trapped vortex combustor

Previous studies on a single-cavity, compact trapped vortex combustor concept showed good flame stability for a wide range of flow conditions. However, achieving good mixing between cavity products and mainstream flow was still a major challenge. In the present study, a passive mixing enhancement strategy of using inclined struts along with a flow guide vane is presented and experimentally tested at atmospheric pressure conditions. Results show excellent mixing and consequently low values of the combustor exit pattern factor in the range of 0.1 and small flame lengths (57 times the main-duct depth). The pressure drop is small in the range of 0.35%, and NOx levels of the order of 12ppm are achieved. The flame stability is excellent, and combustion efficiency is reasonable in the range of 96%. The effectiveness of the proposed strategy is explained on the basis of in-situ OH chemiluminescence images and prior numerical simulations of the resulting complex flow field. The flow guide vane is observed to lead to a counterclockwise cavity vortex, which is conducive to the rise of cavity combustion products along the inclined struts and subsequent mixing with the mainstream flow.

Micro-blast waves using detonation transmission tubing

Micro-blast waves emerging from the open end of a detonation transmission tube were experimentally visualized in this study. A commercially available detonation transmission tube was used (Nonel tube, M/s Dyno Nobel, Sweden), which is a small diameter tube coated with a thin layer of explosive mixture (HMX traces of Al) on its inner side. The typical explosive loading for this tube is of the order of 18 mg/m of tube length. The blast wave was visualized using a high speed digital camera (frame rate 1 MHz) to acquire time-resolved schlieren images of the resulting flow field. The visualization studies were complemented by computational fluid dynamic simulations. An analysis of the schlieren images showed that although the blast wave appears to be spherical, it propagates faster along the tube axis than along a direction perpendicular to the tube axis. Additionally, CFD analysis revealed the presence of a barrel shock and Mach disc, showing structures that are typical of an underexpanded jet. A theory in use for centered large-scale explosions of intermediate strength gave good agreement with the blast trajectory along the tube axis. The energy of these micro-blast waves was found to be J and the average TNT equivalent was found to be . The repeatability in generating these micro-blast waves using the Nonel tube was very good and this opens up the possibility of using this device for studying some of the phenomena associated with muzzle blasts in the near future.

Yanai waves in the western equatorial Indian Ocean

Observations and models have shown the presence of intraseasonal fluctuations in 20-30-day and 10-20-day bands in the equatorial Indian Ocean west of 60 degrees E (WEIO). Their spatial and temporal structures characterize them as Yanai waves, which we label low-frequency (LFYW) and high-frequency (HFYW) Yanai waves, respectively. We explore the dynamics of these intraseasonal signals, using an ocean general circulation model (Modular Ocean Model) and a linear, continuously stratified model. Yanai waves are forced by the meridional wind tau(y) everywhere in the WEIO most strongly during the monsoon seasons. They are forced both directly in the interior ocean and by reflection of the interior response from the western boundary; interference between the interior and boundary responses results in a complex surface pattern that propagates eastward and has nodes. Yanai waves are also forced by instabilities primarily during June/July in a region offshore from the western boundary (52-55 degrees E). At that time, eddies, generated by barotropic instability of the Southern Gyre, are advected southward to the equator. There, they generate a westward-propagating, cross-equatorial flow field, v(eq), with a wave number/frequency spectrum that fits the dispersion relation of a number of Yanai waves, and these waves are efficiently excited. Typically, Yanai waves associated with several baroclinic modes are excited by both wind and eddy forcing; and typically, they superpose to create beams that carry energy vertically and eastward along ray paths. The same processes generate LFYWs and HFYWs, and hence, their responses are similar; differences are traceable to the property that HFYWs have longer wavelengths than LFYWs for each baroclinic mode.

Visualization of coherent structures in turbulent subsonic jet using planar laser induced fluorescence of acetone

In this paper, we present the molecular density distribution measurement in turbulent nitrogen jet (Re approximate to 3 x 10(3)), using acetone as molecular tracer. The tracer was seeded in the nitrogen jet by purging through the liquid acetone at ambient temperature. Planar laser sheet of 266 nm wavelength from frequency quadrupled, Q-switched, Nd:YAG laser was used as an excitation source. Emitted fluorescence images of jet flow field were recorded on CMOS camera. The dependence of planar laser induced fluorescence (PLIF) intensity on acetone vapor density was used to convert PLIF image of nitrogen jet into the density image on pixel-by-pixel basis. Instantaneous quantitative density image of nitrogen jet, seeded with acetone, was obtained. The arrowhead-shaped coherent turbulent structures were observed in the present work. It was found that coherent structures were non-overlapping with separate boundaries. Breaking of coherent structures into turbulence was clearly observed above four times jet width.

Measurement of temperature and pressure on the surface of a blunt cone using FBG sensor in hypersonic wind tunnel

Measurement of temperature and pressure exerted on the leeward surface of a blunt cone specimen has been demonstrated in the present work in a hypersonic wind tunnel using fiber Bragg grating (FBG) sensors. The experiments were conducted on a 30 degrees apex-angle blunt cone with 51 mm base diameter at wind flow speeds of Mach 6.5 and 8.35 in a 300 mm hypersonic wind tunnel of Indian Institute of Science, Bangalore. A special pressure insensitive temperature sensor probe along with the conventional bare FBG sensors was used for explicit temperature and aerodynamic pressure measurement respectively on the leeward surface of the specimen. computational fluid dynamics (CFD) simulation of the flow field around the blunt cone specimen has also been carried out to obtain the temperature and pressure at conditions analogous to experiments. The results obtained from FBG sensors and the CFD simulations are found to be in good agreement with each other.

Effect of side confining walls on the growth rate of compressible mixing layers

Model free simulations are performed to study the effect of the presence of side wall in compressible mixing of two parallel dissimilar gaseous streams with significant temperature difference. The turbulence statistics shows the three dimensional nature of the flow with and without the presence of side walls. The presence of side wall neither makes the flow field two dimensional, nor suppresses three dimensional disturbances. However, the comparison of shear layer growth rate and wall pressures reveal a better match with the two dimensional simulation results. This better match is explained on the basis of formation of oblique structures due to the presence of side walls which also suppress the distribution of momentum in third direction making the pressures to be higher as compared with the case without side walls. (C) 2013 Elsevier Ltd. All rights reserved.

An enthalpy method for modeling eutectic solidification

This paper presents a new micro-scale model for solidification of eutectic alloys. The model is based on the enthalpy method and simulates the growth of adjacent alpha and beta phases from a melt of eutectic composition in a two-dimensional Eulerian framework. The evolution of the two phases is obtained from the solution of volume averaged energy and species transport equations which are formulated using the nodal enthalpy and concentration potential values. The three phases are tracked using the beta-phase fraction and the liquid fraction values in all the computational nodes. Solutal convection flow field in the domain is obtained from the solution of volume-averaged momentum and continuity equations. The governing equations are solved using a coupled explicit-implicit scheme. The model is qualitatively validated with Jackson-Hunt theory. Results show expected eutectic growth pattern and proper species transfer and diffusion field ahead of the interface. Capabilities of the model such as lamella width selection, division of lamella into thinner lamellae and the presence of solutal convection are successfully demonstrated. The present model can potentially be incorporated into the existing framework of enthalpy based micro-scale dendritic solidification models thus leading to an efficient generalized microstructure evolution model. (C) 2014 Elsevier Inc. All rights reserved.

Detonation-driven-shock wave interactions with perforated plates

The study of detonations and their interactions is vital for the understanding of the high-speed flow physics involved and the ultimate goal of controlling their detrimental effects. However, producing safe and repeatable detonations within the laboratory can be quite challenging, leading to the use of computational studies which ultimately require experimental data for their validation. The objective of this study is to examine the induced flow field from the interaction of a shock front and accompanying products of combustion, produced from the detonation taking place within a non-electrical tube lined with explosive material, with porous plates with varying porosities, 0.7-9.7%. State of the art high-speed schlieren photography alongside high-resolution pressure measurements is used to visualise the induced flow field and examine the attenuation effects which occur at different porosities. The detonation tube is placed at different distances from the plates' surface, 0-30 mm, and the pressure at the rear of the plate is recorded and compared. The results indicate that depending on the level of porosity and the Mach number of the precursor shock front secondary reflected and transmitted shock waves are formed through the coalescence of compression waves. With reduced porosity, the plates act almost as a solid surface, therefore the shock propagates faster along its surface.

Shock tunnel experiments on control of shock induced large separation bubble using boundary layer bleed

Shock-Boundary Layer Interaction (SBLI) often occurs in supersonic/hypersonic flow fields. Especially when accompanied by separation (termed strong interaction), the SBLI phenomena largely affect the performance of the systems where they occur, such as scramjet intakes, thus often demanding the control of the interaction. Experiments on the strong interaction between impinging shock wave and boundary layer on a flat plate at Mach 5.96 are carried out in IISc hypersonic shock tunnel HST-2. The experiments are performed at moderate flow total enthalpy of 1.3 MJ/kg and freestream Reynolds number of 4 million/m. The strong shock generated by a wedge (or shock generator) of large angle 30.96 degrees to the freestream is made to impinge on the flat plate at 95 mm (inviscid estimate) from the leading edge, due to which a large separation bubble of length (75 mm) comparable to the distance of shock impingement from the leading edge is generated. The experimental simulation of such large separation bubble with separation occurring close to the leading edge, and its control using boundary layer bleed (suction and tangential blowing) at the location of separation, are demonstrated within the short test time of the shock tunnel (similar to 600 mu s) from time resolved schlieren flow visualizations and surface pressure measurements. By means of suction - with mass flow rate one order less than the mass flow defect in boundary layer - a reduction in separation length by 13.33% was observed. By the injection of an array of (nearly) tangential jets in the direction of mainstream (from the bottom of the plate) at the location of separation - with momentum flow rate one order less than the boundary layer momentum flow defect - 20% reduction in separation length was observed, although the flow field was apparently unsteady. (C) 2014 Elsevier Masson SAS. All rights reserved.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.