SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,477 records · Page 82

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Mixing, Noise and Thrust Benefits Using Corrugated Designs

These projects are directed toward the analysis of several concepts for nozzle and inlet performance improvement and noise reduction from jet exhausts. Currently. The FM&AL also initiates new joint research between the HU/FM&AL, the Hyper-X Program Team at the LaRC, and the Central Institute of Aviation Motors (CIAM), Moscow, Russia in the field of optimization of fuel injection and mixing in air-breathing propulsion systems. The main results of theoretical, numerical simulation and experimental tests obtained in the previous research are in the papers and patents. The goals of the 14U/FM&AL programs are twofold: 1) to improve the working efficiency of the HU/FM&AL team in generating new innovative ideas and in conducting research in the field of fluid dynamics and acoustics, basically for improvement of supersonic and subsonic aircraft engines, and 2) to attract promising minority students to this research and training and, in cooperation with other HU departments, to teach them basic knowledge in Aerodynamics, Gas Dynamics, and Theoretical and Experimental Methods in Aeroacoustics and Computational Fluid Dynamics (CFD). The research at the HU/FM&AL supports reduction schemes associated with the emission of en 'ne pollutants for commercial aircraft and concepts for reduction of 91 observables for military aircraft. These research endeavors relate to the goals of the NASA Strategic Enterprise in Aeronautics concerning the development of environmentally acceptable aircraft. It is in this precise area, where the US aircraft industry, academia, and Government are in great need of trained professionals and which is a high priority goal of the Minority University Research and Education (MLTREP) Program, that the HU/FM&AL can make its most important contribution.

Morgan, Morris H., III↗

NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict three-dimensional ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

NASA NTRS2017

Characterization of a Method for Inverse Heat Conduction Using Real and Simulated Thermocouple Data

It is often impractical to instrument the external surface of high-speed vehicles due to the aerothermodynamic heating. Temperatures can instead be measured internal to the structure using embedded thermocouples, and direct and inverse methods can then be used to estimate temperature and heat flux on the external surface. Two thermocouples embedded at different depths are required to solve direct and inverse problems, and filtering schemes are used to reduce noise in the measured data. Accuracy in the estimated surface temperature and heat flux is dependent on several factors. Factors include the thermocouple location through the thickness of a material, the sensitivity of the surface solution to the error in the specified location of the embedded thermocouples, and the sensitivity to the error in thermocouple data. The effect of these factors on solution accuracy is studied using the methodology discussed in the work of Pizzo, et. al.1 A numerical study is performed to determine if there is an optimal depth at which to embed one thermocouple through the thickness of a material assuming that a second thermocouple is installed on the back face. Solution accuracy will be discussed for a range of embedded thermocouple depths. Moreover, the sensitivity of the surface solution to (a) the error in the specified location of the embedded thermocouple and to (b) the error in the thermocouple data are quantified using numerical simulation, and the results are discussed.

Pizzo, Michelle E.↗

NASA NTRS1988

Numerical Simulation of Subsonic and Transonic Propeller Flow

The numerical simulation of 3-D transonic flow about a system of propeller blades is investigated. In particular, it is shown that the use of helical coordinates significantly simplifies the form of the governing equation when the propeller system is assumed to be surrounded by an irrotational flow field of an inviscid fluid. The unsteady small disturbance equation, valid for lightly loaded blades and expressed in helical coordinates, is derived from the general blade-fixed potential equation, given for an arbitrary coordinate system. The use of a coordinate system which inherently adapts to the mean flow results in a disturbance equation requiring relatively few terms to accurately model the physics of the flow. Furthermore, the helical coordinate system presented here is novel in that it is periodic in the circumferential direction while, simultaneously, maintaining orthogonal properties at the mean blade locations. The periodic characteristic allows a complete cascade of blades to be treated, and the orthogonality property affords straightforward treatment of blade boundary conditions. An ADI numerical scheme is used to compute the solution of the steady flow as an asymptotic limit of an unsteady flow. As an example of the method, solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an 8-bladed cascade. Both high and low advance ratio cases are computed and include a lifting as well as nonlifting cases. The nonlifting solutions obtained are compared to solutions from a Euler code.

Snyder, Aaron↗

NASA NTRS2011

A Deterministic Electron, Photon, Proton and Heavy Ion Radiation Transport Suite for the Study of the Jovian System

