SEARCH · Engineering Papers

Results for “instability criteria”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Random Variables with Moment-Matching Staircase Density Functions

This paper proposes a family of random variables for uncertainty modeling. The variables of interest have a bounded support set, and prescribed values for the first four moments. We present the feasibility conditions for the existence of any of such variables, and propose a class of variables that conforms to such constraints. This class is called staircase because the density of its members is a piecewise constant function. Convex optimization is used to calculate their distributions according to several optimality criteria, including maximal entropy and maximal log-likelihood. The flexibility and efficiency of staircases enable modeling phenomena having a possibly skewed and/or multimodal response at a low computational cost. Furthermore, we provide a means to account for the uncertainty in the distribution caused by estimating staircases from data. These ideas are illustrated by generating empirical staircase predictor models. We consider the case in which the predictor matches the sample moments exactly (a setting applicable to large datasets), as well as the case in which the predictor accounts for the sampling error in such moments (a setting applicable to sparse datasets). A predictor model for the dynamics of an aeroelastic airfoil subject to flutter instability is used as an example. The resulting predictor not only describes the system's response accurately, but also enables carrying out a risk analysis for safe flight.

Luis G. Crespo↗

DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1986

Parametric analytical study of instability-related delamination growth

The effect of various parameters on instability-related delamination growth was studied analytially. The configuration studied consisted of a thick composite laminate with a single through-width delamination located near one surface. Both mechanical and thermal loads were considered. All conclusions were based on the assumption that GI and GII govern delamination growth. An approximate superposition stress analysis was developed which gives closed form expressions for GI and GII. The simplicity of the analysis permitted examination of numerous configurations. Both GI and GII were found to be very sensitive to delamination length and location through the thickness. The magnitude of GI was also very sensitive to initial imperfections, which might be the results of an inclusion of finite thickness in the delamination. Critical loads for delamination growth were calculated based on three growth criteria. Large differences in the predictions highlight the need for a verified mixed-mode delamination growth criterion.

Whitcomb, J. D.↗

DOE OSTI2021

Electronic Origin of T c in Bulk and Monolayer FeSe

FeSe is classed as a Hund’s metal, with a multiplicity of d bands near the Fermi level. Correlations in Hund’s metals mostly originate from the exchange parameter J, which can drive a strong orbital selectivity in the correlations. The Fe-chalcogens are the most strongly correlated of the Fe-based superconductors, with d xy the most correlated orbital. Yet little is understood whether and how such correlations directly affect the superconducting instability in Hund’s systems. By applying a recently developed ab initio theory, we show explicitly the connections between correlations in d xy and the superconducting critical temperature T c . Starting from the ab initio results as a reference, we consider various kinds of excursions in parameter space around the reference to determine what controls T c . We show small excursions in J can cause colossal changes in T c . Additionally we consider changes in hopping by varying the Fe-Se bond length in bulk, in the free standing monolayer M-FeSe, and M-FeSe on a SrTiO 3 substrate (M-FeSe/STO). The twin conditions of proximity of the dxy state to the Fermi energy, and the strength of J emerge as the primary criteria for incoherent spectral response and enhanced single- and two-particle scattering that in turn controls T c . Using c-RPA, we show further that FeSe in monolayer form (M-FeSe) provides a natural mechanism to enhance J. We explain why M-FeSe/STO has a high T c , whereas M-FeSe in isolation should not. Our study opens a paradigm for a unified understanding what controls Tc in bulk, layers, and interfaces of Hund’s metals by hole pocket and electron screening cloud engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1983

Assessment of lubricated contacts: Mechanisms of scuffing and scoring

Scoring and scuffing are defined as two distinct but related forms of failure of hardened ferrous components lubricated by liquids. Experimental observations of these phenomena are described, and criteria for scoring and scuffing are discussed. The mechanisms proposed by various authors to explain these observations are enumerated. The view presented here is that, under conditions yet to be defined, scoring is a gradual breakdown of the lubrication of interacting asperities, such lubrication being boundary or microelastohydrodynamic in nature, or a mixture of the two. The final scuffing stage represents a breakdown of the main elastohydrodynamic system, but this may be triggered by the deterioration in surface topography as a result of scoring. An extension of published theoretical treatment of elastohydrodynamic breakdown is proposed, and a critical experiment to assess the importance of edge effects in band contacts is suggested. The practical implications of the proposed mechanism are reviewed, and recommendations for further work are made. A possible thermal instability mechanism for the breakdown of boundary lubrication is outlined.

