SEARCH · Engineering Papers

Results for “harmonic instability”

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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

NASA NTRSDate not supplied

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

boundary layer transition↗

NASA NTRS2019

Source of the Bursty Bulk Flow Diffuse Aurora: Electrostatic Cyclotron Harmonic and Whistler Waves in the Coupling of Bursty Bulk Flows to Auroral Precipitation

Electron cyclotron harmonic (ECH) and whistler chorus waves are recognized as the two mechanisms responsible for the resonant wave‐particle interactions necessary to precipitate plasma sheet electrons into the ionosphere, producing the diffuse Aurora. Previous work has demonstrated ECH waves dominate electron scattering at L shells >8, while whistler chorus dominates scattering at L shells L < 8. However, we ﬁnd from Time History of Events and Macroscale (THEMIS) Interactions during Substorms observations of fast ﬂows at L = 12 that oblique whistler chorus emissions play the dominant role in scattering electrons. Previous works have identiﬁed whistler‐mode waves within fast ﬂows that are produced by an electron temperature anisotropy Te,⊥/Te,||> 1, consistent with electron betatron acceleration. Here, however, we ﬁnd whistler chorus emissions throughout an interval of fast ﬂows where Te,⊥/Te,||< 1. Parallel electron beams account for the enhanced parallel electron temperature and serve as the instability mechanism for the whistler chorus. The parallel electron beams and associated cigar‐shaped distributions are consistent with Fermi acceleration at dipolarizations in fast ﬂows. We demonstrate that the scattering efﬁciency of the whistler chorus exceeds that of ECH waves, which THEMIS also detects during the fast ﬂows. The obliquity of the whistler waves permits efﬁcient scattering of lower‐energy electrons into the diffuse aurora. We conclude that Fermi acceleration of electrons provides one important free‐energy source for the wave‐particle interactions responsible for coupling plasma sheet electrons into the diffuse aurora during substorm conditions.

Wendel, D. E.↗

DOE OSTI2025

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

On the unique structure of the magnetic fields of Uranus and Neptune

The magnetic fields of Uranus and Neptune, which have comparable dipole, quadrupole, and octupole harmonics, are unique in the present-day solar system, but they resemble the geomagnetic field at the epochs of excursions and reversals known from paleomagnetic data. The precession dynamo model, in which the dominant role in the generation of the planetary magnetic fields is played by external gravitational forces, allows us to propose two scenarios for the formation of the unique topology of the magnetic fields of Uranus and Neptune. In the first case, tidal flows in the 'oceans' of these two planets extend down to the depths where the matter has a noticeable electric conductivity and velocity. A hydromagnetic interaction of the moving conducting fluid with the planetary magnetic field outside the generation region results in the deformation of the field and the deceleration of the motion under the action of the radial magnetic field. In the second case, the deformation of the field facilitates drastic changes in cyclonic cells within the generation region causing instabilities that result in a multi-polar field structure, excursions, and inversions. This paper considers this problem in greater detail by using the Neptune-Triton system as an example.

Dolginov, Sh. SH.↗

DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗

DOE OSTI2020

A regional nuclear conflict would compromise global food security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 °C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30°N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2020

A Regional Nuclear Conflict Would Compromise Global Food Security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 ◦C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30◦N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

Jaegermeyr, Jonas↗

DOE OSTI2023

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗

NASA NTRS1995

Variability of accretion disks surrounding black holes: The role of inertial-acoustic mode instabilities

The global nonlinear time-dependent evolution of the inertial-acoustic mode instability in accretion disks surrounding black holes has been investigated. The viscous stress is assumed to be proportional to the gas pressure only, i.e., tau = alphap(sub g). It is found that an oscillatory nonsteady behavior exists in the inner regions of disks (r is less than 10r(sub g) where r(sub g) is the Schwarzschild radius) for sufficiently large alpha(greater than or approximately equal to 0.2) and for mass accretion rates less than about 0.3 times the Eddington value. The variations of the integrated bolometric luminosity from the disk, Delta L/L, are less than 3%. A power spectrum analysis of these variations reveals a power spectrum which can be fitted to a power-law function of the frequency Pis proportional to f(exp -gamma), with index gamma = 1.4-2.3 and a low-frequency feature at about 4 Hz in one case. In addition, a narrow peak centered at a frequency corresponding to the maximum epicyclic frequency of the disk at approximately 100-130 Hz and its first harmonic is also seen. The low-frequency modulations are remarkably similar to those observed in black hole candidate systems. The possible existence of a scattering corona in the inner region of the disk and/or other processes contributing to the power at high frequencies in the inner region of the accretion disk may make the detection of the high-frequency component difficult.

Chen, Xingming↗

NASA NTRS2010

Harmonic Balance Computations of Fan Aeroelastic Stability

A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study.

Bakhle, Milind A.↗

NASA NTRSDate not supplied

Key Observations and Discoveries from the Vector Electric Field Investigation on the C/NOFS Satellite

