Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 847 records · Page 47

Reynolds averaged Navier-Stokes computations of transonic flows The-state-of-the-art

The state of the art in computer simulations for transonic flowfields requiring solutions for the Navier-Stokes equations is assessed. It is noted that current simulations of transonic flowfields require comparisons with experimental results because the simulations are not free from discretization errors. Cases of turbulence are treated with weighted variables in a time-averaged scheme to yield Reynolds averaged Navier-Stokes equations. The turbulence is modeled in a first-order approach with a Reynolds stress tensor or a second-order approach where the tensor is obtained from the Navier-Stokes equations. The choice of coordinate systems and dependent variables for simulating the flow around airfoils is discussed, with particular attention to curvilinear coordinates. The determinations of boundary conditions is examined, along with numerical methods related to physical phenomena

Mehta, U.↗

Alternative Modal Basis Selection Procedures for Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of the three reduced-order analyses are compared with the results of the computationally taxing simulation in the physical degrees of freedom. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

Dynamic Emulation of NASA Missions for IVandV: A Case Study of JWST and SLS

Software-Only-Simulations are an emerging but quickly developing field of study throughout NASA. The NASA Independent Verification Validation (IVV) Independent Test Capability (ITC) team has been rapidly building a collection of simulators for a wide range of NASA missions. ITC specializes in full end-to-end simulations that enable developers, VV personnel, and operators to test-as-you-fly. In four years, the team has delivered a wide variety of spacecraft simulations ranging from low complexity science missions such as the Global Precipitation Management (GPM) satellite and the Deep Space Climate Observatory (DSCOVR), to the extremely complex missions such as the James Webb Space Telescope (JWST) and Space Launch System (SLS).

Emulation↗

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗

High-Resolution WRF Forecasts of Lightning Threat

Tropical Rainfall Measuring Mission (TRMM)lightning and precipitation observations have confirmed the existence of a robust relationship between lightning flash rates and the amount of large precipitating ice hydrometeors in storms. This relationship is exploited, in conjunction with the capabilities of the Weather Research and Forecast (WRF) model, to forecast the threat of lightning from convective storms using the output fields from the model forecasts. The simulated vertical flux of graupel at -15C is used in this study as a proxy for charge separation processes and their associated lightning risk. Initial experiments using 6-h simulations are conducted for a number of case studies for which three-dimensional lightning validation data from the North Alabama Lightning Mapping Array are available. The WRF has been initialized on a 2 km grid using Eta boundary conditions, Doppler radar radial velocity and reflectivity fields, and METAR and ACARS data. An array of subjective and objective statistical metrics is employed to document the utility of the WRF forecasts. The simulation results are also compared to other more traditional means of forecasting convective storms, such as those based on inspection of the convective available potential energy field.

Goodman, S. J.↗

Simulation of Unsteady Hypersonic Combustion Around Projectiles in an Expansion Tube

The temporal evolution of combustion flowfields established by the interaction between wedge-shaped bodies and explosive hydrogen-oxygen-nitrogen mixtures accelerated to hypersonic speeds in an expansion tube is investigated. The analysis is carried out using a fully implicit, time-accurate, computational fluid dynamics code that we developed recently for solving the Navier-Stokes equations for a chemically reacting gas mixture. The numerical results are compared with experimental data from the Stanford University expansion tube for two different gas mixtures at Mach numbers of 4.2 and 5.2. The experimental work showed that flow unstart occurred for the Mach 4.2 cases. These results are reproduced by our numerical simulations and, more significantly, the causes for unstart are explained. For the Mach 5.2 mixtures, the experiments and numerical simulations both produced stable combustion. However, the computations indicate that in one case the experimental data were obtained during the transient phase of the flow; that is, before steady state had been attained.

Yungster, S.↗

Investigation of Volumetric Sources in Airframe Noise Simulations

