Engineering Papers⌕ Search

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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Real-Gas Effects on Binary Mixing Layers

This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Okong'o, Nora↗
NASA NTRS2007

Analysis of Crossflow Transition Flight Experiment aboard the Pegasus Launch Vehicle

The Pegasus wing-glove flight experiment was designed to provide crossflow transition data at high Mach numbers, specifically to help validate stability based predictions for transition onset in a flight environment. This paper provides an analysis of the flight experiment, with emphasis on computational results for crossflow disturbances and the correlation of disturbance growth factors with in-flight transition locations via the e(sup N) method. Implications of the flight data for attachment line stability are also examined. Analysis of the thermocouple data reveals that transition (from turbulent to laminar flow) was first detected during the ascending flight of the rocket when the free stream Mach number exceeded about 4. Therefore, computations have been performed for flight Mach numbers of 4.13, 4.35, 4.56 and 4.99. Due to continually decreasing unit Reynolds number at higher altitudes, the entire wing-glove boundary layer became laminar at the highest flight Mach number computed above. In contrast, the boundary layer flow over the inboard tile region remained transitional up to and somewhat beyond the time of laminarization over the instrumented glove region. Linear stability predictions confirmed that the tile boundary layer is indeed more unstable to crossflow disturbances than the much colder stainless steel glove boundary layer. The transition locations based on thermocouple data from both the glove and the tile regions are found to correlate with stationary-crossflow N-factors within the range of 7 to 12.4 and with traveling mode N-factors between 7.6 and 14.1. Data from the thermocouples and hot film sensors indicates that transition from turbulent to laminar flow (i.e., laminarization) at a fixed point over the glove is generally completed within a flight time interval of 3 seconds. However, the times at which transition begins and ends as inferred from the hot film sensors are found to differ by about 2 seconds from the corresponding estimates based on the thermocouple data.

Malik, Mujeeb R.↗
DOE OSTI2022

The role of deviatoric stress and dislocations on the α to ω phase transformation in Ti

Under extreme conditions, α-Ti becomes unstable and transforms either into β-Ti at high temperature or into ω-Ti at high pressure. In what concerns the α to ω phase transformation (PT), there has been a wide range of experimentally reported transition pressures from approximately 2 to 15 GPa at room temperature. Deviatoric stresses and internal defects are often assumed to be the root cause of this variation. Here, in this study, these postulates are revisited using both continuum mechanics and molecular dynamics (MD) simulations. First, a simple continuum model, assuming linear elasticity and isotropic plasticity, is developed to describe the effects of applied stress and dislocations on the stability of an ω nucleus in an infinite α domain. Second, a new MD simulation method is developed to generate an ω nucleus in the α domain utilizing the displacement field identified from the topological analysis. Results from MD simulations show that despite the fact that phase diagrams typically delineate the limits between two phases in terms of only P and T, deviatoric stress promotes the α to ω phase transformation by reducing the critical radius above which an ω nucleus is stable. Furthermore, the required deviatoric stress to nucleate and stabilize a nanoscale ω nucleus is likely emanating from the internal stress of defects such as dislocations. The MD-informed micromechanics models are used to identify favorable configurations where dislocations help favor the α to ω transformation. These configurations show that the interaction with a basal or prismatic dislocation reduces the critical radius of a ω nucleus by about 10 or 16 %, respectively. In addition, prismatic edge dislocations are found to promote the growth of ω nucleus when interacting with the ($\bar{1}$$\bar{1}$20) α //(0001) ω interface. Importantly, a simple model of the arrival of dislocation at an ω nucleus suggests that PT does not necessarily require a pile-up to be present but could alternatively be mediated by a constant rapid flow of dislocations.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗
NASA NTRS2011

Data Analysis for the NASA/Boeing Hybrid Laminar Flow Control Crossflow Experiment

The Hybrid-Laminar Flow Control (HLFC) Crossflow Experiment, completed in 1995. generated a large database of boundary layer stability and transition data that was only partially analyzed before data analysis was abruptly ended in the late 1990's. Renewed interest in laminar flow technologies prompted additional data analysis, to integrate all data, including some post-test roughness and porosity measurements. The objective is to gain new insights into the effects of suction on boundary layer stability. A number of challenges were encountered during the data analysis, and their solutions are discussed in detail. They include the effect of the probe vibration, the effect of the time-varying surface temperature on traveling crossflow instabilities, and the effect of the stationary crossflow modes on the approximation of wall location. Despite the low turbulence intensity of the wind tunnel (0.01 to 0.02%), traveling crosflow disturbances were present in the data, in some cases at amplitudes up to 1% of the freestream velocity. However, the data suggests that transition was dominated by stationary crossflow. Traveling crossflow results and stationary data in the presence of suction are compared with linear parabolized stability equations results as a way of testing the quality of the results.

