Method of invariant grid for model reduction of hydrogen combustion

The Method of Invariant Grid (MIG) is a model reduction technique based on the concept of slow invariant manifold (SIM). The MIG approximates the SIM by a set of nodes in the concentration space (invariant grid). In the present work, the MIG is applied to a realistic combustion system: an adiabatic constant volume reactor with H₂-air at stoichiometric proportions. By considering the thermodynamic Lyapunov function of the detailed kinetic system, the notion of the quasi-equilibrium manifold (QEM) is adopted as an initial approximation to the SIM. One- and two-dimensional discrete approximations of the QEM (quasi-equilibrium grids) are constructed and refined via MIG to obtain the corresponding invariant grids. The invariant grids are tabulated and used to integrate the reduced system. Excellent agreement between the reduced and detailed kinetics is demonstrated.     

Slow Motion and Metastability for a Nonlocal Evolution Equation

In this paper we consider a nonlocal evolution equation in one dimension, which describes the dynamics of a ferromagnetic system in the mean field approximation. In the presence of a small magnetic field, it admits two stationary and homogeneous solutions, representing the stable and metastable phases of the physical system. We prove the existence of an invariant, one dimensional manifold connecting the stable and metastable phases. This is the unstable manifold of a distinguished, spatially nonhomogeneous, stationary solution, called the critical droplet.^{(4, 10)} We show that the points on the manifold are droplets longer or shorter than the critical one, and that their motion is very slow in agreement with the theory of metastable patterns. We also obtain a new proof of the existence of the critical droplet, which is supplied with a local uniqueness result.     

An adaptive method for computing invariant manifolds in non-autonomous, three-dimensional dynamical systems

We present a computational method for determining the geometry of a class of three-dimensional invariant manifolds in non-autonomous (aperiodically time-dependent) dynamical systems. The presented approach can be also applied to analyse the geometry of 3D invariant manifolds in three-dimensional, time-dependent fluid flows. The invariance property of such manifolds requires that, at any fixed time, they are given by surfaces in R³. We focus on a class of manifolds whose instantaneous geometry is given by orientable surfaces embedded in R³. The presented technique can be employed, in particular, to compute codimension one (invariant) stable and unstable manifolds of hyperbolic trajectories in 3D non-autonomous dynamical systems which are crucial in the Lagrangian transport analysis. The same approach can also be used to determine evolution of an orientable ‘material surface’ in a fluid flow. These developments represent the first step towards a non-trivial 3D extension of the so-called lobe dynamics — a geometric, invariant-manifold-based framework which has been very successful in the analysis of Lagrangian transport in unsteady, two-dimensional fluid flows. In the developed algorithm, the instantaneous geometry of an invariant manifold is represented by an adaptively evolving triangular mesh with piecewise C² interpolating functions. The method employs an automatic mesh refinement which is coupled with adaptive vertex redistribution. A variant of the advancing front technique is used for remeshing, whenever necessary. Such an approach allows for computationally efficient determination of highly convoluted, evolving geometry of codimension one invariant manifolds in unsteady three-dimensional flows. We show that the developed method is capable of providing detailed information on the evolving Lagrangian flow structure in three dimensions over long periods of time, which is crucial for a meaningful 3D transport analysis.     

The application of chemical reduction methods to a combustion system exhibiting complex dynamics

A number of chemical model reduction techniques have been developed over recent years with a growing range of applications in combustion. The following work demonstrates the application of such reduction techniques for a combustion system describing the oxidation of carbon monoxide + hydrogen in a continuously stirred tank reactor (CSTR) at very low pressure. The system exhibits complex dynamics including oscillatory glow, oscillatory ignition and mixed mode oscillations. It is demonstrated that a range of local reduction methods can be applied to such complex systems, as long as sufficient coverage of the accessed regions of phase space are included in the reduction analysis. The methods include sensitivity analysis, the quasi-steady state approximation (QSSA) and repro-modelling based on the concept of an intrinsic low dimensional manifold (ILDM). The system is qualitatively different from some previous applications of ILDM methods where trajectories tend towards a fixed equilibrium. The underlying dimension of the system is seen to vary throughout selected trajectories with rapid increases occurring over very short time-scales during oscillatory ignition. Nevertheless, a final reduced model of only four variables is developed using fitted orthonormal polynomials describing the system dynamics on a slow manifold. The application serves to demonstrate that the relationship between local reduced model error and global errors can be complex for systems exhibiting complex dynamics, with regions of seemingly small local mapping gradients requiring tighter error control in order to control global errors. This feature may be common in cases where nearby trajectories are seen to diverge within the slow manifold over time.     

