Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

CO2 and CO electron vibrational cooling rates

Electron cooling rates for electron vibrational excitation of CO2 and CO and rotational excitation of CO are calculated. Results of the calculation are expressed as closed form analytic relations that are both convenient for numerical calculations and valid over a wide range of electron and neutral temperatures. It is found that the cooling rates and their temperature dependences differ significantly from other calculations. For CO2, these differences are related mainly to the assumptions about the approach to equilibrium and not to the actual cross-section choices. For electron temperatures below 1000 K and gas temperatures near 300 K (representative of conditions in the Venus atmosphere, where collisional cooling of the electrons dominates), the CO2 cooling rate is as much as an order of magnitude larger than previously considered in ionospheric models.

Porter, H. S.↗
NASA NTRS1991

Use of non-spherical gravity harmonics for relative motion GN/C (category 2)

Detailed analysis of the Automatic Rendezvous and Capture problem indicate a need for three different regions of mathematical description for the GN&C algorithms: (1) multi-vehicle orbital mechanics to the rendezvous interface point, i.e., within 100 nm; (2) relative motion solutions (such as Clohessy-Wiltshire type) from the far-field to the near-field interface, i.e., within 1 nm and; (3) close proximity motion - the near-field motion where the relative differences in the gravitational and orbit inertial accelerations can be neglected from the equations of motion. Limit boundaries to these regions can be precisely defined by further analysis and will be functions of the tracking measurement accuracies and the computer resources available for the solution of the algorithms. This paper analyzes the relative motion in Regions 2 and 3 above and present the derivation and discussion of the general case of non-spherical gravitational perturbed relative motion. Mathematical deviations from the numerically integrated spherical gravity case and solutions from the Clohessy-Wiltshire equations are presented in the analysis. Based upon this preliminary analysis, it is recommended that further efforts be used to assess the relative position and velocity differences in Region 2 due to non-spherical gravity harmonics and that viable GN&C algorithms be developed to include these gravity perturbations (especially the effects of the first gravity harmonic, J2).

Henderson, David M.↗
NASA NTRS1981

Aerodynamic prediction techniques for hypersonic configuration design

An investigation of approximate theoretical techniques for predicting aerodynamic characteristics and surface pressures for relatively slender vehicles at moderate hypersonic speeds was performed. Emphasis was placed on approaches that would be responsive to preliminary configuration design level of effort. Potential theory was examined in detail to meet this objective. Numerical pilot codes were developed for relatively simple three dimensional geometries to evaluate the capability of the approximate equations of motion considered. Results from the computations indicate good agreement with higher order solutions and experimental results for a variety of wing, body, and wing-body shapes for values of the hypersonic similarity parameter M delta approaching one.

Source record↗
NASA NTRS1972

Shock waves and drag in the numerical calculation of isentropic transonic flow

Properties of the shock relations for steady, irrotational, transonic flow are discussed and compared for the full and approximate governing potential in common use. Results from numerical experiments are presented to show that the use of proper finite difference schemes provide realistic solutions and do not introduce spurious shock waves. Analysis also shows that realistic drags can be computed from shock waves that occur in isentropic flow. In analogy to the Oswatitsch drag equation, which relates the drag to entropy production in shock waves, a formula is derived for isentropic flow that relates drag to the momentum gain through an isentropic shock. A more accurate formula for drag, based on entropy production, is also derived, and examples of wave drag evaluation based on these formulas are given and comparisons are made with experimental results.

Steger, J. L.↗
NASA NTRS1983

Nonlinear potential analysis techniques for supersonic-hypersonic configuration design

Approximate nonlinear inviscid theoretical techniques for predicting aerodynamic characteristics and surface pressures for relatively slender vehicles at moderate hypersonic speeds were developed. Emphasis was placed on approaches that would be responsive to preliminary configuration design level of effort. Second order small disturbance and full potential theory was utilized to meet this objective. Numerical pilot codes were developed for relatively general three dimensional geometries to evaluate the capability of the approximate equations of motion considered. Results from the computations indicate good agreement with higher order solutions and experimental results for a variety of wing, body and wing-body shapes for values of the hypersonic similarity parameter M delta approaching one. Case computational times of a minute were achieved for practical aircraft arrangements.

Clever, W. C.↗
NASA NTRS1985

Statistical magnetic anomalies from satellite measurements for geologic analysis

The errors of numerically averaging satellite magnetic anomaly data for geologic analysis are investigated using orbital anomaly simulations of crustal magnetic sources by Gauss-Legendre quadrature integration. These simulations suggest that numerical averaging errors constitute small and relatively minor contributions to the total error-budget of higher orbital estimates (approx. 400 km), whereas for lower orbital estimates the error of averaging may increase substantially. Least-squares collocation is also investigated as an alternative to numerical averaging and found to produce substantially more accurate anomaly estimates as the elevation of prediction is decreased towards the crustal sources.

