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,567 records · Page 87

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Influence of local configuration on the backflow from small rocket thrusters

A numerical study has been made of the low density plume flow from a small thruster motor. For the case studied, the addition of a thin plume shield with a diameter 1.67 times the nozzle exit diameter was found to lead to a reduction in the backflow by several orders of magnitude. Calculations were also made for the three-dimensional flow produced by a pair of motors, in order to determine the increase in the backflow due to the interaction of the plumes. The calculations were made with the direct simulation Monte Carlo (DSMC) method, and one of the purposes of the work was to test a new and faster set of algorithms for this method.

Bird, G. A.↗

NASA NTRS1993

Mapping the stability region of the 3:2 Neptune-Pluto resonance

Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.

Levison, H. F.↗

NASA NTRS2001

Three-Dimensional Upward Flame Spreading in Partial-Gravity Buoyant Flows

Reduced-gravity environments have been used to establish low-speed, purely forced flows for both opposed- and concurrent-flow flame spread studies. Altenkirch's group obtained spacebased experimental results and developed unsteady, two-dimensional numerical simulations of opposed-flow flame spread including gas-phase radiation, primarily away from the flammability limit for thin fuels, but including observations of thick fuel quenching in quiescent environments. T'ien's group contributed some early flame spreading results for thin fuels both in opposed flow and concurrent flow regimes, with more focus on near-limit conditions. T'ien's group also developed two- and three-dimensional numerical simulations of concurrent-flow flame spread incorporating gas-phase radiative models, including predictions of a radiatively-induced quenching limit reached in very low-speed air flows. Radiative quenching has been subsequently observed in other studies of combustion in very low-speed flows including other flame spread investigations, droplet combustion and homogeneous diffusion flames, and is the subject of several contemporary studies reported in this workshop. Using NASA aircraft flying partial-gravity "parabolic" trajectories, flame spreading in purely buoyant, opposed-flow (downward burning) has been studied. These results indicated increases in flame spread rates and enhanced flammability (lower limiting atmospheric oxygen content) as gravity levels were reduced from normal Earth gravity, and were consistent with earlier data obtained by Altenkirch using a centrifuge. In this work, experimental results and a three-dimensional numerical simulation of upward flame spreading in variable partial-gravity environments were obtained including some effects of reduced pressure and variable sample width. The simulation provides physical insight for interpreting the experimental results and shows the intrinsic 3-D nature of buoyant, upward flame spreading. This study is intended to link the evolving understanding of flame spreading in purely-forced flows to the purely-buoyant flow environment, particularly in the concurrent flow regime; provide additional insight into the existence of steady flame spread in concurrent flows; and stimulate direct comparisons between opposed- and concurrent-flow flame spread. Additionally, this effort is intended to provide direct practical understanding applicable to fire protection planning for the habitable facilities in partial gravity environments of anticipated Lunar and Martian explorations.

Sacksteder, Kurt R.↗

NASA NTRS2011

Multi-Dimensional Calibration of Impact Dynamic Models

NASA Langley, under the Subsonic Rotary Wing Program, recently completed two helicopter tests in support of an in-house effort to study crashworthiness. As part of this effort, work is on-going to investigate model calibration approaches and calibration metrics for impact dynamics models. Model calibration of impact dynamics problems has traditionally assessed model adequacy by comparing time histories from analytical predictions to test at only a few critical locations. Although this approach provides for a direct measure of the model predictive capability, overall system behavior is only qualitatively assessed using full vehicle animations. In order to understand the spatial and temporal relationships of impact loads as they migrate throughout the structure, a more quantitative approach is needed. In this work impact shapes derived from simulated time history data are used to recommend sensor placement and to assess model adequacy using time based metrics and orthogonality multi-dimensional metrics. An approach for model calibration is presented that includes metric definitions, uncertainty bounds, parameter sensitivity, and numerical optimization to estimate parameters to reconcile test with analysis. The process is illustrated using simulated experiment data.

Horta, Lucas G.↗

NASA NTRSDate not supplied

