Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,135 records · Page 63

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↗

Vertical Navigation Control Laws and Logic for the Next Generation Air Transportation System

A vertical navigation (VNAV) outer-loop control system was developed to capture and track the vertical path segments of energy-efficient trajectories that are being developed for high-density operations in the evolving Next Generation Air Transportation System (NextGen). The VNAV control system has a speed-on-elevator control mode to pitch the aircraft for tracking a calibrated airspeed (CAS) or Mach number profile and a path control mode for tracking the VNAV altitude profile. Mode control logic was developed for engagement of either the speed or path control modes. The control system will level the aircraft to prevent it from flying through a constraint altitude. A stability analysis was performed that showed that the gain and phase margins of the VNAV control system significantly exceeded the design gain and phase margins. The system performance was assessed using a six-deg-of-freedom non-linear transport aircraft simulation and the performance is illustrated with time-history plots of recorded simulation data.

Hueschen, Richard M.↗

The Theory of Pulsating Flow in Conical Nozzles

A knowledge of the dynamic characteristics of nozzles and orifices is important in many control and stability analyses of engineering devices. It is usual to assume that the instantaneous flow-rate, for a given set of inlet conditions and outlet pressure, is the same as the nontransient value for the same operating conditions. Recently, in connection with the stability analysis of an externally pressurized thrust bearing, the validity of this assumption was questioned. The analysis presented in this paper was undertaken to provide an answer. The present analysis applies to any fluid, liquid, or gas flowing into a simple conical nozzle. The amplitude and phase of the mass-flux response to a sinusoidally time-varying pressure fluctuation at the nozzle exit are determined. An approximate formula is given for these quantities in terms of the nozzle throat area, the solid angle subtended by the cone, the velocity of the fluid at the nozzle throat, the acoustic velocity at the throat, and the frequency of the pressure fluctuation.

Elrod, H. G., Jr.↗

Control Systems with Normalized and Covariance Adaptation by Optimal Control Modification

Disclosed is a novel adaptive control method and system called optimal control modification with normalization and covariance adjustment. The invention addresses specifically to current challenges with adaptive control in these areas: 1) persistent excitation, 2) complex nonlinear input-output mapping, 3) large inputs and persistent learning, and 4) the lack of stability analysis tools for certification. The invention has been subject to many simulations and flight testing. The results substantiate the effectiveness of the invention and demonstrate the technical feasibility for use in modern aircraft flight control systems.

Nguyen, Nhan T.↗

Observations of Large-Amplitude, Parallel, Electrostatic Waves Associated with the Kelvin-Helmholtz Instability by the Magnetospheric Multiscale Mission

On 8 September 2015, the four Magnetospheric Multiscale spacecraft encountered a Kelvin-Helmholtz unstable magnetopause near the dusk flank. The spacecraft observed periodic compressed current sheets, between which the plasma was turbulent. We present observations of large-amplitude (up to 100 mVm) oscillations in the electric field. Because these oscillations are purely parallel to the background magnetic field, electrostatic, and below the ion plasma frequency, they are likely to be ion acoustic-like waves. These waves are observed in a turbulent plasma where multiple particle populations are intermittently mixed, including cold electrons with energies less than 10 eV. Stability analysis suggests a cold electron component is necessary for wave growth.

Wilder, F. D.↗

Correlation of Full-Scale Isolated Proprotor Performance and Loads

A full-scale isolated proprotor test is currently being conducted in the USAF National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel at NASA Ames. The test article is a 3-bladed research rotor derived from the right-hand rotor of the AW609; this rotor was manufactured by Bell Helicopter under contract to NASA. In this paper, this research rotor is referred to as "699". The test, nearly completed, is an integral part of the initial checkout test of the newly developed Tiltrotor Test Rig (TTR), whose purpose is to test advanced, full-scale proprotors in the NFAC. Figure 1 shows the TTR/699 installed in the 40- by 80-Foot test section. The TTR rotor axis is horizontal and the rig rotates in yaw on the wind tunnel turntable for conversion (transition) and helicopter mode testing. To date, a substantial amount of wind tunnel test data has already been acquired. The completed operational conditions include hover, airplane mode (cruise, wind tunnel airspeed V=61 to 267 knots), and the helicopter and conversion conditions (with a comprehensive sweep of the TTR yaw angle ranging, to date, from 90-deg yaw helicopter mode to 30-deg yaw conversion mode, at varying airspeeds). This 699 proprotor performance and loads correlation study uses these newly acquired wind tunnel test data. This paper represents the third analytical study, coming after two earlier analytical studies on the TTR/699; that is, a 2018 paper on pre-test predictions of 699 performance and loads, Ref. 1, and an upcoming January 2019 paper on aeroelastic stability analysis of the TTR/699 installed in the 40- by 80-Foot Wind Tunnel, Ref. 2. Reference 8 will present an overview of the entire TTR/699 test program. For completeness, Ref. 3 addresses the development and initial testing of the TTR. Background information on the TTR effort at NASA Ames can be found at the Aeromechanics website: https://rotorcraft.arc.nasa.gov/Research/Facilities/ttr.html. To the authors' knowledge, the full-scale results presented in this paper are the first of their kind. A literature survey brought up several existing correlation studies, but these were either based on small-scale test data (for example, the studies performed by the University of Maryland) or full-scale aircraft flight test data (for example, flight tests conducted by Bell Helicopter). Separately, the 2009 NASA study involving the JVX rotor is relevant (see Ref.4). The JVX is closely similar to the 699 in size and aerodynamics, and is accordingly a good reference for performance calculations. In Ref. 1 (as mentioned above), pre-test reality checks of the current analytical model were made by comparing JVX and 699 predictions in hover and forward flight (airplane mode).

