Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear normal modes”

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 55 records · Page 3

On the Hilbert-Huang Transform Theoretical Developments

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.

Kizhner, Semion↗

On Certain Theoretical Developments Underlying the Hilbert-Huang Transform

One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,

Kizhner, Semion↗

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗

Initial-value problems for Rossby waves in a shear flow with critical level

The time-dependent evolution of sheared Rossby waves starting from an initial disturbance is studied for the simple case in which the shear is uniform. The uniform-shear assumption allows explicit solutions to be obtained which are useful in addressing the issue of the long-time asymptotic approach to normal modes and in assessing the relative importance of viscosity, nonlinearity and time-dependence in the evolution of Rossby waves in the presence of critical layers.

Tung, K. K.↗

Nonlinear behavior of acoustic waves in combustion chambers

The nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber are considered. A formal framework is provided within which practical problems can be treated with a minimum of effort and expense. The general conservation equations were expanded in two small parameters, one characterizing the mean flow field and one measuring the amplitude of oscillations, and then combined to yield a nonlinear inhomogeneous wave equation. The unsteady pressure and velocity fields were expressed as syntheses of the normal modes of the chamber, but with unknown time-varying amplitudes. This procedure yielded a representation of a general unsteady field as a system of coupled nonlinear oscillators. The system of nonlinear equations was treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. The analysis is applicable to any combustion chamber. The most interesting applications are probably to solid rockets, liquid rockets, or thrust augmentors on jet engines.

Culick, F. E. C.↗

Solar oscillations

The main observational results and related theoretical investigations concerning solar oscillations are reviewed. The normal modes of oscillation of the sun are classified according to their angular eigenvalues and the number of their radial nodes. Observations of excited normal modes are discussed, particularly in relation to five-minute oscillations, nonlinearly coupled unstable g-modes, oscillations of the sun's apparent diameter, and reported oscillations of the solar surface with a period of 2 hr 40 min. Linear stability calculations are briefly examined, and stochastic excitation of p-modes by turbulent convection is considered. The five-minute oscillations are described as the superposition of many excited nonradial p-modes.

Goldreich, P.↗

Gravitational instabilities in protostellar disks

The nonaxisymmetric stability of self-gravitating, geometrically thick accretion disks has been studied for protostellar systems having a wide range of disk-to-central object mass ratios. Global eigenmodes with four distinctly different characters were identified using numerical, nonlinear hydrodynamic techniques. The mode that appears most likely to arise in normal star formation settings, however, resembles the 'eccentric instability' that was identified earlier in thin, nearly Keplerian disks: It presents an open, one-armed spiral pattern that sweeps continuously in a trailing direction through more than 2-pi radians, smoothly connecting the inner and outer edges of the disk, and requires cooperative motion of the point mass for effective amplification. This particular instability promotes the development of a single, self-gravitating clump of material in orbit about the point mass, so its routine appearance in our simulations supports the conjecture that the eccentric instability provides a primary route to the formation of short-period binaries in protostellar systems.

Tohline, J. E.↗

Analysis of nonlinear structures via mode synthesis

An effective procedure for NASTRAN was developed that permits any number of substructures of any size to be synthesized for the purpose of developing normal modes of vibration of the complete structural system. The technique is extended to permit modal transient analysis of the subdivided system. This latter procedure permits the use of NASTRAN's ability to include nonlinear forces in the problem. The five-phase process is accomplished using standard NASTRAN rigid formats with problem-independent alter packages and DMAP sequences.

Gieseke, R. K.↗

The barotropic normal modes in certain shear flows and the traveling waves in the atmosphere

It is shown analytically and numerically that in certain shear flows the linearized nondivergent barotropic vorticity equation has a limited number of neutral normal modes. The latitudinal structures of these shear flows can be expressed as polynomials of the sine of latitude. The first few such shear flows resemble the gross features of the zonal winds in the atmosphere of the earth at different times and altitudes. The spatial structures of the neutral normal modes in these shear flows are spherical harmonics, and, as a consequence, these modes are also the exact solutions of the fully nonlinear equation because the nonlinear interaction term vanishes identically. The spatial structures of the observed 5-, 4-, 2-, and 16-day free traveling waves in the atmosphere are often identified with the spherical harmonics with indices of (m, n) = ( 1, 2), (2, 3), (3, 3), and ( 1, 4), which are known previously as the neutral normal modes of the nondivergent barotropic vorticity equation in a motionless background state. Our results could explain why these free traveling waves can survive the shearing effects of zonal flows that are far different from rest because these spherical harmonics are also normal modes in certain shear flows that resemble the observations of the atmosphere.

Chen, Ping↗

On error detection in the dynamics part of primitive equation models

Rossby-Haurwitz waves have been commonly employed as initial conditions for testing the dynamics part of global atmospheric models. Difficulties arise in connection with the primitive equation (PE) model. Chao and Geller (1982) have suggested that it would be better if normal mode solutions of the linearized primitive equations would be used as initial conditions. However, tests involving normal mode conditions are not adequate for the detection of errors in the nonlinear terms. The present investigation has the objective to review the study conducted by Chao and Geller and to provide tests for the nonlinear terms.

Chao, W. C.↗

Numerical simulation of baroclinic Jovian vortices