Launched in 2008 and operating for 7.5 years, the Air Force Communication /Navigation Outage Forecasting System (C/NOFS) satellite included the Vector Electric Field Investigation (VEFI) designed and built at NASA’s Goddard Space Flight Center. VEFI successfully met its objectives: 1) investigate the role of ambient electric fields in initiating nighttime density depletions and turbulence; 2) determine the quasi-DC electric fields associated with abrupt density depletions, and 3) quantify the spectrum of the irregularities associated with density depletions, providing many key observations and discoveries including: -- Global (low latitude) vector DC electric fields at 16 s/sec revealing variations with local time, longitude, and season between extremely low and moderate solar activity-- Reversed zonal E x B drifts below the F-region ledge at sunset and simultaneous observations of large scale Kelvin-Helmholtz instabilities as seeds of spread-F -- Large expanses of quasi-coherent kilometer-scale vector wave observations (electric field and density) below 450 km and their discovery as a source of scintillations using the C/NOFS GPS -- Reversed zonal DC electric fields and simultaneous observations of afternoon counter electrojet-- Enhanced zonal DC electric fields at sunrise and full vector plasma flow continuity at the terminator-- Vector electric field and density irregularities extending to meter-scales within equatorial plasma depletions -- Intense electric field structure within the equatorial ionosphere at night in the absence of density depletions -- First observations of the westward equatorial electrojet in post-midnight ionosphere and possible association with downward meridional drifts-- Spaceborne Dst observations and implications for ring current asymmetries-- Vector electric and magnetic field Poynting flux within depletions and TIDs-- Measurements of ionospheric reflectance and Alfvenic waves associated with TIDs-- Observations of Alfven resonators -- ULF magnetic field structure within the nightime equatorial ionosphere-- Discovery of Schumman resonances in space and implications of ELF radiowave propagation-- Vector observations of 50-60 Hz powerline radiation without harmonics-- Ion cyclotron resonance absorption lines associated with ELF hiss and the identification of ambient ions-- Parallel electric fields associated with lower hybrid waves driven by thunderstorm lightning-related sferics-- Discovery of Z-mode radiation in the equatorial ionosphere including its associated with density depletions-- Optical lightning detector waveforms and electric field sferics observed up to altitudes of 800 km

Robert F Pfaff↗

DOE OSTI2026

Evidence of nonlinear coupling in the edge harmonic oscillation sustaining quiescent high confinement in a tokamak plasma

In a tokamak plasma, we measure nonlinear coupling between harmonics comprising a saturated edge oscillation, as predicted by nonlinear magnetohydrodynamic simulations. The coherent structure formed by this coupling suppresses bursty edge-localized modes that cause energy and particle loss. We also measure nonlinear coupling of the edge harmonic oscillation to core-resonant magnetic tearing modes. When tearing modes are present, the edge oscillation harmonics decouple, and edge-localized modes return. In this article, we conclude that nonlinear interaction with tearing modes interrupts the formation of the edge harmonic oscillation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2000

Acoustic Receptivity of a Blasius Boundary Layer with 2-D and Oblique Surface Waviness

An experimental investigation was conducted to examine acoustic receptivity and subsequent boundary-layer instability evolution for a Blasius boundary layer formed on a flat plate in the presence of two-dimensional (2-D) and oblique (3-D) surface waviness. The effect of the non-localized surface roughness geometry and acoustic wave amplitude on the receptivity process was explored. The surface roughness had a well defined wavenumber spectrum with fundamental wavenumber k (sub w). A planar downstream traveling acoustic wave was created to temporally excite the flow near the resonance frequency of an unstable eigenmode corresponding to k (sub ts) = k (sub w). The range of acoustic forcing levels, epsilon, and roughness heights, DELTA h, examined resulted in a linear dependence of receptivity coefficients; however, the larger values of the forcing combination epsilon dot DELTA h resulted in subsequent nonlinear development of the Tollmien-Schlichting (T-S) wave. This study provided the first experimental evidence of a marked increase in the receptivity coefficient with increasing obliqueness of the surface waviness in excellent agreement with theory. Detuning of the 2-D and oblique disturbances was investigated by varying the streamwise wall-roughness wavenumber a,, and measuring the T-S response. For the configuration where laminar-to-turbulent breakdown occurred, the breakdown process was found to be dominated by energy at the fundamental and harmonic frequencies, indicative of K-type breakdown.

King, Rudolph A.↗

DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Advanced spectrometric methods for characterizing bio-oils to enable refineries to reduce fuel carbon intensity during co-processing

A promising approach for supplementing petroleum-derived fuels to support reductions in green-house gas emissions is to convert abundant biomass feedstocks into renewable carbon-rich oils using pyrolysis. However, the resultant bio-oils contain various oxygenated compounds that can impart acidity, chemical and thermal instability, and immiscibility with petroleum derived fuels, necessitating further upgrading to derive fuel blendstocks. Co-processing bio-oils and petroleum-derived liquids in existing refineries is a potentially near-term, cost-effective approach for upgrading bio-oils while reducing refinery carbon intensities. However, one cause for hesitation in co-processing bio-oils is limited comprehensive characterization and speciation of the bio-oil components. Advanced analytical techniques are currently under investigation to enable identification of elusive species in bio-oils, enabling researchers to develop strategies to mitigate catalyst deactivation agents and contaminants. This review provides a brief overview of several analytical methods commonly used to analyze bio-oils and their limitations. In addition, advanced techniques currently under development are discussed to further elucidate bio-oil components that may limit its end use. This will help inform the technical and economic feasibility of co-processing bio-oils with petroleum-derived liquids, therefore, improving the overall downstream processes for biofuels blendstock production.

09 BIOMASS FUELS↗