Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eulerian model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

CFD Analysis of Spray Combustion and Radiation in OMV Thrust Chamber

The Variable Thrust Engine (VTE), developed by TRW, for the Orbit Maneuvering Vehicle (OMV) uses a hypergolic propellant combination of Monomethyl Hydrazine (MMH) and Nitrogen Tetroxide (NTO) as fuel and oxidizer, respectively. The propellants are pressure fed into the combustion chamber through a single pintle injection element. The performance of this engine is dependent on the pintle geometry and a number of complex physical phenomena and their mutual interactions. The most important among these are (1) atomization of the liquid jets into fine droplets; (2) the motion of these droplets in the gas field; (3) vaporization of the droplets (4) turbulent mixing of the fuel and oxidizer; and (5) hypergolic reaction between MMH and NTO. Each of the above phenomena by itself poses a considerable challenge to the technical community. In a reactive flow field of the kind occurring inside the VTE, the mutual interactions between these physical processes tend to further complicate the analysis. The objective of this work is to develop a comprehensive mathematical modeling methodology to analyze the flow field within the VTE. Using this model, the effect of flow parameters on various physical processes such as atomization, spray dynamics, combustion, and radiation is studied. This information can then be used to optimize design parameters and thus improve the performance of the engine. The REFLEQS CFD Code is used for solving the fluid dynamic equations. The spray dynamics is modeled using the Eulerian-Lagrangian approach. The discrete ordinate method with 12 ordinate directions is used to predict the radiative heat transfer in the OMV combustion chamber, nozzle, and the heat shield. The hypergolic reaction between MMH and NTO is predicted using an equilibrium chemistry model with 13 species. The results indicate that mixing and combustion is very sensitive to the droplet size. Smaller droplets evaporate faster than bigger droplets, leading to a well mixed zone in the combustion chamber. The radiative heat flux at combustion chamber and nozzle walls are an order of negligible less than the conductive heat flux. Simulations performed with the heat shield show that a negligible amount of fluid is entrained into the heat shield region. However, the heat shield is shown to be effective in protecting the OMV structure surrounding the engine from the radiated heat.

Giridharan, M. G.↗
NASA NTRS1979

Mathematical model investigation of long-term transport of ocean-dumped sewage sludge related to remote sensing

An existing, three-dimensional, Eulerian-Lagrangian finite-difference model was modified and used to examine the transport processes of dumped sewage sludge in the New York Bight. Both in situ and laboratory data were utilized in an attempt to approximate model inputs such as mean current speed, horizontal diffusion coefficients, particle size distributions, and specific gravities. The results presented are a quantitative description of the fate of a negatively buoyant sewage sludge plume resulting from continuous and instantaneous barge releases. Concentrations of the sludge near the surface were compared qualitatively with those remotely sensed. Laboratory study was performed to investigate the behavior of sewage sludge dumping in various ambient density conditions.

Kuo, C. Y.↗
NASA NTRS1978

Investigation of turbulent diffusion in the extreme lower atmosphere

Turbulent diffusion in the extreme lower layer of the atmosphere (up to 5 m) has been investigated. Turbulent flow was simulated under dry, stable and calm conditions by means of a 3.04 m diameter fan installed at a field site situated on flat grassland. The ambient wind was continuously monitored by means of a cup anemometer placed outside the wake, and the temperature distribution was measured by four thermometers placed on an 18 m tower, also outside the wake. Balloons and red smoke were used to visualize the wake flow and investigate the predominant sizes of turbulent eddies and their streamwise behavior. The mean and turbulent velocities along the turbulence line were measured using an array of hot-wire anomometers. Results provide substantial verification of a recently proposed model (Koper and Sadeh, 1975; Koper et al., 1978) relating the Lagrangian to the Eulerian turbulent velocity autocorrelation. In this model the Lagrangian autocorrelation is given by a domain integral over a set of ordinary Eulerian autocorrelations acquired simultaneously at all points within the flow field in question, which is viewed as a turbulence 'box'.

Koper, C. A., Jr.↗
NASA NTRS1974

Study of effects of uncertainties of comet and asteroid encounter and contact guidance requirements. Part 2: Tumbling problem studies

