Explicit filtering LES for turbulent non-premixed combustion

Two models for large eddy simulation of turbulent reacting flow in homogeneous turbulence were studied. The sub-grid stress arising out of non-linearities of the Navier-Stokes equations were modeled using an explicit filtering approach. A filtered mass density function (FMDF) approach was used for closure of the sub-grid scalar fluctuations. A posteriori calculations, when compared with the results from the direct numerical simulation, indicate that the explicit filtering is adequate in representing the effect of sub-grid stress on the filtered velocity field in the absence of reaction. Discrepancies arise when reactions occur, but the FMDF approach suffices to account for sub-grid scale fluctuations of the reacting scalars, accurately.

Regional Flood Frequency Analysis using Two-level Clustering Approach

Clustering techniques are used in regional flood frequency analysis (RFFA) to partition watersheds into natural groups or regions with similar hydrologic responses. The linear Kohonen's self‐organizing feature map (SOFM) has been applied as a clustering technique for RFFA in several recent studies. However, it is seldom possible to interpret clusters from the output of an SOFM, irrespective of its size and dimensionality. In this study, we demonstrate that SOFMs may, however, serve as a useful precursor to clustering algorithms. We present a two‐level. SOFM‐based clustering approach to form regions for FFA. In the first level, the SOFM is used to form a two‐dimensional feature map. In the second level, the output nodes of SOFM are clustered using Fuzzy c‐means algorithm to form regions. The optimal number of regions is based on fuzzy cluster validation measures. Effectiveness of the proposed approach in forming homogeneous regions for FFA is illustrated through application to data from watersheds in Indiana, USA. Results show that the performance of the proposed approach to form regions is better than that based on classical SOFM.

Multicomponent solids of lamotrigine with some selected coformers and their characterization by thermoanalytical, spectroscopic and X-ray diffraction methods

The present study investigates the structural and pharmaceutical properties of different multicomponent crystalline forms of lamotrigine (LTG) with some pharmaceutically acceptable coformers viz. nicotinamide (1), acetamide (2), acetic acid (3), 4-hydroxy-benzoic acid (4) and saccharin (5). The structurally homogeneous phases were characterized in the solid state by DSC/TGA, FT-IR and XRD (powder and single crystal structure analysis) as well as in the solution phase. Forms 1 and 2 were found to be cocrystal hydrate and cocrystal, respectively, while in forms 3, 4 and 5, proton transfer was observed from coformer to drug. The enthalpy of formation of multicomponent crystals from their components was determined from the enthalpy of solution of the cocrystals and the components separately. Higher exothermic values of the enthalpy of formation for molecular complexes 3, 4 and 5 suggest these to be more stable than 1 and 2. The solubility was measured in water as well as in phosphate buffers of varying pH. The salt solvate 3 exhibited the highest solubility of the drug in water as well as in buffers over the pH range 7-3 while the cocrystal hydrate 1 showed the maximum solubility in a buffer of pH 2. A significant lowering of the dosage profile of LTG was observed for 1, 3 and 5 in the animal activity studies on mice.

Thermal Weight Functions and Stress Intensity Factors for Bonded Dissimilar Media Using Body Analogy

In this study, an analytical method is presented for the computation of thermal weight functions in two dimensional bi-material elastic bodies containing a crack at the interface and subjected to thermal loads using body analogy method. The thermal weight functions are derived for two problems of infinite bonded dissimilar media, one with a semi-infinite crack and the other with a finite crack along the interface. The derived thermal weight functions are shown to reduce to the already known expressions of thermal weight functions available in the literature for the respective homogeneous elastic body. Using these thermal weight functions, the stress intensity factors are computed for the above interface crack problems when subjected to an instantaneous heat source.

Calculation of the front part of the sonic boom signature for a Maneuvering aerofopil

The sonic boom at a large distance from its source consists of a leading shock, a trailing shock and a one parameter family of nonlinear wavefronts in between these shocks. A new ray theoretical method using a shock ray theory and a weakly nonlinear lay theory has been used to obtain the shock fronts and wavefronts respectively, for a maneuvering aerofoil in a homogeneous medium. This method introduces a one parameter family of Cauchy problems to calculate the shock and wave fronts emerging from the surface of the aerofoil. These problems are solved numerically to obtain the leading shock front and the nonlinear wavefronts emerging from the front portion of the aerofoil.

New Insights into Selective Heterogeneous Nucleation of Metal Nanoparticles on Oxides by Microwave-Assisted Reduction: Rapid Synthesis of High-Activity Supported Catalysts

