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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Evolution of Astrodynamics for Multibody Environments, Numerical Precision and Dynamical Understanding

The innovative nature of dynamical systems theory and its applications in trajectory design enabled many recent missions by opening up new transfer and capture opportunities to meet mission requirements, reduce fuel costs, and address various launch constraints and contingencies. With a combination of proven numerical methods and an advanced understanding of the dynamics, an ‘experienced’ user can rapidly design multibody and cislunar missions such as Lunar IceCube and Artemis/Themis and relatively quickly conceive trajectory concepts for Sun-Earth mission orbits, such as the Space Weather Follow On-L1 Lissajous orbit—one that can be achieved and maintained.

Dynamical Systems↗
NASA NTRSDate not supplied

A Novel Machine Learning Method for Surface PM2.5 Estimations from Geostationary Satellites

Particulate matter (PM) with a diameter of less or equal to 2.5 μm, known as PM , affects human health as it penetrates the respiratory system. The Environmental Protection Agency (EPA) measures the atmospheric concentration of PM using air quality monitors stationed throughout the Continental United States (CONUS). Such measurements are points on a spatial domain and therefore, might not be representative of the air quality at nearby areas considering that the composition of the atmosphere is highly variable from place to place. Satellite based AOD permits a spatially uniform means of estimating PM and new geostationary satellites provide high temporal and spatial resolution estimation of AOD. However, the concentration of PM is non-linearly dependent on other atmospheric parameters that include relative humidity, temperature, and height of the planetary boundary layer. This information may be estimated at similar spatial and temporal resolutions as AOD from numerical modeling such as from the National Oceanic and Atmospheric Administration’s (NOAA) High Resolution Rapid Refresh (HRRR) model which resolves near real-time atmospheric conditions over the CONUS. The estimation of PM concentration is a multi-parametric problem that considers the effect of temporal dependencies among the different parameters. Deep learning approaches are appropriate for such complex estimation problems as they intrinsically capture relations among multiple non-linear parameters. This study compares deep-learning methods to traditional regression analysis to demonstrate the capabilities of these methods in predicting PM2.5 concentrations. Additionally, a novel ensemble learning approach is employed to identify scientific processes that could further improve the estimation of PM concentration. Utilizing Long Short-Term Memory (LSTM) neural networks, which are suitable for multivariate time series estimation problems as they are capable of learning long-term dependencies, individual models are created for each EPA station and trained on the aforementioned dataset collocated over each station. Individual station models are merged if the model's performance is improved by reducing the root mean squared error (RMSE) metric. This ensemble training method ultimately reduces the RMSE value. Evaluation of these results provide insights into physical processes and related observable parameters that may contribute to PM concentrations. Identified parameters evaluated to be statistically different between the merged and unmerged models are expected to improve overall performance. These new parameters are then utilized for reevaluation of the deep learning methods with an extreme gradient boosting model with an RMSE of 5.5 providing the best results.

George Priftis↗
NASA NTRS1998

Study of Interesting Solidification Phenomena on the Ground and in Space (MEPHISTO)

Real-time Seebeck voltage variations in a Sn-Bi melt during directional solidification in the MEPHISTO spaceflight experiment flown on the USMP-3 mission, can be correlated with well characterized thruster firings and an Orbiter Main System (OMS) burn. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-crystal interface. This allowed us to make a direct comparison of numerical simulations with the experimentally obtained Seebeck signals. Based on the results of preflight and real-time computations, several well-defined thruster firing events were programmed to occur at specific times during the experiment. In particular, we simulated the effects of the thruster firings on melt and crystal composition in a directionally solidifying Sn-Bi alloy. The relative accelerations produced by the firings were simulated by impulsive accelerations of the same magnitude, duration and orientation as the requested firings. A comparison of the simulation results with the Seebeck signal indicates that there is a good agreement between the two. This unique opportunity allows us, for the first time, to quantitatively characterize actual g-jitter effects on an actual crystal growth experiment and to properly calibrate our models of g-jitter effects on crystal growth.

Favier, J.-J.↗
NASA NTRS1999

