Assessment of conditional moment closure models of turbulent autoignition using DNS data

A systematic assessment of the submodels of conditional moment closure (CMC) formalism for the autoignition problem is carried out using direct numerical simulation (DNS) data. An initially non-premixed, n-heptane/air system, subjected to a three-dimensional, homogeneous, isotropic, and decaying turbulence, is considered. Two kinetic schemes, (1) a one-step and (2) a reduced four-step reaction mechanism, are considered for chemistry An alternative formulation is developed for closure of the mean chemical source term , based on the condition that the instantaneous fluctuation of excess temperature is small. With this model, it is shown that the CMC equations describe the autoignition process all the way up to near the equilibrium limit. The effect of second-order terms (namely, conditional variance of temperature excess sigma(2) and conditional correlations of species q(ij)) in modeling is examined. Comparison with DNS data shows that sigma(2) has little effect on the predicted conditional mean temperature evolution, if the average conditional scalar dissipation rate is properly modeled. Using DNS data, a correction factor is introduced in the modeling of nonlinear terms to include the effect of species fluctuations. Computations including such a correction factor show that the species conditional correlations q(ij) have little effect on model predictions with a one-step reaction, but those q(ij) involving intermediate species are found to be crucial when four-step reduced kinetics is considered. The "most reactive mixture fraction" is found to vary with time when a four-step kinetics is considered. First-order CMC results are found to be qualitatively wrong if the conditional mean scalar dissipation rate is not modeled properly. The autoignition delay time predicted by the CMC model compares excellently with DNS results and shows a trend similar to experimental data over a range of initial temperatures.

Autoignition in a non-premixed medium: DNS studies on the effects of three-dimensional turbulence

Direct numerical simulation (DNS) results of autoignition in anon-premixed medium under an isotropic, homogeneous, and decaying turbulence are presented. The initial mixture consists of segregated fuel parcels randomly distributed within warm air, and the entire medium is subjected to a three-dimensional turbulence. Chemical kinetics is modeled by a four-step reduced reaction mechanism for autoignition of n-heptane/air mixture. Thus, this work overcomes the principal limitations of a previous contribution of the authors on two-dimensional DNS of autoignition with a one-step reaction model. Specific attention is focused on the differences in the effects of two- and three-dimensional turbulence on autoignition characteristics. The three-dimensional results show that ignition spots are most likely to originate at locations jointly corresponding to the most reactive mixture fraction and low scalar dissipation rate. Further, these ignition spots are found to originate at locations corresponding to the core of local vortical structures, and after ignition, the burning gases move toward the vortex periphery Such a movement is explained as caused by the cyclostrophic imbalance developed when the local gas density is variable. These results lead to the conclusion that the local ignition-zone structure does not conform to the classical stretched flamelet description. Parametric studies show that the ignition delay time decreases with an increase in turbulence intensity. Hence, these three-dimensional simulation results resolve the discrepancy between trends in experimental data and predictions from DNSs of two-dimensional turbulence. This qualitative difference between DNS results from three- and two-dimensional simulations is discussed and attributed to the effect of vortex stretching that is present in the former, but not in the latter.

Surfaces of Bounded Mean Curvature In Riemannian Manifolds

Consider a sequence of closed, orientable surfaces of fixed genus g in a Riemannian manifold M with uniform upper bounds on the norm of mean curvature and area. We show that on passing to a subsequence, we can choose parametrisations of the surfaces by inclusion maps from a fixed surface of the same genus so that the distance functions corresponding to the pullback metrics converge to a pseudo-metric and the inclusion maps converge to a Lipschitz map. We show further that the limiting pseudo-metric has fractal dimension two. As a corollary, we obtain a purely geometric result. Namely, we show that bounds on the mean curvature, area and genus of a surface F subset of M, together with bounds on the geometry of M, give an upper bound on the diameter of F. Our proof is modelled on Gromov's compactness theorem for J-holomorphic curves.

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic-perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near-tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T-stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T-stress becomes more negative irrespective of lattice orientation. Also, the near-tip field for a range of constraint levels can be characterized by two-parameters such as K-T or J-Q as in isotropic plastic solids.

Almost quarter-pinched Kähler metrics and Chern numbers