Dyson, A.↗

NASA NTRS1993

Solar Nebula Magnetohydrodynamic Dynamos: Kinematic Theory, Dynamical Constraints, and Magnetic Transport of Angular Momentum

A hydromagnetic dynamo provides the best mechanism for contemporaneously producing magnetic fields in a turbulent solar nebula. We investigate the solar nebula in the framework of a steady-state accretion disk model and establish the criteria for a viable nebular dynamo. We have found that typically a magnetic gap exists in the nebula, the region where the degree of ionization is too small for the magnetic field to couple to the gas. The location and width of this gap depend on the particular model; the supposition is that gaps cover different parts of the nebula at different evolutionary stages. We have found, from several dynamical constraints, that the generated magnetic field is likely to saturate at a strength equal to equipartition with the kinetic energy of turbulence. Maxwell stress arising from a large-scale magnetic field may significantly influence nebular structure, and Maxwell stress due to small-scale fields can actually dominate other stresses in the inner parts of the nebula. We also argue that the bulk of nebular gas, within the scale height from the midplane, is stable against Balbus-Hawley instability.

Stepinski, Tomasz F.↗

NASA NTRS2010

A Parametric Study on Using Active Debris Removal to Stabilize the Future LEO Debris Environment

Recent analyses of the instability of the orbital debris population in the low Earth orbit (LEO) region and the collision between Iridium 33 and Cosmos 2251 have reignited the interest in using active debris removal (ADR) to remediate the environment. There are; however, monumental technical, resources, operational, legal, and political challenges in making economically viable ADR a reality. Before a consensus on the need for ADR can be reached, a careful analysis of the effectiveness of ADR must be conducted. The goal is to demonstrate the feasibility of using ADR to preserve the future environment and to guide its implementation to maximize the benefit-cost ratio. This paper describes a comprehensive sensitivity study on using ADR to stabilize the future LEO debris environment. The NASA long-term, orbital debris evolutionary model, LEGEND, is used to quantify the effects of many key parameters. These parameters include (1) the starting epoch of ADR implementation, (2) various target selection criteria, (3) the benefits of collision avoidance maneuvers, (4) the consequence of targeting specific inclination or altitude regimes, (5) the consequence of targeting specific classes of vehicles, and (6) the timescale of removal. Additional analyses on the importance of postmission disposal and how future launches might affect the requirements to stabilize the environment are also included.

Liou, J.C.↗

DOE OSTI2021

Beam Envelope Stability in an Advanced Linear Induction Accelerator

A new linear induction accelerator (LIA) is under development for multipulse flash radiography. Because it has substantially more cells than present LIAs, higher magnetic focusing fields are needed to suppress beam breakup (BBU). It is, therefore, more susceptible to the parametric beam envelope instability (PEI), which has an instability threshold that has usually been typified by the vacuum phase advance per cell exceeding some large fraction of π . Here we derive a threshold criterion for PEI that depends not only on the magnetic field, but also on the beam space charge and emittance. A tune designed to suppress BBU in Scorpius is shown to be stable to the PEI according to this criterion, and also by the lack of emittance growth in particle-in-cell (PIC) code simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1966

Storable Propellant Combustion Instability Program at Lewis Research Center