Eppink, Jenna L.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRS1973

The finite element method in shell stability analysis.

A development of the finite element method for thin shell instability analysis is presented, covering three principal aspects: (1) representation of shell geometry, (2) representation of element behavior, and (3) algorithmic tools for solution of the large-order systems of nonlinear algebraic equations which characterize various phases of shell instability. Two shell elements are described, an arbitrary quadrilateral and a triangle, and numerical results are presented for two widely-employed comparison problems for linear (stable) analysis. Two shell problems which include instability effects are also solved.-

Gallagher, R. H.↗
NASA NTRS1995

Euler/Navier-Stokes flow computations on flexible configurations for stability analysis

Longitudinal dynamic stability derivatives required for design of aircraft are computed by using the state-of-the-art numerical methods for wing-body configurations. The flow is modeled using the Euler/Navier-Stokes equations with turbulence models and solved using an efficient finite-difference scheme suitable for patched structured grids. Computations are made at a flow regime that is beyond the limits of the current linear methods mostly used for computing stability derivatives. Flow conditions include shockwaves and viscous dominated vortical flows. Effect of Mach number and angle-of-attack on stability derivatives are demonstrated for a typical wing-body configuration. For the same configuration the effects of wing flexibility on the magnitude and phase angles of stability derivatives are also demonstrated.

Guruswamy, G.↗
DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗
NASA NTRS1995

Euler/Navier-Stokes Flow Computations on Flexible Configurations for Stability Analysis

Abstract Longitudinal dynamic stability derivatives required for design of aircraft are computed by using the state-of-the-art numerical methods for wing-body configurations. The flow is modeled using the Euler/Navier-Stokes equations with turbulence models and solved using an efficient finite-difference scheme suitable for patched structured grids. Computations are made at a flow regime that is beyond the limits of the current linear methods mostly used for computing stability derivatives. Flow conditions include shock-waves and viscous dominated vortical flows. Effect of Mach number and angle-of-attack on stability derivatives are demonstrated for a typical wing-body configuration. For the same configuration the effects of wing flexibility on the magnitude and phase angles of stability derivatives are also demonstrated.

Guruswamy, Guru P.↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring nonlinear Rashba effect and spin Hall conductivity in Janus MXenes W 2 ⁢CO ⁢𝑋 (𝑋=S, Se, Te)

Rashba spin-orbit coupling (RSOC) facilitates spin manipulation without relying on an external magnetic field, opening up exciting possibilities for advanced spintronic devices. In this paper, we examine the effects of crystal momentum (𝑘) nonlinearity and anisotropy on the conventional Rashba effect, with a particular focus on their impact on the spin Hall conductivity (SHC) in a newly predicted family of 2D Janus materials, W 2 ⁢CO⁢𝑋 (𝑋 =S, Se, Te). Using first-principles density functional theory calculations, we confirm the dynamical and mechanical stability of the studied 2D materials. Strikingly, this materials family exhibits pronounced nonlinear Rashba spin splitting at the Γ point of Brillouin zone near the Fermi level, which cannot be adequately described by the linear-𝑘 Rashba model. Therefore, third-order momentum contributions (𝑘 3 ) must be incorporated into the Rashba Hamiltonian. Our analysis reveals that among the studied systems, W 2 ⁢COS exhibits the highest 𝑘 3 contribution of −45.9 eV Å 3 , despite having the lowest linear Rashba constant. Here, a detailed analysis of electronic structure reveals topologically nontrivial behavior in these 2D materials, yielding sizable SHC that is primarily governed by the nonlinear Rashba effect. Notably, these materials also exhibit large spin Hall angle (0.018–2.5 at E 𝐹 ), which is comparable to that of in bulk topological insulators like Bi 2 ⁢Se 3 and Bi 2 ⁢Te 3 , and surpassing those in narrow bandgap bulk semiconductors GeTe and SnTe, as well as heavy metals such as Pt. Sizable SHC, large spin Hall angles, and the ability to tune SHC via electric fields without altering the topological properties, rooted in the crystal field splitting, underscore the potential of these materials for spintronic applications.

Electronic structure↗
NASA NTRS1981

Advanced stability analysis for laminar flow control