Given n is an element of Z(+) and epsilon > 0, we prove that there exists delta = delta(epsilon, n) > 0 such that the following holds: If (M(n),g) is a compact Kahler n-manifold whose sectional curvatures K satisfy -1 -delta <= K <= -1/4 and c(I)(M), c(J)(M) are any two Chern numbers of M, then |c(I)(M)/c(J)(M) - c(I)(0)/c(J)(0)| < epsilon, where c(I)(0), c(J)(0) are the corresponding characteristic numbers of a complex hyperbolic space form. It follows that the Mostow-Siu surfaces and the threefolds of Deraux do not admit Kahler metrics with pinching close to 1/4.

A note on surface water waves for finite depth in the presence of an ice-cover

A class of I boundary value problems involving propagation of two-dimensional surface water waves, associated with water of uniform finite depth, against a plane vertical wave maker is investigated under the assumption that the surface is covered by a thin sheet of ice. It is assumed that the ice-cover behaves like a thin isotropic elastic plate. Then the problems under consideration lead to those of solving the two-dimensional Laplace equation in a semi-infinite strip, under Neumann boundary conditions on the vertical boundary as well as on one of the horizontal boundaries, representing the bottom of the fluid region, and a condition involving upto fifth order derivatives of the unknown function on the top horizontal ice-covered boundary, along with the two appropriate edge-conditions, at the ice-covered corner, ensuring the uniqueness of the solutions. The mixed boundary value problems are solved completely, by exploiting the regularity property of the Fourier cosine transform.

The effect of T-stress on plane strain dynamic crack growth in elastic-plastic materials

In this paper, the effects of T -stress on steady, dynamic crack growth in an elastic-plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J (2) flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T -stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T -stress for both ductile and cleavage mode of crack growth.

Geometrically Nonlinear Dynamics of Composite Four-Bar Mechanisms

This work intends to demonstrate the importance of geometrically nonlinear crosssectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically nonlinear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and nonlinear 1-D analyses along the four beam reference curves. For thin rectangular cross-sections considered here, the 2-D cross-sectional nonlinearity is overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the nonlinear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the nonlinear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses, more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the nonlinear, flexible fourbar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we shall attempt to identify and investigate a few problems where the cross-sectional nonlinearities are significant. This will be carried out by varying stacking sequences and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form nonlinear beam stiffness matrix. Numerical examples will be presented and results from this analysis will be compared with those available in the literature, for linear cross-sectional analysis and isotropic materials as special cases.

Generalized asymptotic expansions for the wavenumbers in infinite flexible in vacuo orthotropic cylindrical shells

Analytical expressions are found for the wavenumbers and resonance frequencies in flexible, orthotropic shells using the asymptotic methods. These expressions are valid for arbitrary circumferential orders n. The Donnell-Mushtari shell theory is used to model the dynamics of the cylindrical shell. Initially, an in vacuo cylindrical isotropic shell is considered and expressions for all the wavenumbers (bending, near-field bending, longitudinal and torsional) are found. Subsequently, defining a suitable orthotropy parameter epsilon, the problem of wave propagation in an orthotropic shell is posed as a perturbation on the corresponding problem for an isotropic shell. Asymptotic expressions for the wavenumbers in the in vacuo orthotropic shell are then obtained by treating epsilon as an expansion parameter. In both cases (isotropy and orthotropy), a frequency-scaling parameter (eta) and Poisson's ratio (nu) are used to find elegant expansions in the different frequency regimes. The asymptotic expansions are compared with numerical solutions in each of the cases and the match is found to be good. The main contribution of this work lies in the extension of the existing literature by developing closed-form expressions for wavenumbers with arbitrary circumferential orders n in the case of both, isotropic and orthotropic shells. Finally, we present natural frequency expressions in finite shells (isotropic and orthotropic) for the axisymmetric mode and compare them with numerical and ANSYS results. Here also, the comparison is found to be good. (C) 2011 Elsevier Ltd. All rights reserved.

A family of Runge-Kutta based explicit methods for rotational dynamics

The present paper develops a family of explicit algorithms for rotational dynamics and presents their comparison with several existing methods. For rotational motion the configuration space is a non-linear manifold, not a Euclidean vector space. As a consequence the rotation vector and its time derivatives correspond to different tangent spaces of rotation manifold at different time instants. This renders the usual integration algorithms for Euclidean space inapplicable for rotation. In the present algorithms this problem is circumvented by relating the equation of motion to a particular tangent space. It has been accomplished with the help of already existing relation between rotation increments which belongs to two different tangent spaces. The suggested method could in principle make any integration algorithm on Euclidean space, applicable to rotation. However, the present paper is restricted only within explicit Runge-Kutta enabled to handle rotation. The algorithms developed here are explicit and hence computationally cheaper than implicit methods. Moreover, they appear to have much higher local accuracy and hence accurate in predicting any constants of motion for reasonably longer time. The numerical results for solutions as well as constants of motion, indicate superior performance by most of our algorithms, when compared to some of the currently known algorithms, namely ALGO-C1, STW, LIEMID[EA], MCG, SUBCYC-M.

