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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Turbine blade forced response prediction using FREPS

This paper describes a software system called FREPS (Forced REsponse Prediction System) that integrates structural dynamic, steady and unsteady aerodynamic analyses to efficiently predict the forced response dynamic stresses in axial flow turbomachinery blades due to aerodynamic and mechanical excitations. A flutter analysis capability is also incorporated into the system. The FREPS system performs aeroelastic analysis by modeling the motion of the blade in terms of its normal modes. The structural dynamic analysis is performed by a finite element code such as MSC/NASTRAN. The steady aerodynamic analysis is based on nonlinear potential theory and the unsteady aerodynamic analyses is based on the linearization of the non-uniform potential flow mean. The program description and presentation of the capabilities are reported herein. The effectiveness of the FREPS package is demonstrated on the High Pressure Oxygen Turbopump turbine of the Space Shuttle Main Engine. Both flutter and forced response analyses are performed and typical results are illustrated.

Murthy, Durbha, V.↗
NASA NTRS1994

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗
DOE OSTI2026

Impact of n > 1 neoclassical tearing modes on fast ions, plasma rotation, and the onset of disruptive tearing modes in DIII-D ITER baseline scenario

We present integrated TRANSP-kick analyses quantitatively characterizing how n > 1 magnetic islands redistribute fast ions (FIs) and modify momentum balance in low-torque DIII-D H-mode scenarios operated with the ITER normalized parameter set and shape (ITER baseline scenario, or IBS). In this plasma scenario, disruptive neoclassical tearing modes (NTMs) are most commonly seeded by nonlinear 3-wave coupling when the differential rotation between the q = 1 and q = 2 surfaces ($δf_{1, 2}$) approaches zero. Analysis of a DIII-D IBS database of shots unstable to the 2/1 NTM demonstrates that the flattening rotation profile is correlated with n > 1 NTM amplitude growth, with a strong drop in $δf_{1, 2}$ occurring when the radial magnetic perturbation amplitude at the rational surface ($δB$ r s ) is estimated to be about 20 G. The interpretive TRANSP-kick simulations show that the experimentally observed reduction in core rotation is due to the fast-ion kicks. Nonlinearity in momentum losses versus NTM amplitude is found to be correlated with overlap of FI island chains, occurring when $δB$ r s is about 35 G, in rough agreement with the experimental results. This agreement suggests the primary mechanism for the conductive momentum loss in the core is due to chaotic FI orbits and FI losses developing from overlapping resonant FI island chains.

ITER baseline scenario↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
NASA NTRS1987

Second-overtone models of RR Lyrae stars

Linear models of RR Lyrae stars pulsating in the second overtone indicate that this mode is self-excited only for rather low luminosities if the masses are normal. An unstable second-overtone model has been selected from the linear survey and has been followed up to limiting amplitude by using a nonlinear hydrodynamical computer program. Deep splitting of the main light peak in the theoretical light curve occurs but does not resemble any observed feature. Moreover, the distribution of periods of RR Lyrae stars in globular clusters can be accounted for adequately by fundamental-mode and first-overtone pulsation. It is concluded that second-overtone pulsators probably do not exist among RR Lyrae stars.

Stothers, Richard B.↗
NASA NTRS1990

Rapid near-optimal trajectory generation and guidance law development for single-stage-to-orbit airbreathing vehicles