Study of Interesting Solidification Phenomena on the Ground and in Space (MEPHISTO)

Real-time Seebeck voltage variations in a Sn-Bi melt during directional solidification in the MEPHISTO spaceflight experiment flown on the USMP-3 mission, have been correlated with well-characterized thruster firings and an Orbiter Main System (OMS) burn. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-crystal interface. This allowed us to make a direct comparison of numerical simulations with the experimentally obtained Seebeck signals. Based on the results of preflight and real-time computations, several well-defined thruster firing events were programmed to occur at specific times during the experiment. In particular, we simulated the effects of the thruster firings on melt and crystal composition in a directionally solidifying Sn-Bi alloy. The relative accelerations produced by the firings were simulated by impulsive accelerations of the same magnitude, duration and orientation as the requested firings. A comparison of the simulation results with the Seebeck signal indicates that there is a good agreement between the two. This unique opportunity allows us to make the first quantitative characterization of actual g-jitter effects on an actual crystal growth experiment and to calibrate our models of g-jitter effects on crystal growth.

Alexander, J. Iwan D.↗
NASA NTRS1980

The Finite Analytic Method for steady and unsteady heat transfer problems

A new numerical method called the Finite Analytical Method for solving partial differential equations is introduced. The basic idea of the finite analytic method is the incorporation of the local analytic solution in obtaining the numerical solution of the problem. The finite analytical method first divides the total region of the problem into small subregions in which local analytic solutions are obtained. Then an algebraic equation is derived from the local analytic solution for each subregion relating an interior nodal value at a point P in the subregion to its neighboring nodal values. The assembly of all the local analytic solutions thus provides the finite-analytic numerical solution of the problem. In this paper the finite analytic method is illustrated in solving steady and unsteady heat transfer problems.

Chen, C.-J.↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
NASA NTRS2013

Towards Bridging the Gaps in Holistic Transition Prediction via Numerical Simulations

The economic and environmental benefits of laminar flow technology via reduced fuel burn of subsonic and supersonic aircraft cannot be realized without minimizing the uncertainty in drag prediction in general and transition prediction in particular. Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper provides a summary of selected research activities targeting the current gaps in high-fidelity transition prediction, specifically those related to the receptivity and laminar breakdown phases of crossflow induced transition in a subsonic swept-wing boundary layer. The results of direct numerical simulations are used to obtain an enhanced understanding of the laminar breakdown region as well as to validate reduced order prediction methods.

Choudhari, Meelan M.↗
NASA NTRS1996

Analysis of Residual Acceleration Effects on Transport and Segregation During Directional Solidification of Tin-Bismuth in the MEPHISTO Furnace Facility

The objective of this work is to approach the problem of determining the transport conditions (and effects of residual acceleration) during the plane-front directional solidification of a tin-bismuth alloy under low gravity conditions. The work involves using a combination of 2- and 3-D numerical models, scaling analyses, ID models and the results of ground-based and low-gravity experiments. The latter are to be conducted during the MEPHISTO experiment scheduled for USMP-3 in early 1996. The models will be used to predict the response of the transport conditions and consequent solute segregation in directionally solidifying tin-bismuth melt. Real-time Seebeck voltage variations across a Sn-Bi melt during directional solidification in MEPHISTO on USMP-1 show a distinct variation which can be correlated with thruster firings. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-solid interface. This allows a direct comparison of numerical simulations with the Seebeck signals obtained on USMP-1. The effects of such accelerations on composition for a directionally solidifying Sn-Bi alloy have been simulated numerically. USMP-1 acceleration data was used to assist in our choice of acceleration magnitude and orientation. The results show good agreement with experimental observations. The USMP-3 experiments took place earlier this year (February 22 through March 6). There were several differences between the USMP-3 experiments as compared to USMP-1. Firstly a more concentrated alloy was solidified and, secondly, Primary Reaction Control System thruster burns were requested at particular times during four separate growth runs. This allowed us to monitor the response Seebeck response under well-characterized growth conditions. In addition, we carried out simulations during the experiment in order to interpret the Seebeck signal. Preliminary results are described here.