Five classes of problems are addressed: (1) the extension of the SALLY stability analysis code to the full eighth order compressible stability equations for three dimensional boundary layer; (2) a comparison of methods for prediction of transition using SALLY for incompressible flows; (3) a study of instability and transition in rotating disk flows in which the effects of Coriolis forces and streamline curvature are included; (4) a new linear three dimensional instability mechanism that predicts Reynolds numbers for transition to turbulence in planar shear flows in good agreement with experiment; and (5) a study of the stability of finite amplitude disturbances in axisymmetric pipe flow showing the stability of this flow to all nonlinear axisymmetric disturbances.

Orszag, S. A.↗
NASA NTRS1981

Analytical techniques for the analysis of stall/spin flight test data

Analytical techniques for the analysis of stall/spin flight test data are reviewed by discussing (1) certain special flight instrumentation issues, (2) the mathematical modeling techniques, and (3) the analysis of post stall and spinning flight of general aviation airplanes. The angles of attack, sideslip, roll, pitch, and yaw are derived from measurements of angular velocity and linear acceleration. The key to the success of this approach is to simultaneously estimate both the biases in the instrumentation and the initial conditions. Techniques for determining stability derivatives from flight data are applied to angles of attack too high for stabilized flight. This practice greatly expands the range over which aerodynamic characteristics can be determined from flight test. Nonlinear terms in certain aerodynamic functions are shown to be valid by comparing them with the trends of results at different angles of attack. A very old technique of studying spins is extended and applied to some modern light airplanes. Airplanes for which the wing provides the dominant moments during spins, offer the possibility of linking spin characteristics to longitudinal data.

Taylor, L. W., Jr.↗
NASA NTRS1994

The kinetic response of a stochastic plasma to low frequency perturbations

Following suggestion that substorm breakup might be caused by an interchange or ballooning instability, several magnetohydrodynamic (MHD) and gyro-kinetic stability analysis have been performed for plasma sheet magnetic field geometries. However, the stochastic ion dynamics in the highly stressed, thin high-beta near earth plasma sheet violate the locality requirements of MHD and the invarience of the magnetic moment, mu, required by gyro-kinetic theories. In this paper we develop a new linear Vlasov kinetic theory (for low frequency modes omega less than omega(sub b), omega(sub b) being the bounce frequency) which includes the dynamics of stochastic ions.

Hurricane, Omar A.↗
NASA NTRS1996

Aeroservoelastic Modeling and Validation of a Thrust-Vectoring F/A-18 Aircraft

An F/A-18 aircraft was modified to perform flight research at high angles of attack (AOA) using thrust vectoring and advanced control law concepts for agility and performance enhancement and to provide a testbed for the computational fluid dynamics community. Aeroservoelastic (ASE) characteristics had changed considerably from the baseline F/A-18 aircraft because of structural and flight control system amendments, so analyses and flight tests were performed to verify structural stability at high AOA. Detailed actuator models that consider the physical, electrical, and mechanical elements of actuation and its installation on the airframe were employed in the analysis to accurately model the coupled dynamics of the airframe, actuators, and control surfaces. This report describes the ASE modeling procedure, ground test validation, flight test clearance, and test data analysis for the reconfigured F/A-18 aircraft. Multivariable ASE stability margins are calculated from flight data and compared to analytical margins. Because this thrust-vectoring configuration uses exhaust vanes to vector the thrust, the modeling issues are nearly identical for modem multi-axis nozzle configurations. This report correlates analysis results with flight test data and makes observations concerning the application of the linear predictions to thrust-vectoring and high-AOA flight.

Brenner, Martin J.↗
NASA NTRS1997

PLATSIM: A Simulation and Analysis Package for Large-Order Flexible Systems

The software package PLATSIM provides efficient time and frequency domain analysis of large-order generic space platforms. PLATSIM can perform open-loop analysis or closed-loop analysis with linear or nonlinear control system models. PLATSIM exploits the particular form of sparsity of the plant matrices for very efficient linear and nonlinear time domain analysis, as well as frequency domain analysis. A new, original algorithm for the efficient computation of open-loop and closed-loop frequency response functions for large-order systems has been developed and is implemented within the package. Furthermore, a novel and efficient jitter analysis routine which determines jitter and stability values from time simulations in a very efficient manner has been developed and is incorporated in the PLATSIM package. In the time domain analysis, PLATSIM simulates the response of the space platform to disturbances and calculates the jitter and stability values from the response time histories. In the frequency domain analysis, PLATSIM calculates frequency response function matrices and provides the corresponding Bode plots. The PLATSIM software package is written in MATLAB script language. A graphical user interface is developed in the package to provide convenient access to its various features.

Maghami, Peiman G.↗