Minimal curvature trajectories: Riemannian geometry concepts for slow manifold computation in chemical kinetics

In dissipative ordinary differential equation systems different time scales cause anisotropic phase volume contraction along solution trajectories. Model reduction methods exploit this for simplifying chemical kinetics via a time scale separation into fast and slow modes. The aim is to approximate the system dynamics with a dimension-reduced model after eliminating the fast modes by enslaving them to the slow ones via computation of a slow attracting manifold. We present a novel method for computing approximations of such manifolds using trajectory-based optimization. We discuss Riemannian geometry concepts as a basis for suitable optimization criteria characterizing trajectories near slow attracting manifolds and thus provide insight into fundamental geometric properties of multiple time scale chemical kinetics. The optimization criteria correspond to a suitable mathematical formulation of “minimal relaxation” of chemical forces along reaction trajectories under given constraints. We present various geometrically motivated criteria and the results of their application to four test case reaction mechanisms serving as examples. We demonstrate that accurate numerical approximations of slow invariant manifolds can be obtained.     

f(R,L m ) gravity

We generalize the f(R) type gravity models by assuming that the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R and of the matter Lagrangian L _m . We obtain the gravitational field equations in the metric formalism, as well as the equations of motion for test particles, which follow from the covariant divergence of the energy-momentum tensor. The equations of motion for test particles can also be derived from a variational principle in the particular case in which the Lagrangian density of the matter is an arbitrary function of the energy density of the matter only. Generally, the motion is non-geodesic, and it takes place in the presence of an extra force orthogonal to the four-velocity. The Newtonian limit of the equation of motion is also considered, and a procedure for obtaining the energy-momentum tensor of the matter is presented. The gravitational field equations and the equations of motion for a particular model in which the action of the gravitational field has an exponential dependence on the standard general relativistic Hilbert–Einstein Lagrange density are also derived.     

One-Dimensional Structures Behind Twisted and Untwisted SuperYang–Mills Theory

We give a one-dimensional interpretation of the four-dimensional twisted N = 1 superYang–Mills theory on a Kähler manifold by performing an appropriate dimensional reduction. We prove the existence of a 6-generator superalgebra, which does not possess any invariant Lagrangian but contains two different subalgebras that determine the twisted and untwisted formulations of the N = 1 superYang–Mills theory.     

基于Foliation条件的离散动力系统二维流形计算

研究离散动力系统双曲不动点的二维流形计算,利用不变流形轨道上Jacobian矩阵能够传递导数这一特殊性质,提出一种新的一维流形计算方法,通过预测-校正两个步骤迅速确定流形上新网格点,避免重复计算,并简化精度控制条件.在此基础上,将基于流形面Foliation条件进行推广,推广后的Foliation条件能够控制二维流形上的一维子流形的增长速度,从而实现二维流形在各个方向上的均匀增长.此外,算法可以同时用于二维稳定和不稳定流形的计算.以超混沌三维Hénon映射和具有蝶形吸引子的Lorenz系统为例验证了算法的有效性.     

Probability distribution for the number of cycles between successive regime transitions for the Lorenz model

The Lorenz model has been widely used for exploring many real world problems. In this paper we obtain, with the help of an invariant manifold technique, the return map for the maximum value of the variable x of the model and use this return map to derive the simple, empirically obtained, regime transition rules for forecasting regime changes and length in the new regime for the model. The probability distribution for number of cycles between successive regime transitions of the Lorenz model may be of interest in many disciplines. We apply the Perron-Frobenius algorithm over the return map to estimate the probability distribution for the number of cycles between successive regime transitions. These probabilities are also estimated for the forced Lorenz model, which is a conceptual model to explore the effects of sea surface temperature on seasonal rainfall.     

