Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

Phase diagram of a hard-sphere system in a quenched random potential: A numerical study

We report numerical results for the phase diagram in the density-disorder plane of a hard-sphere system in the presence of quenched, random, pinning disorder. Local minima of a discretized version of the Ramakrishnan-Yussouff free energy functional are located numerically and their relative stability is studied as a function of the density and the strength of disorder. Regions in the phase diagram corresponding to liquid, glassy, and nearly crystalline states are mapped out, and the nature of the transitions is determined. The liquid to glass transition changes from first to second order as the strength of the disorder is increased. For weak disorder, the system undergoes a first-order crystallization transition as the density is increased. Beyond a critical value of the disorder strength, this transition is replaced by a continuous glass transition. Our numerical results are compared with those of analytical work on the same system. Implications of our results for the field-temperature phase diagram of type-II superconductors are discussed.

A numerical study of the role of the vertical structure of vorticity during tropical cyclone genesis

An eight-level axisymmetric model with simple parameterizations for clouds and the atmospheric boundary layer was developed to examine the evolution of vortices that are precursors to tropical cyclones. The effect of vertical distributions of vorticity, especially that arising from a merger of mid-level vortices, was studied by us to provide support for a new vortex-merger theory of tropical cyclone genesis. The basic model was validated with the analytical results available for the spin-down of axisymmetric vortices. With the inclusion of the cloud and boundary layer parameterizations, the evolution of deep vortices into hurricanes and the subsequent decay are simulated quite well. The effects of several parameters such as the initial vortex strength, radius of maximum winds, sea-surface temperature and latitude (Coriolis parameter) on the evolution were examined. A new finding is the manner in which mid-level vortices of the same strength decay and how, on simulated merger of these mid-level vortices, the resulting vortex amplifies to hurricane strength in a realistic time frame. The importance of sea-surface temperature on the evolution of full vortices was studied and explained. Also it was found that the strength of the surface vortex determines the time taken by the deep vortex to amplify to hurricane strength.

Calculation of Coulomb interaction strengths for 3d transition metals and actinides

Coulomb interaction strengths (Udd and Uff) have been calculated from Hartree-Fock-Slater atomic calculations for 3d transition and 5f actinide elements, respectively. By decomposing the different contributions to the response (screening) to the 3d charge fluctuation, we show that a substantial reduction in Udd arises due to the relaxation of the 3d charge distribution itself. This, combined with the screening due to the response of the 4s charge density, is shown to provide a very compact screening charge comparable to the metallic case, explaining the success of the atomic calculations for estimating U even in the metals. A pronounced dependence of Udd (or Uff) on the number of electrons nd (nf) or the electronic configuration is also shown here.

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

Sliding Modes for the Simulation of Mechanical and Electrical Systems Defined by Differential-Algebraic Equations

We offer a technique, motivated by feedback control and specifically sliding mode control, for the simulation of differential-algebraic equations (DAEs) that describe common engineering systems such as constrained multibody mechanical structures and electric networks. Our algorithm exploits the basic results from sliding mode control theory to establish a simulation environment that then requires only the most primitive of numerical solvers. We circumvent the most important requisite for the conventionalsimulation of DAEs: the calculation of a set of consistent initial conditions. Our algorithm, which relies on the enforcement and occurrence of sliding mode, will ensure that the algebraic equation is satisfied by the dynamic system even for inconsistent initial conditions and for all time thereafter. [DOI:10.1115/1.4001904]

Spatially dependent correlation functions in the Anderson model

We report results from a first principles calculation of spatially dependent correlation functions around a magnetic impurity in metals described by the nondegenerate Anderson model. Our computations are based on a combination of perturbative scaling theory and numerical renormalization group methods. Results for the conduction election charge density around the impurity and correlation functions involving the conduction electron and impurity charge and spin densities will be presented. The behavior in various regimes including the mixed valent regime will be explored.

Correlation of Oxygen Storage Capacity and Structural Distortion in Transition-Metal-, Noble-Metal-, and Rare-Earth-Ion-Substituted CeO2 from First Principles Calculation