We examine the evolution of baroclinic vortices in a time-dependent, nonlinear numerical model of a Jovian atmosphere. The model uses a normal-mode expansion in the vertical, using the barotropic and first two baroclinic modes. Results for the stability of baroclinic vortices on an f plane in the absence of a mean zonal flow are similar to results of Earth vortex models, although the presence of a fluid interior on the Jovian planets shifts the stability boundaries to smaller length scales. The presence of a barotropic mean zonal flow in the interior stabilizes vortices against instability and significantly modifies the finite amplitude form of baroclinic instabilities. The effect of a zonal flow on a form of barotropic instability produces periodic oscillations in the latitude and longitude of the vortex as observed at the level of the cloud tops. This instability may explain some, but not all, observations of longitudinal oscillations of vortices on the outer planets. Oscillations in aspect ratio and orientation of stable vortices in a zonal shear flow are observed in this baroclinic model, as in simpler two-dimensional models. Such oscillations are also observed in the atmospheres of Jupiter and Neptune. The meridional propagation and decay of vortices on a beta plane is inhibited by the presence of a mean zonal flow. The direction of propagation of a vortex relative to the mean zonal flow depends upon the sign of the meridional potential vorticity gradient; combined with observations of vortex drift rates, this may provide a constraint on model assumption for the flow in the deep interior of the Jovian planets.

Achterberg, Richard K.↗

Dynamical criterion for a marginally unstable, quasi-linear behavior in a two-layer model

A two-layer quasi-geostrophic flow forced by meridional variations in heating can be in regimes ranging from radiative equilibrium to forced geostrophic turbulence. Between these extremes is a regime where the time-mean (zonal) flow is marginally unstable. Using scaling arguments, it is concluded that such a marginally unstable state should occur when a certain parameter, measuring the strength of wave-wave interactions relative to the beta effect and advection by the thermal wind, is small. Numerical simulations support this proposal. A transition from the marginally unstable regime to a more nonlinear regime is then examined through numerical simulations with different radiative forcings. It is found that transition is not caused by secondary instability of waves in the marginally unstable regime. Instead, the time-mean flow can support a number of marginally unstable normal modes. These normal modes interact with each other, and if they are of sufficient amplitude, the flow enters a more nonlinear regime.

Ebisuzaki, W.↗

Disruption avoidance via radio frequency current condensation in magnetic islands produced by off-normal events

This article discusses the use of radio frequency (RF) current drive to stabilize large islands, focusing on nonlinear effects that appear when relatively high powers are used to stabilize large islands. We are interested in developing a capability to stabilize large islands via RF driven currents to avoid the need for mitigation to the extent possible. As tokamaks are designed and built with increasing levels of stored energy in the plasma, disruptions become increasingly dangerous. It has been reported that 95% of the disruptions in the Joint European Torus tokamak with the ITER-like wall are preceded by the growth of large locked islands. These large islands are mostly produced by off-normal events other than neoclassical tearing modes. This paper presents theory and modeling for a nonlinear “RF current condensation” effect that can concentrate the RF driven current near the center of a large island, thereby increasing the efficiency of the stabilization. A nonlinear shadowing effect can hinder the stabilization of islands if the aiming of the ray trajectories does not properly consider the nonlinear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Plasma Emission in Magnetized Plasmas

The recent Parker Solar Probe observations of type III radio bursts show that the effects of the finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of the background magnetic field on the plasma-emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of the ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyrofrequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave–particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second-harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less clear. The underlying nonlinear-mode coupling processes indicate that to properly understand and interpret the simulation results requires sophisticated analyses involving interactions among magnetized plasma normal modes, including the two transverse modes of the magneto-active plasma, namely, the extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex-mode coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.

79 ASTRONOMY AND ASTROPHYSICS↗

Differential stiffness effects

Differential stiffness as developed in NASTRAN is a linear change in stiffness caused by applied loads. Examples of differential stiffness are the stiffening effects of gravity forces in a pendulum, centrifugal forces in rotor blades and pressure loading of shell structures. In cases wherein this stiffness caused by a load is destabilizing, the differential stiffness concept lends itself to nonlinear structural analysis. Rigid Formats 4 (static analysis with differential stiffness) and 13 (normal modes with differential stiffness) are specifically designed to account for such stiffness changes. How pressure loading may be treated in these rigid formats is clarified. This clarification results from modal correlation of Ground Vibration Test (GVT) results from the empty and pressurized Filament Wound Case (FWC) quarter-scale Space Shuttle solid rocket booster (QSSRB). A sketch of the QSSRB cantilevered to the floor by its external tank attachments is shown.

Ujihara, B. H.↗

Instabilities of two-dimensional inviscid compressible vortices

An investigation of the linear stability and subsequent nonlinear evolution and acoustic radiation of a planar inviscid compressible vortex is presented. The effects of the entropy gradient are investigated, and for the particular entropy profile chosen, the positive average entropy in the vortex core is destabilizing, while the opposite is true for the negative average entropy gradient. Finite-difference methods are used to study the linear initial value problem. These methods are found to be capable of accurately computing the frequencies and weak growth rates of the normal modes. When the initial condition consists of random perturbations, the long-time behavior is found to correspond to the most unstable normal mode in all cases. The numerical solution of the Euler equations is used to study the nonlinear development of an elliptic-mode perturbation.

Chan, W. M.↗

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Chemistry↗

Nonlinear Pressurization and Modal Analysis Procedure for Dynamic Modeling of Inflatable Structures

An introduction and set of guidelines for finite element dynamic modeling of nonrigidized inflatable structures is provided. A two-step approach is presented, involving 1) nonlinear static pressurization of the structure and updating of the stiffness matrix and 2) hear normal modes analysis using the updated stiffness. Advantages of this approach are that it provides physical realism in modeling of pressure stiffening, and it maintains the analytical convenience of a standard bear eigensolution once the stiffness has been modified. Demonstration of the approach is accomplished through the creation and test verification of an inflated cylinder model using a large commercial finite element code. Good frequency and mode shape comparisons are obtained with test data and previous modeling efforts, verifying the accuracy of the technique. Problems encountered in the application of the approach, as well as their solutions, are discussed in detail.

Smalley, Kurt B.↗