Goyal, H. K.↗
NASA NTRS2003

First Results from a Hardware-in-the-Loop Demonstration of Closed-Loop Autonomous Formation Flying

A closed-loop system for the demonstration of autonomous satellite formation flying technologies using hardware-in-the-loop has been developed. Making use of a GPS signal simulator with a dual radio frequency outlet, the system includes two GPS space receivers as well as a powerful onboard navigation processor dedicated to the GPS-based guidance, navigation, and control of a satellite formation in real-time. The closed-loop system allows realistic simulations of autonomous formation flying scenarios, enabling research in the fields of tracking and orbit control strategies for a wide range of applications. The autonomous closed-loop formation acquisition and keeping strategy is based on Lyapunov's direct control method as applied to the standard set of Keplerian elements. This approach not only assures global and asymptotic stability of the control but also maintains valuable physical insight into the applied control vectors. Furthermore, the approach can account for system uncertainties and effectively avoids a computationally expensive solution of the two point boundary problem, which renders the concept particularly attractive for implementation in onboard processors. A guidance law has been developed which strictly separates the relative from the absolute motion, thus avoiding the numerical integration of a target trajectory in the onboard processor. Moreover, upon using precise kinematic relative GPS solutions, a dynamical modeling or filtering is avoided which provides for an efficient implementation of the process on an onboard processor. A sample formation flying scenario has been created aiming at the autonomous transition of a Low Earth Orbit satellite formation from an initial along-track separation of 800 m to a target distance of 100 m. Assuming a low-thrust actuator which may be accommodated on a small satellite, a typical control accuracy of less than 5 m has been achieved which proves the applicability of autonomous formation flying techniques to formations of satellites as close as 50 m.

Gill, E.↗
NASA NTRS1978

NASA aviation safety reporting system

A sample of reports relating to operations during winter weather is presented. Several reports involving problems of judgment and decisionmaking have been selected from the numerous reports representative of this area. Problems related to aeronautical charts are discussed in a number of reports. An analytic study of reports involving potential conflicts in the immediate vicinity of uncontrolled airports was performed; the results are discussed in this report. It was found that in three-fourths of 127 such conflicts, neither pilot, or only one of the pilots, was communicating position and intentions on the appropriate frequency. The importance of providing aural transfer of information, as a backup to the visual see and avoid mode of information transfer is discussed. It was also found that a large fraction of pilots involved in potential conflicts on final approach had executed straight-in approaches, rather than the recommended traffic pattern entries, prior to the conflicts. A selection of alert bulletins and responses to them by various segments of the aviation community is presented.

Source record↗
NASA NTRS1986

Statistics of voids in hierarchical universes

As one alternative to the N-point galaxy correlation function statistics, the distribution of holes or the probability that a volume of given size and shape be empty of galaxies can be considered. The probability of voids resulting from a variety of hierarchical patterns of clustering is considered, and these are compared with the results of numerical simulations and with observations. A scaling relation required by the hierarchical pattern of higher order correlation functions is seen to be obeyed in the simulations, and the numerical results show a clear difference between neutrino models and cold-particle models; voids are more likely in neutrino universes. Observational data do not yet distinguish but are close to being able to distinguish between models.

Fry, J. N.↗
NASA NTRS2012

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗
NASA NTRS1977

Application of the method of integral relations to laminar boundary layers in three dimensions

The method of integral relations is extended to general three-dimensional compressible laminar boundary layer flows. The transformation employed to transform the basic three-dimensional compressible boundary layer equations into quasi-incompressible form is an extension of the Howarth transformation. The resulting system of differential equations is integrated numerically by the method of integral relations as proposed by Dorodnitsyn. To demonstrate the accuracy of the method, it is applied to calculation of the parabolic flow over a flat plate and the boundary flow over an infinite yawed cylinder, for which solutions are known. It is then applied to the flow over a flat plate disturbed by a cylinder normal to the plate, for which a finite-difference solution is available for comparison. It is finally applied to calculating the crossflow velocity variation for supersonic flow over a swept wedge.

Holt, M.↗
NASA NTRS2011

NASA Water Resources Program

