SEARCH · Engineering Papers

Results for “Linear stability analysis”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Small-angle stability analysis of a linear control system for a high power communication satellite

A small angle stability analysis is presented for one particular configuration of a high power communication satellite having a linear control system. Both the central body and the solar array are treated as rigid bodies. The control system studied consists of three-axis control of the central body and one-axis control of the solar array rotation relative to the central body. The results yield preliminary indications of the relation of stability to satellite inertias and control gains.

Omalley, T. A.↗

NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗

DOE OSTI2023

Insights into ENSO Diversity from an Intermediate Coupled Model. Part I: Uniqueness and Sensitivity of the ENSO Mode

Abstract The basic dynamics of the spatiotemporal diversity for El Niño–Southern Oscillation (ENSO) has been the subject of extensive research and, while several hypotheses have been proposed, remains elusive. One promising line of studies suggests that the observed eastern Pacific (EP) and central Pacific (CP) ENSO may originate from two coexisting leading ENSO modes. We show that the coexistence of unstable EP-like and CP-like modes in these studies arises from contaminated linear stability analysis due to unnoticed numerical scheme caveats. In this two-part study, we further investigate the dynamics of ENSO diversity within a Cane–Zebiak-type model. We first revisit the linear stability issue to demonstrate that only one ENSO-like linear leading mode exists under realistic climate conditions. This single leading ENSO mode can be linked to either a coupled recharge-oscillator (RO) mode favored by the thermocline feedback or a wave-oscillator (WO) mode favored by the zonal advective feedback at the weak air–sea coupling end. Strong competition between the RO and WO modes for their prominence in shaping this ENSO mode into a generalized RO mode makes it sensitive to moderate changes in these two key feedbacks. Modulations of climate conditions yield corresponding modulations in spatial pattern, amplitude, and period associated with this ENSO mode. However, the ENSO behavior undergoing this linear climate condition modulations alone does not seem consistent with the observed ENSO diversity, suggesting the inadequacy of linear dynamics in explaining ENSO diversity. A nonlinear mechanism for ENSO diversity will be proposed and discussed in Part II.

Geng, Licheng↗

NASA NTRS1989

Experimental transition and boundary-layer stability analysis for a slotted swept laminar flow control airfoil

A swept, supercritical laminar flow control (LFC) airfoil designated NASA SCLFC(1)-0513F was tested at subsonic and transonic speeds in the NASA Langley eight-foot Transonic Pressure Tunnel. This paper examines Tollmien-Schlichting and crossflow disturbance amplification for this airfoil using the linear stability method. The design methodology using linear stability analysis is evaluated and the results of the incompressible and compressible methods are compared. Experimental data on the swept, supercritical LFC airfoil and reference wind tunnel and flight results are used to correlate and evaluate the N-factor method for transition prediction over a speed range M(infinity) from zero to one.

Harvey, William D.↗

DOE OSTI2024

Prediction of pellet mass thresholds for ELM triggering in low-collisionality, ITER-like discharges

Abstract In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2004

Cross Flow Effects on Glaze Ice Roughness Formation

The present study examines the impact of large-scale cross flow on the creation of ice roughness elements on the leading edge of a swept wing under glaze icing conditions. A three-dimensional triple-deck structure is developed to describe the local interaction of a 3 D air boundary layer with ice sheets and liquid films. A linear stability analysis is presented here. It is found that, as the sweep angle increases, the local icing instabilities enhance and the most linearly unstable modes are strictly three dimensional.

Tsao, Jen-Ching↗

DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗

NASA NTRS1993

Wave-turbulence dynamics in the stably stratified boundary layer

New data obtained at the Boulder Atmospheric Observatory are analyzed to obtain separation of wave, turbulence, and mean field necessary for a complete treatment of wave-turbulence interaction. The data were compared with a linear stability analysis of the background atmospheric state, showing good agreement between measured wave parameters (such as wavelength, period, and vector phase velocity) and the eigenvalues of the linear solution. The analysis of the budgets of wave heat flux and temperature variance revealed the essential role of wave-turbulence interaction in maintaining a large amplitude temperature wave and countergradient heat flux. A mechanism for the maintenance of turbulence by waves in strongly stratified boundary layers is described, which emphasizes that the time-mean Richardson number is an irrelevant parameter at such times.

Einaudi, F.↗

NASA NTRS1974

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.↗

DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2003

Stability Analysis of Flow Induced by the Traveling Magnetic Field

Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or.crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Mazuruk, Konstantin↗

NASA NTRS2003

Stability Analysis of Flow Induced by the Traveling Magnetic Field

Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Mazuruk, Konstantin↗

DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗

DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗

NASA NTRS1998

Dendritic Growth of Hard-Sphere Crystals

Recent observations of the disorder-order transition for colloidal hard spheres under microgravity revealed dendritic crystallites roughly 1-2 mm in size for samples in the coexistence region of the phase diagram. Order-of-magnitude estimates rationalize the absence of large or dendritic crystals under normal gravity and their stability to annealing in microgravity. A linear stability analysis of the Ackerson and Schaetzel model for crystallization of hard spheres establishes the domain of instability for diffusion-limited growth at small supersaturations. The relationship between hard-sphere and molecular crystal growth is established and exploited to relate the predicted linear instability to the well-developed dendrites observed.

Russel, W. B.↗

NASA NTRS2001

Hypersonic Boundary Layer Instability Over a Corner

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Balakumar, Ponnampalam↗

DOE OSTI2020

Stabilizing effect of enhanced resistivity on peeling-ballooning instabilities on EAST

Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.

EAST tokamak↗

DOE OSTI2020

Wobbly discs – corrugations seen in the dust lanes of edge-on galaxies

ABSTRACT We report the detection of small-scale bending waves, also known as corrugations, in the dust lanes of five nearby edge-on disc galaxies. This phenomenon, where the disc mid-plane bends to become wavy, just as in warps but on a smaller scale, is seen here for the first time, in the dust lanes running across the discs. Because they are seen in absorption, this feature must be present in the dust disc in the outskirts of these galaxies. We enhance the visibility of these features using unsharp masking, trace the dust mid-plane across the disc, measure the corrugation amplitude by eye and the corrugation wavelength using Fourier analysis. The corrugation amplitude is found to be in the range of 70–300pc and the wavelengths lie between 1 and 5 kpc. In this limited sample, we find that the amplitude of the corrugations tends to be larger for lower mass galaxies, whereas the wavelength of corrugation does not seem to depend on the mass of host galaxies. Linear stability analysis is performed to find out the dynamical state of these dust discs. Based on WKB analysis, we find that the dust corrugations in about half of our sample are stable. Further analysis, on a larger sample would be useful to strengthen the above results.

Narayan, Chaitra A.↗