SEARCH · Engineering Papers

Results for “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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Numerical Simulations of Single Flow Element in a Nuclear Thermal Thrust Chamber

The objective of this effort is to develop an efficient and accurate computational methodology to predict both detailed and global thermo-fluid environments of a single now element in a hypothetical solid-core nuclear thermal thrust chamber assembly, Several numerical and multi-physics thermo-fluid models, such as chemical reactions, turbulence, conjugate heat transfer, porosity, and power generation, were incorporated into an unstructured-grid, pressure-based computational fluid dynamics solver. The numerical simulations of a single now element provide a detailed thermo-fluid environment for thermal stress estimation and insight for possible occurrence of mid-section corrosion. In addition, detailed conjugate heat transfer simulations were employed to develop the porosity models for efficient pressure drop and thermal load calculations.

Cheng, Gary↗

NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

DOE OSTI2022

Supersaturation Variability from Scalar Mixing: Evaluation of a New Subgrid-Scale Model Using Direct Numerical Simulations of Turbulent Rayleigh–Bénard Convection

We report that supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars - water vapor and temperature - leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1997

Rocket Engine Numerical Simulator (RENS)

Work is being done at three universities to help today's NASA engineers use the knowledge and experience of their Apolloera predecessors in designing liquid rocket engines. Ground-breaking work is being done in important subject areas to create a prototype of the most important functions for the Rocket Engine Numerical Simulator (RENS). The goal of RENS is to develop an interactive, realtime application that engineers can utilize for comprehensive preliminary propulsion system design functions. RENS will employ computer science and artificial intelligence research in knowledge acquisition, computer code parallelization and objectification, expert system architecture design, and object-oriented programming. In 1995, a 3year grant from the NASA Lewis Research Center was awarded to Dr. Douglas Moreman and Dr. John Dyer of Southern University at Baton Rouge, Louisiana, to begin acquiring knowledge in liquid rocket propulsion systems. Resources of the University of West Florida in Pensacola were enlisted to begin the process of enlisting knowledge from senior NASA engineers who are recognized experts in liquid rocket engine propulsion systems. Dr. John Coffey of the University of West Florida is utilizing his expertise in interviewing and concept mapping techniques to encode, classify, and integrate information obtained through personal interviews. The expertise extracted from the NASA engineers has been put into concept maps with supporting textual, audio, graphic, and video material. A fundamental concept map was delivered by the end of the first year of work and the development of maps containing increasing amounts of information is continuing. Find out more information about this work at the Southern University/University of West Florida. In 1996, the Southern University/University of West Florida team conducted a 4day group interview with a panel of five experts to discuss failures of the RL10 rocket engine in conjunction with the Centaur launch vehicle. The discussion was recorded on video and audio tape. Transcriptions of the entire proceedings and an abbreviated video presentation of the discussion highlights are under development. Also in 1996, two additional 3year grants were awarded to conduct parallel efforts that would complement the work being done by Southern University and the University of West Florida. Dr. Prem Bhalla of Jackson State University in Jackson, Mississippi, is developing the architectural framework for RENS. By employing the Rose Rational language and Booch Object Oriented Programming (OOP) technology, Dr. Bhalla is developing the basic structure of RENS by identifying and encoding propulsion system components, their individual characteristics, and cross-functionality and dependencies. Dr. Ruknet Cezzar of Hampton University, located in Hampton, Virginia, began working on the parallelization and objectification of rocket engine analysis and design codes. Dr. Cezzar will use the Turbo C++ OOP language to translate important liquid rocket engine computer codes from FORTRAN and permit their inclusion into the RENS framework being developed at Jackson State University. The Southern University/University of West Florida grant was extended by 1 year to coordinate the conclusion of all three efforts in 1999.

Davidian, Kenneth O.↗

NASA NTRS2018

Numerical Simulation of AFP Nip Point Temperature Prediction for Complex Geometries

Material placement at the ideal nip point temperature over complex surfaces with uniformity across the width of the compaction rollers results in optimized part properties for Automated Fiber Placement (AFP) processes. However, current AFP systems utilize heat control models and methodologies, based on multiple process parameters such as feed-rate and orientation, that are mostly open-loop. Here, infrared (IR) heater input is calibrated as a function of process parameters during machine qualification. This work presents a numerical simulation to predict arrayed-infrared (AIR) emitter radiation onto a substrate that includes view factor implementation, IR radiative heat flow calculation, energy rate balance, and a transient heat transfer model. The purpose of this numerical model is to predict nip point temperature on complex surfaces, serving as a baseline for a new arrayed-infrared (AIR) thermoset heater to improve AFP process control. It is anticipated that this simulation will accurately control the temperature for high-speed AFP layup of complex geometries. An anticipated result of an AIR heater system is that material calibration and testing will be reduced as temperature is instantaneously monitored and controlled. Therefore, temperature across the roller width will be uniform during placement of complex parts, independent of their geometry.