Isolated Proprotor Performance↗

Periodic Orbit-Attitude Solutions in the Planar Elliptic Restricted Three-Body Problem

The pitch motion of a spacecraft in the planar elliptic restricted three-body system is studied. Previous studies laid the foundation for spacecraft stability analysis with a small perturbation to the zero pitch motion. In this study, a cell mapping approach that combines analytical and numerical techniques is used to study the global behavior of the full nonlinear spacecraft attitude in which coupling between orbital dynamics and attitude occurs. The effect of gravity gradient torques, orbital eccentricity, and the spacecraft configuration at different Lagrangian points is analyzed. Multiple-period periodic solutions and invariant surfaces are presented for different cases. Reference trajectories around the Lagrangian points are also considered to study coupled dynamics.

Koh, Dayung↗

MRAC with Time-Varying Reference Model from Convex Combination

An extension of model reference adaptive control is presented that accommodates use of a time-varying reference model. Specifically, the reference model is taken to be a time-varying convex combination of two linear, time-invariant models. The design is intended to act as a way to smoothly transition between two different reference models without resorting to a scheduled switch. It also provides the ability to use an interpolated reference model when the plant is operating between design points. The time variation of the combination must satisfy some requirements to ensure stability but is otherwise user choice. Subject to these requirements, bounded tracking error behavior is demonstrated via Lyapunov stability analysis for the single-input, single-output, output feedback case. Tracking error convergence is asymptotic when time variation ceases. The proposed design is demonstrated in simulation of a numerical model.

Hashemi, Kelley E.↗

Assessment of Transition Modeling Capability in OVERFLOW with Emphasis on Swept-Wing Configurations

In preparation for comparisons with data obtained from the recently concluded experiments in the National Transonic Facility at the NASA Langley Research Center on the new common research model with natural laminar flow (CRM-NLF), an assessment of the current transition modeling capability in the NASA OVERFLOW 2.2o code has been carried out. A combination of the available experimental data and linear stability analysis is used to evaluate the accuracy and robustness of these models for selected swept-wing type configurations, with significant crossflow. An additional goal for this work involves providing a comparative assessment of the relevant transition models in the context of a single flow solver and identifying model limitations as well as the potential for future improvements that would help strengthen the physical basis of such transition models. Included in this investigation is an assessment of the sensitivities of the underlying transition models to grid resolution (wall-normal, as well as streamwise and spanwise) and the values of extra input parameters such as the level of surface roughness and freestream turbulence variables. The flow configurations targeted in this assessment include the NASA NLF(2)-0415 swept-wing configuration, the sickle-shaped wing introduced by the Technical University of Braunschweig, and the wing-body configuration of the CRM model from the fourth and fifth AIAA CFD Drag Prediction Workshops.

Freestream velocity↗

Assessment and Improvement of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon the correlations that were derived from a limited number of low-speed experiments and do not account for all of the transition mechanisms and/or their variation with the significant flow parameters. Available data from a recent experiment in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.2o code for a swept wing configuration at transonic cruise conditions. Specifically, the OVERFLOW solutions are used together with detailed stability analysis of the boundary layer flow over the new Common Research Model with Natural Laminar Flow (CRM-NLF) to evaluate the accuracy and the robustness of the transport-equation-based transition models, with the goal of proposing improvements that would help to strengthen the physical basis of these models for the important class of flows involving the combined effects of crossflow and flow compressibility. Results highlight the significant underprediction of the laminar flow extent within the inboard region of the wing, wherein the onset of transition may be attributed to a gradual amplification of Tollmien-Schlichting instabilities.

Boundary layer transition↗

Stability and targeting in Dawn's final orbit

The Dawn spacecraft conducted two extended missions at Ceres following the completion of the primary mission in June 2016. The final orbit of the second extended mission was designed to have a 35-km periapsis altitude for 10x higher resolution science. The mission ended in this orbit when the spacecraft ran out of attitude control propellant. In this paper, we describe the final orbit and discuss the challenges of flying this low at Ceres. We also include our stability analysis showing the spacecraft will remain in orbit for more than 20 years, as stipulated by the planetary protection requirements.

