SEARCH · Engineering Papers

Results for “Numerical stability”

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 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS1994

(abstract) Short Time Period Variations in Jupiter's Synchrotron Radiation

The long term time variability of Jupiter's synchrotron radiation on yearly time scales has been established for some time. For many years, theorists have speculated about the effects variations in the solar wind, solar flux, Io, the Io torus, and Jupiter's magnetic field have on the ultra-relativistic electron population responsible for the emission. Early observational results suggested the additional possibility of a short term time variability, on timescales of days to weeks. In 1989 a program designed to investigate the existence of short term time variability using the 85 foot Hat Creek radio telescope operating at 1400 MHz was initiated. The availability of a dedicated telescope provided the opportunity, for the first time, to obtain numerous observations over the full Jupiter rotation period. These and future observations will enable two important studies, characterization and confirmation of possible short term variations, and the investigation of the stability of Jupiter's synchrotron emission beaming curve. Analysis of Hat Creek observations and early results from the Maryland Point Naval research Laboratory will be presented.

telescope magnetic field↗

NASA NTRS2006

AIRS: Improving Weather Forecasting and Providing New Data on Greenhouse Gases

This paper discusses the performance of AIRS and examines how it is meeting its operational and research objectives based on the experience of more than 2 yr with AIRS data. We describe the science background and the performance of AIRS in terms of the accuracy and stability of its observed spectral radiances. We examine the validation of the retrieved temperature and water vapor profiles against collocated operational radiosondes, and then we assess the impact thereof on numerical weather forecasting of the assimilation of the AIRS spectra and the retrieved temperature. We close the paper with a discussion on the retrieval of several minor tropospheric constituents from AIRS spectra.

remote sensing↗

NASA NTRS1991

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗

NASA NTRS2018

Exploration of Double Clad Fibers for Increased Stability of Bidirectional Free Space Optical Links

Bidirectional, high data rate, low size, weight, and power (SWaP), and low cost free space optical links are needed for space communication applications to send and receive large volumes of data. We are exploring design strategies for optical transceivers to reduce SWaP and cost through increased misalignment tolerance (pointing requirement reduction) and sharing the optical transmit and receive paths (imposing optical symmetry). In applications where the detector is fiber coupled, the fiber numerical aperture is the main driver of the pointing accuracy requirement. Unfortunately, reducing the pointing requirement by increasing the fiber numerical aperture symmetrically causes instability in received power over small environmental changes. This paper explores double clad fibers as a solution to both reduce the power instability and increase the pointing accuracy tolerance. Double clad fibers can transmit a Gaussian beam from a single mode fiber and receive in a multi-mode aperture. Results show that double clad fiber have an improved misalignment tolerance and a higher stability for small changes in temperature when compared to single mode fibers and multimode fibers. Also, double clad fibers are shown to match the performance of an asymmetrical link design with a single mode transmit fiber and a multimode receive fiber.

Double Clad fibers↗

NASA NTRS2018

Exploration of Double Clad Fibers for Increased Stability of Bidirectional Free Space Optical Links

Double Clad fibers↗

NASA NTRS1996

The Effects of Acoustic Treatment on Pressure Disturbances From a Supersonic Jet in a Circular Duct

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined circular duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions at each axial location. These values are used to determine the growth rates and phase velocities of the instability waves and the near field pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, and wall impedance, and the effects of these changes on non-axisymmetric modes. The primary results indicate that the effects of the duct wall on stability characteristics diminish as the outer flow increases and as the jet azimuthal mode number increases. Also, wall reflections are reduced when using a finite impedance boundary condition at the wall; but in addition, reflections are reduced and growth rates diminished by keeping the imaginary part of the impedance negative when using the negative exponential for the harmonic dependence.

Dahl, Milo D.↗

NASA NTRS1993

On the spectral stability of time integration algorithms for a class of constrained dynamics problems

Incomplete field formulations have recently been the subject of intense research because of their potential in coupled analysis of independently modeled substructures, adaptive refinement, domain decomposition, and parallel processing. This paper discusses the design and analysis of time-integration algorithms for these formulations and emphasizes the treatment of their inter-subdomain constraint equations. These constraints are shown to introduce a destabilizing effect in the dynamic system that can be analyzed by investigating the behavior of the time-integration algorithm at infinite and zero frequencies. Three different approaches for constructing penalty-free unconditionally stable second-order accurate solution procedures for this class of hybrid formulations are presented, discussed and illustrated with numerical examples. The theoretical results presented in this paper also apply to a large family of nonlinear multibody dynamics formulations. Some of the algorithms outlined herein are important alternatives to the popular technique consisting of transforming differential/algebraic equations into ordinary differential equations via the introduction of a stabilization term that depends on arbitrary constants and that influences the computed so1ution.

Farhat, Charbel↗