Oxygen storage/release (OSC) capacity is an important feature common to all three-way catalysts to combat harmful exhaust emissions. To understand the mechanism of improved OSC for doped CeO2, we undertook the structural investigation by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), H-2-TPR (temperature-programmed hydrogen reduction) and density functional theoretical (DFT) calculations of transition-metal-, noble-metal-, and rare-earth (RE)-ion-substituted ceria. In this report, we present the relationship between the OSC and structural changes induced by the dopant ion in CeO2. Transition metal and noble metal ion substitution in ceria greatly enhances the reducibility of Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu, Pd, Pt, Ru), whereas rare-earth-ion-substituted Ce(1-x)A(x)O(2-delta) (A = La, Y) have very little effect in improving the OSC. Our simulated optimized structure shows deviation in cation oxygen bond length from ideal bond length of 2.34 angstrom (for CeO2). For example, our theoretical calculation for Ce28Mn4O62 structure shows that Mn-O bonds are in 4 + 2 coordination with average bond lengths of 2.0 and 3.06 angstrom respectively. Although the four short Mn-O bond lengths spans the bond distance region of Mn2O3, the other two Mn-O bonds are moved to longer distances. The dopant transition and noble metal ions also affects Ce coordination shell and results in the formation of longer Ce-O bonds as well. Thus longer cation oxygen bonds for both dopant and host ions results in enhanced synergistic reduction of the solid solution. With Pd ion substitution in Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu) further enhancement in OSC is observed in H-2-TPR. This effect is reflected in our model calculations by the presence of still longer bonds compared to the model without Pd ion doping. The synergistic effect is therefore due to enhanced reducibility of both dopant and host ion induced due to structural distortion of fluorite lattice in presence of dopant ion. For RE ions (RE = Y, La), our calculations show very little deviation of bonds lengths from ideal fluorite structure. The absence of longer Y-O/La-O and Ce-O bonds make the structure much less susceptible to reduction.

Shell-model calculation of defect energies in alkali halides employing crystal-independent interionic potential parameters

Shell model calculation of defect energies in alkali halides have been carried out using the ion-dependent, crystal-independent potential parameters of Sangster and Atwood (1978). Results indicate that appreciable differences exist between barrier heights for migration of cations and anions. While barrier heights for cations are generally lower than for anions in alkali halides of NaCl structure, the opposite is true in alkali halides of CsCl structure.

Theoretical Treatment of Interaction Parameters in Multicomponent Metallic Solutions

A new theoretical equation for interaction parameter in multicomponent metallic solutions is developed using the pseudopotential formalism coupled with the free energy of the hard sphere system. The approximate expression for the pseudopotential term is given in terms of the heat of solution at infinite dilution, to allow easy evaluation of the interaction parameter in various multicomponent systems. This theory has been applied to 23 non-ferrous alloys based on Pb, Sn, Bi and indium. Comparison with the results of previous theoretical calculations using only the hard sphere model suggests that the inclusion of the pseudopotential term yields a quantitatively more correct prediction of interaction parameters in multicomponent metallic solutions. Numerical calculations were also made for 320 Fe-base solutions relevant to steelmaking and the agreement between calculation and experimental data appears reasonable, with 90% reliability in predicting the correct sign.

A numerical survey of relativistic rotating neutron star structures using the Hartle-Thorne formalism

A broad numerical survey of relativistic rotating neutron star structures was compiled using an exhaustive list of presently available equation of state models for neutron star matter. The structure parameters (spherical deformations in mass and radii, the moment of inertia and quadrupole moment, oblateness, and free precession) are calculated using the formalism proposed by Hartle and Thorne (1968). The results are discussed in relation to the relevant observational information. Binary pulsar data and X-ray burst sources provide information on the bulk properties of neutron stars, enabling the derivation of constraints that can be put on the structure of neutron stars and equation of state models.

Variational calculation for the spin-(1/2 Heisenberg antiferromagnet on a square lattice

The ground-state properties of the spin-(1/2 Heisenberg antiferromagnet on a square lattice are studied by using a simple variational wave function that interpolates continuously between the Néel state and short-range resonating-valence-bond states. Exact calculations of the variational energy for small systems show that the state with the lowest energy has long-range antiferromagnetic order. The staggered magnetization in this state is approximately 70% of its maximum possible value. The variational estimate of the ground-state energy is substantially lower than the value obtained for the nearest-neighbor resonating-valence-bond wave function.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.