Engineering Papers⌕ Search

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.

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

At least 1,189 records · Page 66

NASA aerodynamics program

The annual accomplishments is reviewed for the Aerodynamics Division during FY 1991. The program includes both fundamental and applied research directed at the full spectrum of aerospace vehicles, from rotorcraft to planetary entry probes. A comprehensive review is presented of the following aerodynamics elements: computational methods and applications; CFD validation; transition and turbulence physics; numerical aerodynamic simulation; test techniques and instrumentation; configuration aerodynamics; aeroacoustics; aerothermodynamics; hypersonics; subsonics; fighter/attack aircraft and rotorcraft.

Williams, Louis J.↗

Overview Presentation

This report provides an overview presentation of the 2000 NPSS (Numerical Propulsion System Simulation) Review and Planning Meeting. Topics include: 1) a background of the program; 2) 1999 Industry Feedback; 3) FY00 Status, including resource distribution and major accomplishments; 4) FY01 Major Milestones; and 5) Future direction for the program. Specifically, simulation environment/production software and NPSS CORBA Security Development are discussed.

Lytle, John↗

Neutralizer Hollow Cathode Simulations and Validation with Experiments

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially-ionized gas in these devices. The results for the plasma are compared directly with Langmuir probe measurements. The computed keeper voltages are also compared with the observed values. Wherever model inputs and/or specific physics of the cathode discharge are uncertain, additional sensitivity calculations have been performed to quantify the uncertainties.

partially ionized gases↗

Investigation of Unsteady Flow Interaction Between an Ultra-Compact Inlet and a Transonic Fan

In the study presented, unsteady flow interaction between an ultra-compact inlet and a transonic fan stage is investigated. Future combat aircraft engines require ultra-compact inlet ducts as part of an integrated, advanced propulsion system to improve air vehicle capability and effectiveness to meet future mission needs. The main purpose of the current study is to advance the understanding of the flow interaction between a modern ultra-compact inlet and a transonic fan for future design applications. Many experimental/ analytical studies have been reported on the aerodynamics of compact inlets in aircraft engines. On the other hand, very few studies have been reported on the effects of flow distortion from these inlets on the performance of the following fan/compressor stages. The primary goal of the study presented is to investigate how flow interaction between an ultra-compact inlet and a transonic compressor influence the operating margin of the compressor. Both Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) approaches are used to calculate the unsteady flow field, and the numerical results are used to study the flow interaction. The present study indicates that stall inception of the following compressor stage is affected directly based on how the distortion pattern evolves before it interacts with the fan/compressor face. For the present compressor, the stall initiates at the tip section with clean inlet flow and distortion pattern away from the casing itself seems to have limited impacts on the stall inception of the compressor. A counter-rotating swirl, which is generated due to flow separation inside the s-shaped compact duct, generates an increased flow angle near the blade tip. This increased flow angle near the rotor tip due to the secondary flow from the counter-rotating vortices is the primary reason for the reduced compressor stall margin.

Compressor↗

Three-dimensional incompressible Navier-Stokes simulations of slender-wing vortices

The application of an implicit flux-difference splitting scheme to leading-edge vortex flows is described. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. Second-order spatial accuracy is achieved by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Good agreement between computed results and experimental data has been obtained for three different low-aspect-ratio wings.

Hsu, Chung-Hao↗

Interrogating Solar Flare Loop Models with IRIS Observations 2: Plasma Properties, Energy Transport, and Future Directions

During solar flares a tremendous amount of magnetic energy is released and transported through the Sun’s atmosphere and out into the heliosphere. Despite over a century of study, many unresolved questions surrounding solar flares are still present. Among those are how does the solar plasma respond to flare energy deposition, and what are the important physical processes that transport that energy from the release site in the corona through the transition region and chromosphere? Attacking these questions requires the concert of advanced numerical simulations and high spatial-, temporal-, and spectral-resolution observations. While flares are 3D phenomenon, simulating the NLTE flaring chromosphere in 3D and performing parameter studies of 3D models is largely outwith our current computational capabilities. We instead rely on state-of-the-art 1D field-aligned simulations to study the physical processes that govern flares. Over the last decade, data from the Interface Region Imaging Spectrograph (IRIS) have provided the crucial observations with which we can critically interrogate the predictions of those flare loop models. Here in Paper 2 of a two-part review of IRIS and flare loop models, I discuss how forward modelling flares can help us understand the observations from IRIS, and how IRIS can reveal where our models do well and where we are likely missing important processes, focussing in particular on the plasma properties, energy transport mechanisms, and future directions of flare modelling.

solar flares↗

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗

Interaction Between Tropical Convection and its Embedding Environment: An Energetics Analysis of a 2-D Cloud Resolving Simulation

The phase relation between the perturbation kinetic energy (K') associated with the tropical convection and the horizontal-mean moist available potential energy (bar-P) associated with environmental conditions is investigated by an energetics analysis of a numerical experiment. This experiment is performed using a 2-D cloud resolving model forced by the TOGA-COARE derived vertical velocity. The imposed upward motion leads to a decrease of bar-P directly through the associated vertical advective cooling, and to an increase of K' directly through cloud related processes, feeding the convection. The maximum K' and its maximum growth rate lags and leads, respectively, the maximum imposed large-scale upward motion by about 1-2 hours, indicating that convection is phase locked with large-scale forcing. The dominant life cycle of the simulated convection is about 9 hours, whereas the time scales of the imposed large-scale forcing are longer than the diurnal cycle. In the convective events, maximum growth of K' leads maximum decay of the perturbation moist available potential energy (P') by about 3 hours through vertical heat transport by perturbation circulation, and perturbation cloud heating. Maximum decay of P' leads maximum decay of bar-P by about one hour through the perturbation radiative, processes, the horizontal-mean cloud heating, and the large-scale vertical advective cooling. Therefore, maximum gain of K' occurs about 4-5 hours before maximum decay of bar-P.

Li, Xiaofan↗

CFD 2030 Grand Challenge: CFD-in-the-Loop Monte Carlo Flight Simulation for Space Vehicle Design

Flight qualification of space vehicles is markedly different from those typically employed for aircraft. The concept of an extensive flight test campaign for a space vehicle does not exist, and vehicle designers must look to alternative techniques for demonstrating robust and reliable performance of their vehicles prior to operational flight. A space vehicle may undergo only a handful of flight tests in its development cycle, with each flight representing a drastically different flight phase or flight configuration. For instance, NASA’s Space Launch System (SLS) launch vehicle and Orion spacecraft will only see a total of four flight demonstrations before flying a crew on its first operational mission, and each flight demonstrates a unique vehicle configuration and/or set of flight conditions. The SLS will be flown only one time before it becomes operational (Artemis 1). The Orion spacecraft Crew Module (CM) will have been tested twice, once on a Delta IV launch vehicle (Exploration Flight Test 1) and once as a fully integrated system with the SLS launch vehicle (Artemis 1). The Orion Launch abort system will have been tested twice, once in a pad abort scenario (Pad Abort 1) and once in an inflight abort scenario (Ascent Abort 2) on a modified Peacekeeper booster. Both of these latter tests involve only a boiler plate CM, not a functional Orion spacecraft. Thus, unlike aircraft, there is very little opportunity for engineers to assess and evaluate their preflight predictions. Instead, space vehicle designers rely on Monte Carlo flight simulations with detailed dispersions of predicted nominal flight behavior to determine how robust their design is to errors and uncertainties in the flight conditions their vehicle may encounter. These Monte Carlo analyses entail thousands of trajectory simulations to demonstrate that the vehicle can meet design requirements at a specified level of reliability. From an aerodynamics and aerothermodynamics perspective, these trajectory simulations are fueled by an extensive aerodynamic database that covers the complete range of expected flight conditions, vehicle configurations, and flight attitudes expected in a given mission. Today, these databases amount to a table of engineering parameters that can be quickly interrogated by the trajectory simulator. The aerodynamic and aerothermodynamic databases are assembled via a series of ground tests, empirical and analytical analysis, physics-based computational analysis, applicable past flight performance data, and in some cases, engineering judgment. These databases generally take years to assemble for a new space vehicle system and in the case of SLS/Orion, over a decade of test and analysis have been expended to develop the extensive databases required to cover the myriad of configurations and potential flight conditions required for the system. Recently, it has been proposed that Computational Fluid Dynamic (CFD) and computing capability may be reaching a point where it is foreseeable that CFD could be integrated directly into the production trajectory simulation tools used to design NASA’s space vehicles. To demonstrate this, NASA has embarked on two demonstrations of this type of capability, one where six degree of freedom flight trajectory simulation equations are embedded in an existing CFD solver and another where a production CFD solver is loosely coupled with a production trajectory simulation tool. These efforts represent an initial demonstration of a future approach to flight trajectory simulation, but they are a far cry from the capability required to perform a full-up CFD-in-the-loop Monte Carlo trajectory simulation. Therefore, this represents a viable grand challenge for computational methods addressing space vehicle design and development. The final paper/presentation will discuss the many hurdles, beyond simply raw computational power, to realizing this grand challenge and how they map directly to the CFD Vision 2030 ojectives. Among these are the wide range of flight conditions, including accelerating/decelerating flight, encountered by a space vehicle during launch and/or entry. The vehicle can also encounter numerous configuration changes, some of which can be quite drastic, during the course of its flight, so robust, automated geometry modeling, grid generation, and adaptation will play a huge role in reaching this goal. Multiply this by 1000’s of trajectory simulations occurring simultaneously in a given Monte Carlo analysis, and the problem readily scales to absorb virtually any size of supercomputer envisioned today. The concept of CFD-in-the-loop Monte Carlo trajectory simulation poses a formidable challenge for emerging and future computing systems, and it has the potential to shave years off the development cycle for aerodynamic and aerothermodynamic performance predictions as compared to today’s space vehicle design approach.

CFD 2030↗

A spacecraft-borne gradiometer mission analysis

Numerical simulations were performed to obtain the orbit- and attitude-determination requirements of a spacecraft-borne gradiometer mission. Results demonstrated that position determination of 300 meters in the along-track and cross-track directions and 50 meters in the radial direction are mission requirements. The optimal orientation of the gradiometer sensing plane is achieved when the spin vector elevation is 0 degrees. The attitude-determination requirements are 5 degrees resolution for spin-vector azimuth and 0.2 degree resolution for spin-vector elevation. When these requirements are met, 3-degree gravity anomalies can be recovered globally with an accuracy of 0.025/mm/sq s (2.5 mgals). The Appendix documents the mathematical procedures for estimating detailed gravity fields from gradiometer data.

Argentiero, P.↗

Efficient numerical simulation of an electrothermal de-icer pad

In this paper, a new approach to calculate the transient thermal behavior of an iced electrothermal de-icer pad was developed. The method of splines was used to obtain the temperature distribution within the layered pad. Splines were used in order to create a tridiagonal system of equations that could be directly solved by Gauss elimination. The Stefan problem was solved using the enthalpy method along with a recent implicit technique. Only one to three iterations were needed to locate the melt front during any time step. Computational times were shown to be greatly reduced over those of an existing one dimensional procedure without any reduction in accuracy; the curent technique was more than 10 times faster.

Roelke, R. J.↗

Analysis and visualization of complex unsteady three-dimensional flows

Flow field animation is the natural choice as a tool in the analysis of the numerical simulations of complex unsteady three-dimensional flows. The PLOT4D extension of the widely used PLOT3D code to allow the interactive animation of a broad range of flow variables was developed and is presented. To allow direct comparison with unsteady experimental smoke and dye flow visualization, the code STREAKER was developed to produce time accurate streaklines. Considerations regarding the development of PLOT4D and STREAKER, and example results are presented.

Van Dalsem, William R.↗

NASA's supercomputing experience

A brief overview of NASA's recent experience in supercomputing is presented from two perspectives: early systems development and advanced supercomputing applications. NASA's role in supercomputing systems development is illustrated by discussion of activities carried out by the Numerical Aerodynamical Simulation Program. Current capabilities in advanced technology applications are illustrated with examples in turbulence physics, aerodynamics, aerothermodynamics, chemistry, and structural mechanics. Capabilities in science applications are illustrated by examples in astrophysics and atmospheric modeling. Future directions and NASA's new High Performance Computing Program are briefly discussed.

Bailey, F. Ron↗

Mean Line Pump Flow Model in Rocket Engine System Simulation

A mean line pump flow modeling method has been developed to provide a fast capability for modeling turbopumps of rocket engines. Based on this method, a mean line pump flow code PUMPA has been written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The pump code can model axial flow inducers, mixed-flow and centrifugal pumps. The code can model multistage pumps in series. The code features rapid input setup and computer run time, and is an effective analysis and conceptual design tool. The map generation capability of the code provides the map information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of the code permit parametric design space exploration of candidate pump configurations and provide pump performance data for engine system evaluation. The PUMPA code has been integrated with the Numerical Propulsion System Simulation (NPSS) code and an expander rocket engine system has been simulated. The mean line pump flow code runs as an integral part of the NPSS rocket engine system simulation and provides key pump performance information directly to the system model at all operating conditions.

Veres, Joseph P.↗

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗

Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes

Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.

Zimmerman, Michael I.↗

Updates to Simulation of a Single-Element Lean-Direct Injection Combustor Using Arbitary Polyhedral Meshes

This paper summarizes the procedures of (1) generating control volumes anchored at the nodes of a mesh; and (2) generating staggered control volumes via mesh reconstructions, in terms of either mesh realignment or mesh refinement, as well as presents sample results from their applications to the numerical solution of a single-element LDI combustor using a releasable edition of the National Combustion Code (NCC).

Power↗