With increasing population pressure and water usage coupled with climate variability and change, water issues are being reported by numerous groups as the most critical environmental problems facing us in the 21st century. Competitive uses and the prevalence of river basins and aquifers that extend across boundaries engender political tensions between communities, stakeholders and countries. In addition to the numerous water availability issues, water quality related problems are seriously affecting human health and our environment. The potential crises and conflicts especially arise when water is competed among multiple uses. For example, urban areas, environmental and recreational uses, agriculture, and energy production compete for scarce resources, not only in the Western U.S. but throughout much of the U.S. and also in numerous parts of the world. Mitigating these conflicts and meeting water demands and needs requires using existing water resources more efficiently. The NASA Water Resources Program Element works to use NASA products and technology to address these critical water issues. The primary goal of the Water Resources is to facilitate application of NASA Earth science products as a routine use in integrated water resources management for the sustainable use of water. This also includes the extreme events of drought and floods and the adaptation to the impacts from climate change. NASA satellite and Earth system observations of water and related data provide a huge volume of valuable data in both near-real-time and extended back nearly 50 years about the Earth's land surface conditions such as precipitation, snow, soil moisture, water levels, land cover type, vegetation type, and health. NASA Water Resources Program works closely to use NASA and Earth science data with other U.S. government agencies, universities, and non-profit and private sector organizations both domestically and internationally. The NASA Water Resources Program organizes its projects under five functional themes. I) Streamflow and Flood Forecasting 2) Water Supply and Irrigation (includes evapotranspiration) 3) Drought 4) Water Quality 5) Climate and Water Resources. To maximize this activity NASA Water Resources Program works closely with other government agencies (e.g., the National Oceanic and Atmospheric Administration (NOAA); the U.S. Department of Agriculture (USDA); the U.S. Geological Survey (USGS); the Environmental Protection Agency (EPA), USAID, the Air Force Weather Agency (AFWA)), universities, non-profit national and international organizations, and the private sector. The NASA Water Resources program currently is funding 21 active projects under the functional themes (http://wmp.gsfc.nasa.gov & http://science.nasa.gov/earth-science/applied-sciences/).

Toll, David L.↗
NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗
NASA NTRS2025

Robust Trajectory Optimization for NRHO Rendezvous Using SPICE Kernel Relative Motion

In this paper, robust optimization is performed on trajectory correction maneuvers during the lunar lander return phase of an Artemis mission, treating the trajectory from one hour after low lunar orbit departure to arrival in the vicinity of the lunar Gateway as a relative motion problem. To enable rapid stochastic optimization techniques requiring many candidate trajectories, SPICE kernel relative motion as implemented by the Quadratic Interpolated State Transition (QIST) system is used as the underlying dynamics propagation. The optimization is performed with a genetic optimizer using linear covariance (LinCov) software in a simplified operational context, taking into account the availability of navigation sensors with varying measurement models, ranges, and accuracies. No numerical integration is used, since the relative motion around Gateway is fully characterized with the a priori computation of the QIST coefficients. Maneuver placements are computed to optimize the minimum 3σ delta-v of the trajectory, the position dispersion at a target point, and a convex combination of these two metrics. An order of magnitude runtime improvement is provided over legacy methods with less than 10% error introduced. All QIST results are shown to be in-family with legacy methods. The tradespace for optimal delta-v design is found to range from 77.0 to 93.9 m/s, while the range of optimal dispersion is between 1.4 and 11.7 km.

Relative Motion↗
NASA NTRS1980

Numerical simulation of advection fog formation on multi-disperse aerosols due to combustion-related pollutants

The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.

Hung, R. J.↗
NASA NTRS1991

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗
NASA NTRSDate not supplied

In-Depth Analysis of Subject Functional Performance and Subjective Feedback Data for Evaluating Argos Simulation Quality

Mechanical offloading systems such as the Active Response Gravity Offload System (ARGOS) at NASA Johnson Space Center (JSC) are used to simulate various partial gravity environments. ARGOS offloads pressurized suited subjects through a gimbaling pivot point system attached to a volumetric Portable Life Support System (PLSS). The pivot point can be configured to multiple locations with respect to the human-suit system center of gravity (CG). Previous work has tested and documented the interaction between different pivot point locations and measured CG. Observations indicated that small changes in pivot point location substantially affect the stability and functional performance (simulation quality) of a subject. This led to the development of a standardized gimbal assessment protocol to assess the functional performance of a pressurized suited subject as a function of the pivot point location, using the data obtained from a variety of planetary suit mobility tasks. As a result, an expansive repository of CG-related task metrics was generated from numerous ARGOS spacesuit test events, for different gimbal configurations and different test subjects. Based on the evaluation metrics and observed trends, several pivot point locations were iteratively identified to determine an “appropriate/optimal” configuration. This study thus describes the in-depth analysis of the CG evaluation metrics to better quantify trends and determine if specific factors are strongly associated with appropriate/optimal configurations. Quantitative and qualitative variables relating to subject task performance, subjective feedback, and anthropometry will be evaluated via probabilistic methods. The results from this study are expected to improve our understanding of optimal ARGOS gimbal settings, which will better inform the gimbal configuration identification process and improve simulation quality for extravehicular activity (EVA) training.

Joseph Yao↗