Hybrid methods for the prediction of airframe noise involve a simulation of the near field flow that is used as input to an acoustic propagation formula. The acoustic formulations discussed herein are those based on the Ffowcs Williams and Hawkings equation. Some questions have arisen in the published literature in regard to an apparently significant dependence of radiated noise predictions on the location of the integration surface used in the solution of the Ffowcs Williams and Hawkings equation. These differences in radiated noise levels are most pronounced between solid-body surface integrals and off-body, permeable surface integrals. Such differences suggest that either a non-negligible volumetric source is contributing to the total radiation or the input flow simulation is suspect. The focus of the current work is the issue of internal consistency of the flow calculations that are currently used as input to airframe noise predictions. The case study for this research is a computer simulation for a three-element, high-lift wing profile during landing conditions. The noise radiated from this flow is predicted by a two-dimensional, frequency-domain formulation of the Ffowcs Williams and Hawkings equation. Radiated sound from volumetric sources is assessed by comparison of a permeable surface integration with the sum of a solid-body surface integral and a volume integral. The separate noise predictions are found in good agreement.

Casper, Jay H.↗

Simplified Worst-Case Analysis of Complex Systems

Statistical method avoids excessive computer time. Technique involves simplified Monte Carlo simulation of system with randomly chosen parameters and comparison of tolerance extremes of several of worst-case situations found. Resulting combination of tolerance extremes then used in detailed analysis - one that makes use of full complex nonlinear model capable of accurate simulations.

Pond, C. L.↗

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗

Spacecraft Trajectory Analysis and Mission Planning Simulation (STAMPS) Software

STAMPS simulates either three- or six-degree-of-freedom cases for all spacecraft flight phases using translated HAL flight software or generic GN&C models. Single or multiple trajectories can be simulated for use in optimization and dispersion analysis. It includes math models for the vehicle and environment, and currently features a "C" version of shuttle onboard flight software. The STAMPS software is used for mission planning and analysis within ascent/descent, rendezvous, proximity operations, and navigation flight design areas.

Puckett, Nancy↗

Study of Graphite Ablation at Arc-Jet Conditions using Finite-Rate and Equilibrium Chemistry Models

Moving at hypersonic speeds through the atmosphere generates immense amount of heat that triggers thermo-chemical non-equilibrium effects in the flow approaching the vehicle and leads to surface chemical reactions that cause ablation and consume the material mass. This complex behavior and interaction of the aerothermal environment with a graphite surface is com- putationally modeled in the present work using a fully-coupled simulation between flow solver CHAMPS NBS-Cart and material response solver KATS-MR. In this study, we account for the thermo-chemical non-equilibrium effects occurring in the flow and use an 11-species air-carbon re- action model to simulate the ablation of the material. The reactions at the surface such oxidation, nitridation and recombination are modeled with a newly developed air-carbon reaction model. The coupled simulation is applied to study two experimental cases performed at IHF arc-jet facility at NASA Ames research center. The obtained results are compared against the measured data of surface recession, as well as surface and in-depth temperatures. Additionally, uncoupled sim- ulations, performed with equilibrium chemistry at the wall, point out the substantial differences between the coupled and uncoupled solutions, leading to higher deviation of the predicted temper- atures. Overall, we show excellent agreement between the coupled simulation predictions and the experimental data in all measured quantities.

EDL↗

An automated procedure for developing hybrid computer simulations of turbofan engines

A systematic, computer-aided, self-documenting methodology for developing hybrid computer simulations of turbofan engines is presented. The methodology makes use of a host program that can run on a large digital computer and a machine-dependent target (hybrid) program. The host program performs all of the calculations and date manipulations needed to transform user-supplied engine design information to a form suitable for the hybrid computer. The host program also trims the self contained engine model to match specified design point information. A test case is described and comparisons between hybrid simulation and specified engine performance data are presented.

Szuch, J. R.↗

Vibration-dissociation coupling in nonequilibrium flows

The purpose of the study is to evaluate all important physical phenomena contributing to the vibration-dissociation problem by developing a numerical technique for the simulation of the vibrational relaxation and dissociation/recombination of N(2) in an isothermal box. A set of coupled master equations representing each vibrational level of N(2) is solved, and both V-V and V-T exchanges between molecules and atoms are included, along with the effects of dissociation and recombustion from/to each level. It is observed that the inverse range parameter has a strong impact on the vibrational exchange rates and relaxation. In a heating case, the lower vibrational levels relax through a series of equilibrium distributions, while in a cooling case, a population inversion is formed in the upper levels. Other cases of heating and cooling will be simulated, and together with this work, used to develop a matrix of nonequilibrium relaxation runs.