General problems associated with on-board trajectory optimization, propulsion system cycle selection, and with the synthesis of guidance laws were addressed for an ascent to low-earth-orbit of an air-breathing single-stage-to-orbit vehicle. The NASA Generic Hypersonic Aerodynamic Model Example and the Langley Accelerator aerodynamic sets were acquired and implemented. Work related to the development of purely analytic aerodynamic models was also performed at a low level. A generic model of a multi-mode propulsion system was developed that includes turbojet, ramjet, scramjet, and rocket engine cycles. Provisions were made in the dynamic model for a component of thrust normal to the flight path. Computational results, which characterize the nonlinear sensitivity of scramjet performance to changes in vehicle angle of attack, were obtained and incorporated into the engine model. Additional trajectory constraints were introduced: maximum dynamic pressure; maximum aerodynamic heating rate per unit area; angle of attack and lift limits; and limits on acceleration both along and normal to the flight path. The remainder of the effort focused on required modifications to a previously derived algorithm when the model complexity cited above was added. In particular, analytic switching conditions were derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another for two cases: the case in which engine cycle operations can overlap, and the case in which engine cycle operations are mutually exclusive. The resulting guidance algorithm was implemented in software and exercised extensively. It was found that the approximations associated with the assumed time scale separation employed in this work are reasonable except over the Mach range from roughly 5 to 8. This phenomenon is due to the very large thrust capability of scramjets in this Mach regime when sized to meet the requirement for ascent to orbit. By accounting for flight path angle and flight path angle rate in construction of the flight path over this Mach range, the resulting algorithm provides the means for rapid near-optimal trajectory generation and propulsion cycle selection over the entire Mach range from take-off to orbit.

Calise, A. J.↗
NASA NTRS1983

Large-amplitude hydromagnetic waves in collisionless relativistic plasma - Exact solution for the fast-mode magnetoacoustic wave

An exact nonlinear solution is found to the relativistic kinetic and electrodynamic equations (in their hydromagnetic limit) that describes the large-amplitude fast-mode magnetoacoustic wave propagating normal to the magnetic field in a collisionless, previously uniform plasma. It is pointed out that a wave of this kind will be generated by transverse compression of any collisionless plasma. The solution is in essence independent of the detailed form of the particle momentum distribution functions. The solution is obtained, in part, through the method of characteristics; the wave exhibits the familiar properties of steepening and shock formation. A detailed analysis is given of the ultrarelativistic limit of this wave.

Barnes, A.↗
DOE OSTI2021

Instability Issue of Paralleled Dies in an SiC Power Module in Solid-State Circuit Breaker Applications

Paralleled dies in a power module could have instability issues during high current switching transients. Here, the instability is caused by the differential-mode oscillation among paralleled MOSFETs. Conventional analyses of paralleled MOSFETs’ stability are normally limited to a single operating point, which ignores the influences of the switching trajectory and nonlinear device parameters on stability. This article reveals that the switching trajectory can significantly influence parallel stability. The analysis is improved by solving eigenvalues of state-space modeling system matrices of all operating points that the switching trajectory goes through considering nonlinear device parameters. Higher voltage and current stresses result in greater real parts of complex eigenvalues, which explains why the paralleled MOSFETs are more unstable with higher voltage and current stresses. To improve stability in solid-state circuit breaker applications, we propose a method to manipulate the switching trajectory to avoid the unstable region where the conventional hard switching trajectory normally goes through. Experimental results show that the turn- off current capability can be increased from ~five times of rated current with the gate oscillation using the conventional turn- off trajectory to ~ten times of rated current without the gate oscillation using the optimal turn- off trajectory.

42 ENGINEERING↗
NASA NTRS1986

Circulation regime-dependent nonlinear interactions during northern hemisphere winter

Data of 500 mb geopotential heights, obtained in Januaries from 1963-1977 are used to study the relationship between the interaction of different scales of motion and the circulation regime characterizing the large-scale pattern. Geostrophic and 500 mb height streamfunctions and energy and enstrophy spectra are computed. Differences in the height variance, kinetic energy, and enstrophy spectra of the large positive anomaly events are compared to normal circulation. Systematic differences in the nonlinear cascading are examined. The data reveal that transient eddies perform a regime-dependent function in the interactions with atmospheric circulation modes on the scale of the persistent anomalies.

Brown, P. S.↗
NASA NTRS2014

DNS of Laminar-Turbulent Transition in Swept-Wing Boundary Layers

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.

Duan, L.↗
NASA NTRS1993

Evaluation of MARC for the analysis of rotating composite blades

