Engineering Papers⌕ Search

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.

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

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Liquefied Bleed for Stability and Efficiency of High Speed Inlets

A mission analysis code was developed to perform a trade study on the effectiveness of liquefying bleed for the inlet of the first stage of a TSTO vehicle. By liquefying bleed, the vehicle weight (TOGW) could be reduced by 7 to 23%. Numerous simplifying assumptions were made and lessons were learned. Increased accuracy in future analyses can be achieved by: Including a higher fidelity model to capture the effect of rescaling (variable vehicle TOGW). Refining specific thrust and impulse models ( T m a and Isp) to preserve fuel-to-air ratio. Implementing LH2 for T m a and Isp. Correlating baseline design to other mission analyses and correcting vehicle design elements. Implementing angle-of-attack effects on inlet characteristics. Refining aerodynamic performance (to improve L/D ratio at higher Mach numbers). Examining the benefit with partial cooling or densification of the bleed air stream. Incorporating higher fidelity weight estimates for the liquefied bleed system (heat exchange and liquid storage versus bleed duct weights) could be added when more fully developed. Adding trim drag or 6-degree-of-freedom trajectory analysis for higher fidelity. Investigating vehicle optimization for each of the bleed configurations.

Saunders, J. David↗
NASA NTRS2019

Stability and Spin-Orbit Resonance Analysis of Low Altitude Martian Orbits

Orbit stability has been thoughtfully studied in various celestial bodies. In particular the focus has been placed on the Moon, due to the unstable natural perturbations of its gravity field. The increasing interest in Mars orbiters brings the question of the likelihood of natural decay in low altitude regimes. This paper studies the change in shape of low altitude Mars orbits by carrying out large sets of numerical high fidelity simulations. Results showed that various configurations of the orbital elements gave perturbations that result-ed in unstable orbits. The paper also studies the potential causes of the observed unstable regions. First by taking a close look at zonal and tesseral harmonics to find the implications of Mars mass concentrations of the used gravity fields, and second by computing theoretical spin-orbit resonances to study their implications in the stability at low altitudes.

gravity perturbations↗
NASA NTRS1972

Analytic prediction of airplane equilibrium spin characteristics

The nonlinear equations of motion are solved algebraically for conditions for which an airplane is in an equilibrium spin. Constrained minimization techniques are employed in obtaining the solution. Linear characteristics of the airplane about the equilibrium points are also presented and their significance in identifying the stability characteristics of the equilibrium points is discussed. Computer time requirements are small making the method appear potentially applicable in airplane design. Results are obtained for several configurations and are compared with other analytic-numerical methods employed in spin prediction. Correlation with experimental results is discussed for one configuration for which a rather extensive data base was available. A need is indicated for higher Reynolds number data taken under conditions which more accurately simulate a spin.

Adams, W. M., Jr.↗
NASA NTRS1978

Ride quality flight testing

The ride quality experienced by passengers is a function of airframe rigid-body, elastic dynamic responses, autopilot, and stability augmentation system control inputs. A frequency response method has been developed to select sinusoidal elevator input time histories yielding vertical load factor distributions, within a given limit, as a function of fuselage station. The numerical technique is illustrated by applying two-degree-of-freedom short-period and first symmetric mode equations of motion to a B-1 aircraft at Mach 0.85 during sea level flight conditions.

Swaim, R. L.↗
NASA NTRS1979

A fast, conservative algorithm for solving the transonic full-potential equation

A fast, fully implicit approximate factorization (AF) algorithm designed to solve the conservative transonic full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., 'rotated differencing'), and thereby greatly enhances the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line overrelaxation algorithms.

Holst, T. L.↗
NASA NTRS1979

Perturbation solutions of combustion instability problems

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Googerdy, A.↗
NASA NTRS1983

A diagnostic study of the forcing of the Ferrel cell by eddies, with latent heat effects included

A diagnostic study of the forcing of the Ferrel cell by eddy fluxes in the Northern Hemisphere is carried out. The quasi-geostrophic omega equation, and Oort and Rasmusson's (1971) data set, are used. The effects of condensation associated with the large scale motions are introduced to the omega equation by using the quasi-geostrophic moisture conservation equation. Thus, the dry static stability is replaced by a moist static stability, and the forcing of the Ferrel cell by eddy latent heat fluxes as well as sensible heat and momentum fluxes is included. Both effects tend to enhance the forcing of the Ferrel cell. The numerical analysis indicates that the effects are small in January, but in July the maximum vertical velocities are enhanced by about 30 percent.

Salustri, G.↗
NASA NTRS1985

Calibration stability of some hot-cathode ion gauges