Microwave-based methods are widely employed to synthesize metal nanoparticles on various substrates. However, the detailed mechanism of formation of such hybrids has not been addressed. In this paper, we describe the thermodynamic and kinetic aspects of reduction of metal salts by ethylene glycol under microwave heating conditions. On the basis of this analysis, we identify the temperatures above which the reduction of the metal salt is thermodynamically favorable and temperatures above which the rates of homogeneous nucleation of the metal and the heterogeneous nucleation of the metal on supports are favored. We delineate different conditions which favor the heterogeneous nucleation of the metal on the supports over homogeneous nucleation in the solvent medium based on the dielectric loss parameters of the solvent and the support and the metal/solvent and metal/support interfacial energies. Contrary to current understanding, we show that metal particles can be selectively formed on the substrate even under situations where the temperature of the substrate Is lower than that of the surrounding medium. The catalytic activity of the Pt/CeO(2) and Pt/TiO(2) hybrids synthesized by this method for H(2) combustion reaction shows that complete conversion is achieved at temperatures as low as 100 degrees C with Pt-CeO(2) catalyst and at 50 degrees C with Pt-TiO(2) catalyst. Our method thus opens up possibilities for rational synthesis of high-activity supported catalysts using a fast microwave-based reduction method.

Thermodynamic properties of oxygen in RE-O (RE = Gd, Tb, Dy, Er) solid solutions

The oxygen potentials of four rare-earth metal – oxygen (RE–O: RE=Gd, Dy, Tb, Er) solid solutions have been measured by equilibration with yttrium – oxygen (Y–O) and titanium – oxygen (Ti–O) solid solutions. Rare-earth metal, yttrium and titanium samples were immersed in calcium-saturated CaCl2 melt at temperatures between 1093 and 1233 K. Homogeneous oxygen potential was established in the metallic samples through the fused salt, which contains some dissolved CaO. The metallic samples were analyzed for oxygen after quenching. The oxygen potentials of RE–O solid solutions were determined using either Y–O or Ti–O solid solution as the reference. This method enabled reliable measurement of extremely low oxygen potentials at high temperature (circa pO2=10−48 atm at 1173 K). It was found that the oxygen affinity of the metals decreases in the order: Y>Er>Dy>Tb>Gd>Ti. Values for the standard Gibbs energy of solution of oxygen in RE metals obtained in this study, permit assessment of the extent of deoxidation that can be achieved with various purification techniques. It may be possible to achieve an oxygen level of 10 mass ppm using an electrochemical deoxidation method.

Trapped fermions in a synthetic non-Abelian gauge field

On increasing the coupling strength (lambda) of a non-Abelian gauge field that induces a generalized Rashba spin-orbit interaction, the topology of the Fermi surface of a homogeneous gas of noninteracting fermions of density rho similar to k(F)(3) undergoes a change at a critical value, lambda(T) approximate to k(F) [Phys. Rev. B 84, 014512 ( 2011)]. In this paper we analyze how this phenomenon affects the size and shape of a cloud of spin-1/2 fermions trapped in a harmonic potential such as those used in cold atom experiments. We develop an adiabatic formulation, including the concomitant Pancharatnam-Berry phase effects, for the one-particle states in the presence of a trapping potential and the gauge field, obtaining approximate analytical formulas for the energy levels for some high symmetry gauge field configurations of interest. An analysis based on the local density approximation reveals that, for a given number of particles, the cloud shrinks in a characteristic fashion with increasing.. We explain the physical origins of this effect by a study of the stress tensor of the system. For an isotropic harmonic trap, the local density approximation predicts a spherical cloud even for anisotropic gauge field configurations. We show, via a calculation of the cloud shape using exact eigenstates, that for certain gauge field configurations there is a systematic and observable anisotropy in the cloud shape that increases with increasing gauge coupling lambda. The reasons for this anisotropy are explained using the analytical energy levels obtained via the adiabatic approximation. These results should be useful in the design of cold atom experiments with fermions in non-Abelian gauge fields. An important spin-off of our adiabatic formulation is that it reveals exciting possibilities for the cold-atom realization of interesting condensed matter Hamiltonians by using a non-Abelian gauge field in conjunction with another potential. In particular, we show that the use of a spherical non-Abelian gauge field with a harmonic trapping potential produces a monopole field giving rise to a spherical geometry quantum Hall-like Hamiltonian in the momentum representation.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

Structural, optical, and electrical properties of bulk single crystals of InAsxSb(1-x) grown by rotatory Bridgman method