NASA NTRS1958

Estimation of Directional Stability Derivatives at Small Angles and Subsonic and Supersonic Speeds

Methods are presented for estimating the directional stability derivative increments contributed by the stabilizing surfaces of subsonic and supersonic aircraft. These methods are strictly applicable at zero angle of attack and small angles of sideslip. The procedure of totaling the incremental coefficients to obtain an estimation of the total empennage side-force and yawing-moment coefficient derivatives is also shown, together with numerical examples. A correlation is presented between estimated and experimental incremental coefficients which indicates that the methods of this report generally estimate the increment of side force gained by the addition of a panel to within +/-10 percent of the experimental value while the yawing-moment increment is generally estimated to within +/-20 percent. This is true for both subsonic and supersonic Mach numbers. An example application of the methods to one of the problems in directional stability, that of minimizing the effect of Mach number on the side-force coefficient derivative of the empennage, is discussed.

Goodwin, Frederick K.↗

NASA NTRS2020

Convolutional Neural Network for Transition Modeling Based on Linear Stability Theory

Transition prediction is an important aspect of aerodynamic design because of its impact on skin friction and potential coupling with flow separation characteristics. Traditionally, the modeling of transition has relied on correlation-based empirical formulas based on integral quantities such as the shape factor of the boundary layer. However, in many applications of computational fluid dynamics, the shape factor is not straightforwardly available or not well-defined. We propose using the complete velocity profile along with other quantities (e.g., frequency, Reynolds number) to predict the perturbation amplification factor. While this can be achieved with regression models based on a classical fully connected neural network, such a model can be computationally more demanding. We propose a novel convolutional neural network inspired by the underlying physics as described by the stability equations. Specifically, convolutional layers are first used to extract integral quantities from the velocity profiles, and then fully connected layers are used to map the extracted integral quantities, along with frequency and Reynolds number, to the output (amplification ratio). Numerical tests on classical boundary layers clearly demonstrate the merits of the proposed method. More importantly, we demonstrate that, for Tollmien-Schlichting instabilities in two-dimensional, low-speed boundary layers, the proposed network encodes information in the boundary layer profiles into an integral quantity that is strongly correlated to a well-known, physically defined parameter – the shape factor.

Laminar-turbulent transition↗

NASA NTRS2019

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

NASA NTRS1973

Computer simulation for time optimal or energy optimal attitude control of spin-stabilized spacecraft.

An original technique for determining the optimal magnetic torque strategy for control of the attitude of spin stabilized spacecraft is presented. By employing Lagrange multipliers and the Calculus of Variations, optimal control equations are derived which define minimum time and minimum energy attitude maneuvers. Computer program algorithms to numerically solve these optimal control equations are also described. The performance of this technique is compared with a commonly employed planning method.

Woolley, R. D.↗

NASA NTRS1978

ORACLS: A system for linear-quadratic-Gaussian control law design

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

Armstrong, E. S.↗

NASA NTRS1982

Stability of boundary layers with porous suction strips: Experiment and theory

Low turbulence tunnel experiments on the stability and transition of 2 D boundary layers on flat plates with and without suction are described. A number of general suction cases are discussed. Test results showed that the maximum stabilization occurred when the suction was moved toward the Branch I neutral point. An analytical study of the stability of two dimensional, incompressible boundary layer flows over plates with suction through porous strips was performed. The mean flow was calculated using linearized triple deck, closed form solutions. The stability results of the triple deck theory are shown to be in good agreement with those of the interacting boundary layers. An analytical optimization scheme for the suction configuration was developd. Numerical calculations were performed corresponding to the experimental configurations. In each case, the theory correctly predicts the experimental results.

Reynolds, G. A.↗

NASA NTRS2019

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

NASA NTRSDate not supplied

Interpretation of Vehicle Tumbling Predictions from 6-DOF Entry and Descent Simulation

Blunt body entry vehicles are subject to dynamic instability during terminal descent. This often manifests as limit cycle oscillations in total angle of attack, but can diverge into tumbling behavior under certain conditions. For the Mars Sample Return Earth Entry Vehicle (MSR EEV), there is a constraint on the orientation of the sample tubes so backward impact is impermissible. Past missions have chosen to deploy parachutes to preclude tumbling, but active events after release of MSR EEV have been ground-ruled out with the intent to maximize system reliability. Prevention of tumbling during subsonic descent is a design driver for MSR EEV. During preliminary design of the MSR EEV, six degree-of-freedom numerical simulations indicated an unacceptably high probability of tumbling for a 60degree sphere-cone forebody geometry, which necessitated a design change. Decreasing the forebody angle was expected to improve dynamic stability, but would also adversely impact mass, aeroheating, and manufacturing risk. Hence there was strong motivation to understand the physical drivers for the onset of tumbling, and to determine: (a) whether the causes of tumbling are representative of physically realizable vehicle entry configurations and (b) what changes can be made to existing design and analysis practices to ensure a stable vehicle.