The suitability of the MARC code for the analysis of rotating composite blades was evaluated using a four-task process. A nonlinear displacement analysis and subsequent eigenvalue analysis were performed on a rotating spring mass system to ensure that displacement-dependent centrifugal forces were accounted for in the eigenvalue analysis. Normal modes analyses were conducted on isotropic plates with various degrees of twist to evaluate MARC's ability to handle blade twist. Normal modes analyses were conducted on flat composite plates to validate the newly developed coupled COBSTRAN-MARC methodology. Finally, normal modes analyses were conducted on four composite propfan blades that were designed, analyzed, and fabricated at NASA Lewis Research Center. Results were compared with experimental data. The research documented herein presents MARC as a viable tool for the analysis of rotating composite blades.

Bartos, Karen F.↗
DOE OSTI2025

Flat nonlinear optics with intersubband polaritonic metasurfaces

Nonlinear intersubband polaritonic metasurfaces produce some of the strongest second- and third-order nonlinear optical responses reported for condensed matter systems at infrared frequencies. These metasurfaces are fabricated as two-dimensional arrays of nanoresonators from multi-quantum-well semiconductor heterostructures, designed to produce strong nonlinear responses associated with intersubband transitions. By optimally coupling the optical modes of the nanoresonators to vertically polarized intersubband transitions in semiconductor heterostructures, one can boost the nonlinear response associated with intersubband transitions, make intersubband transitions interact with free-space radiation at normal incidence, and hence produce optically thin flat nonlinear optical elements compatible with free-space optical setups. As a result of the strong nonlinear response in these systems, significant nonlinear conversion efficiencies (>0.1 %) can be attained in deeply subwavelength optical films using modest pumping intensities of only 10–100 kW/cm 2 . Subwavelength metasurface thickness relaxes phase-matching constraints limiting the operation of bulk nonlinear crystals. Furthermore, the amplitude and phase of the nonlinear optical response in intersubband polaritonic metasurfaces can be tailored for a specific pump wavelength and a nonlinear process of interest through the co-optimization of quantum engineering of electron states in semiconductor heterostructures and photonic engineering of the metasurface nanoresonators design. Additionally, an applied voltage can dynamically control the amplitude and phase of the nonlinear optical response at a nanoresonator level. Here, we review the current state of the art in this rapidly expanding field, focusing on nonlinear processes supporting second-harmonic generation, saturable absorption, and optical power limiting.

flat optics↗
NASA NTRS1998

Study of the effects of MeV Ag, Cu, Au, and Sn implantation on the optical properties of LiNbO3

The authors present the results of characterization of linear absorption and nonlinear refractive index of Au, Ag, Cu and Sn ion implantation into LiNbO3. Ag was implanted at 1.5 MeV to fluences of 2 to 17 x 17(exp 16)/sq cm at room temperature. Au and Cu were implanted to fluences of 5 to 20 x 10(exp 16)/sq cm at an energy of 2.0 MeV. Sn was implanted to a fluence of 1.6 x 10(exp 17)/sq cm at 160 kV. Optical absorption spectrometry indicated an absorption peak for the Au implanted samples after heat treatment at 1,000 C at approx. 620 nm. The Ag implanted samples absorption peaks shifted from approx. 450 nm before heat treatment to 550 nm after 500 C for 1h. Heat treatment at 800 C returned the Ag implanted crystals to a clear state. Cu nanocluster absorption peaks disappears at 500 C. No Sn clusters were observed by optical absorption or XRD. The size of the Ag and Au clusters as a function of heat treatment were determined from the absorption peaks. The Ag clusters did not change appreciably in size with heat treatment. The Au clusters increased from 3 to 9 nm diameter upon heat treatment at 1000 C. TEM analysis performed on a Au implanted crystal indicated the formation of Au nanocrystals with facets normal to the c-axis. Measurements of the nonlinear refractive indices were carried out using the Z-scan method with a tunable dye laser pumped by a frequency doubled mode-locked Nd:YAG laser. The dye laser had a 4.5 ps pulse duration time and 76 MHz pulse repetition rate (575 nm).