An algebraic formulation of exact force-, moment-isotropy in spatial parallel manipulators

In this paper, we present an algebraic method to study and design spatial parallel manipulators that demonstrate isotropy in the force and moment distributions.We use the force and moment transformation matrices separately,and derive conditions for their isotropy individually as well as in combination. The isotropy conditions are derived in closed-form in terms of the invariants of the quadratic forms associated with these matrices. The formulation has been applied to a class of Stewart platform manipulators. We obtain multi-parameter families of isotropic manipulator analytically. In addition to computing the isotropic configurations of an existing manipulator,we demonstrate a procedure for designing the manipulator for isotropy at a given configuration.

Anomalous temperature dependence of phonons and photoluminescence bands in pyrochlore Er(2)Ti(2)O(7): signatures of structural deformation at 130 K

We report a Raman study of single crystal pyrochlore Er(2)Ti(2)O(7) as a function of temperature from 12 to 300 K. In addition to the phonons, various photoluminescence (PL) lines of Er(3+) in the visible range are also observed. Our Raman data show an anomalous red-shift of two phonons (one at similar to 200 cm(-1) and another at similar to 520 cm(-1)) upon cooling from room temperature which is attributed to phonon-phonon anharmonic interactions. However, the phonons at similar to 310, 330, and 690 cm(-1) initially show a blue-shift upon cooling from room temperature down to about 130 K, followed by a red-shift, indicating a structural deformation at similar to 130 K. The intensities of the PL bands associated with the transitions between the various levels of the ground state manifold ((4)I(15/2)) and the (2)H(11/2) as well as (4)S(3/2) excited state manifolds of Er(3+) show a change at similar to 130 K. Moreover, the temperature dependence of the peak position of the two PL bands shows a change in their slope (d(omega)/d(T)) at similar to 130 K, thus further strengthening the proposal of a structural deformation. The temperature dependence of the peak positions of the PL bands has been analyzed using the theory of optical dephasing in crystals.

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.

Development of novel multiple quantummethodologies for the analyses of complex protonNMR spectra of scalar coupled spins

One of the significant advancements in Nuclear Magnetic Resonance spectroscopy (NMR) in combating the problem of spectral complexity for deriving the structure and conformational information is the incorporation of additional dimension and to spread the information content in a two dimensional space. This approach together with the manipulation of the dynamics of nuclear spins permitted the designing of appropriate pulse sequences leading to the evolution of diverse multidimensional NMR experiments. The desired spectral information can now be extracted in a simplified and an orchestrated manner. The indirect detection of multiple quantum (MQ) NMR frequencies is a step in this direction. The MQ technique has been extensively used in the study of molecules aligned in liquid crystalline media to reduce spectral complexity and to determine molecular geometries. Unlike in dipolar coupled systems, the size of the network of scalar coupled spins is not big in isotropic solutions and the MQ 1H detection is not routinely employed,although there are specific examples of spin topology filtering. In this brief review, we discuss our recent studies on the development and application of multiple quantum correlation and resolved techniques for the analyses of proton NMR spectra of scalar coupled spins.

Molecular Dynamics of Thermotropic Liquid Crystals: Anomalous relaxation dynamics of calamitic and discotic liquid crystals

Recent optical kerr effect (OKE) studies have demonstrated that orientational relaxation of rod-like nematogens exhibits temporal power law decay at intermediate times not only near the isotropic–nematic (I–N) phase boundary but also in the nematic phase. Such behaviour has drawn an intriguing analogy with supercooled liquids. We have investigated both collective and single-particle orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I–N phase boundary, the system behaves remarkably like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system. A model liquid crystal, consisting of disk-like molecules, has also been investigated in molecular dynamics simulations for orientational relaxation along two isobars starting from the high temperature isotropic phase. The isobars have been so chosen that the phase sequence isotropic (I)–nematic (N)–columnar (C) appears upon cooling along one of them and the sequence isotropic (I)–columnar(C) along the other. While the orientational relaxation in the isotropic phase near the I–N phase transition shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I–C phase boundary. The origin of the power law decay in the single-particle second-rank orientational time correlation function (OTCF) is traced to the growth of the orientational pair distribution functions near the I–N phase boundary. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational OTCF follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.