SEARCH · Engineering Papers
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.
The beam-foil time-of-flight method has been used to investigate radiative lifetimes and transition rates involving allowed intrashell transitions within the L shell of highly ionized sulfur. The results for these transitions, which can be particularly correlation-sensitive, are compared with current calculations based upon multiconfigurational models.
The efficiency characteristics of finite element and finite difference approximations for the steady-state solution of the Navier-Stokes equations are examined. The finite element method discussed is a standard Galerkin formulation of the incompressible, steady-state Navier-Stokes equations. The finite difference formulation uses simple centered differences that are O(delta x-squared). Operation counts indicate that a rapidly converging Newton-Raphson-Kantorovitch iteration scheme is generally preferable over a Picard method. A split NOS Picard iterative algorithm for the finite difference method was most efficient.
A method for determining the optimum inlet geometry of a liquid rocket engine swirl injector includes obtaining a throttleable level phase value, volume flow rate, chamber pressure, liquid propellant density, inlet injector pressure, desired target spray angle and desired target optimum delta pressure value between an inlet and a chamber for a plurality of engine stages. The tangential inlet area for each throttleable stage is calculated. The correlation between the tangential inlet areas and delta pressure values is used to calculate the spring displacement and variable inlet geometry. An injector designed using the method includes a plurality of geometrically calculated tangential inlets in an injection tube; an injection tube cap with a plurality of inlet slots slidably engages the injection tube. A pressure differential across the injector element causes the cap to slide along the injection tube and variably align the inlet slots with the tangential inlets.
An analytic theory is presented in which the classical slender wing theory is modified to account for the combined effects of large angle of attack and nonsonic Mach number on the unsteady aerodynamics. The computed results agree well with available static and dynamic experimental data for slender delta wings in the freestream Mach number range between 0 and 2.8. The method was extended to compute the unsteady aerodynamics of the space shuttle orbiter by defining an equivalent slender wing using static experimental data. The results obtained in this manner are in good agreement with dynamic experimental results for the freestream Mach number range between 0.3 and 1.2.
Approximation concepts and dual method algorithms are combined to create a new method for minimum weight design of structural systems. Approximation concepts convert the basic mathematical programming statement of the structural synthesis problem into a sequence of explicit primal problems of separable form. These problems are solved by constructing explicit dual functions, which are maximized subject to nonnegativity constraints. The dual method is successfully extended to deal with pure discrete and mixed continuous-discrete design variable problems. The power of the method presented is illustrated with numerical results for example problems, including a thin delta wing with fiber composite skins.
An analysis of the lift augmentation due to a thin jet of air issuing from a slot along the leading edge of a delta wing is presented. The problem is treated with an extension of the method of Brown and Michael, representing the separated flow on the lee side of the wing by a pair of concentrated vortices and corresponding feeding sheets. It is assumed that the jet is not affected by Coanda forces. The analysis produces qualitative agreement with experiments.
The effects of F content on the Tc of superconducting YBa2Cu3F(x)O(6.5+delta-0.5x) (x = 0.0165-1.65) prepared by the solid-state reaction method (Bandsahl and Sandkuhl, 1988) are investigated experimentally. The results are presented in tables and graphs and characterized. Tc is found to increase steadily from 90.8 to 93.4 K as x is increased from 0.0165 to 0.066, but to decrease as x is further increased (Tc = 90.3 K at x = 1.65). The structure of the samples with x = 0.066 or lower is found to be single-phase perovskite.
The three-dimensional vortical viscous flow past a sideslipping delta wing at a high angle of attack is numerically investigated. A computational method for calculating low-speed viscous flowfields is developed. The time-iterative method uses an implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. An algebraic turbulence model is implemented to account for the eddy viscosity. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. Computed results are compared with experimental data.
Spacecraft such as Magellan and Mars Observer will use earth-based interferometric delay-rate tracking to meet their navigation requirements. Examples of methods used to implement these measurements are Differenced One-Way Doppler, delta-Differenced One-Way Doppler, and Two-Way Minus Three-Way Doppler. These measurements have a greater capability to observe certain orbital elements than that possessed by two-way Doppler, which has been the sole data type used in all U.S. planetary orbiter missions to date. In this paper, an approximate analytic model is developed for delay-rate measurements, as a function of classical orbital elements. The resulting expressions are used to show how the information content of these data types varies with orbit size, shape, and orientation. Comparison is made with the information content of Doppler data, using orbital element sets derived from the Magellan and Mars Observer missions. Results indicate that navigation accuracy improvements obtained by augmenting Doppler data with delay-rate measurements are greatest for spacecraft in low-altitude circular orbits; but decrease steadily for orbits of progressively greater size and eccentricity. The differences in Doppler-only and Doppler plus delay-rate orbit determination accuracies were also found to be highly dependent on Doppler measurement accuracy.
This paper presents the status and on-going development of a new model for the prediction of aft-radiated turbofan broadband noise. The proposed model is developed using data from two full-scale flight tests, the Boeing Quiet Technology Demonstrator 2 (QTD2) and the NASA/Boeing Propulsion Airframe Aeroacoustics and Aircraft System Noise (PAA&ASN) flight tests. The proposed model accounts for an observed coupling between directivity and fan relative tip Mach number. Additional features such as improved spectral and directivity characteristics are discussed. Significant improvements over prior methods implemented in the NASA Aircraft Noise Prediction Program are demonstrated, and predicted deltas to measured levels are reduced from +10/-20 dB to +/-3 dB or better at the most relevant frequencies. This is expected to greatly improve the ability of NASA to predict the noise of both current and future aircraft concepts.
A method is presented which finds Mars parking orbits which allow tangential periapsis burns at both arrival and departure. This method accounts for the actual geometry at both arrival and departure between the hyperbolic asymptotes and the orbital plane, along with the precession effects caused by the oblateness of Mars. Thus, realistic Delta-V values (and hence initial low-earth orbit masses) are obtained for these orbits. The results obtained from the present method compare very well with a trajectory integration program while only requiring CPU times of about one minute. Therefore, due to the computational efficiency and accuracy, the present method would be an ideal tool to use in preliminary mission design, since it provides the opportunity to incorporate realistic Mars parking orbits effects.
The present report concerns a method of computing the velocity and pressure distributions on bodies of revolution in axially symmetrical flow in the subsonic range. The differential equation for the velocity potential Phi of a compressible fluid motion is linearized tn the conventional manner, and then put in the form Delta(Phi) = 0 by affine transformation. The quantity Phi represents the velocity potential of a fictitious incompressible flow, for which a constant superposition of sources by sections is secured by a method patterned after von Karman which must comply with the boundary condition delta(phi)/delta(n) = 0 at the originally specified contour. This requirement yields for the "pseudo-stream function" psi a differential equation which must be fulfilled for as many points on the contour as source lengths are assumed. In this manner, the problem of defining the still unknown source intensities is reduced to the solution of an inhomogeneous equation system. The pressure distribution is then determined with the aid of Bernoulli's equation and adiabatic equation of state. Lastly, the pressure distributions in compressible and incompressible medium are compared on a model problem.
A quasitransient in-ground effect method is developed to study the effect of vertical landing on a launch vehicle base-heating environment. This computational methodology is based on a three-dimensional, pressure-based, viscous flow, chemically reacting, computational fluid dynamics formulation. Important in-ground base-flow physics such as the fountain-jet formation, plume growth, air entrainment, and plume afterburning are captured with the present methodology. Convective and radiative base-heat fluxes are computed for comparison with those of a flight test. The influence of the laminar Prandtl number on the convective heat flux is included in this study. A radiative direction-dependency test is conducted using both the discrete ordinate and finite volume methods. Treatment of the plume afterburning is found to be very important for accurate prediction of the base-heat fluxes. Convective and radiative base-heat fluxes predicted by the model using a finite rate chemistry option compared reasonably well with flight-test data.
An on-line compliance-based method to account for residual stress effects in stress-intensity factor and fatigue crack growth property determinations has been evaluated. Residual stress intensity factor results determined from specimens containing friction stir weld induced residual stresses are presented, and the on-line method results were found to be in excellent agreement with residual stress-intensity factor data obtained using the cut compliance method. Variable stress-intensity factor tests were designed to demonstrate that a simple superposition model, summing the applied stress-intensity factor with the residual stress-intensity factor, can be used to determine the total crack-tip stress-intensity factor. Finite element, VCCT (virtual crack closure technique), and J-integral analysis methods have been used to characterize weld-induced residual stress using thermal expansion/contraction in the form of an equivalent delta T (change in local temperature during welding) to simulate the welding process. This equivalent delta T was established and applied to analyze different specimen configurations to predict residual stress distributions and associated residual stress-intensity factor values. The predictions were found to agree well with experimental results obtained using the crack- and cut-compliance methods.
A structured singular value (mu) analysis method of computing flutter margins has robust stability of a linear aeroelastic model with uncertainty operators (Delta). Flight data is used to update the uncertainty operators to accurately account for errors in the computed model and the observed range of aircraft dynamics of the aircraft under test caused by time-varying aircraft parameters, nonlinearities, and flight anomalies, such as test nonrepeatability. This mu-based approach computes predict flutter margins that are worst case with respect to the modeling uncertainty for use in determining when the aircraft is approaching a flutter condition and defining an expanded safe flight envelope for the aircraft that is accepted with more confidence than traditional methods that do not update the analysis algorithm with flight data by introducing mu as a flutter margin parameter that presents several advantages over tracking damping trends as a measure of a tendency to instability from available flight data.
The electronic structure and the hole concentrations in the high Tc superconductor HgBa2CuO(4+delta) (delta = O, 1) has been investigated by employing a first principles full potential self-consistent LMTO method with the local density functional theory. The scalar relativistic effects have been considered. The hole concentrations of the Cu-d and O-p(x,y) orbitals are seen to be larger for the HgBaCuO5 system than those of the HgBaCuO4 solid. However, the van Hove singularity (vHs) induced Cu-d and O-p peak which is seen to lie comparatively away and above the Fermi level in the delta = 1 system shifts towards the Fermi level in the delta = 0 system. Thus, the superconducting behavior appears to originate from the occurrence of the vHs peak at the Fermi level. The Fermi surface nesting area in the delta = 0 compound is seen to be larger than in the delta = 1 compound. The calculation reveals that the increase in pressure on the crystal enhances the hole concentrations but without showing any optimum value, On the other hand, the vHs peak approaches to-wards the Fermi level with pressure and crosses the Fermi surface near V/Vo approximately equals 0.625 (V and Vo are the crystal volumes at high and normal pressures, respectively). Our calculated value of the bulk modulus equal to 0.626 Mbar predicts the occurrence of this crossover at about 24 GPa which is in complete agreement with the experimental value. At this pressure the compound has maximum nesting area and self-doped behavior.
The current primary method of deepspace navigation is the NASA Deep Space Network (DSN). High-performance navigation is achieved using Delta Differential One-Way Range techniques that utilize simultaneous observations from multiple DSN sites, and incorporate observations of quasars near the line-of-sight to a spacecraft in order to improve the range and angle measurement accuracies. Over the past four decades, x-ray astronomers have identified a number of xray pulsars with pulsed emissions having stabilities comparable to atomic clocks. The x-ray pulsar-based navigation and time determination (XNAV) system uses phase measurements from these sources to establish autonomously the position of the detector, and thus the spacecraft, relative to a known reference frame, much as the Global Positioning System (GPS) uses phase measurements from radio signals from several satellites to establish the position of the user relative to an Earth-centered fixed frame of reference. While a GPS receiver uses an antenna to detect the radio signals, XNAV uses a detector array to capture the individual xray photons from the x-ray pulsars. The navigation solution relies on detailed xray source models, signal processing, navigation and timing algorithms, and analytical tools that form the basis of an autonomous XNAV system. Through previous XNAV development efforts, some techniques have been established to utilize a pulsar pulse time-of-arrival (TOA) measurement to correct a position estimate. One well-studied approach, based upon Kalman filter methods, optimally adjusts a dynamic orbit propagation solution based upon the offset in measured and predicted pulse TOA. In this delta position estimator scheme, previously estimated values of spacecraft position and velocity are utilized from an onboard orbit propagator. Using these estimated values, the detected arrival times at the spacecraft of pulses from a pulsar are compared to the predicted arrival times defined by the pulsar s pulse timing model. A discrepancy provides an estimate of the spacecraft position offset, since an error in position will relate to the measured time offset of a pulse along the line of sight to the pulsar. XNAV researchers have been developing additional enhanced approaches to process the photon TOAs to arrive at an estimate of spacecraft position, including those using maximum-likelihood estimation, digital phase locked loops, and "single photon processing" schemes that utilize all available time data associated with each photon. Using pulsars from separate, non-coplanar locations provides range and range-rate measurements in each pulsar s direction. Combining these different pulsar measurements solves for offsets in position and velocity in three dimensions, and provides accurate overall navigation for deep space vehicles.
The charge, energy spectra, and isotopic composition of hydrogen through nitrogen have been studied with the GSFC/UNH cosmic ray experiments on board the Pioneer 10 and 11 spacecraft. Good charge and mass resolution have been obtained by utilizing three parameter analysis (i.e., two measurements of delta E and a measurement of the residual kinetic energy of stopping particles). This method allows isotopic resolution for the High Energy Telescope (HET) which operates over a typical energy range of 30-56-MeV/nuc for He. The observations will be over solar quiet times from September 1976-September 1977. The resolution, methods of analysis, spectra, and isotopic composition will be presented and discussed where the statistical significance of H to C isotopes is meaningful.