Data are presented on the stability of calibration of some hot-cathode ion gauges: a conventional triode used in the range of 0.0002-0.2 Pa, and Bayard-Alpert-type gauges used in the range of 0.0002-0.02 Pa. Numerical magnitudes are given of the systematic corrections that must be applied to the gauge indication in order to achieve the highest accuracy of pressure indication. The relative sensitivities of these gauges to Ar, N2, Ne, and He are also given.

Warshawsky, I.↗
NASA NTRS1985

Measurements of distributed polarized radio sources from spinning spacecraft - Effect of a tilted axial antenna ISEE-3 application and results

An analysis is presented of the system response of a satellite receiver-antenna system to locate a radio source when the satellite is tilted on its axis. The satellite is spin stabilized but experiences a tilt due to either a mechanical misalignment or a shift in the electrical axis caused by parasitic currents in other spacecraft structures. The shorter the antenna, the more significant the effects. Numerical techniques are developed for obtaining the Stokes parameters and the angular parameters of a uniform conical source sensed by a linear antenna in order to derive the average power response of a synthesized dipole to a point on a distributed polarized source. Relative gains are calculated along the antenna at different angles to the source. The techniques are applied to sample ISEE-3 satellite data for Type III solar radio bursts which were sensed by an axial and an equatorial antenna. The two antennas permit localization of the source and quantification of the polarization and angular extent of the source. The resulting high precision in calculations of all three source parameters commends use of the model in analyses of data from the planned ULYSSES mission.

Fainberg, J.↗
NASA NTRS1988

Ion-ion waves in the auroral region - Wave excitation and ion heating

The properties of the ion-ion mode which is excited in plasmas when two or more cold ion beams are streaming relative to one another are investigated assuming a warm electron distribution at rest (the model consistent with particle distributions observed in auroral plasma cavities). Numerical solutions are derived for the generalized electrostatic dispersion equation for parallel propagation and for oblique propagation of ion waves. It is shown that the relative ion temperature and concentration have significant effects on the stability of the accelerated plasma and the expected ion heating. Finally, a relationship between ion drift and thermal speed is derived using the marginal stability of ion-ion waves for nonzero ion temperature; the relationship was found to compare favorably with DE-1 ion observations at the high-altitude boundary of the auroral cavity.

Dusenbery, P. B.↗
NASA NTRS2000

An Indirect Mixed-Sensitivity Approach to Microgravity Vibration Isolation: The Exploitation of Kinematic Coupling In Frequency-weighting Design-Filter Selections

Many space-science experiments need an active isolation system to provide them with the requisite microgravity environment. The isolation systems planned for use with the International Space Station have been appropriately modeled using relative position relative velocity, and acceleration states. In theory, frequency-weighting design filters can be applied to these state-space models, in order to develop optimal H2 or mixed-norm controllers with desired stability and performance characteristics. In practice, however, the kinematic coupling among the various states can lead, through the associated frequency-weighting-filters, to conflicting demands on the Riccati design "machinery." The results can be numerically ill-conditioned regulator and estimator Riccati equations and/or reduced intuition in the design process. In addition, kinematic coupling can result in a redundancy in the demands imposed by the frequency weights. Failure properly to account for this type of coupling can lead to an unnecessary increase in controller dimensionality and, in turn, controller complexity. This paper suggests a rational approach to the assignment of frequency-weighting design filters, in the presence of the kinematic coupling among states that exists in the microgravity vibration isolation problem.

Hampton, R. David↗
NASA NTRS2000

An Indirect Mixed-Sensitivity Approach to Microgravity Vibration Isolation: The Exploitation of Kinematic Coupling In Frequency-Weighting Design-Filter Selections

Many space science experiments need an active isolation system to provide them with the requisite microgravity environment. The isolation systems planned for use with the International Space Station have been appropriately modeled using relative position, relative velocity, and acceleration states. In theory, frequency design filters can be applied to these state-space models, in order to develop optimal H, or mixed-norm controllers with desired stability- and performance characteristics. In practice. however, the kinematic coupling among the various states can lead, through the associated frequency-weighting-filters, to conflicting demands on the Riccati design "machinery." The results can be numerically ill-conditioned regulator and estimator Riccati equations and/or reduced intuition in the design process. In addition, kinematic coupling can result in a redundancy in the demands imposed by the frequency weights. Failure properly to account for this type of coupling can lead to an unnecessary increase in controller dimensionality and, in turn, controller complexity. This paper suggests a rational approach to the assignment of frequency weighting design filters, in the presence of the kinematic coupling among states that exists in the microgravity vibration isolation problem.

Hampton, R. David↗
NASA NTRS2003

Airborne Measurements of Formaldehyde Employing a Tunable Diode Laser Absorption Spectrometer During TRACE-P