Current Test Capabilities of NASA Langley’s Arc-Heated Scramjet Test Facility

The Arc-Heated Scramjet Test Facility at NASA Langley Research Center has demonstrated itself repeatedly as a key contributor to fundamental hypersonic research and scramjet engine technology development in the Mach 4.7 to 8 range. Since its inception in the 1970’s, the facility has hosted numerous scramjet engine configurations that explored the performance and operability of airframe-integrated hydrogen-fueled scramjet flow path concepts, successfully maturing both analysis tools and design technique and methodology for hypersonic propulsion. Contributions to major programs include extensive tests of the hydrogen-fueled NASP and Hyper-X scramjet engines, and the hydrocarbon-fueled HIFiRE Flight 2 scramjet flow path, in addition to fundamental in-house research. The facility capabilities have expanded to include direct-connect flow path support, gaseous hydrocarbon fuels, dynamic trajectory simulation and extension of flight simulation down to Mach 2.2. The low-Mach upgrades include increased air flow rate and controls capability and two new low-Mach nozzles (Mach 2.2 and 3.5). The facility recently completed verification of the new operational range through successful envelope expansion tests with the Mach 2.2 nozzle. With these new modifications, the operating envelope for the facility is now Mach 2.2 – 8.0(+), with a demonstrated dynamic pressure simulation of 2000 psf at the low end, and margin available for even more. Lastly, a continuous run time of over 5 minutes has been demonstrated, indicating that run time is limited only by test hardware or back-pressure requirements. This paper will provide a detailed description of these previously unreported upgrades and expanded simulation capabilities.

Hypersonic↗

NASA NTRS1992

Simulation of space measurements of vegetation canopy bidirectional reflectance factors

A vegetation canopy receives both monodirectional and scattered radiation which is then reflected back into the atmosphere according to its bidirectional reflectance distribution function (BRDF). Remote measurements are typically the angular and spectral distributions of radiations exiting the atmosphere. The evaluation of vegetation canopy bidirectional reflectance factors requires a numerical solution of the canopy radiative transfer equation. A horizontally homogeneous vegetation canopy, of finite physical depth, filled densely with small leaves and bounded by a flat Lambertian ground surface is considered. The canopy is illuminated spatially uniformly by monodirectional radiation. The scattered intensity distribution emerging at the top of the canopy in all directions in the upper hemisphere is evaluated. The magnitude, sign and angular distribution of atmospheric effects are wavelength specific and depend on the composition of the atmosphere. The net atmospheric effect is positive (negative) at the red (near-infrared) wavelength due to strong scattering (absorption) in the atmosphere. The angular distribution of these effects is bowl-shaped. Thus, off-nadir remote observations, especially at shorter wavelengths, are most affected by atmospheric perturbations.

Myneni, R. B.↗

NASA NTRS2002

Dual-Mode Combustion

The Department of Mechanical and Aerospace Engineering at the University of Virginia has conducted an investigation of the mixing and combustion processes in a hydrogen fueled dual-mode scramjet combustor. The experiment essentially consisted of the "direct connect" continuous operation of a Mach 2 rectangular combustor with a single unswept ramp fuel injector. The stagnation enthalpy of the test flow simulated a flight Mach number of 5. Measurements were obtained using conventional wall instrumentation and laser based diagnostics. These diagnostics included, pressure and wall temperature measurements, Fuel Plume Imaging (FPI) and Particle Image Velocimetry (PIV). A schematic of the combustor configuration and a summary of the measurements obtained are presented. The experimental work at UVa was parallel by Computational Fluid Dynamics (CFD) work at NASA Langley. The numerical and experiment results are compared in this document.

Goyne, Christopher P.↗

NASA NTRS2015

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Turblent boundary layer↗

NASA NTRS1984

Some features of auroral electric fields as seen in 2D numerical simulations