Alexander J. Iwan D.↗
NASA NTRS1996

Refraction and Shielding of Noise in Non-Axisymmetric Jets

This paper examines the shielding effect of the mean flow and refraction of sound in non-axisymmetric jets. A general three-dimensional ray-acoustic approach is applied. The methodology is independent of the exit geometry and may account for jet spreading and transverse as well as streamwise flow gradients. We assume that noise is dominated by small-scale turbulence. The source correlation terms, as described by the acoustic analogy approach, are simplified and a model is proposed that relates the source strength to 7/2 power of turbulence kinetic energy. Local characteristics of the source such as its strength, time- or length-scale, convection velocity and characteristic frequency are inferred from the mean flow considerations. Compressible Navier Stokes equations are solved with a k-e turbulence model. Numerical predictions are presented for a Mach 1.5, aspect ratio 2:1 elliptic jet. The predicted sound pressure level directivity demonstrates favorable agreement with reported data, indicating a relative quiet zone on the side of the major axis of the elliptic jet.

Khavaran, Abbas↗
NASA NTRS1998

One Year Report for SAMS and OARE on STS-73/USML-2

The Second United States Microgravity Laboratory (USML-2) payload flew on the orbiter Columbia on mission STS-73 from October 20 to November 5, 1995. The USML-2 payload on STS-73 was dedicated to microgravity experiments. Two accelerometer systems managed by the NASA Lewis Research Center (LeRC) flew to support these experiments, namely the Orbital Acceleration Research Experiment (OARE) and the Space Acceleration Measurements System (SAMS). OARE downlinked real-time quasi-steady acceleration data, which were provided to the investigators. The SAMS recorded higher frequency data onboard for post-mission analysis. The Principal Investigator Microgravity Services (PIMS) project at NASA LeRC supports principal investigators of microgravity experiments as they evaluate the effects of varying acceleration levels on their experiments. A summary report was prepared by PIMS to furnish interested experiment investigators with a guide to evaluate the acceleration environment during STS-73, and as a means of identifying areas which require further study. The summary report provides an overview of the STS-73 mission, describes the accelerometer systems flown on this mission, discusses some specific analyses of the accelerometer data in relation to the various activities which occurred during the mission, and presents plots resulting from these analyses as a snapshot of the environment during the mission. Numerous activities occurred during the STS-73 mission that are of interest to the low-gravity community. Specific activities of interest during this mission were crew exercise, payload bay door motion, Glovebox fan operations, water dumps, Ku band antenna activity, orbital maneuvering system, and primary reaction control system firings, and attitude changes. The low-gravity environment related to these activities is discussed in the summary report.

Hakimzadeh, Roshanak↗
NASA NTRS1972

The influence of laser ranging on selenodetic control.

In this paper, a mathematical model is presented which can be used to perform an adjustment using laser distances to improve coordinates of points of the earth and on the moon. The observation equations are given and the orientation of the earth and the moon is defined in terms of three Eulerian angles. Parameters related to the orientation of the two bodies are the six initial conditions in each case, for the numerical integration of the three angles and their time derivatives. Numerical experiments were performed for the purpose of investigating the expected accuracies of the various parameters in the adjustment model under differing conditions. The results of the experiments indicate that the statistics for both the coordinates of points on earth and on the moon, and for the orientation parameters of the two bodies, are favorable.

Mueller, I. I.↗
NASA NTRS1972

Optimum performance of static propellers and rotors.

A criterion for the optimum performance of static propellers and hovering helicopter rotors is developed. Numerical results are presented for the optimum radial distributions of circulation and inflow at the propeller disk and relating the optimum power coefficient and the figure of merit to the thrust coefficient. It is shown that the present theory, which fully accounts for the effect of slipstream rotation, predicts optimum distributions of circulation and inflow that differ significantly from those based on approximate methods.

Wu, J. C.↗
NASA NTRS1979

Delay compensation - Its effect in reducing sampling errors in Fourier spectroscopy

