Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 37 records · Page 2

Direct simulations of chemically reacting turbulent mixing layers, part 2

The results of direct numerical simulations of chemically reacting turbulent mixing layers are presented. This is an extension of earlier work to a more detailed study of previous three dimensional simulations of cold reacting flows plus the development, validation, and use of codes to simulate chemically reacting shear layers with heat release. Additional analysis of earlier simulations showed good agreement with self similarity theory and laboratory data. Simulations with a two dimensional code including the effects of heat release showed that the rate of chemical product formation, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release. Subsequent three dimensional simulations showed similar behavior, in agreement with laboratory observations. Baroclinic torques and thermal expansion in the mixing layer were found to produce changes in the flame vortex structure that act to diffuse the pairing vortices, resulting in a net reduction in vorticity. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers were shown to result from vorticity generation by baroclinic torques.

Metcalfe, Ralph W.↗

Double Compton scatter telescope calibration

Calibration techniques for a medium energy gamma ray telescope are described. Gain calibration using Compton edge spectra involves comparisons of pulse height spectra with spectra simulated by a Monte Carlo computer code which includes Compton scattering and pair production, plural scattering and variable energy resolution, and cell size. The telescope considered comprises 56 cells of liquid scintillator in four size groups, with a total liquid volume of 325 l; each cell has its own photomultiplier tube. Energy and angular resolutions and the PMT gain calibration procedure are verified with double scatter data for monoenergetic gamma rays at a known location. Detection probabilities for any cell combination in the two telescope arrays are calculated per steradian as a function of the scattering for a number of different energies with a Van de Graaff accelerator.

Dayton, B.↗

Cascade model of gamma-ray bursts: Power-law and annihilation-line components

If, in a neutron star magnetosphere, an electron is accelerated to an energy of 10 to the 11th or 12th power eV by an electric field parallel to the magnetic field, motion of the electron along the curved field line leads to a cascade of gamma rays and electron-positron pairs. This process is believed to occur in radio pulsars and gamma ray burst sources. Results are presented from numerical simulations of the radiation and photon annihilation pair production processes, using a computer code previously developed for the study of radio pulsars. A range of values of initial energy of a primary electron was considered along with initial injection position, and magnetic dipole moment of the neutron star. The resulting spectra was found to exhibit complex forms that are typically power law over a substantial range of photon energy, and typically include a dip in the spectrum near the electron gyro-frequency at the injection point. The results of a number of models are compared with data for the 5 Mar., 1979 gamma ray burst. A good fit was found to the gamma ray part of the spectrum, including the equivalent width of the annihilation line.

Harding, A. K.↗

High speed commercial transport fuels considerations and research needs

NASA is currently evaluating the potential of incorporating High Speed Civil Transport (HSCT) aircraft in the commercial fleet in the beginning of the 21st century. NASA sponsored HSCT enabling studies currently underway with airframers and engine manufacturers, are addressing a broad range of technical, environmental, economic, and related issues. Supersonic cruise speeds for these aircraft were originally focused in the Mach 2 to 5 range. At these flight speeds, both jet fuels and liquid methane were considered potential fuel candidates. For the year 2000 to 2010, cruise Mach numbers of 2 to 3+ are projected for aircraft fuel with thermally stable liquid jet fuels. For 2015 and beyond, liquid methane fueled aircraft cruising at Mach numbers of 4+ may be viable candidates. Operation at supersonic speeds will be much more severe than those encountered at subsonic flight. One of the most critical problems is the potential deterioration of the fuel due to the high temperature environment. HSCT fuels will not only be required to provide the energy necessary for flight, but will also be subject to aerodynamic heating and, will be required to serve as the primary heat sink for cooling the engine and airframe. To define fuel problems for high speed flight, a fuels workshop was conducted at NASA Lewis Research Center. The purpose of the workshop was to gather experts on aviation fuels, airframe fuel systems, airport infrastructure, and combustion systems to discuss high speed fuel alternatives, fuel supply scenarios, increased thermal stability approaches and measurements, safety considerations, and to provide directional guidance for future R and D efforts. Subsequent follow-up studies defined airport infrastructure impacts of high speed fuel candidates. The results of these activities are summarized. In addition, an initial case study using modified in-house refinery simulation model Gordian code (1) is briefly discussed. This code can be used to simulate different types of refineries, emphasizing jet fuel production and relative cost factors.

Lee, C. M.↗

Cascade model of gamma-ray bursts