Kennedy, Brian↗

Design Exploration of a Transonic Cruise Slotted Airfoil

A knowledge-based aerodynamic design method, CDISC, has been leveraged for the computational design exploration of the cruise slotted airfoil as a drag-saving technology for single aisle, transonic transport aircraft. Aerodynamic predictions were generated using the NASA USM3D Reynolds-averaged Navier-Stokes flow solver, and laminar flow assessments were conducted using the NASA BLSTA3D boundary layer profile solver paired with the LASTRAC stability analysis and transition prediction software. The CDISC design method was used to parametrically vary several slotted-airfoil design variables to understand their impact on aerodynamic performance. In addition to establishing best practices for slotted airfoil design, it was observed that the historically noted skin-friction penalty for cruise slotted airfoils was attributed to the low-Reynolds-number boundary layer of the flap. At cruise, this skin-friction penalty negated any potential decrease in pressure drag enabled by the cruise slotted airfoil architecture. This observation led to the development of an Aft-Laminar Multi-element Airfoil concept that uses airfoil shaping to promote natural laminar flow on the flap. At cruise conditions, the concept is predicted to offer a 2.7% reduction in sectional drag relative to a fully turbulent supercritical airfoil. Off-design analyses showed that laminar flow could be maintained with near-cruise variations in angle of attack and Mach number, resulting in sustained aerodynamic efficiency improvements and a delay in the drag rise Mach number. A low-speed, high-lift analysis at takeoff conditions predicted a maximum lift coefficient of 2.4 could be achieved by use of a variable-camber leading edge on the main element and simple deflection of the flap. These computationally-predicted benefits have motivated future research toward an aft-laminar cruise slotted wing design as an incremental step toward the application of natural laminar flow technology on transonic commercial transports.

Brett R Hiller↗

On the origin & thermal stability of Arrokoth's and Pluto's ice

In this paper we discuss in a thermodynamic, geologically empirical way the long-term nature of the stable majority ices that could be present in Kuiper Belt object (KBO) 2014 MU69 (also called Arrokoth; hereafter “MU(69)”) after its 4.6 Gyr residence in the Edgeworth-Kuiper belt (EKB) as a cold classical object. We compare the upper bounds for the gas production rate (~10^(24) molecules/s) measured by the New Horizons (NH) spacecraft flyby on 01 Jan 2019 to estimates for the outgassing flux rates from a suite of common cometary and KBO ices at the average ~ 40 K sunlit surface temperature of MU69, but do not find the upper limit very constraining except for the most volatile of species (e.g. CO, N2, CH4). More constraining is the stability versus sublimation into vacuum requirement over Myr to Gyr, and from this we find only 3 common ices that are truly refractory: HCN, CH3OH, and H2O (in order of increasing stability), while NH3 and H2CO ices are marginally stable and may be removed by any positive temperature excursions in the EKB, as produced every 10^(8)–10^(9) years by nearby supernovae and passing O/B stars. To date the NH team has reported the presence of abundant CH3OH and H2O on MU69's surface (Stern et al., 2019; Grundy et al., 2020). NH3 has been searched for, but not found. We predict that future absorption feature detections, if any are ever derived from higher signal-to-noise ratio spectra, will be due to an HCN or poly-H2CO based species. Consideration of the conditions present in the EKB region during the formation era of MU69 lead us to state that it is highly likely that it “formed in the dark”, in an optically thick mid-plane, unable to see the nascent, variable, highly luminous Young Stellar Object (YSO)/TTauri Sun, and that KBOs contain HCN and CH3OH ice phases in addition to the H2O ice phases found in their short period (SP) comet descendants. Finally, when we apply our ice thermal stability analysis to bodies/populations related to MU69, we find that methanol ice is likely ubiquitous in the outer solar system; that if Pluto isn't a fully differentiated body, then it must have gained its hypervolatile ices from proto-planetary disk (PPD) sources in the first few Myr of the solar system's existence; and that hypervolatile rich, highly primordial comet C/2016 R2 was placed onto an Oort Cloud orbit on a similar few Myr timescale.

Pluto↗

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary layer transition↗

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗

Laminar-Turbulent Transition in High Mach-Number Flows

Accurate prediction of laminar-turbulent transition is critical for the design and characterization of high-speed flight systems. Presently, laminar-turbulent transition prediction in computational fluid dynamics (CFD) is accomplished using PDE-based models coupled to the governing equations of the flow system. Current industrial models suffer significant performance issues when analyzing high-speed systems because they were developed for lower speed flows. Current research focuses on multiple solutions to this problem. The coupling of physics-based stability analysis tools which produce quality transition predictions for high-speed flows to NASA’s OVERFLOW CFD solver is being pursued, as well as the implementation and improvement of emerging models developed specifically for high-speed flow.

High-Speed Flow↗