This final report summarizes the progress achieved over the entire 3-year proposal period including two extensions spanning 1 year. These activities include: 1) Preparation for and participation in the NASA 2001 TRACE-P campaign using our airborne tunable diode laser system to acquire measurements of formaldehyde (CH2O); 2) Comprehensive data analysis and data submittal to the NASA archive; 3) Follow up data interpretation working with NASA modelers to place our ambient CH2O measurements into a broader photochemical context; 4) Publication of numerous JGR papers using this data; 5) Extensive follow up laboratory tests on the selectivity and efficiency of our CH20 scrubbing system; and 6) An extensive follow up effort to assess and study the mechanical stability of our entire optical system, particularly the multipass absorption cell, with aircraft changes in cabin pressure.

Fried, Alan↗
NASA NTRS2009

Addition of Improved Shock-Capturing Schemes to OVERFLOW 2.1

Existing approximate Riemann solvers do not perform well when the grid is not aligned with strong shocks in the flow field. Three new approximate Riemann algorithms are investigated to improve solution accuracy and stability in the vicinity of strong shocks. The new algorithms are compared to the existing upwind algorithms in OVERFLOW 2.1. The new algorithms use a multidimensional pressure gradient based switch to transition to a more numerically dissipative algorithm in the vicinity of strong shocks. One new algorithm also attempts to artificially thicken captured shocks in order to alleviate the errors in the solution introduced by "stair-stepping" of the shock resulting from the approximate Riemann solver. This algorithm performed well for all the example cases and produced results that were almost insensitive to the alignment of the grid and the shock.

Burning, Pieter G.↗
NASA NTRS2010

Receptivity and Transition of Supersonic Boundary Layers Over Swept Wings

The receptivity, stability, and transition of three-dimensional supersonic boundary layers over (1) a swept cylinder, (2) a swept wing with a sharp leading edge, and (3) a swept wing with a blunt leading edge are numerically investigated for a free-stream Mach number of 3. These computations are compared to an earlier experimental and computational study performed by Archambaud et al.1 The steady flow fields with and without roughness elements are obtained by solving the full Navier-Stokes equations. The N-factors computed in this study at the transition onset locations reported in Ref. 1 for flow over the swept cylinder are approximately 16.5 for traveling crossflow disturbances and 9 for stationary disturbances. The N-factors for the traveling crossflow are high based on our past experiences. However, they are comparatively smaller than those reported by Archambaud et al., who found N-factor values in the range of 20 to 25 for traveling disturbances and 13 to 20 for stationary disturbances. Similarly, the N-factors computed in this study for the traveling and stationary disturbances for the flow over the sharp wing are approximately 7 and 2.5, respectively, and for the flow over the blunt wing are 6.5 and 4.8, respectively. Using the envelope method, Archambaud et al. obtained values of approximately 8.0 and 4.0 for the sharp wing case and 16.0 and 12.0 for the blunt wing case.

Balakumar, Ponnampalam↗
NASA NTRSDate not supplied

Logical shadow tomography: Efficient estimation of error-mitigated observables

We introduce a technique to estimate error-mitigated expectation values on noisy quantum computers. Our technique performs shadow tomography on a logical state to produce a memory-efficient classical reconstruction of the noisy density matrix. Using efficient classical post-processing, one can mitigate errors by projecting into the codespace as in subspace expansion and taking powers of the density matrix as in virtual distillation. Relative to subspace expansion which requires Ω (2^((n-1)k) samples to estimate a Pauli observable with an [[n; k]] stabilizer code, our technique requires only Ө(2^k) samples. Relative to virtual distillation, our technique can compute powers of the density matrix without implementing additional copies of quantum states the quantum computer. We present numerical results using logical states encoded with up to sixty physical qubits and show fast convergence to error-free expectation values with only 10^5 samples under 1% depolarizing noise.

quantum computing↗
NASA NTRS1999

Numerical Modeling of Crystal of ZnSe by Physical Vapor Transport - Towards a more Comprehensive Formulations

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of PVT growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided. The effect of gravity on the process from numerical computations shows subtle effects although experimental evidence from vertically and horizontally grown samples show dramatic evidence of gravitational effects. The shortcomings of the problem formulation will be discussed and a framework will be provided leading up towards a more comprehensive model of PVT systems.

Ramachandran, N.↗
NASA NTRS1973

Numerical analysis of stiffened shells of revolution. Volume 1: Theory manual for STARS-2S, 2B, 2V digital computer programs

The theoretical analysis background for the STARS-2 (shell theory automated for rotational structures) program is presented. The theory involved in the axisymmetric nonlinear and unsymmetric linear static analyses, and the stability and vibrations (including critical rotation speed) analyses involving axisymmetric prestress are discussed. The theory for nonlinear static, stability, and vibrations analyses, involving shells with unsymmetric loadings are included.

Svalbonas, V.↗