If, in a neutron star magnetosphere, an electron is accelerated to an energy of 10 to the 11th or 12th power eV by an electric field parallel to the magnetic field, motion of the electron along the curved field line leads to a cascade of gamma rays and electron-positron pairs. This process is believed to occur in radio pulsars and gamma ray burst sources. Results are presented from numerical simulations of the radiation and photon annihilation pair production processes, using a computer code previously developed for the study of radio pulsars. A range of values of initial energy of a primary electron was considered along with initial injection position, and magnetic dipole moment of the neutron star. The resulting spectra was found to exhibit complex forms that are typically power law over a substantial range of photon energy, and typically include a dip in the spectrum near the electron gyro-frequency at the injection point. The results of a number of models are compared with data for the 5 Mar., 1979 gamma ray burst. A good fit was found to the gamma ray part of the spectrum, including the equivalent width of the annihilation line.

Sturrock, P. A.↗

Investigation of three-dimensional separation at wing/body junctions in supersonic flows using TVD MacCormack's scheme

A Total Variation Diminishing (TVD) MacCormack scheme is used to study the three-dimensional separation at wing/body junctions in supersonic flows. The new scheme is robust, resolves discontinuities well with-out any numerical oscillations present in the flow field. Moreover, this new scheme does not contain any problem dependent parameters to be adjusted and many production codes employing MacCormack algorithm can be easily updated to good effect. Numerical simulation carried out in laminar supersonic junction flows using the new scheme yields improved prediction for the separation location and the axial velocity profiles in the separated flow region.

Lakshmanan, B.↗

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗

Comparative study of MacCormack and TVD MacCormack schemes for three-dimensional separation at wing/body junctions in supersonic flows

A robust, discontinuity-resolving TVD MacCormack scheme containing no dependent parameters requiring adjustment is presently used to investigate the 3D separation of wing/body junction flows at supersonic speeds. Many production codes employing MacCormack schemes can be adapted to use this method. A numerical simulation of laminar supersonic junction flow is found to yield improved separation location predictions, as well as the axial velocity profiles in the separated flow region.

Lakshmanan, Balakrishnan↗

Ares I Flight Control System Overview

This paper describes the control challenges posed by the Ares I vehicle, the flight control system design and performance analyses used to test and verify the design. The major challenges in developing the control system are structural dynamics, dynamic effects from the powerful first stage booster, aerodynamics, first stage separation and large uncertainties in the dynamic models for all these. Classical control techniques were employed using innovative methods for structural mode filter design and an anti-drift feature to compensate for translational and rotational disturbances. This design was coded into an integrated vehicle flight simulation and tested by Monte Carlo methods. The product of this effort is a linear, robust controller design that is easy to implement, verify and test.

Hall, Charles↗

A Preliminary Study on the Toxic Combustion Products Testing of Polymers Used in High-Pressure Oxygen Systems

One likely cause of polymer ignition in a high-pressure oxygen system is adiabatic-compression heating of polymers caused by pneumatic impact. Oxidative _ pyrolysis or combustion of polymers in a high-pressure oxygen system could generate toxic gases. This paper reports the preliminary results of toxic combustion product testing of selected polymers in a pneumatic-impact test system. Five polymers commonly used in high-pressure oxygen systems, Nylon 6/6, polychlorotrifluoroethylene (CTFE), polytetrafluoroethylene (PTFE), fluoroelastomer (Viton(TradeMark) A), and nitrile rubber (Buna N), were tested in a pneumatic-impact test system at 2500- or 3500-psia oxygen pressure. The polymers were ignited and burned, then combustion products were collected in a stainless-steel sample bottle and analyzed by GC/MS/IRD, GC/FID, and GC/Methanizer/FID. The results of adiabatic-compression tests show that combustion of hydrocarbon polymers, nitrogen-containing polymers, and halogenated polymers in high-pressure oxygen systems are relatively complete. Toxicity of the combustion product gas is presumably much lower than the combustion product gas generated from ambient-pressure oxygen (or air) environments. The NASA-Lewis equilibrium code was used to determine the composition of combustion product gas generated from a simulated, adiabatic-compression test of nine polymers. The results are presented and discussed.

Hshieh, Fu-Yu↗

Comparison between Hayabusa 2 Spectral Measurements and Simulations

This paper compares recent Hayabusa 2 spectral measurements with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the coupled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator, based on the limited available published information. The comparison between the simulations and measurements focuses on two atomic nitrogen lines and two atomic oxygen lines, as well as the CN Violet band system. For the atomic lines, the measurements and simulations agree within 25% over most of the trajectory. This excellent agreement is unprecedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield/radiation modeling. To provide a link between these observed radiation measurements and the radiative heating to a vehicle surface, the recently developed flowfield property binning approach is used to identify flowfield properties that provide the dominant emission contribution to the measured spectrum. These identified flowfield properties are then shown to match those for the surface radiative heating to the currently developed Mars Sample Return (MSR) Earth Entry System (EES) at a specific trajectory point and surface location. This matching of flowfield properties indicates that this EES radiative heating and Hayabusa 2 observed radiation simulations are equivalent radiation problems. Therefore, the excellent agreement between Hayabusa 2 measurements and simulations may be leveraged to inform the radiation heating margin for the MSR EES.