A deterministic suite of radiation transport codes, developed at NASA Langley Research Center (LaRC), which describe the transport of electrons, photons, protons, and heavy ions in condensed media is used to simulate exposures from spectral distributions typical of electrons, protons and carbon-oxygen-sulfur (C-O-S) trapped heavy ions in the Jovian radiation environment. The particle transport suite consists of a coupled electron and photon deterministic transport algorithm (CEPTRN) and a coupled light particle and heavy ion deterministic transport algorithm (HZETRN). The primary purpose for the development of the transport suite is to provide a means for the spacecraft design community to rapidly perform numerous repetitive calculations essential for electron, proton and heavy ion radiation exposure assessments in complex space structures. In this paper, the radiation environment of the Galilean satellite Europa is used as a representative boundary condition to show the capabilities of the transport suite. While the transport suite can directly access the output electron spectra of the Jovian environment as generated by the Jet Propulsion Laboratory (JPL) Galileo Interim Radiation Electron (GIRE) model of 2003; for the sake of relevance to the upcoming Europa Jupiter System Mission (EJSM), the 105 days at Europa mission fluence energy spectra provided by JPL is used to produce the corresponding dose-depth curve in silicon behind an aluminum shield of 100 mils ( 0.7 g/sq cm). The transport suite can also accept ray-traced thickness files from a computer-aided design (CAD) package and calculate the total ionizing dose (TID) at a specific target point. In that regard, using a low-fidelity CAD model of the Galileo probe, the transport suite was verified by comparing with Monte Carlo (MC) simulations for orbits JOI--J35 of the Galileo extended mission (1996-2001). For the upcoming EJSM mission with a potential launch date of 2020, the transport suite is used to compute the traditional aluminum-silicon dose-depth calculation as a standard shield-target combination output, as well as the shielding response of high charge (Z) shields such as tantalum (Ta). Finally, a shield optimization algorithm is used to guide the instrument designer with the choice of graded-Z shield analysis.

Norman, Ryan B.↗

NASA NTRS1988

Liquid rocket boosters for the Space Shuttle - An outline of trade issues and analysis techniques

Numerous liquid rocket booster configurations were studied to replace the Solid Rocket Boosters on the current Space Shuttle. While improved safety, reliability, performance and economy of operation were sought, the designs also sought to minimize volumetric growth to reduce aerodynamic loads and structural modifications. Aero loads could be significant as a function of booster diameter if extrapolated directly from solid rocket characteristics. Procedures to meet these goals with analytical booster sizing and simple dynamic simulations are shown. Engine chamber pressure proved to be a good independent variable for implications of pressure-fed and pump-fed engines, fuel and material selection. For some configurations, ascent trajectory modifications were suggested: longer vertical rise, lofting, heads-up orientation, and higher first and second stage target altitudes. Since engine manufacture would dominate booster costs, booster simplification of reuse are discussed based on cost estimates.

Kelly, W. D.↗

NASA NTRS1995

The 3D widgets for exploratory scientific visualization

Computational fluid dynamics (CFD) techniques are used to simulate flows of fluids like air or water around such objects as airplanes and automobiles. These techniques usually generate very large amounts of numerical data which are difficult to understand without using graphical scientific visualization techniques. There are a number of commercial scientific visualization applications available today which allow scientists to control visualization tools via textual and/or 2D user interfaces. However, these user interfaces are often difficult to use. We believe that 3D direct-manipulation techniques for interactively controlling visualization tools will provide opportunities for powerful and useful interfaces with which scientists can more effectively explore their datasets. A few systems have been developed which use these techniques. In this paper, we will present a variety of 3D interaction techniques for manipulating parameters of visualization tools used to explore CFD datasets, and discuss in detail various techniques for positioning tools in a 3D scene.

Herndon, Kenneth P.↗

NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗

NASA NTRS2023

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Chris Ostrom↗

NASA NTRS1988

Monte Carlo studies of liquid semiconductor surfaces - Si and Ge

The liquid-vapor interface of Si and Ge, and of Si doped with impurities is studied using the empirical Stillinger-Weber two- and three-body potentials. The surface tension of the pure elements is calculated by use of a direct evaluation of the free energy required to create the surface. For Si, both the surface tension and its temperature derivative are in good agreement with experiment. To within numerical accuracy, the free surfaces of both Si and Ge have a monotonically decreasing density, with a 10-90 percent surface width of about 2.2 A in both cases. When large and small model impurities are introduced into pure liquid Si, they are found in the simulations to migrate, respectively, towards the surface and away from the surface. This behavior is consistent with the interpretation that impurities with atoms larger than Si tend to lower the surface tension.

Wang, Z. Q.↗

NASA NTRS2011

Atomic Data for Neutron-capture Elements I. Photoionization and Recombination Properties of Low-charge Selenium Ions