An experimental program at Lewis Research Center investigates acoustical-mode combustion instability in liquid propellant rockets. One phase of this program is concerned with nitrogen tetroxide and a 50-50 fuel blend of hydrazine-UDMH. Effects on combustor stability and performance by variations in injection velocities, impingement angle and distance, amd thrust per element of triplet injectors were studied in a 10.77-inch- diameter cylindrical combustor at chamber pressures of 100 and 300 psia. Nominal thrust levels were 6700 and 20 000 pounds. Stability rating was accomplished by the use of various size RDX explosive charges to generate tangential pressure disturbances. The injection velocity effect correlated by using a ratio V(sub o)/V(sub f). Stability correlation with V(sub o)/V(sub f) resulted in a sharply humped curve with maximum stability at a ratio of 1.2. Performance decreased with increase in V(sub o)/V(sub f) after a maximum at 0.7, which agreed fairly well with the "uniform mixture ratio distribution" criteria of the Jet Propulsion Laboratory. Variation of impingement angle from 38 deg to 120 deg indicated that performance correlated with the absolute velocity ratio V(sub o)/V(sub f), but stability correlation appeared to be best with a velocity ratio using the radial component of the fuel injection velocity V(sub o)/V(sub fr). Impingement distances from 0.5 to 1.0 inch appeared to have only slight effect on stability and performance. Investigation of injector thrust per element for V(sub o)/V(sub fr) ranging from 0.8 to 2.8 yielded convex maxima curves. Acoustic liners were also tested to improve a liner design computer program, to obtain design criteria, and to develop flight-type liners. A marginally stable and spontaneously unstable combustor were used to evaluate linear configurations. All liners were effective in the marginally stable configuration, and only liners with large theoretical absorptivity were successful in the spontaneously unstable combustor.

Tabata, William K.↗

DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

DOE OSTI2022

Database-wide hazard modelling of the onset of DIII-D tearing modes with field features

The rate of onset (hazard) of tearing modes is modelled probabilistically using statistical learning algorithms. Axisymmetric energy-density equilibrium fields are taken as raw high-dimensional input features which are reduced with principal component analysis. Signal processing of non-axisymmetric magnetics fluctuation array data provides the target information from which to learn. Model selection, visualization and calibration assessment procedures are detailed. Here, the analysis is deployed at large scale across the DIII-D tokamak database. Standard model selection criteria suggest that the energy-density post-processed feature is a better choice for modelling the onset rate compared to the non-processed equilibrium reconstruction solution. Two example applications of the learned rate function are demonstrated: (i) proximity-to-onset discharge monitoring and (ii) database analysis showing an (expected) observational global trend that the general hazard increases as a plasma performance metric increases. An important connection between the hazard function and its use as a conditional probability generator is reviewed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2010

Effects of Temperature on the Performance and Stability of Recent COTS Silicon Oscillators

Silicon oscillators have lately emerged to serve as potential replacement for crystal and ceramic resonators to provide timing and clock signals in electronic systems. These semiconductor-based devices, including those that are based on MEMS technology, are reported to be resistant to vibration and shock (an important criteria for systems to be deployed in space), immune to EMI, consume very low current, require few or no external components, and cover a wide range of frequency for analog and digital circuits. In this work, the performance of five recently-developed COTS silicon oscillator chips from different manufacturers was determined within a temperature range that extended beyond the individual specified range of operation. In addition, restart capability at extreme temperatures, i.e. power switched on while the device was soaking at extreme (hot or cold) temperature, and the effects of thermal cycling under a wide temperature range on the operation of these silicon oscillators were also investigated. Performance characterization of each oscillator was obtained in terms of its output frequency, duty cycle, rise and fall times, and supply current at specific test temperatures. The five different oscillators tested operated beyond their specified temperature region, with some displaying excellent stability throughout the whole test temperature range. Others experienced some instability at certain temperature test points as evidenced by fluctuation in the output frequency. Recovery from temperature-induced changes took place when excessive temperatures were removed. It should also be pointed out that all oscillators were able to restart at the extreme test temperatures and to withstand the limited thermal cycling without undergoing any significant changes in their characteristics. In addition, no physical damage was observed in the packaging material of any of these silicon oscillators due to extreme temperature exposure and thermal cycling. It is recommended that additional and more comprehensive testing under long term cycling be carried out to fully establish the reliability of these devices and to determine their suitability for use in space exploration missions under extreme temperature conditions.

Patterson, Richard L.↗

NASA NTRS2018