The problem of determining the rotational motion of a tumbling celestial body of the asteroid type using spacecraft-acquired data is addressed. The rotational motion of the body is modeled by free-Eulerian motion of a triaxial, rigid body and its translational motion with respect, to a nonrotating, observing spacecraft, which is not thrusting, is assumed to be uniform during the time observations are made. The mathematical details which form the basis for a digital simulation of the motion and observations are presented. Two algorithms for determining the motion from observations for the special case of uniform rotational motion are given.

Cochran, J. E., Jr.↗
NASA NTRS1977

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Pavish, D. L.↗
NASA NTRS1992

Hypervelocity impact simulation for micrometeorite and debris shield design

A new capability has been developed for direct computer simulation of hypervelocity impacts on multi-plate orbital debris shields, for combinations of low shield thickness and wide shield spacing which place extreme demands on conventional Eulerian analysis techniques. The modeling methodology represents a novel approach to debris cloud dynamics simulation, a problem of long term interest in the design of space structures. Software implementation of the modeling methodology provides a new design tool for engineering analysis of proposed orbital debris protection systems.

Fahrenthold, Eric P.↗
NASA NTRS1994

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods for the analysis of complete jet engines. We have entitled this program by applying the two dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition, and solution capabilities were successfully tested. We then focused attention on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion that results from these structural displacements. This is treated by a new arbitrary Lagrangian-Eulerian (ALE) technique that models the fluid mesh motion as that of a fictitious mass-spring network. New partitioned analysis procedures to treat this coupled three-component problem are developed. These procedures involved delayed corrections and subcycling. Preliminary results on the stability, accuracy, and MPP computational efficiency are reported.

Felippa, C. A.↗
NASA NTRS2015

An Internal Thermal Environment Model of an Aluminized Solid Rocket Motor with Experimental Validation

Due to the severity of the internal solid rocket motor (SRM) environment, very few direct measurements of that environment exist; therefore, the appearance of such data provides a unique opportunity to assess current thermal/fluid modeling capabilities. As part of a previous study of SRM internal insulation performance, the internal thermal environment of a laboratory-scale SRM featuring aluminized propellant was characterized with two types of custom heat-flux calorimeters: one that measured the total heat flux to a graphite slab within the SRM chamber and another that measured the thermal radiation flux. Therefore, in the current study, a thermal/fluid model of this lab-scale SRM was constructed using ANSYS Fluent to predict not only the flow field structure within the SRM and the convective heat transfer to the interior walls, but also the resulting dispersion of alumina droplets and the radiative heat transfer to the interior walls. The dispersion of alumina droplets within the SRM chamber was determined by employing the Lagrangian discrete phase model that was fully coupled to the Eulerian gas-phase flow. The P1-approximation was engaged to model the radiative heat transfer through the SRM chamber where the radiative contributions of the gas phase were ignored and the aggregate radiative properties of the alumina dispersion were computed from the radiative properties of its individual constituent droplets, which were sourced from literature. The convective and radiative heat fluxes computed from the thermal/fluid model were then compared with those measured in the lab-scale SRM test firings and the modeling approach evaluated.

Martin, Heath T.↗
NASA NTRS2000

Age-of-Air, Tape Recorder, and Vertical Transport Schemes

A numerical-analytic investigation of the impacts of vertical transport schemes on the model simulated age-of-air and the so-called 'tape recorder' will be presented using an idealized 1-D column transport model as well as a more realistic 3-D dynamical model. By comparing to the 'exact' solutions of 'age-of-air' and the 'tape recorder' obtainable in the 1-D setting, useful insight is gained on the impacts of numerical diffusion and dispersion of numerical schemes used in global models. Advantages and disadvantages of Eulerian, semi-Lagrangian, and Lagrangian transport schemes will be discussed. Vertical resolution requirement for numerical schemes as well as observing systems for capturing the fine details of the 'tape recorder' or any upward propagating wave-like structures can potentially be derived from the 1-D analytic model.

Lin, S.-J.↗
NASA NTRS1980

Instrumented aircraft verification of clear-air radar detection of low-level wind shear