Radially-homogeneous and single-phase InAsxSb(1−x) crystals, up to 5.0 at. % As concentration, have been grown using the rotatory Bridgman method. Single crystallinity has been confirmed by x-ray and electron diffraction studies. Infrared transmission spectra show a continuous decrease in optical energy gap with the increase of arsenic content in InSb. The measured values of mobility and carrier density at room temperature (for x = .05) are 5.6×104 cm2/V s and 2.04×1016 cm−3, respectively.

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.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Diffuse optical tomographic imager using a single light source

Near-infrared diffuse optical tomography (DOT) technique has the capability of providing good quantitative reconstruction of tissue absorption and scattering properties with additional inputs such as input and output modulation depths and correction for the photon leakage. We have calculated the two-dimensional (2D) input modulation depth from three-dimensional (3D) diffusion to model the 2D diffusion of photons. The photon leakage when light traverses from phantom to the fiber tip is estimated using a solid angle model. The experiments are carried for single (5 and 6 mm) as well as multiple inhomogeneities (6 and 8 mm) with higher absorption coefficient in a homogeneous phantom. Diffusion equation for photon transport is solved using finite element method and Jacobian is modeled for reconstructing the optical parameters. We study the development and performance of DOT system using modulated single light source and multiple detectors. The dual source methods are reported to have better reconstruction capabilities to resolve and localize single as well as multiple inhomogeneities because of its superior noise rejection capability. However, an experimental setup with dual sources is much more difficult to implement because of adjustment of two out of phase identical light probes symmetrically on either side of the detector during scanning time. Our work shows that with a relatively simpler system with a single source, the results are better in terms of resolution and localization. The experiments are carried out with 5 and 6 mm inhomogeneities separately and 6 and 8 mm inhomogeneities both together with absorption coefficient almost three times as that of the background. The results show that our experimental single source system with additional inputs such as 2D input/output modulation depth and air fiber interface correction is capable of detecting 5 and 6 mm inhomogeneities separately and can identify the size difference of multiple inhomogeneities such as 6 and 8 mm. The localization error is zero. The recovered absorption coefficient is 93% of inhomogeneity that we have embedded in experimental phantom.

Influence of substrate temperature on the properties of oxygen‐ion‐assisted deposited CeO2 films

Substrate temperature and ion bombardment during deposition have been observed to modify significantly the optical and structural properties of dielectric thin films. Single‐layer films of CeO2 have been deposited by electron beam evaporation with simultaneous oxygen‐ion bombardment using a Kaufman broad beam ion source and maintaining the substrates at elevated temperature. A systematic study has been made on the influence of (a) substrate temperature in the range ambient to 300 °C, (b) ion energy in the range 300–700 eV, and (c) ion current density 100–220 μA/cm2 on optical properties such as refractive index, extinction coefficient, inhomogeneity, packing density, and structural properties. The refractive index increased with in increase in substrate temperature: ion energy up to 600 eV and ion current density. Homogeneous, absorption free and high index (2.48) films have been obtained at 600 eV, 220 μA/cm2 and at substrate temperature of 300 °C. The packing density of the films was observed to be unity for the same deposition conditions. Substrate temperature with simultaneous ion bombardment modified the structure of the films from highly ordered to fine grain structure.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 1. DNS and fluctuating force model

The fluctuating force model is developed and applied to the turbulent flow of a gas-particle suspension in a channel in the limit of high Stokes number, where the particle relaxation time is large compared to the fluid correlation time, and low particle Reynolds number where the Stokes drag law can be used to describe the interaction between the particles and fluid. In contrast to the Couette flow, the fluid velocity variances in the different directions in the channel are highly non-homogeneous, and they exhibit significant variation across the channel. First, we analyse the fluctuating particle velocity and acceleration distributions at different locations across the channel. The distributions are found to be non-Gaussian near the centre of the channel, and they exhibit significant skewness and flatness. However, acceleration distributions are closer to Gaussian at locations away from the channel centre, especially in regions where the variances of the fluid velocity fluctuations are at a maximum. The time correlations for the fluid velocity fluctuations and particle acceleration fluctuations are evaluated, and it is found that the time correlation of the particle acceleration fluctuations is close to the time correlations of the fluid velocity in a `moving Eulerian' reference, moving with the mean fluid velocity. The variances of the fluctuating force distributions in the Langevin simulations are determined from the time correlations of the fluid velocity fluctuations and the results are compared with direct numerical simulations. Quantitative agreement between the two simulations are obtained provided the particle viscous relaxation time is at least five times larger than the fluid integral time.