Christopher O Johnston↗

Comparison between Hayabusa 2 Spectral Measurements and Simulations

Recent Hayabusa 2 spectral measurements are compared with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the cou- pled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator based on available published information. The comparison between the simulations and measurements focuses on three atomic nitrogen lines, two atomic oxygen lines, one atomic hydrogen line, and the CN Violet band system. Parametric uncertainties are evaluated for these features using a recently developed capability in LAURA/HARA. These uncertainties capture the impact of flowfield and radiation modeling uncertainties on the observed radiation simulations, and include the impact of coupled radiation and ablation product emission. The resulting parametric uncertainties range from 50 to over 100% for the atomic lines, and 80% for CN Violet. Uncertainties in flowfield kinetics and non-Boltzmann modeling provide the dominant contribution to these uncertainties. For the atomic nitrogen and oxygen lines, com- parisons between the nominal simulations and measurements show agreement well within these simulation uncertainty bounds, with the measurements within 20% of the nominal simulation over most of the trajectory, and within 10% at peak emission. This excellent agreement is un- precedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield and radiation modeling. For the ablation product emission from the 656 nm atomic hydrogen line, the agreement between measurements and simulations is within ±20% for all trajectory points except one outlier. This good comparison is surprising considering this species depends on the ablation rate and is shown to have strong wake contribution, which both contribute to the large ≈±70% parametric uncertainty. For the ablation product emission from CN Violet, the simulations over-predict the measurements by over a factor of two early in the trajectory, with the comparison improving later in the trajectory. This disagreement for CN Violet is consistent with the ≈±80% parametric uncertainty bounds evaluated for the simulations. To utilize the excellent agreement between the atomic line simulations and measurements in the assessment of the radiative heating margin for a future flight vehicle, the flowfield property binning (FPB) approach is applied. Using the Mars Sample Return (MSR) Earth Entry System (EES) as an example, the combination of the excellent Hayabusa comparisons and the FPB analysis justify nearly a 10% reduction in the radiative heating margin from an assumed baseline margin of 30%. This direct quantitative link between observed radiation measurements and the radia- tive heating margin of future flight vehicles demonstrates the value of these observed radiation measurements.

Shock-Layer Radiative Heating↗

Computer simulation of auroral kilometric radiation

An investigation is carried out of the linear amplification and nonlinear saturation of electromagnetic waves associated with two components of electrons consisting of a cold Maxwellian background and a weakly relativistic population of electrons possessing a loss-cone distribution. The goal is to understand the mechanism responsible for the production of auroral kilometric radiation (AKR). A relativistic 1-2/2 dimensional electromagnetic particle simulation code is used in studying an initial value problem modeling the AKR source region. The simulation makes it possible to follow the growth of all modes of radiation, including the ordinary mode (O-mode), slow extraordinary mode (Z-mode), whistler mode and fast extraordinary mode (X-mode), past the point of saturation.

Wagner, J. S.↗

Using the Diurnal Variability in GeoNEX TOA Reflectances for Earth Monitoring

Observations from the third-generation geostationary satellite instruments (GOES 16/17 ABI, Himawari 8/9 AHI, and etc.) have spatial resolution and spectral band configurations comparable to flagship LEO sensors (e.g., MODIS/VIIRS). More importantly, these data are acquired at very high temporal resolution, faithfully recording the variations of the full disk of Earth at every 5-10 minutes. They thus provide unique information about Earth’s atmosphere and surface. In order to explore the unique information content of geostationary data, this study systematically analyzes the diurnal variability in the GeoNEX L1G TOA reflectance products and compares them to simulated results by state-of-the-art radiative transfer codes. Our results show that • The smoothness of the TOA reflectance diurnal cycle provides a convenient and reliable way to identify stable atmospheric conditions and filter out passing clouds/shadows. • The diurnal variability of the blue band (0.47µm) reflectance is regulated mainly by atmospheric optical conditions over a majority of land cover types. As such, the diurnal variability of the blue band data allows us to retrieve AOD without invoking the use of spectral band ratios (SRC) as in previous algorithms. • In comparison, the diurnal variability of the short-wave infrared band (2.2µm) BRFs is mainly regulated by surface reflectance and the sun-target-satellite geometry. This information allows us to test and, if suitable, retrieve surface BRDF parameters. • Spectral band ratios, especially those between the 2.2µm and 0.47µm bands, are not constant but vary by locations and sun-target-satellite geometries. Our analysis clearly demonstrates that the information provided in high-frequent geostationary observations is unique and complementary to LEO sensors. Therefore, a synergy of GEO and LEO (and other) sensors has the great potential to improve existing remote sensing models and algorithms for better Earth monitoring.

Diurnal Variability↗