An approximate formula is derived for the spectrum ghosts caused by periodic drive speed variations in a Michelson interferometer. The solution represents the case of fringe-controlled sampling and is applicable when the reference fringes are delayed to compensate for the delay introduced by the electrical filter in the signal channel. Numerical results are worked out for several common low-pass filters. It is shown that the maximum relative ghost amplitude over the range of frequencies corresponding to the lower half of the filter band is typically 20 times smaller than the relative zero-to-peak velocity error, when delayed sampling is used. In the lowest quarter of the filter band it is more than 100 times smaller than the relative velocity error. These values are ten and forty times smaller, respectively, than they would be without delay compensation if the filter is a 6-pole Butterworth.

Zachor, A. S.↗
NASA NTRS1986

Turbulent relaxation to a force-free field-reversed state

The evolution of nonequilibrium initial conditions of an incompressible magnetohydrodynamic Z pinch is described by a three-dimensional, pseudospectral numerical code. Magnetohydrodynamic turbulence develops in the resistive, nonviscous magnetofluid, resulting in the selective decay of the energy relative to the magnetic helicity, at Lundquist numbers of only a few hundred. An interior force-free region grows with time and achieves spontaneous reversal of the toroidal magnetic field at the wall, without the necessity of an external electric field.

Dahlburg, J. P.↗
NASA NTRS1986

Detailed near-wake flowfield surveys with comparison to an Euler method of an aspect ratio 4 rectangular wing

An experimental investigation of the flowfield in the near-wake of an aspect ratio 4 rectangular wing was conducted, providing a complete detailed set of data for use in the validation of computational methods. An angle of attack of 8 degrees and two Reynolds numbers 530,000 and 391,000 were investigated using pitot and six-hole probes. In addition, two types of flow visualization were employed. The data presented includes contours of total pressure, mean velocity, flow angularity, and vorticity distribution data at five chordwise stations of the near-wake ranging from 0.167 to 5.00 chord lengths aft of the trailing edge. The experimental results were compared to the predicted results of a 2-D Euler numerical method. The results predicted by an Euler method failed to accurately define the flowfield. Tangential velocities remained relatively constant over the range of X/C considered though increased in angle of attack and Reynolds number did bring about corresponding increases. Axial velocities also increased with angle of attack and Reynolds number but showed greater sensitivity to increases in X/C. Graphic displays and contours of the total pressure data indicate that roll-up of the wing tip vortex is essentially complete one and one half chords downstream of the trailing edge.

Klinge, M. D.↗
NASA NTRS1988

Fatigue testing and damage development in continuous fiber reinforced metal matrix composites

A general overview of the fatigue behavior of metal matrix composites (MMC) is presented. The first objective is to present experimental procedures and techniques for conducting a meaningful fatigue test to detect and quantify fatigue damage in MMC. These techniques include interpretation of stress-strain responses, acid etching of the matrix, edge replicas of the specimen under load, radiography, and micrographs of the failure surfaces. In addition, the paper will show how stiffness loss in continuous fiber reinforced metal matrix composites can be a useful parameter for detecting fatigue damage initiation and accumulation. Second, numerous examples of how fatigue damage can initiate and grow in various MMC are given. Depending on the relative fatigue behavior of the fiber and matrix, and the interface properties, the failure modes of MMC can be grouped into four categories: (1) matrix dominated, (2) fiber dominated, (3) self-similar damage growth, and (4) fiber/matrix interfacial failures. These four types of damage will be discussed and illustrated by examples with the emphasis on the fatigue of unnotched laminates.

Johnson, W. S.↗
NASA NTRS1988

Focal region fields of distorted reflectors

The problem of the focal region fields scattered by an arbitrary surface reflector under uniform plane wave illumination is solved. The physical optics (PO) approximation is used to calculate the current induced on the reflector. The surface of the reflector is described by a number of triangular domain-wise 5th degree bivariate polynomials. A 2-dimensional Gaussian quadrature is employed to numerically evaluate the integral expressions of the scattered fields. No Freshnel or Fraunhofer zone approximations are made. The relation of the focal fields problem to surface compensation techniques and other applications are mentioned. Several examples of distorted parabolic reflectors are presented. The computer code developed is included, together with instructions on its usage.

Buris, N. E.↗