We present multi-configuration Breit-Pauli AUTOSTRUCTURE calculations of distorted-wave photoionization (PI) cross sections. and total and partial final-state resolved radiative recombination (RR) and dielectronic recombination (DR) rate coefficients for the first six ions of the trans-iron element Se. These calculations were motivated by the recent detection of Se emission lines in a large number of planetary nebulae. Se is a potentially useful tracer of neutron-capture nucleosynthesis. but accurate determinations of its abundance in photoionized nebulae have been hindered by the lack of atomic data governing its ionization balance. Our calculations were carried out in intermediate coupling with semi re1ativistic radial wavefunctions. PI and recombination data were determined for levels within the ground configuration of each ion, and experimental PI cross-section measurements were used to benchmark our results. For DR, we allowed (Delta)n = 0 core excitations, which are important at photoionized plasma temperatures. We find that DR is the dominant recombination process for each of these Se ions at temperatures representative of photoionized nebulae (approx.10(exp 4) K). In order to estimate the uncertainties of these data, we compared results from three different configuration-interaction expansions for each ion, and also tested the sensitivity of the results to the radial scaling factors in the structure calculations. We find that the internal uncertainties are typically 30-50% for the direct PI cross sections and approx.10% for the computed RR rate coefficients, while those for low-temperature DR can be considerably larger (from 15-30% up to two orders of magnitude) due to the unknown energies of near-threshold autoionization resonances. These data are available at the CDS, and fitting coefficients to the total RR and DR rate coefficients are presented. The results are suitable for incorporation into photoionization codes used to numerically simulate astrophysical nebulae, and will enable robust determinations of nebular Se abundances.

Sterling, N. C.↗

NASA NTRS1981

Daylight imagery of ocean surface waves for wave spectra

A simple second-order model is developed to describe the transformation between surface wave slope and radiance incident on an optical sensor. The model includes azimuth as well as elevation variations in the sky-radiance distribution and their effects on the imaging of wave slopes approaching the sensor azimuth from various directions. The model also includes upwelling radiance. Comparison between the exact transformation from wave slope to radiance and the second-order model suggests that they yield quantitatively similar results. The model is used in physical interpretations of how wave slopes are imaged and in indicating a way of optimizing the linearity and contrast of the transfer function relating wave slope to radiance for various sensor geometries, thereby minimizing sources of error. One-dimensional numerical simulations of the imaging of an analytically generated wave-slope profile demonstrate the utility of this technique in obtaining slope spectra.

Monaldo, F. M.↗

NASA NTRS1987

Numerical study of finite-rate supersonic combustion using parabolized equations

A set of partial differential equations, describing the two-dimensional supersonic chemically-reacting flow of the hydrogen-air system, is formulated such that the equations are parabolic in the streamwise direction. A fully-implicit fully-coupled finite-difference algorithm is used to develop a computer code which solves the governing equations by marching in the streamwise direction. The combustion process is modeled by a two-step finite-rate chemistry whereas turbulence is simulated by an algebraic turbulence model. Results of two calculations of internal supersonic reacting flow show fairly good agreement with the results obtained by the more costly full Navier-Stokes procedure.

Chitsomboon, T.↗

NASA NTRS2009

Application of TRMM PR and TMI Measurements to Assess Cloud Microphysical Schemes in the MM5 Model for a Winter Storm

Seen by the human eye, precipitation particles are commonly drops of rain, flakes of snow, or lumps of hail that reach the ground. Remote sensors and numerical models usually deal with information about large collections of rain, snow, and hail (or graupel --also called soft hail ) in a volume of air. Therefore, the size and number of the precipitation particles and how particles interact, evolve, and fall within the volume of air need to be represented using physical laws and mathematical tools, which are often implemented as cloud and precipitation microphysical parameterizations in numerical models. To account for the complexity of the precipitation physical processes, scientists have developed various types of such schemes in models. The accuracy of numerical weather forecasting may vary dramatically when different types of these schemes are employed. Therefore, systematic evaluations of cloud and precipitation schemes are of great importance for improvement of weather forecasts. This study is one such endeavor; it pursues quantitative assessment of all the available cloud and precipitation microphysical schemes in a weather model (MM5) through comparison with the observations obtained by National Aeronautics and Space Administration (NASA) s and Japan Aerospace Exploration Agency (JAXA) s Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) and microwave imager (TMI). When satellite sensors (like PR or TMI) detect information from precipitation particles, they cannot directly observe the microphysical quantities (e.g., water species phase, density, size, and amount etc.). Instead, they tell how much radiation is absorbed by rain, reflected away from the sensor by snow or graupel, or reflected back to the satellite. On the other hand, the microphysical quantities in the model are usually well represented in microphysical schemes and can be converted to radiative properties that can be directly compared to the corresponding PR and TMI observations. This study employs this method to evaluate the accuracy of the simulated radiative properties by the MM5 model with different microphysical schemes. It is found that the representations of particle density, size, and mass in the different schemes in the MM5 model determine the model s performance when predicting a winter storm over the eastern Pacific Ocean. Schemes lacking moderate density particles (i.e. graupel), with snow flakes that are too large, or with excessive mass of snow or graupel lead to degraded prediction of the radiative properties as observed by the TRMM satellite. This study demonstrates the uniqueness of the combination of both an active microwave sensor (PR) and passive microwave sensor (TMI) onboard TRMM on assessing the accuracy of numerical weather forecasting. It improves our understanding of the physical and radiative properties of different types of precipitation particles and provides suggestions for better representation of cloud and precipitation processes in numerical models. It would, ultimately, contribute to answering questions like "Why did it not rain when the forecast says it would?"