Physical Processes for Driving Ionospheric Outflows in Global Simulations

We review and assess the importance of processes thought to drive ionospheric outflows, linking them as appropriate to the solar wind and interplanetary magnetic field, and to the spatial and temporal distribution of their magnetospheric internal responses. These begin with the diffuse effects of photoionization and thermal equilibrium of the ionospheric topside, enhancing Jeans' escape, with ambipolar diffusion and acceleration. Auroral outflows begin with dayside reconnexion and resultant field-aligned currents and driven convection. These produce plasmaspheric plumes, collisional heating and wave-particle interactions, centrifugal acceleration, and auroral acceleration by parallel electric fields, including enhanced ambipolar fields from electron heating by precipitating particles. Observations and simulations show that solar wind energy dissipation into the atmosphere is concentrated by the geomagnetic field into auroral regions with an amplification factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Internal plasmas thus enable electromagnetic driving via coupling to the plasma, neutral gas and by extension, the entire body " We assess the Importance of each of these processes in terms of local escape flux production as well as global outflow, and suggest methods for their implementation within multispecies global simulation codes. We complete 'he survey with an assessment of outstanding obstacles to this objective.

Moore, Thomas Earle↗

Mechanisms of Ionospheric Mass Ejection

Ionospheric outflows are directly responsive to solar wind disturbances, particularly in the dayside auroral cusp or cleft regions. Inputs of both electromagnetic energy (Poynting flux) and kinetic energy (particle precipitation) are closely correlated with these outflows. We assess the importance of processes thought to drive ionospheric outflows. These begin with the diffuse effects of photoionization and thermal equilibrium of the ionospheric topside, enhancing Jeans' escape, with ambipolar diffusion and acceleration. Auroral outflows begin with dayside reconnexion and resultant field-aligned currents and driven convection. These produce plasmaspheric plumes, collisional heating and wave-particle interactions, centrifugal acceleration, and auroral acceleration by parallel electric fields, including enhanced ambipolar fields from electron heating by precipitation particles. Solar wind energy dissipation is concentrated by the geomagnetic field into auroral regions with an amplification factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Internal plasmas thus enable electromagnetic driving via coupling to the plasma and neutral gas. We assess the importance of each of these processes in terms of local escape flux production as well as global outflow, and suggest methods for their implementation within multi-species global simulation codes. We conclude by assessing outstanding obstacles to this objective.

Moore, Thomas Earle↗

Modelling the Water and Carbon Dioxide Production Rates of Comet 67P/Churyumov-Gerasimenko

The European Space Agency Rosetta/Philae mission to Comet 67P/Churyumov–Gerasimenko in 2014–2016 is the most complete and diverse investigation of a comet carried out thus far. Yet, many physical and chemical properties of the comet remain uncertain or unknown, and cometary activity is still not a well–understood phenomenon. We here attempt to place constraints on the nucleus abundances and sublimation front depths of H 2 O and CO 2 ice, and to reconstruct how the nucleus evolved throughout the perihelion passage. We employ the thermophysical modelling code ‘Numerical Icy Minor Body evolUtion Simulator’, or NIMBUS, to search for conditions under which the observed H 2 O and CO 2 production rates are simultaneously reproduced before and after perihelion. We find that the refractories to water–ice mass ratio of relatively pristine nucleus material is μ ≈ 1, that airfall material has μ ≈ 2, and that the molar abundance of CO 2 relative H 2 O is near 30 per cent. The dust mantle thickness is typically ≲ 2 cm. The average CO 2 sublimation front depths near aphelion were ~ 3.8 m and ~ 1.9 m on the northern and southern hemispheres, respectively, but varied substantially with time. We propose that airfall material is subjected to substantial fragmentation and pulverisation due to thermal fatigue during the aphelion passage. Sub–surface compaction of material due to CO 2 activity near perihelion seems to have reduced the diffusivity in a measurable way.

67P/Churyumov–Gerasimenko↗

Comparison of DAC and MONACO DSMC Codes with Flat Plate Simulation

Various implementations of the direct simulation Monte Carlo (DSMC) method exist in academia, government and industry. By comparing implementations, deficiencies and merits of each can be discovered. This document reports comparisons between DSMC Analysis Code (DAC) and MONACO. DAC is NASA's standard DSMC production code and MONACO is a research DSMC code developed in academia. These codes have various differences; in particular, they employ distinct computational grid definitions. In this study, DAC and MONACO are compared by having each simulate a blunted flat plate wind tunnel test, using an identical volume mesh. Simulation expense and DSMC metrics are compared. In addition, flow results are compared with available laboratory data. Overall, this study revealed that both codes, excluding grid adaptation, performed similarly. For parallel processing, DAC was generally more efficient. As expected, code accuracy was mainly dependent on physical models employed.

Padilla, Jose F.↗