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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CFD investigation of low-attrition air-reactor designs for the NETL chemical-looping combustion reactor

The flow dynamics and attrition characteristics of three air reactor designs are investigated using an Eulerian-Eulerian kinetic theory model. The hydrodynamics of these reactor designs are analyzed to identify the flow regimes in which the particles evolve. The investigated designs which differ in the configuration of the secondary flow ports, exhibit bubbling and turbulent fluidization regimes within similar spatial regions. The most substantial differences in the hydrodynamics was in the turbulent fluidized regions near the secondary inlets. The mechanical stressing environments undergone by the particles were analyzed using the collision energy spectra, reconstructed using the local flow field properties (i.e. granular temperature and particle number density). The spectra show that the increase in the specific collision power is associated with the regions near the secondary injector ports, which are in a turbulent fluidization regime. A simple analysis relying upon a partitioning of the energy spectrum, showed that attrition rate may be reduced by about 34% over the original design.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Lagrangian dynamics of particle transport in oral and nasal breathing

Here, we present a large-eddy simulation (LES) of saliva particle transport during normal human breathing through the nose and mouth. The flow of the air–saliva mixture is modeled using an Eulerian LES that is coupled with a Lagrangian particle tracking module to obtain trajectories of saliva particles in a room with stagnant air conditions. The coupled Eulerian–Lagrangian simulation yields novel insights into the intricate dynamics of Lagrangian coherent structures (LCS) and fundamental material lines that emerge from the saliva particles' trajectories during several breathing cycles. Further, we systematically compare the quantitative LCS diagnostics of mouth breathing with those of mouth and nose normal breathing. Analyzing the simulation results of human breathing from the mouth and nose, we show that, soon after the first breathing cycle, saliva particles form a series of roll-up vortex rings that propagate forward. The forward propagation of these vortex rings leads to the formation of an asymmetrical primary forefront vortex. The individual vortex rings continuously propagate forward, merging with the forefront vortex, and ascending along the limb of the leading vortex.

42 ENGINEERING↗
DOE OSTI2023

Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design

Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Impact ejecta on the moon

The partitioning of energy and the distribution of the resultant ejecta on the moon is numerically modeled using a Eulerian finite difference grid. The impact of an iron meteoroid at 15 km/sec on a gabbroic anorthosite lunar crust is examined. The high speed impact induced flow is described over the entire hydrodynamic regime from a time where the peak pressures are 6 Mbar until the stresses everywhere in the flow are linearly elastic, and less than 5 kbar. Shock-induced polymorphic phase changes, (plagioclase and pyroxene to hollandite and perovskite), and the subsequent reversion to low pressure phases are demonstrated to enhance shock wave attenuation. A rate-dependent equation of state is used for describing the hysteretic effect of the phase change. Ballistic equations for a spherical planet are then applied to material with net velocity away from the moon.

Okeefe, J. D.↗
NASA NTRS1988

A model-based diagnostic study of the rapid development phase of the Presidents' Day cyclone

A model simulation of the rapid development phaser of the Presidents' Day cyclone on February 19, 1979, is studied. Model trajectories and Eulerian analyses show that, during this phase, three airstreams converge into the cyclogenetic region. The evolution of the PJ-trough system and its associated tropopause fold is described, including the simulated surface low and the vertical profile of mass divergence which contributes to the decreasing sea-level pressure. The various terms of the vorticity equation are evaluated to identify the processes contributing to the increase in absolute vorticity. The nature and origin of airstreams entering the low-level circulation associated with rapid cyclogenesis are examined from a Lagrangian perspective.

Whitaker, Jeffrey S.↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1981

A diagnostic model for studying daytime urban air quality trends

A single cell Eulerian photochemical air quality simulation model was developed and validated for selected days of the 1976 St. Louis Regional Air Pollution Study (RAPS) data sets; parameterizations of variables in the model and validation studies using the model are discussed. Good agreement was obtained between measured and modeled concentrations of NO, CO, and NO2 for all days simulated. The maximum concentration of O3 was also predicted well. Predicted species concentrations were relatively insensitive to small variations in CO and NOx emissions and to the concentrations of species which are entrained as the mixed layer rises.

Brewer, D. A.↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2013

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗
NASA NTRS2006

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗
DOE OSTI2020

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗
NASA NTRS1992

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗
DOE OSTI2021

Study the thermal radiation effects in gas-solid flows with gray and non-gray P1 models implemented in MFiX

Thermal radiation is a dominant mode of heat transfer in combustion/gasification, packed/circulating bed reactors, and energy storage/conversion devices. In hydrocarbon fuel combustion, the high absorption and emission of product gases like CO 2 and H 2 O significantly affect the heat transfer characteristics. It is essential to understand the radiative energy propagation in such applications. In the present study, we developed a framework to implement the P-1 radiation model in MFiX for Eulerian-Eulerian (Two fluid-TFM) and Eulerian-Lagrangian (Discrete element-DEM and Particle in cell-PIC) models. The P-1 is a simplification of the spherical harmonics method. The framework has been verified on various flue gas conditions. The verified solver results are validated against the benchmark results available in the literature. The verified and validated gray and non-gray weighted sum of gray gas models is evaluated for dry and wet flue gas conditions. Further, they are also used to study gas-solid flows in a fluidized bed/riser.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
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.↗