Results of 2D plasma simulations are presented and related to auroral observations. The formation of V-shaped potentials is studied with a 2 1/2 dimensional electrostatic particle-in-cell code for a magnetized plasma. It is shown that amplitudes for perpendicular electric fields are larger than for parallel electric fields, and for Te less than 100 eV, the amplitudes are comparable to the electric fields associated with the electrostatic shocks observed from the S3-3 satellite. The excitation of electrostatic ion-cyclotron EIC waves which occurs in the region below the parallel potential drop is discussed. In auroral plasmas EIC waves are observed above the V-shaped double layers in association with ion beams and field-aligned currents. The results also show that oppositely directed electric fields in the center and at the edges of the simulation region produce oppositely directed currents. Precipitating auroral ions in association with electron inverted-V events are seen by the DMSP-F6 satellite.

Thiemann, H.↗

NASA NTRS1991

GEOSIM: A numerical model for geophysical fluid flow simulation

A numerical model which simulates geophysical fluid flow in a wide range of problems is described in detail, and comparisons of some of the model's results are made with previous experimental and numerical studies. The model is based upon the Boussinesq Navier-Stokes equations in spherical coordinates, which can be reduced to a cylindrical system when latitudinal walls are used near the pole and the ratio of latitudinal length to the radius of the sphere is small. The equations are approximated by finite differences in the meridional plane and spectral decomposition in the azimuthal direction. The user can specify a variety of boundary and initial conditions, and there are five different spectral truncation options. The results of five validation cases are presented: (1) the transition between axisymmetric flow and baroclinic wave flow in the side heated annulus; (2) the steady baroclinic wave of the side heated annulus; (3) the wave amplitude vacillation of the side heated annulus; (4) transition to baroclinic wave flow in a bottom heated annulus; and (5) the Spacelab Geophysical Fluid Flow Cell (spherical) experiment.

Butler, Karen A.↗

NASA NTRS1992

Simulation of spatially evolving turbulence and the applicability of Taylor's hypothesis in compressible flow

For the numerical simulation of inhomogeneous turbulent flows, a method is developed for generating stochastic inflow boundary conditions with a prescribed power spectrum. Turbulence statistics from spatial simulations using this method with a low fluctuation Mach number are in excellent agreement with the experimental data, which validates the procedure. Turbulence statistics from spatial simulations are also compared to those from temporal simulations using Taylor's hypothesis. Statistics such as turbulence intensity, vorticity, and velocity derivative skewness compare favorably with the temporal simulation. However, the statistics of dilatation show a significant departure from those obtained in the temporal simulation. To directly check the applicability of Taylor's hypothesis, space-time correlations of fluctuations in velocity, vorticity, and dilatation are investigated. Convection velocities based on vorticity and velocity fluctuations are computed as functions of the spatial and temporal separations. The profile of the space-time correlation of dilatation fluctuations is explained via a wave propagation model.

Lee, Sangsan↗

NASA NTRS1994

The nature of the source of CO in Comet P/Halley

We show by numerical simulations that the CO distribution in the coma of Comet P/Halley as measured by neutral gas mass spectrometer aboard Giotto spacecraft (P. Eberhardt et al., 1987) can be explained by an extended jet source originating from vent 1 in the rotational model of M. J. S. Belton et al. (1991). This is an alternative interpretation for the origin of the CO distribution where Eberhardt et al. (1987) explained the observations in terms of a spherically symmetric extended source function. We still find nearly 1/3 of CO is due to CO sublimating directly at the nucleus, possibly from the vent itself. Approximately 50% of CO that originates in the extended jet source may be formed via photolysis of H2CO while the rest can be due to CO trapped in grain mantels and other C=O-bearing molecules. The overall production rates of CO in the two models differ by less than 29%. The prediction of the jet source model compares satisfactorily with the distribution of CO observed in the coma by the IUE satellite but unfortunately does not allow one to discriminate between the two models.

Samarasinha, Nalin H.↗

NASA NTRS2011

Hybrid CFD/CAA Modeling for Liftoff Acoustic Predictions