Outer Solar System Material in Inner Solar System Regolith Breccias

There is excellent evidence that a dynamical instability in the early solar system led to gravitational interactions between the giant planets and trans-Neptunian objects (TNOs). Giant planetary migration triggered by the instability dispersed a disk of primordial TNOs and created a number of small body reservoirs (e.g. the Kuiper Belt, scattered disk, irregular satellites, and the Jupiter/Neptune Trojan populations). It also injected numerous bodies into the main asteroid belt, where modeling shows they can successfully reproduce the observed P and D-type asteroid populations. During the injection era and after implantation, some of these “main belt TNOs” would have collided with S-class asteroids. Some of this material may have survived as a component of asteroid regolith breccias. Thus, we have been searching for evidence of these impact events in the form of carbonaceous xenoliths in brecciated ordinary chondrites. These xenoliths would have experienced a wide range of impact velocities, and therefore we should expect to see everything between relatively unaltered material to completely shock-melted lithologies. This material might also be different from the carbonaceous chondrites that represent standard C-complex asteroids. A goal of this research is to define useful criteria for distinguishing between these two classes of materials, including O, Cr, N and C isotopes, petrographic characteristics, and chronology.

Zolensky, M.↗

DOE OSTI2023

Use of intrinsic hysteresis for the active control of internal transport barriers in magnetically confined fusion plasmas

In magnetically confined fusion devices, control of internal transport barriers (ITBs) is important both to enhance and suppress the turbulent transport to improve confinement control. Barrier control should allow for the improvement of confinement to aid in achieving the needed fusion criteria while also permitting the degradation of confinement to control profiles and clean the device by moving out the impurities accumulated near the core. In this work, we present a novel control scenario that takes advantage of the hysteresis intrinsic to transport barriers to easily cycle through enhanced and degraded confinement regimes. The control scenario is illustrated using a five-field simplified transport model for an ITB using typical parameters of a neutral beam injection-heated DIII-D tokamak discharge. Pellets and ion cyclotron resonance frequency power are used as control knobs for this active control scenario. These knobs adequately modify at will the local gradients and, therefore, the growth rates and shearing rates, allowing for easy and efficient control of the barrier by taking advantage of the barrier hysteresis. The result is a control cycle that could be operated with a relatively small amount of power in high performance regimes which, nowadays, typically require large power to control. It may also have advantages to avoid, or at least ameliorate, the appearance of magnetohydrodynamic instabilities in the barrier region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Co-Occurring Atmospheric Features and Their Contributions to Precipitation Extremes

Object-based identification algorithms for atmospheric features are commonly utilized to attribute global precipitation. This study employs a systematic approach to examine feature co-occurrences and their relationships to mean and extreme precipitation. Four features are identified using existing data sets for atmospheric rivers (ARs), mesoscale convective systems (MCSs), low-pressure systems (LPSs), and fronts (FTs). Often, a single atmospheric phenomenon satisfies the criteria set by multiple feature identification algorithms, yielding an association between precipitation and multiple features. Over the extra-tropics, the number of features attributed to a single event typically increases with precipitation intensity. Over two-thirds of the precipitation is from co-occurring features, with a considerable fraction related to AR-FT co-occurrences. Over the tropics, about one-quarter of precipitation is associated with co-occurring features, with LPS-MCS co-occurrences contributing substantially in monsoon regions. MCSs are the leading single-feature contributors over tropical land and oceans. In the extra-tropics, FTs, ARs, and their co-occurrences account for over half of the total precipitation over oceans. AR-FT-MCS and FT-MCS co-occurrences contribute to extremes (precipitation exceeding the 95th percentile) over both oceans (over 30%) and land (over 20%). Any combination of features involving MCSs shows a larger contribution to high percentiles of precipitation intensity. A case analysis indicates that AR-FT-MCS co-occurrences exhibit convective instability and deep vertical motion, suggesting that the feature trackers and reanalysis are capturing physics relevant to both convective and frontal systems. The results here emphasize the need for simultaneous identifications of multiple features when attributing precipitation to atmospheric phenomena.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

NASA NTRS2019

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