Xia, Kaishu↗

NASA NTRS1991

Experiment and numerical simulation of RF heating in the Tandem Mirror plasma propulsion device

Recent results of experimental and theoretical studies of ICRF heating of the plasma in the Tendem Mirror rocket are presented. The radial and axial profile of the magnetic field of the wave have been measured, and the data agree with the results from numerical simulation. One very important new finding is that the wave damped when it approached the resonance plane. This is a strong indication of beach heating effect and that RF power is absorbed by the ions in the plasma as expected. This power absorption phenomenon was also further confirmed by the refined analytical study of the wave propagation in a slab model. The electron density in the central cell has been measured by a microwave interferometer. The existence of two types of discharge similar to those of H-alpha emission have also been observed from the density measurements. A preliminary design of a flight system for a 100-day Mars transit has been undertaken. The specific weight for a 10 MW rocket is only 0.04 kg/kW.

Yang, T. F.↗

NASA NTRS1989

A mesoscale numerical simulation of an index for use in the forecasting of the potential severity of wet microbursts

This paper demonstrates the feasibility of using the wet microburst potential intensity index developed by Proctor (1989) for predicting the onset of a wet microburst. After describing the index and the principles involved in its derivation, the version 4.0 of the Mesoscale Atmospheric Simulation System (MASS) model described by Karyampudi et al. (1988) was used to illustrate how the index would perform as a forecast product. An analysis is presented of the use of the microburst potential intensity index in the numerical simulation of the DFR aircraft accident on August 2, 1985 in a downburst environment.

Zack, John W.↗

DOE OSTI2017

Numerical simulation of vortex-induced motion of a deep-draft paired-column semi-submersible offshore platform

Understanding and predicting vortex-induced motion (VIM) of offshore systems for deep seawater applications is crucial to improve the system safety and integrity. In this paper, we report on experimental tow-tank measurements and numerical simulations of VIM of a deep-draft offshore platform, specifically Paired-Column Semisubmersible (PC-Semi). The study is carried out in model scale (1:54), at subcritical flow regime with Re~10 4 . Motion of the floating structure has three degrees of freedom: in-line, cross-flow, and yaw. Large periodic cross-flow motion is measured for headings 0°, 11.25°, and 22.5°, for reduced velocities (U r ) between 5 and 10. Considerably smaller cross-flow amplitude is recorded at 45° heading across the U r range considered. An extensive sensitivity study is performed using computational fluid dynamics (CFD) to capture the transient displacement history of VIM (in-line, cross-flow, and yaw motion components). Amplitude and period of cross-flow (transverse) motion are obtained from statistical analysis of VIM time history and subsequently used as the validation criterion between the CFD simulation and the model tests. Satisfactory agreement between the CFD results and tow-tank measurements is achieved with a Delayed Detached Eddy Simulation – Shear Stress Transport (DDES-SST) formulation. Finally, this work provides experimental results and serves as a practical starting point to set up a CFD problem to estimate amplitude and period of cross-flow VIM motion for offshore engineering applications.

42 ENGINEERING↗

DOE OSTI2022

Numerical simulations of noisy quantum circuits for computational chemistry

The opportunities afforded by near-term quantum computers to calculate the ground-state properties of small molecules depend on the structure of the computational ansatz as well as the errors induced by device noise. Here we investigate the behavior of these noisy quantum circuits using numerical simulations to estimate the accuracy and fidelity of the prepared quantum states relative to the ground truth obtained by conventional means. We implement several different types of ansatz circuits derived from unitary coupled cluster theory for the purposes of estimating the ground-state energy of sodium hydride using the variational quantum eigensolver algorithm. We show how relative error in the energy and the fidelity scale with the levels of gate-based noise, the internuclear configuration, the ansatz circuit depth, and the parameter optimization methods.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Numerical simulations challenged on the prediction of massive subhalo abundance in galaxy clusters: the case of Abell 2142