A technique is described in which the longitudinal wind component at airports, measured by Doppler radar, is put through a Fourier transformation to produce an energy density spectrum. The result is multiplied by the aircraft response characteristics to produce quantitative data on aircraft performance in a wind shear situation. The Doppler radar measurements are made along the intended flight path. Two case studies, along a 3 deg approach path, are described from measurements in two modes: Lagrangian, where predictions were computed one range gate ahead of the plane, and Eulerian, where samples were measured instantaneously along the flight path. Wind data obtained every second were interpolated to fit a numerical model by using a cubic spline. Eulerian data were superior to Lagrangian data in terms of ranging, and nearly equal for wind speeds. Improvements in current 75% accuracy are noted to be possible with shrouded antennas to reduce sidelobes and use of Doppler radar with higher resolution to eliminate effects of small scale disturbances within the pulse volume.

Mccarthy, J.↗
NASA NTRS1988

Mesoscale spiral vortex embedded within a Lake Michigan snow squall band - High resolution satellite observations and numerical model simulations

It is known that Great Lakes snow squall convection occurs in a variety of different modes depending on various factors such as air-water temperature contrast, boundary-layer wind shear, and geostrophic wind direction. An exceptional and often neglected source of data for mesoscale cloud studies is the ultrahigh resolution multispectral data produced by Landsat satellites. On October 19, 1972, a clearly defined spiral vortex was noted in a Landsat-1 image near the southern end of Lake Michigan during an exceptionally early cold air outbreak over a still very warm lake. In a numerical simulation using a three-dimensional Eulerian hydrostatic primitive equation mesoscale model with an initially uniform wind field, a definite analog to the observed vortex was generated. This suggests that intense surface heating can be a principal cause in the development of a low-level mesoscale vortex.

Lyons, Walter A.↗
NASA NTRS1993

CFD flowfield simulation of Delta Launch Vehicles in a power-on configuration

This paper summarizes recent work at McDonnell Douglas Aerospace (MDA) to develop and validate computational fluid dynamic (CFD) simulations of under expanded rocket plume external flowfields for multibody expendable launch vehicles (ELVs). Multi engine reacting gas flowfield predictions of ELV base pressures are needed to define vehicle base drag and base heating rates for sizing external nozzle and base region insulation thicknesses. Previous ELV design programs used expensive multibody power-on wind tunnel tests that employed chamber/nozzle injected high pressure cold or hot-air. Base heating and pressure measurements were belatedly made during the first flights of past ELV's to correct estimates from semi-empirical engineering models or scale model tests. Presently, CFD methods for use in ELV design are being jointly developed at the Space Transportation Division (MDA-STD) and New Aircraft Missiles Division (MDA-NAMD). An explicit three dimensional, zonal, finite-volume, full Navier-Stokes (FNS) solver with finite rate hydrocarbon/air and aluminum combustion kinetics was developed to accurately compute ELV power-on flowfields. Mississippi State University's GENIE++ general purpose interactive grid generation code was chosen to create zonal, finite volume viscous grids. Axisymmetric, time dependent, turbulent CFD simulations of a Delta DSV-2A vehicle with a MB-3 liquid main engine burning RJ-1/LOX were first completed. Hydrocarbon chemical kinetics and a k-epsilon turbulence model were employed and predictions were validated with flight measurements of base pressure and temperature. Zonal internal/external grids were created for a Delta DSV-2C vehicle with a MB-3 and three Castor-1 solid motors burning and a Delta-2 with an RS-27 main engine (LOX/RP-1) and 9 GEM's attached/6 burning. Cold air, time dependent FNS calculations were performed for DSV-2C during 1992. Single phase simulations that employ finite rate hydrocarbon and aluminum (solid fuel) combustion chemistry are currently in progress. Reliable and efficient Eulerian algorithms are needed to model two phase (solid-gas) momentum and energy transfer mechanisms for solid motor fuel combustion products.

Pavish, D. L.↗
NASA NTRS1984

A Lagrangian method for a two layer shallow water model

One of the biggest drawbacks of Eulerian methods is the presence of the nonlinear convective terms which leads to inaccurate representations of advection and fronts. An alternative and conceptually simpler approach is the Lagrangian approach in which the fluid particles themselves are tracked and equations are derived based on local spatial interaction. In this formulation the nonlinear convective terms do not appear. The main drawback with the Lagrangian formulation, however, is that the local spatial interactions are time dependent. That is, at each time level, one must know the current neighbors of a given particle to accurately compute the forces acting on the particle. The Voroni Mesh which is continuously deforming and always linking nearest neighbors is needed. The effectiveness of a simple Lagrangian method for simulating the onset of a Kelvin-Helmholtz instability in a 2 layer fluid is demonstrated.

