Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear manifolds”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Self- and Air-Broadened Line Shapes in the 2v3 P and R Branches of 12CH4

In this paper we report line shape parameters of 12CH4 for several hundred 2V(sub 3) transitions in the spectral regions 5891-5996 cm( exp -1) (P branch) and 6015-6115 cm(exp -1) (R branch). Air- and self-broadening coefficients were measured as a function of temperature; line mixing via off-diagonal relaxation matrix element coefficients was also obtained for 47 transition pairs. In total, nearly 1517 positions and intensities were retrieved, but many transitions were too weak for the line shape study. For this analysis, we used 25 high-resolution (0.0056 and 0.0067 cm(ex[ -1) and high signal-to-noise (S/N) spectra of high-purity 12CH4 and the same high-purity 12CH4 broadened by dry air recorded at different sample temperatures between 130 K and 295 K with the Bruker IFS 125HR Fourier transform spectrometer at JPL. Three different absorption cells were used (1) a White cell set to a path length of 13.09 m for room temperature data, (2) a single-pass 0.2038 m long coolable cell (for self-broadening) and (3) a multipass cell with 20.941 m total path coolable Herriott cell (for air-broadening). In total there were 13 spectra with pure 12CH4 (0.27-599 Torr) and 12 air-broadened spectra with total sample pressures of 80-805 Torr and volume mixing ratios (VMR) of methane between 0.18 and 1.0. An interactive multispectrum nonlinear least-squares technique was employed to fit the individual P10-P1 and R0-R10 manifolds in all the spectra simultaneously. Results obtained from the present analysis are compared to other recent measurements.

Devi, V. Malathy↗

Exploiting Power Flow Manifold to Solve AC Optimal Power Flow

AC optimal power flow has proven difficult to solve with interior point methods on GPUs. This is largely due to challenging linear algebra problems that current state of the art massively parallel linear solvers struggle with. However, the advent of Riemannian optimization techniques and the fact that the power flow equations form a smooth manifold present an alternative approach. In this talk, we present the basics of Riemannian optimization techniques in which optimization is done directly on a manifold. Then we present computational results showing that Riemannian techniques are capable of producing solutions of comparable quality as interior point methods.

AC optimal power flow↗

Topological analysis of chaotic orbits: Revisiting Hyperion

There is emerging interest in the possibility of chaotic evolution in astrophysical systems. To mention just one example, recent well-sampled ground-based observations of the Saturian satellite Hyperion strongly suggest that it is exhibiting chaotic behavior. We present a general technique, the method of close returns, for the analysis of data from astronomical objects believed to be exhibiting chaotic motion. The method is based on the extraction of pieces of the evolution that exhibit nearly periodic behavior-episodes during which the object stays near in phase space to some unstable periodic orbit. Such orbits generally act as skeletal features, tracing the topological organization of the manifold on which the chaotic dynamics takes place. This method does not require data sets as lengthy as other nonlinear analysis techniques do and is therefore well suited to many astronomical observing programs. Well sampled data covering between twenty and forty characteristic periods of the system have been found to be sufficient for the application of this technique. Additional strengths of this method are its robustness in the presence of noise and the ability for a user to clearly distinguish between periodic, random, and chaotic behavior by inspection of the resulting two-dimensional image. As an example of its power, we analyze close returns in a numerically generated data set, based on a model for Hyperion extensively studied in the literature, corresponding to nightly observations of the satellite. We show that with a small data set, embedded unstable periodic orbits can be extracted and that these orbits can be responsible for nearly periodic behavior lasting a substantial fraction of the observing run.

Boyd, Patricia T.↗

ssys

`ssys` is a Python toolkit for exact algebraic transformation of ordinary differential equation (ODE) models into canonical S-system or Generalized Mass Action (GMA) form. Given a model in Antimony or SBML format, `ssys` produces a mathematically equivalent representation. The transformation introduces auxiliary variables as needed to decompose a broad class of nonlinearities into products of power-law terms. The recast is exact: the original and transformed systems have identical dynamics on the invariant constraint manifold defined by auxiliary variable definitions, given consistent initial conditions.

Hlavacek, William [Los Alamos National Laboratory]↗

On-shell correlators and color-kinematics duality in curved symmetric spacetimes

We define a perturbatively calculable quantity — the on-shell correlator — which furnishes a unified description of particle dynamics in curved spacetime. Specializing to the case of flat and anti-de Sitter space, on-shell correlators coincide precisely with on-shell scattering amplitudes and boundary correlators, respectively. Remarkably, we find that symmetric manifolds admit a generalization of on-shell kinematics in which the corresponding momenta are literally the isometry generators of the spacetime acting on the external kinematic data. These isometric momenta are intrinsically non-commutative but exhibit on-shell conditions that are identical to those of flat space, thus providing a common language for computing and representing on-shell correlators which is agnostic about the underlying geometry. Afterwards, we compute tree-level on-shell correlators for biadjoint scalar (BAS) theory and the nonlinear sigma model (NLSM) and learn that color-kinematics duality is manifested at the level of fields under a mapping of the color algebra to the algebra of gauged isometries on the spacetime manifold. Last but not least, we present a field theoretic derivation of the fundamental BCJ relations for on-shell correlators following from the existence of certain conserved currents in BAS theory and the NLSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Differential Flatness and Cooperative Tracking in the Lorenz System

In this paper the control of the Lorenz system for both stabilization and tracking problems is studied via feedback linearization and differential flatness. By using the Rayleigh number as the control, only variable physically tunable, a barrier in the controllability of the system is incidentally imposed. This is reflected in the appearance of a singularity in the state transformation. Composite controllers that overcome this difficulty are designed and evaluated. The transition through the manifold defined by such a singularity is achieved by inducing a chaotic response within a boundary layer that contains it. Outside this region, a conventional feedback nonlinear control is applied. In this fashion, the authority of the control is enlarged to the whole. state space and the need for high control efforts is mitigated. In addition, the differential parametrization of the problem is used to track nonlinear functions of one state variable (single tracking) as well as several state variables (cooperative tracking). Control tasks that lead to integrable and non-integrable differential equations for the nominal flat output in steady-state are considered. In particular, a novel numerical strategy to deal with the non-integrable case is proposed. Numerical results validate very well the control design.

Crespo, Luis G.↗

Data-Driven Modeling and Correction of Vehicle Dynamics

We develop a data-driven framework for learning and correcting nonautonomous vehicle dynamics. Physics-based vehicle models are often simplified for tractability and therefore exhibit inherent model-form uncertainty, motivating the need for data-driven correction. Moreover, nonautonomous dynamics are governed by time-dependent control inputs, which pose challenges in learning predictive models directly from temporal snapshot data. To address these, we reformulate the vehicle dynamics via a local parameterization of the time-dependent inputs, yielding a modified system composed ofa sequence of local parametric dynamical systems. Here, we approximate these parametric systems using two complementary approaches. First, we employ the dimension reduction and interpolation in parameter space (DRIPS) methodology to construct efficient linear surrogate models, equipped with lifted observable spaces and manifold-based operator interpolation. This enables data-efficient learning of vehicle models whose dynamics admit accurate linear representations in the lifted spaces. Second, for more strongly nonlinear systems, we employ flow map learning (FML), a deep neural network (DNN) approach that approximates the parametric evolution map without requiring special treatment of nonlinearities. We further extend FML with a transfer-learning-based model correction procedure, enabling the correction of misspecified prior models using only a sparse set of high-fidelity or experimental measurements, without assuming a prescribed form for the correction term. Through a suite of numerical experiments on unicycle, simplified bicycle, and slip-based bicycle models, we demonstrate that DRIPS offers robust and highly data-efficient learning of nonautonomous vehicle dynamics, while FML provides expressive nonlinear modeling and effective correction of model-form errors under severe data scarcity.

data-driven modeling↗

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗

Super-resonant four-photon collinear laser frequency multiplication in plasma

Resonant four-photon scattering could nearly double frequencies of intense laser pulses in plasma. However, transverse slippage between pulses presents a technological challenge, while collinear four-photon scattering is forbidden for classical light dispersion in plasma. Nonlinear renormalization of intense laser pulses can enable collinear four-photon resonance. However, such a very intensity-sensitive resonance is difficult to maintain for evolving pulses. Remarkably, there is a lower-dimensionality submanifold of the resonant four-photon manifold where the evolving pulses stay in resonance. Furthermore this could enable an all-optical frequency doubling of mildly relativistic-intense laser pulses in collinear geometry, advantageously free of the transverse slippage challenges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perturbed nonlinear differential equations

The existence of a solution defined for all t and possessing a type of boundedness property is established for the perturbed nonlinear system y = f(t,y) + F(t,y). The unperturbed system x = f(t,x) has a dichotomy in which some solutions exist and are well behaved as t increases to infinity, and some solution exists and are well behaved as t decreases to minus infinity. A similar study is made for a perturbed nonlinear differential equation defined on a half line, R+, and the existence of a family of solutions with special boundedness properties is established. The ideas are applied to integral manifolds.

Proctor, T. G.↗

Covariance and the hierarchy of frame bundles

This is an essay on the general concept of covariance, and its connection with the structure of the nested set of higher frame bundles over a differentiable manifold. Examples of covariant geometric objects include not only linear tensor fields, densities and forms, but affinity fields, sectors and sector forms, higher order frame fields, etc., often having nonlinear transformation rules and Lie derivatives. The intrinsic, or invariant, sets of forms that arise on frame bundles satisfy the graded Cartan-Maurer structure equations of an infinite Lie algebra. Reduction of these gives invariant structure equations for Lie pseudogroups, and for G-structures of various orders. Some new results are introduced for prolongation of structure equations, and for treatment of Riemannian geometry with higher-order moving frames. The use of invariant form equations for nonlinear field physics is implicitly advocated.

Estabrook, Frank B.↗

Estimation for bilinear stochastic systems

Three techniques for the solution of bilinear estimation problems are presented. First, finite dimensional optimal nonlinear estimators are presented for certain bilinear systems evolving on solvable and nilpotent lie groups. Then the use of harmonic analysis for estimation problems evolving on spheres and other compact manifolds is investigated. Finally, an approximate estimation technique utilizing cumulants is discussed.

Willsky, A. S.↗

Robust and real-time control of magnetic bearings for space engines

Currently, NASA Lewis Research Center is developing magnetic bearings for Space Shuttle Main Engine (SSME) turbopumps. The control algorithms which have been used are based on either the proportional-intergral-derivative control (PID) approach or the linear quadratic (LQ) state space approach. These approaches lead to an acceptable performance only when the system model is accurately known, which is seldom true in practice. For example, the rotor eccentricity, which is a major source of vibration at high speeds, cannot be predicted accurately. Furthermore, the dynamics of a rotor shaft, which must be treated as a flexible system to model the elastic rotor shaft, is infinite dimensional in theory and the controller can only be developed on the basis of a finite number of modes. Therefore, the development of the control system is further complicated by the possibility of closed loop system instability because of residual or uncontrolled modes, the so called spillover problem. Consequently, novel control algorithms for magnetic bearings are being developed to be robust to inevitable parametric uncertainties, external disturbances, spillover phenomenon and noise. Also, as pointed out earlier, magnetic bearings must exhibit good performance at a speed over 30,000 rpm. This implies that the sampling period available for the design of a digital control system has to be of the order of 0.5 milli-seconds. Therefore, feedback coefficients and other required controller parameters have to be computed off-line so that the on-line computational burden is extremely small. The development of the robust and real-time control algorithms is based on the sliding mode control theory. In this method, a dynamic system is made to move along a manifold of sliding hyperplanes to the origin of the state space. The number of sliding hyperplanes equals that of actuators. The sliding mode controller has two parts; linear state feedback and nonlinear terms. The nonlinear terms guarantee that the systems would reach the intersection of all sliding hyperplanes and remain on it when bounds on the errors in the system parameters and external disturbances are known. The linear part of the control drives the system to the origin of state space. Another important feature is that the controller parameter can be computed off-line. Consequently, on-line computational burden is small.

Sinha, Alok↗

Perturbed nonlinear differential equations

For perturbed nonlinear systems, a norm, other than the supremum norm, is introduced on some spaces of continuous functions. This makes possible the study of new types of behavior. A study is presented on a perturbed nonlinear differential equation defined on a half line, and the existence of a family of solutions with special boundedness properties is established. The ideas developed are applied to the study of integral manifolds, and examples are given.

Proctor, T. G.↗

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

We report recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry– Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extended harmonic mapping connects the equations in classical, statistical, fluid, quantum physics and general relativity

Abstract One potential pathway to find an ultimate rule governing our universe is to hunt for a connection among the fundamental equations in physics. Recently, Ren et al. reported that the harmonic maps with potential introduced by Duan, named extended harmonic mapping (EHM), connect the equations of general relativity, chaos and quantum mechanics via a universal geodesic equation. The equation, expressed as Euler–Lagrange equations on the Riemannian manifold, was obtained from the principle of least action. Here, we further demonstrate that more than ten fundamental equations, including that of classical mechanics, fluid physics, statistical physics, astrophysics, quantum physics and general relativity, can be connected by the same universal geodesic equation. The connection sketches a family tree of the physics equations, and their intrinsic connections reflect an alternative ultimate rule of our universe, i.e. , the principle of least action on a Finsler manifold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