In this Letter we compare the abundance of the member galaxies of a rich, nearby (z = 0.09) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant (≳7 sigma confidence level) smaller amount of massive (circular velocity above 200 km s{sup −1}) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Clustering instabilities in sedimenting fluid–solid systems: critical assessment of kinetic-theory-based predictions using direct numerical simulation data

In this paper the quantitative and qualitative ability of a kinetic-theory-based two-fluid model (KT-TFM) is assessed in a state of fully periodic sedimentation (fluidization), with a focus on statistically steady, unstable (clustered) states. The accuracy of KT-TFM predictions is evaluated via direct comparison to direct numerical simulation (DNS) data. The KT-TFM and DNS results span a rather wide parameter space: mean-flow Reynolds numbers on the order of 1 and 10, mean solid volume fractions from 0.1 to 0.4, solid-to-fluid density ratios from 10 to 1000 and elastic and moderately inelastic (restitution coefficient of 0.9) conditions. Data from both KT-TFM and DNS display a rich variety of statistically steady yet unstable structures (clusters). Instantaneous snapshots of KT-TFM and DNS demonstrate remarkable qualitative agreement. This qualitative agreement is quantified by calculating the critical density ratio at which the structure transitions from a chaotic, dynamic state to a regular, plug-flow state, with good overall comparisons. Further quantitative assessments of mean and fluctuating velocities show good agreement at high density ratios but weaker agreement at intermediate to low density ratios depending on the mean-flow Reynolds numbers and solid fractions. Deviations of the KT-TFM results from the DNS data were traced to a breakdown in one of the underlying assumptions of the kinetic theory derivation: high thermal Stokes number. Surprisingly, however, even though the low Knudsen number assumption, also associated with the kinetic theory derivation, is violated throughout most of the parameter space, it does not seem to affect the good quantitative accuracy of KT-TFM simulations.

42 ENGINEERING↗

NASA NTRS1992

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗

NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

NASA NTRS2020

High Ice Water Content in Tropical Cyclones during NASA/FAA Radar Flight Campaigns with Comparison to Numerical Simulations

Results from two recent flight campaigns into High Ice Water Content (HIWC) are summarized. HIWC was detected with airborne radar and measured via microphysical probes mounted on NASA’s DC-8. Results from the campaign demonstrate detectability of HIWC using airborne radar. In the cases examined, tropical storms are more proficient at producing large areas of HIWC than hurricanes. Three-dimensional, numerical simulations of tropical cyclones encountered during the campaign are also presented and are compared to measurements in order to understand conditions associated with HIWC. Within Hurricane Lane (2018), the highest ice water contents were found in regenerating thermal plumes located in the eye wall and in feeder bands. Within Tropical Storm Danny (2015), large areas of ice water content greater than 1 g m-3 were fed by broad-areas of convection occurring downwind from the center of circulation.

Aviation hazards↗

DOE OSTI2019

Numerical simulations of yield-based sooting tendencies of aromatic fuels using ReaxFF molecular dynamics

In this work, we present the first ReaxFF Molecular Dynamics (MD) simulations to quantitatively predict the sooting tendencies of various fuels. The specific sooting tendency metric used in this work is the Yield Sooting Index (YSI), which quantifies the effects of fuel molecular structure on soot yield. YSI has been experimentally measured and numerically simulated using computational fluid dynamics for a large range of fuels, but there is no existing reactive MD framework for YSI simulations. To adopt the experimental YSI concept, a multi-stage simulation procedure is designed using ReaxFF. As a proof-of-concept, toluene and phenol are selected as test fuels, since both have relatively well-understood reaction pathways. The ReaxFF YSI simulations are shown to capture key reaction events for both fuels selected that are consistent with existing chemical kinetic understanding. Toluene is shown to mostly retain its original aromatic ring structure and directly grow to larger aromatic compounds with multiple rings. On the other hand, the aromatic growth process from phenol is accompanied by carbon-loss reactions with CO release. In addition, a quantitative YSI formulation is also derived in ReaxFF and the ReaxFF-predicted YSI values are compared with measurement data and reasonably good agreement is achieved. The results reported in this work demonstrates that the ReaxFF-based framework can potentially be used to quantitatively predict the relative sooting tendencies, especially for fuels with unknown or poorly-known chemistry, to understand their sooting properties and in search for soot relevant reaction pathways from these fuels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1989

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗

NASA NTRS2008

Gas Core Reactor Numerical Simulation Using a Coupled MHD-MCNP Model

Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.

Kazeminezhad, F.↗

NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