Augenbaum, J. M.↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
NASA NTRS2001

Impact Cratering Calculations

This research is computational /theoretical and complements the Caltech experimental program. We have developed an understanding of the basic physical processes and produced computational models and implemented these into Eulerian and Lagrangian finite element codes. The key issues we have addressed include the conditions required for: faulting (strain localization), elastic moduli weakening, dynamic weakening (layering elastic instabilities and fluidization), bulking (creation of porosity at zero pressure) and compaction of pores, frictional melting (creation of pseudotachylytes), partial and selective devolatilization of materials (e.g. CaCO3, water/ice mixtures), and debris flows.

Ahrens, Thomas J.↗
NASA NTRS2010

Analysis of Flame Deflector Spray Nozzles in Rocket Engine Test Stands

The development of a unified tightly coupled multi-phase computational framework is described for the analysis and design of cooling spray nozzle configurations on the flame deflector in rocket engine test stands. An Eulerian formulation is used to model the disperse phase and is coupled to the gas-phase equations through momentum and heat transfer as well as phase change. The phase change formulation is modeled according to a modified form of the Hertz-Knudsen equation. Various simple test cases are presented to verify the validity of the numerical framework. The ability of the methodology to accurately predict the temperature load on the flame deflector is demonstrated though application to an actual sub-scale test facility. The CFD simulation was able to reproduce the result of the test-firing, showing that the spray nozzle configuration provided insufficient amount of cooling.

Sachdev, Jai S.↗
NASA NTRS2020

Winter Arctic Sea Ice Thickness from ICESat-2 Freeboards

NASA’s ICESat-2 mission was launched in September 2018 with the primary goal of monitoring our rapidly changing polar regions. The sole instrument onboard, the Advanced Topographic Laser Altimeter System (ATLAS), is now providing routine, very high-resolution, surface elevation data across the globe, including the Arctic and Southern Oceans. In this study we demonstrate our new processing chain for converting the along-track ICESat-2 sea ice freeboard product (ATL10) into sea ice thickness, focusing our initial efforts on the Arctic Ocean. For this conversion we primarily make use of snow depth and density data from the NASA Eulerian Snow on Sea Ice Model (NESOSIM). The coarse resolution (~100 km) snow data are redistributed onto the high resolution (~30 - 100 m) ATL10 freeboards using relationships obtained from snow depth and freeboard data collected by NASA’s Operation IceBridge mission. We present regional sea ice thickness distributions and highlight their seasonal evolution through our first winter-season of data collection. We include ice thickness uncertainty estimates, while also acknowledging the limitations of these estimates. We generate a gridded monthly thickness product and compare this with various monthly sea ice thickness estimates obtained from ESA’s CryoSat-2 satellite mission, with ICE 33 Sat-2 showing consistently lower thicknesses. Finally, we compare our February/March 2019 thickness estimates to ICESat February/March (19th February – 21st March) 2008 ice thickness estimates using the same input assumptions, which show a ~0.40 m or ~21% thinning across an inner Arctic Ocean domain in this 11-year time period.

Alek A Petty↗
NASA NTRS1989

Disk-instability model for soft-X-ray transients containing black holes

Consideration is given to the time-dependent behavior of the accretion disks in low-mass X-ray binaries as a model for soft X-ray transients. The thermal stability and vertical structure properties of the disk are studied in the range between 10 to the 7th and 10 to the 8th cm. An Eulerian time-dependent thin disk model is used to show that, for relevant accretion rates, the disk suffers a thermal instability leading to intermittent accretion onto the central compact object. The disk instability model reproduces the basic features of outbursts of soft X-ray transients. Also, the physical state of quiescence of the binary A0620-00 is examined, focusing on the possible existence of a black hole. It is suggested that the modulation of the accretion flow in the inner disk may be responsible for the bimodal behavior observed in the spectra of black hole candidates such as Cyg X-1, GX 339-4, and A0620-00.

Mineshige, Shin↗