Landrum, D. Brian↗

An automated procedure for developing hybrid computer simulations of turbofan engines

This paper offers a systematic, computer-aided, self-documenting methodology for developing hybrid computer simulations of turbofan engines. The methodology that is presented makes use of a host program that can run on a large digital computer and a machine-dependent target (hybrid) program. The host program performs all of the calculations and data manipulations that are needed to transform user-supplied engine design information to a form suitable for the hybrid computer. The host program also trims the self-contained engine model to match specified design point information. A test case is described and comparisons between hybrid simulation and specified engine performance data are presented.

Szuch, J. R.↗

Numerically simulated two-dimensional auroral double layers

A magnetized 2 1/2-dimensional particle-in-cell system which is periodic in one direction and bounded by reservoirs of Maxwellian plasma in the other is used to numerically simulate electrostatic plasma double layers. For the cases of both oblique and two-dimensional double layers, the present results indicate periodic instability, Debye length rather than gyroradii scaling, and low frequency electrostatic turbulence together with electron beam-excited electrostatatic electron-cyclotron waves. Estimates are given for the thickness of auroral doule layers, as well as the separations within multiple auroral arcs. Attention is given to the temporal modulation of accelerated beams, and the possibilities for ion precipitation and ion conic production by the double layer are hypothesized. Simulations which include the atmospheric backscattering of electrons imply the action of an ionospheric sheath which accelerates ionospheric ions upward.

Borovsky, J. E.↗

Altimetry data and the elastic stress tensor of subduction zones

The maximum shear stress (mss) field due to mass anomalies is estimated in the Apennines, the Kermadec-Tonga Trench, and the Rio Grande Rift areas and the results for each area are compared to observed seismicity. A maximum mss of 420 bar was calculated in the Kermadec-Tonga Trench region at a depth of 28 km. Two additional zones with more than 300 bar mss were also observed in the Kermadec-Tonga Trench study. Comparison of the calculated mss field with the observed seismicity in the Kermadec-Tonga showed two zones of well correlated activity. The Rio Grande Rift results showed a maximum mss of 700 bar occurring east of the rift and at a depth of 6 km. Recorded seismicity in the region was primarily constrained to a depth of approximately 5 km, correlating well to the results of the stress calculations. Two areas of high mss are found in the Apennine region: 120 bar at a depth of 55 km, and 149 bar at the surface. Seismic events observed in the Apennine area compare favorably with the mss field calculated, exhibiting two zones of activity. The case of loading by seamounts and icecaps are also simulated. Results for this study show that the mss reaches a maximum of about 1/3 that of the applied surface stress for both cases, and is located at a depth related to the diameter of the surface mass anomaly.

Caputo, Michele↗

Numerical Simulation of Screech Tones from Supersonic Jets: Physics and Prediction

The objectives of this project are to: (1) perform a numerical simulation of the jet screech phenomenon; and (2) use the data of the simulations to obtain a better understanding of the physics of jet screech. The original grant period was for three years. This was extended at no cost for an extra year to allow the principal investigator time to publish the results. We would like to report that our research work and results (supported by this grant) have fulfilled both objectives of the grant. The following is a summary of the important accomplishments: (1) We have now demonstrated that it is possible to perform accurate numerical simulations of the jet screech phenomenon. Both the axisymmetric case and the fully three-dimensional case were carried out successfully. It is worthwhile to note that this is the first time the screech tone phenomenon has been successfully simulated numerically; (2) All four screech modes were reproduced in the simulation. The computed screech frequencies and intensities were in good agreement with the NASA Langley Research Center data; (3) The staging phenomenon was reproduced in the simulation; (4) The effects of nozzle lip thickness and jet temperature were studied. Simulated tone frequencies at various nozzle lip thickness and jet temperature were found to agree well with experiments; (5) The simulated data were used to explain, for the first time, why there are two axisymmetric screech modes and two helical/flapping screech modes; (6) The simulated data were used to show that when two tones are observed, they co-exist rather than switching from one mode to the other, back and forth, as some previous investigators have suggested; and (7) Some resources of the grant were used to support the development of new computational aeroacoustics (CAA) methodology. (Our screech tone simulations have benefited because of the availability of these improved methods.)

Tam, Christopher K. W.↗