Han, Mei↗

NASA NTRS1979

Application of a microprocessor controlled lidar to tropospheric ozone measurements

A microprocessor controlled lidar system under construction at the NASA Goddard Space Flight Center is described and the problems in making space-based measurements of tropospheric ozone are considered. The differential absorption lidar using a dual wavelength, pulsed CO2 laser and direct detection receiver can significantly improve the existing global data base on tropospheric ozone burden. Sensitivity to tropospheric ozone can be obtained in the spaceborne version of lidar by selecting laser lines located in the wings of the target zone lines. Simulation studies using various laser line pairs in the P-branch of the CO2 9.4 micron band show that the ozone burden retrieval may be weighted to particular altitude regions. These simulation studies are based on numerical integration of differences in absorption coefficient at the two selected laser wavelengths using the AFGL absorption line parameter compilations, U.S. Standard Atmosphere ozone profile, and laser software. Simulation, data collection and reduction are performed by a microprocessor subsystem of the CO2 lidar.

Stewart, R. W.↗

NASA NTRS1984

Effect of a buried-wire gage on the separation bubble Numerical study

A numerically simulated buried-wire separation gage is investigated with emphasis on its effect on the separation bubble. The conjugated problem of a supersonic, time-dependent, two-dimensional flowfield above a conductive solid wall with an embedded heat source is solved using implicit finite difference algorithms. Steady-state and transient cases were computed for different locations of the heat source within the bubble. Results show that by using a steady heat source, the flow direction near the wall can be detected, without distorting the flowfield, only if the source is located in regions where the bubble is thick (i.e., not too close to the separation). The flow direction near separation can be detected by using a temperature pulse at the solid/fluid interface with insignificant distortion of the flowfield.

Degani, D.↗

NASA NTRS1993

Damage detection using experimentally measured mass and stiffness matrices

A method is presented for locating physical damage or change in a structure using experimentally measured mass and stiffness matrices. The approach uses a recently developed algorithm for transforming a state-space realization into a second order structural model with physical displacements as the generalized coordinates. This is accomplished by first rotating a state-space model of the identified structural dynamics into modal coordinates and approximating the mass normalized modal vectors for the output measurement set. Next, the physical mass, damping and stiffness matrices are synthesized directly from the measured modal parameters. This yields experimental mass and stiffness matrices for the structure without the use of a finite element model or a numerical search. The computed mass and stiffness are asymptotically equivalent to a static condensation of the global physical coordinate model. Techniques for solving the inverse connectivity problem are then developed whereby it is possible to assess the stiffness in a region of the structure bounded by several sensors. Applications to both simulated data and experimental data are used to discuss the effectiveness of the approach.

Peterson, L. D.↗

NASA NTRS1993

The cosmic Mach number - Direct comparisons of observations and models

A robust method is used to derive the cosmic Mach number, M, from existing peculiar velocity data sets. For the IR Tully-Fisher distances of spirals in the Local Supercluster of Aaronson et al. (1979, 1982), M of 1.03 and a characteristic distance from the Local Group of 1464 km/s are found. For the elliptical galaxy sample of Faber et al. (1987, 1988), M of 0.57, with a characteristic distance of 2572 km/s is found. These results are compared with Monte Carlo simulations of the observational realizations drawn from numerical simulations of the universe based on various scenarios, including the standard cold dark matter (CDM) scenario. The effect of velocity bias on the derived Mach number is found to be small. Only 5 percent of the CDM simulations have M as large as or larger than that derived for the Aaronson et al. sample.

Strauss, Michael A.↗