This paper presents development efforts at the NASA Marshall Space flight Center to establish a hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) simulation system for launch vehicle liftoff acoustics environment analysis. Acoustic prediction engineering tools based on empirical jet acoustic strength and directivity models or scaled historical measurements are of limited value in efforts to proactively design and optimize launch vehicles and launch facility configurations for liftoff acoustics. CFD based modeling approaches are now able to capture the important details of vehicle specific plume flow environment, identifY the noise generation sources, and allow assessment of the influence of launch pad geometric details and sound mitigation measures such as water injection. However, CFD methodologies are numerically too dissipative to accurately capture the propagation of the acoustic waves in the large CFD models. The hybrid CFD/CAA approach combines the high-fidelity CFD analysis capable of identifYing the acoustic sources with a fast and efficient Boundary Element Method (BEM) that accurately propagates the acoustic field from the source locations. The BEM approach was chosen for its ability to properly account for reflections and scattering of acoustic waves from launch pad structures. The paper will present an overview of the technology components of the CFD/CAA framework and discuss plans for demonstration and validation against test data.

Strutzenberg, Louise L.↗

NASA NTRS2014

Turbulent Concentration of MM-Size Particles in the Protoplanetary Nebula: Scaled-Dependent Multiplier Functions

The initial accretion of primitive bodies (asteroids and TNOs) from freely-floating nebula particles remains problematic. Here we focus on the asteroids where constituent particle (read "chondrule") sizes are observationally known; similar arguments will hold for TNOs, but the constituent particles in those regions will be smaller, or will be fluffy aggregates, and are unobserved. Traditional growth-bysticking models encounter a formidable "meter-size barrier" [1] (or even a mm-cm-size barrier [2]) in turbulent nebulae, while nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids [3]. Even if growth by sticking could somehow breach the meter size barrier, other obstacles are encountered through the 1-10km size range [4]. Another clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids [5]; scenarios leading directly from independent nebula particulates to this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios [6-8]. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. The typical sizes of planetesimals and the rate of their formation [7,8] are determined by a statistical model with properties inferred from large numerical simulations of turbulence [9]. Nebula turbulence can be described by its Reynolds number Re = L/eta sup(4/3), where L = ETA alpha sup (1/2) the largest eddy scale, H is the nebula gas vertical scale height, and α the nebula turbulent viscosity parameter, and η is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time tη. In the nebula, Re is far larger than any numerical simulation can handle, so some physical model is needed to extend the results of numerical simulations to nebula conditions.

protoplanetary nebula↗

NASA NTRS2019

Plasma simulations for the assessment of pole erosion in the Magnetically Shielded Miniature Hall Thruster (MaSMi)

Numerical simulations have been conducted to assess the erosion rates at the pole covers of the magnetically shielded miniature (MaSMi) Hall thruster. MaSMi is part of the Ascendant Sub-kW Transcelestial Electric Propulsion System (ASTRAEUS), whose design requirements demand a xenon throughput of 100 kg, thus enabling deep space missions using SmallSats. The results of our simulations show that the proposed thickness of the graphite pole covers will not be completely eroded at end of life. We find that the erosion rates are larger near the corners of the covers, which are directly exposed to the plasma in the acceleration channel.

Chaplin, Vernon H.↗

NASA NTRSDate not supplied

Assessment of Retrieved GMI Emissivity Over Land, Snow and Sea Ice in the GEOS System