Williams, E. K.↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

A real-time guidance algorithm for aerospace plane optimal ascent to low earth orbit

Problems of onboard trajectory optimization and synthesis of suitable guidance laws for ascent to low Earth orbit of an air-breathing, single-stage-to-orbit vehicle are addressed. A multimode propulsion system is assumed which incorporates turbojet, ramjet, Scramjet, and rocket engines. An algorithm for generating fuel-optimal climb profiles is presented. This algorithm results from the application of the minimum principle to a low-order dynamic model that includes angle-of-attack effects and the normal component of thrust. Maximum dynamic pressure and maximum aerodynamic heating rate constraints are considered. Switching conditions are derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another. A nonlinear transformation technique is employed to derived a feedback controller for tracking the computed trajectory. Numerical results illustrate the nature of the resulting fuel-optimal climb paths.

Calise, A. J.↗
DOE OSTI2023

Gyrokinetic simulations of core turbulence and thermal transport in the high-β P discharge on EAST

The properties of core turbulence and thermal transport are investigated for EAST high-β P (β P ~ 3.1) plasmas with dominant electron heating (Te/Ti > 1) via gyrokinetic simulation with the NLT code. Linear simulations identify that the electrostatic η e -driven trapped electron mode (η e -TEM) dominates in the core region (ρ< 0.7) and the ion temperature gradient (ITG) mode dominates in the out region (ρ≥ 0.7), consistent with the linear threshold analysis of micro-instabilities. Sensitivity analysis shows that the normalized electron density gradient (R/L ne ) and ITG (R/L Ti ) are two effective parameters to stabilize TEM instability. Nonlinear simulations are also carried out and compared with the experimental results, which show that the electron thermal internal transport barrier (ITB) in EAST high-β P plasma is determined by the TEM-induced turbulence. A higher zonal flow shearing rate is observed in the ITB region (0.2 < ρ < 0.34), which can regulate energy transport induced by TEM turbulence and facilitate the formation of e-ITB. A plausible positive feedback mechanism to mitigate turbulence transport and improve energy confinement via enhanced ion heating is proposed for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A beamline for ultrafast extreme ultraviolet magneto-optical spectroscopy in reflection near the shot noise limit

High harmonic generation (HHG) makes it possible to measure spin and charge dynamics in materials on femtosecond to attosecond timescales. However, the extreme nonlinear nature of the high harmonic process means that intensity fluctuations can limit measurement sensitivity. Here we present a noise-canceled, tabletop high harmonic beamline for time-resolved reflection mode spectroscopy of magnetic materials. We use a reference spectrometer to independently normalize the intensity fluctuations of each harmonic order and eliminate long term drift, allowing us to make spectroscopic measurements near the shot noise limit. Here, these improvements allow us to significantly reduce the integration time required for high signal-to-noise (SNR) measurements of element-specific spin dynamics. Looking forward, improvements in the HHG flux, optical coatings, and grating design can further reduce the acquisition time for high SNR measurements by 1–2 orders of magnitude, enabling dramatically improved sensitivity to spin, charge, and phonon dynamics in magnetic materials.

47 OTHER INSTRUMENTATION↗
NASA NTRS1981

Dynamic stability of a rotor blade using finite element analysis

The aeroelastic stability of flap bending, lead-lag bending, and torsion of a helicopter rotor blade in hover is examined using a finite element formulation based on the principle of virtual work. Quasi-steady two-dimensional airfoil theory is used to obtain the aerodynamic loads. The rotor blade is discretized into beam elements, each with ten modal degrees of freedom. The resulting nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The flutter solution is calculated assuming blade motion to be a small perturbation about the steady solution. The normal mode method based on the coupled rotating natural modes about the steady deflections is used to reduce the number of equations in the flutter eigenanalysis. Results are presented for hingeless and articulated rotor blade configurations.

Sivaneri, N. T.↗