Quantum Gravity Effect on Neutrino Oscillation

The quantum gravity may have strong consequence for neutrino oscillation phenemomenon over a large distance.We found a significant modification of neutrino oscillation due to quantum gravity effects. Quantum gravity (Planck scale effects) leads to an effective S U(2) _L ×U(1) invariant dimension-5 Lagrangian involving, neutrino and Higgs fields. On symmetry breaking, this operator gives rise to correction to the neutrino masses and mixing. The gravitational interaction (M _X =M _{p l} ) demands that the element of this perturbation matrix should be independent of flavor indices. In this paper, we study the quantum gravity effects on neutrino oscillation, namely modified dispersion relation for neutrino oscillations parameter.     

Where do inertial particles go in fluid flows?

We derive a general reduced-order equation for the asymptotic motion of finite-size particles in unsteady fluid flows. Our inertial equation is a small perturbation of passive fluid advection on a globally attracting slow manifold. Among other things, the inertial equation implies that particle clustering locations in two-dimensional steady flows can be described rigorously by the Q parameter, i.e., by one-half of the squared difference of the vorticity and the rate of strain. Use of the inertial equation also enables us to solve the numerically ill-posed problem of source inversion, i.e., locating initial positions from a current particle distribution. We illustrate these results on inertial particle motion in the Jung-Tél-Ziemniak model of vortex shedding behind a cylinder in crossflow.     

Quantum Gravity on Neutrino Oscillation Length

The oscillation length of neutrino oscillation could be discussed in the frame work of quantum gravity. Quantum gravity (Planck scale effects) leads to an effective SU(2) _L ×U(1) invariant dimension-5 Lagrangian involving, neutrino and Higgs fields. On symmetry breaking, this operator gives rise to correction to the neutrino masses and mixing. We compute the neutrino oscillation length due to Planck scale effects. The gravitational interaction (M _X =M _pl ) demands that the element of this perturbation matrix should be independent of flavor indices. In this paper, we study the quantum gravity effects on neutrino oscillation length, namely modified dispersion relation for neutrino oscillation phases.     

Neutrino Oscillation Phase Shift from Quantum Gravity

The phase shift of neutrino oscillation could be discussed in the frame work of quantum gravity. Quantum gravity (Planck scale effects) leads to an effective SU(2) _L ×U(1) invariant dimension-5 Lagrangian involving, neutrino and Higgs fields. On symmetry breaking, this operator gives rise to correction to the neutrino masses and mixing. We compute the neutrino oscillation phase due to Planck scale effects. The gravitational interaction (M _X =M _pl ) demands that the element of this perturbation matrix should be independent of flavor indices. In this paper, we study the quantum gravity effects on neutrino oscillation phases, namely modified dispersion relation for neutrino oscillation phases.     

Classical mechanics on the GL(n, ℝ) group and Euler-Calogero-Sutherland model

Relations between free motion on the GL ⁺(n, ?) group manifold and the dynamics of an n-particle system with spin degrees of freedom on a line interacting with a pairwise 1/sinh² x “potential” (Euler-Calogero-Sutherland model) are discussed within a Hamiltonian reduction. Two kinds of reductions of the degrees of freedom are considered: that which is due to continuous invariance and that which is due to discrete symmetry. It is shown that, upon projecting onto the corresponding invariant manifolds, the resulting Hamiltonian system represents the Euler-Calogero-Sutherland model in both cases.     

Effect of Majorana Phases in Neutrino Oscillation

Quantum gravity (Planck scale effects) lead to an effective SU(2) _L ×U(1) invariant dimension-5 Lagrangian involving neutrino and Higgs fields. On symmetry breaking, this operator gives rise to correction to the above masses and mixing. The gravitational interaction M _X =M _pl , we find that for degenerate neutrino mass spectrum, it is shown that the Majorana phase of the neutrino mixing matrix can effects in neutrino oscillation probability.     

