SEARCH · Search NASA
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.
Indirect method considers effect of nonlinearity and obtains approximate probability densities of certain characteristics, from these densities various statistical properties are calculated. Direct method determines response moments directly. Neither method is restricted by system motion with regard to whether or not it is stationary.
The nonlinear evolution of marginally unstable wave packets in rotating pipe flow is studied. These flows depend on two control parameters, which may be taken to be the axial Reynolds number R and a Rossby number, q. Marginal stability is realized on a curve in the (R, q)-plane, and the entire marginal stability boundary is explored. As the flow passes through any point on the marginal stability curve, it undergoes a supercritical Hopf bifurcation and the steady base flow is replaced by a traveling wave. The envelope of the wave system is governed by a complex Ginzburg-Landau equation. The Ginzburg-Landau equation admits Stokes waves, which correspond to standing modulations of the linear traveling wavetrain, as well as traveling wave modulations of the linear wavetrain. Bands of wavenumbers are identified in which the nonlinear modulated waves are subject to a sideband instability.
Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.
This paper provides an experimental and theoretical investigation of azimuthal instabilities in a lean, premixed combustor. It presents data from a high Reynolds number, industrial-scale facility with multiple fuel stages that enable variation of radial and azimuthal symmetries. The combustor exhibits a variety of behaviors – including standing, spinning, mixed waves, and intermittent behaviors, as well as periodic limit cycles and quasi-periodic oscillations – depending upon varying power and azimuthal symmetry. A second key objective of this paper is to present a weakly nonlinear framework that is used to understand and interpret these results, showing the influence of eigenmode degeneracy, background noise, and azimuthal non-uniformities in mean quantities such as temperature and flame-acoustic coupling. In particular, this analysis shows the importance of the linear instability growth/damping rate and azimuthal non-uniformities in controlling the relative prominence of standing, spinning, or mixed waves at the limit cycle, as well as the preference in the spinning direction. It also shows the conditions under which single frequency or quasi-periodic oscillations are dominant, and how these are more fundamentally influenced by azimuthal non-uniformities and background noise.
Herein we study the nonlinear evolution of very oblique small-scale Alfvén waves with k ⊥ d i ≳1. At these scales, the waves become significantly compressive, unlike in magnetohydrodynamics, due to the Hall term in the equations. We demonstrate that when frequencies are small compared to the ion gyrofrequency and amplitudes are small compared to unity, no new nonlinear interaction appears due to the Hall term alone at the lowest non-trivial order, even when k ⊥ d i ~1. However, at the second non-trivial order, we discover that the Hall physics leads to a slow but resonant nonlinear interaction between co-propagating Alfvén waves, an inherently three-dimensional effect. Including the effects of finite temperature, finite frequency, and electron inertia, the two-fluid Alfvén wave also becomes dispersive once one or more of k ⊥ ρ s , k ⊥ d e , or k ∥ d i becomes significant: for oblique waves at low β as studied here, this can be at a much smaller scale than d i . We show that the timescale for one-dimensional steepening of two-fluid Alfvén waves is only significant at these smaller dispersive scales, and also derive an expression for the amplitude of driven harmonics of a primary wave. Importantly, both new effects are absent in gyrokinetics and other commonly used reduced two-fluid models. Our calculations have relevance for the interpretation of laboratory Alfvén wave experiments, as well as shedding light on the physics of turbulence in the solar corona and inner solar wind, where the dominant nonlinear interaction between counter-propagating waves is suppressed, allowing these new effects to become important.
Extensive linear and nonlinear simulations have been carried out to investigate the energetic particle-driven fishbone instability in spherical tokamak plasmas with weakly reversed q profile and the q{sub min} slightly above unity. The global kinetic-MHD hybrid code M3D-K is used. Numerical results show that a fishbone instability is excited by energetic beam ions preferentially at higher q{sub min} values, consistent with the observed appearance of the fishbone before the “long-lived mode” in MAST and NSTX experiments. In contrast, at lower q{sub min} values, the fishbone tends to be stable. In this case, the beam ion effects are strongly stabilizing for the non-resonant kink mode. Nonlinear simulations show that the fishbone saturates with strong downward frequency chirping as well as radial flattening of the beam ion distribution. An (m, n) = (2, 1) magnetic island is found to be driven nonlinearly by the fishbone instability, which could provide a trigger for the (2, 1) neoclassical tearing mode sometimes observed after the fishbone instability in NSTX.
We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (quadratic gravity) in the spherically-symmetric sector. The formulation relies on (i) the harmonic gauge to cast the evolution system into quasilinear form (ii) the Cartoon method to reduce to spherical symmetry in keeping with the harmonic gauge, and (iii) order reduction to first order (in time) by means of introducing auxiliary variables. The well posedness of the respective initial-value problem is numerically confirmed by evolving randomly perturbed flat-space and black-hole initial data. In conclusion, our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives.
In this paper we propose a new Koopman operator approach for nonlinear model decomposition using Koopman gramians. We introduce Koopman operators and discuss the notion of a input-Koopman operator. We show how input-Koopman operators can be used to cast a nonlinear system into the classical state-space form, as well as conditions under which input and state observable functions are well separated. We extend an existing method of dynamic mode decomposition for learning Koopman operators from data known as deep dynamic mode decomposition, for systems with control or disturbances. We illustrate the accuracy of the method on learning an input-state separable Koopman operator for an example system, even when the underlying system exhibits mixed state-input terms. Finally, we introduce a nonlinear decomposition algorithm, based on Koopman gramians, that maximizes internal subsystem observability and disturbance rejection from unwanted noise in other subsystems. To relieve the NP-hardness of the decomposition problem, we derive a relaxation based on Koopman gramians and multi-way partitioning. We illustrate our algorithm on a swing dynamics model for a IEEE 39 bus and discuss our findings.
We propose an ensemble score filter (EnSF) for solving high-dimensional nonlinear filtering problems with superior accuracy. A major drawback of existing filtering methods, e.g., particle filters or ensemble Kalman filters, is the low accuracy in handling high-dimensional and highly nonlinear problems. EnSF attacks this challenge by exploiting the score-based diffusion model, defined in a pseudo-temporal domain, to characterizing the evolution of the filtering density. EnSF stores the information of the recursively updated filtering density function in the score function, in stead of storing the information in a set of finite Monte Carlo samples (used in particle filters and ensemble Kalman filters). Unlike existing diffusion models that train neural networks to approximate the score function, we develop a training-free score estimation that uses mini-batch-based Monte Carlo estimator to directly approximate the score function at any pseudo-spatial-temporal location, which provides sufficient accuracy in solving high-dimensional nonlinear problems as well as saves tremendous amount of time spent on training neural networks. Another essential aspect of EnSF is its analytical update step, gradually incorporating data information into the score function, which is crucial in mitigating the degeneracy issue faced when dealing with very high-dimensional nonlinear filtering problems. High-dimensional Lorenz systems are used to demonstrate the performance of our method. EnSF provides surprisingly impressive performance in reliably tracking extremely high-dimensional Lorenz systems (up to 1,000,000 dimension) with highly nonlinear observation processes, which is a well-known challenging problem for existing filtering methods.
A simple nonlinear, axisymmetric, shallow-water numerical model has been used to study the asymmetry in the neutral flow between the dusk and dawn sides of the auroral oval. The results indicate that the Coriolis force and the curvature terms are nearly in balance on the evening side and require only a small pressure gradient to effect adjustment. The result is smaller neutral velocities near dawn and larger velocities near dusk than would be the case for a linearized treatment. A consequence is that more gravity wave energy is produced on the morning side than on the evening side.
The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.
The cycling combustion instabilities in a diesel engine have been analyzed based on chaos theory. The objective was to investigate the dynamical characteristics of combustion in diesel engine. In this study, experiments were performed under the entire operating range of a diesel engine (the engine speed was changed from 600 to 1400 rpm and the engine load rate was from 0% to 100%), and acquired real-time series of in-cylinder combustion pressure using a piezoelectric transducer installed on the cylinder head. Several methods were applied to identify and quantitatively analyze the combustion process complexity in the diesel engine including delay-coordinate embedding, recurrence plot (RP), Recurrence Quantification Analysis, correlation dimension (CD), and the largest Lyapunov exponent (LLE) estimation. The results show that the combustion process exhibits some determinism. If LLE is positive, then the combustion system has a fractal dimension and CD is no more than 1.6 and within the diesel engine operating range. We have concluded that the combustion system of diesel engine is a low-dimensional chaotic system and the maximum values of CD and LLE occur at the lowest engine speed and load. This means that combustion system is more complex and sensitive to initial conditions and that poor combustion quality leads to the decrease of fuel economy and the increase of exhaust emissions.
No abstract prepared.
We report on the measurement of the instantaneous Kerr nonlinearity and the retarded alignment of air molecules CO{sub 2}, N{sub 2}, and O{sub 2} triggered by an intense, lambda-cubic terahertz pulse, a diffraction- and transform-limited single-cycle pulse. The strong-field, impulsive low-frequency excitation (3.9 THz) leads to field-free alignment dynamics of these molecules thanks to the terahertz-induced transient dipole moments in the otherwise non-polar molecules. The strong coupling to the terahertz electric transient results in the excitation of coherent large amplitude long-living rotational states at room temperature and ambient pressure. Beyond fundamental investigations of nonlinear properties in gases, our results suggest a route towards field-free molecular alignment at laser intensity well below the ionization threshold.
The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.
Curvature and diamagnetic effects are both known to have an influence on tearing mode dynamics. In this paper, we investigate the impact of these effects on the nonlinear stability and saturation of a (2, 1) island using non-linear two-fluid MHD simulations and we apply our results to Tore Supra experiments, where its behavior is not well understood from the single fluid MHD model. Simulations show that a metastable state induced by diamagnetic effect exists for this mode and that it also produces a reduction of the saturated island size, in presence of toroidal curvature. The mode is found to be nonlinearly destabilized by a seed island and it saturates at a macroscopic level causing a significant confinement degradation. The interpretation of dual states, with either no island on q = 2 or a large one, observed on discharges with high non inductive current source on Tore Supra, is revisited.
Motivations (Ensure a Higher Level of Confidence in the Predictability & Reliability of Numerical Simulation for Multiscale Complex Nonlinear Fluid Problems) - The last two decades have been an era when computation is ahead of analysis & when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems & non-traditional fields (Especially when computations offer the ONLY way of generating this type of data limited simulations). - At present some of the numerical uncertainties can be explained and minimized by traditional numerical analysis and standard CFD practices. However, such practices, usually based on linearized analysis, MIGHT NOT be sufficient for strongly nonlinear and/or stiff problems. - We need a good understanding of the nonlinear behavior of numerical schemes being used as an integral part of code verification, validation and certification.
Recent developments in the study of integrable Hamiltonian systems have led to nonlinear accelerator lattice designs with two transverse invariants. These lattices may drastically improve the performance of high-power machines, providing wide tune spreads and Landau damping to protect the beam from instabilities, while preserving dynamic aperture. To test the feasibility of these concepts, the Integrable Optics Test Accelerator (IOTA) is being designed and built at Fermilab. One way to obtain a nonlinear integrable lattice is by using the fields generated by a magnetically confined electron beam (electron lens) overlapping with the circulating beam. The parameters of the required device are similar to the ones of existing electron lenses. We present theory, numerical simulations, and first design studies of electron lenses for nonlinear integrable optics.