Entry Vehicle↗

NASA NTRS2023

Implementation, Realization and an Effective Solver of Two-Equation Turbulence Models

Currently, when the Reynolds-Averaged Navier-Stokes (RANS) equations are solved using turbulence modeling, most often the one-equation model of Spalart and Allmaras is used. Then, it is only necessary to solve the RANS equations in conjunction with a single transport equation for modeling turbulence. For this model, considerable assessment and analysis has been performed, allowing the possibility of a reliable solution method for an eddy viscosity required to compute the Reynolds stresses in the RANS equations. Such evaluation along with analysis has not been achieved to realize similar performance with two-equation models of the k-w type. The primary objective of this paper is to present and discuss the components of an effective numerical algorithm for solving the RANS equations and the two transport equations of k-w type turbulence models. All the important details of the turbulence model as actually implemented are given, which is sometimes not done in various papers considering such modeling. The viability and effectiveness of this solution algorithm are demonstrated by solving both two-dimensional and three-dimensional aerodynamic flows. In all applications, a linear rate of convergence without oscillations or other evidence of unstable behavior is observed. This behavior is also particularly true when the proposed algorithm is applied to systematically renewed mesh sequences, which is generally not observed with algorithms solving more than one transport equation. Thus, numerical integration errors are systematically reduced, allowing for a significantly more reliable assessment of the effectiveness of the turbulence model. Additionally, in this paper, analysis of the solution algorithm, including linear stability, is also performed for a particular flow problem.

Turbulence↗

NASA NTRS2024

Space Algae-2 Ground and Lunar Analog Studies in Preparation for Long-Duration Propagation of Cyanobacteria in Spaceflight

There are numerous applications for microalgae in spaceflight missions and on Earth, such as oxygen production, carbon dioxide removal, nutrition, wastewater processing, and biofuel production. Space Algae-2 aims to test the genetic stability of Arthrospira platensis, commonly known as spirulina, during six-months of continuous culture on the International Space Station. Long-duration exposure to ionizing radiation and microgravity may impact growth, nutrient composition, and genetic stability. The high protein, vitamin, antioxidant content, and radiation resistance make spirulina a promising candidate for bioregenerative life support systems. A concept of operations was developed to grow and harvest algal biomass in space. Preflight testing experiments optimized conditions for an extended growth period in a gas permeable bioreactor bag. Preflight and post-harvest storage methods were developed in addition to a novel cryopreservation method. After sample return, multi-omics analysis will be conducted to determine the mutation rate, gene expression, and the protein and metabolite profile. The concept of operations for Space Algae-2 was tested during a lunar mission simulation within a semi-controlled environment. During a six-day lunar analog mission at the Hawai’i Space Exploration Analog and Simulation (HI-SEAS) A. platensis was successfully grown using flight-like hardware. The cyanobacteria were harvested and used to supplement bread as an example of spirulina biomass utilization. Overall, the data collected from Space Algae-2 will inform potential bioengineering of spirulina for space and terrestrial applications.

Algae↗

NASA NTRSDate not supplied

Space Algae-2 Ground and Lunar Analog Studies in Preparation for Long-Duration Propagation of Cyanobacteria in Spaceflight

There are numerous applications for microalgae in spaceflight missions and on Earth, such as, oxygen production, carbon dioxide removal, nutrition, wastewater processing, and biofuel production. Space Algae-2 aims to test the genetic stability of Arthrospira platensis, commonly known as spirulina, during six-months of continuous culture in spaceflight on the International Space Station. Long-duration exposure to ionizing radiation and microgravity may impact growth, nutrient composition, and genetic stability. The high protein, vitamin, antioxidant content, and radiation resistance make spirulina a promising candidate for bioregenerative life support systems during long-duration missions. A concept of operations was developed to grow and harvest algal biomass in space. Preflight testing experiments were conducted to optimize conditions for an extended growth period in a gas permeable bioreactor bag. Preflight and post-harvest storage methods were developed in addition to a novel cryopreservation method. After sample return, multi-omics analysis will be conducted to determine the mutation rate, gene expression, and protein and metabolite profile. The concept of operations for Space Algae-2 was tested at HI-SEAS (Hawai’i Space Exploration Analog and Simulation) during a six-day lunar analog mission (EMMIHS23, EuroMoonMars, International MoonBase Alliance, HI-SEAS, 2023). A. platensis was grown in the semi-controlled environment using flight-like hardware and solar powered LED lights. Then, the biomass was harvested and used to supplement bread as an example of A. platensis utilization. Overall, the data collected from Space Algae-2 will inform potential bioengineering of spirulina for space and terrestrial applications.

Algae↗