Renormalization-Group Method for Reduction of Evolution Equations; Invariant Manifolds and Envelopes

The renormalization group (RG) method as a powerful tool for reduction of evolution equations is formulated in terms of the notion of invariant manifolds. We start with derivation of an exact RG equation which is analogous to the Wilsonian RG equations in statistical physics and quantum field theory. It is clarified that the perturbative RG method constructs invariant manifolds successively as the initial value of evolution equations, thereby the meaning to set t₀=t is naturally understood where t₀ is the arbitrary initial time. We show that the integral constants in the unperturbative solution constitutes natural coordinates of the invariant manifold when the linear operator A in the evolution equation is semi-simple, i.e., diagonalizable; when A is not semi-simple and has a Jordan cell, a slight modification is necessary because the dimension of the invariant manifold is increased by the perturbation. The RG equation determines the slow motion of the would-be integral constants in the unperturbative solution on the invariant manifold. We present the mechanical procedure to construct the perturbative solutions hence the initial values with which the RG equation gives meaningful results. The underlying structure of the reduction by the RG method as formulated in the present work turns out to completely fit to the universal one elucidated by Kuramoto some years ago. We indicate that the reduction procedure of evolution equations has a good correspondence with the renormalization procedure in quantum field theory; the counter part of the universal structure of reduction elucidated by Kuramoto may be Polchinski's theorem for renormalizable field theories. We apply the method to interface dynamics such as kink–anti-kink and soliton–soliton interactions in the latter of which a linear operator having a Jordan-cell structure appears.     

Reduction of Vaisman structures in complex and quaternionic geometry

We consider locally conformal Kähler geometry as an equivariant (homothetic) Kähler geometry: a locally conformal Kähler manifold is, up to equivalence, a pair $(K,\Gamma )$, where $K$ is a Kähler manifold and $\Gamma$ is a discrete Lie group of biholomorphic homotheties acting freely and properly discontinuously. We define a new invariant of a locally conformal Kähler manifold $(K,\Gamma )$ as the rank of a natural quotient of $\Gamma$, and prove its invariance under reduction. This equivariant point of view leads to a proof that locally conformal Kähler reduction of compact Vaisman manifolds produces Vaisman manifolds and is equivalent to a Sasakian reduction. Moreover, we define locally conformal hyperKähler reduction as an equivariant version of hyperKähler reduction and in the compact case we show its equivalence with 3-Sasakian reduction. Finally, we show that locally conformal hyperKähler reduction induces hyperKähler with torsion (HKT) reduction of the associated HKT structure and the two reductions are compatible, even though not every HKT reduction comes from a locally conformal hyperKähler reduction.     

Invariant manifold methods for metabolic model reduction

After the decay of transients, the behavior of a set of differential equations modeling a chemical or biochemical system generally rests on a low-dimensional surface which is an invariant manifold of the flow. If an equation for such a manifold can be obtained, the model has effectively been reduced to a smaller system of differential equations. Using perturbation methods, we show that the distinction between rapidly decaying and long-lived (slow) modes has a rigorous basis. We show how equations for attracting invariant (slow) manifolds can be constructed by a geometric approach based on functional equations derived directly from the differential equations. We apply these methods to two simple metabolic models. (c) 2001 American Institute of Physics.     

Gauged Kaluza-Klein AdS pseudo-supergravity

We obtain the pseudo-supergravity extension of the D-dimensional Kaluza-Klein theory, which is the circle reduction of pure gravity in D+1 dimensions. The fermionic partners are pseudo-gravitino and pseudo-dilatino. The full Lagrangian is invariant under the pseudo-supersymmetric transformation, up to quadratic order in fermion fields. We find that the theory possesses a U(1) global symmetry that can be gauged so that all the fermions are charged under the Kaluza-Klein vector. The gauging process generates a scalar potential that has a maximum, leading to the AdS vacuum. Whist the highest dimension for gauged AdS supergravity is seven, our gauged AdS pseudo-supergravities can exist in arbitrary dimensions.