Measurements from microwave sounders and imagers provide a valuable source of information including atmospheric temperature and water vapor in Numerical Weather Prediction (NWP) systems that assimilate these observations directly over water surfaces (oceans and other large water bodies). In a recent decadal survey, targeted observables in the Planetary Boundary Layer (PBL) were cited as a key need for future observations (NASEM, 2018). Microwave observations which sense in the PBL are currently available, however, utilizing surface-sensitive microwave observations for atmospheric data assimilation remains a challenge over land, snow and sea ice. This is in part due to the inability of surface emissivity models used by NWP data assimilation systems to simulate observations with sufficient accuracy. The GEOS-ADAS (Todling and el Akkraoui, 2018) which utilizes the Community Radiative Transfer Model (CRTM) (Han, 2006; Chen 2009) is no exception. The ECMWF system has retrieved instantaneous surface emissivity from surface-sensitive channels for SSMI/S and MHS radiance observations, and apply these estimates to the closest channels higher in frequency (Baordo and Geer 2016) in the calculation of simulated radiances. This approach currently is also being tested in the GEOS-ADAS for AMSU-A and ATMS radiances (Zhu et al. 2021). No or minimal emissivity spectral variability has been assumed in the above-mentioned studies. Recently, work by Munchak et al., 2020 (hereby referred to as M2020) provided a new database for emissivity over land, snow and sea ice retrieved from the NASA Global Precipitation Mission (GPM). Compared with Tool to Estimate Land Surface Emissivities at Microwave (TELSEM2; Wang et al., 2017), M2020 provides emissivities for more frequencies(i.e., 10.7 GHz V/H). Moreover, this database is unique in that it utilizes both active and passive data to retrieve surface emissivity and normalized radar cross section. While the emissivity values may be useful for other sensors, they are most applicable to the GPM Microwave Imager (GMI). In this work the GEOS-ADAS is modified to utilize emissivity values from Munchak et al, 2020 in place of values used by CRTM. Presently, only GMI radiances over ocean are used in the operational GEOS-ADAS. This study will focus on the GMI radiances over land, snow, and ice, as a first attempt to evaluate GMI radiances over these non-water surface types. Two cases are then presented, one with one week of observation minus background departures using the modified GEOS-ADAS, and one utilizing the original GEOS-ADAS. It should be noted that the surface emissivity models in CRTM are not state of the art and are scheduled to be replaced by the Community Surface Emissivity Module (CSEM; Chen and Weng, 2016). Simulations using default CRTM emissivity values are used merely as reference comparing against M2020, and is not a thorough comparison against other more state of the art modules such as CSEM.

GMI↗

NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗

NASA NTRS2011

Evaluating the Impacts of NASA/SPoRT Daily Greenness Vegetation Fraction on Land Surface Model and Numerical Weather Forecasts

The NASA Short-term Prediction Research and Transition (SPoRT) Center develops new products and techniques that can be used in operational meteorology. The majority of these products are derived from NASA polar-orbiting satellite imagery from the Earth Observing System (EOS) platforms. One such product is a Greenness Vegetation Fraction (GVF) dataset, which is produced from Moderate Resolution Imaging Spectroradiometer (MODIS) data aboard the NASA EOS Aqua and Terra satellites. NASA SPoRT began generating daily real-time GVF composites at 1-km resolution over the Continental United States (CONUS) on 1 June 2010. The purpose of this study is to compare the National Centers for Environmental Prediction (NCEP) climatology GVF product (currently used in operational weather models) to the SPoRT-MODIS GVF during June to October 2010. The NASA Land Information System (LIS) was employed to study the impacts of the new SPoRT-MODIS GVF dataset on land surface models apart from a full numerical weather prediction (NWP) model. For the 2010 warm season, the SPoRT GVF in the western portion of the CONUS was generally higher than the NCEP climatology. The eastern CONUS GVF had variations both above and below the climatology during the period of study. These variations in GVF led to direct impacts on the rates of heating and evaporation from the land surface. The second phase of the project is to examine the impacts of the SPoRT GVF dataset on NWP using the Weather Research and Forecasting (WRF) model. Two separate WRF model simulations were made for individual severe weather case days using the NCEP GVF (control) and SPoRT GVF (experimental), with all other model parameters remaining the same. Based on the sensitivity results in these case studies, regions with higher GVF in the SPoRT model runs had higher evapotranspiration and lower direct surface heating, which typically resulted in lower (higher) predicted 2-m temperatures (2-m dewpoint temperatures). The opposite was true for areas with lower GVF in the SPoRT model runs. These differences in the heating and evaporation rates produced subtle yet quantifiable differences in the simulated convective precipitation systems for the selected severe weather case examined.

Bell, Jordan R.↗