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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Overview of the Predictive Simulation Capability Element of the Plume Surface Interaction Project

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 GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element 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 PSC Element, 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 describes the PSC Element of the PSI Project as well as providing descriptions of recent accomplishments and remaining work. The overall structure of the PSC Element is broken down into four areas of focus. The first area is the Prediction of Plume Flow in low pressure environments. The development approach taken is the augmentation of the existing production-mode computational fluid dynamics (CFD) tool Loci/Chem, with targeted extensions necessary to accurately model rarefied conditions found in both Martian and Lunar applications. Production readiness and validation of predictive capability are the major objectives of this task. The second area is the effect of mixed continuum/rarefied flow on crater development and ejecta sheets. A new CFD application, Loci/GGFS (Gas Granular Flow Solver), is being developed which implements an Eulerian/Eulerian two-phase model of gas- and soil-phases in order to simulate the soil erosion, crater formation, and soil ejecta transport in a fully coupled simulation. This task seeks to verify Loci/GGFS is production-ready as well as perform validation studies to determine the degree of predictive capability achieved by Loci/GGFS. The third area is focused on the details and extension of particle phase modeling of soil. In this task, Discrete Element Modeling (DEM) techniques are used to perform direct simulations of complex soil particles under the action of forcing similar to that to be cause by PSI. The simulation results are then used to construct closures to the Eulerian model of the soil phase used by Loci/GGFS. The fourth area is gas-particle interaction modeling. In this task, experiments are being conducted as well as detailed simulation results are being studied to further understand the complexities of gas-particle interactions in dilute, intermediate, and high soil volume fraction regimes. Improved models of particle drag and the particle turbulent kinetic energy (PTKE) resulting from the interaction of gas flows within particle clouds are the objective of this task.

Jeff West↗
NASA NTRS2014

On the TFNS Subgrid Models for Liquid-Fueled Turbulent Combustion

This paper describes the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models used to emulate the major processes occurring in the turbulence-chemistry interaction. These two subgrid models are termed as LEM-like model and EUPDF-like model (Eulerian probability density function), respectively. Two-phase turbulent combustion in a single-element lean-direct-injection (LDI) combustor is calculated by employing the TFNS/LEM-like approach as well as the TFNS/EUPDF-like approach. Results obtained from the TFNS approach employing these two different subgrid models are compared with each other, along with the experimental data, followed by more detailed comparison between the results of an updated calculation using the TFNS/LEM-like model and the experimental data.

Liu, Nan-Suey↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Persistent meanders and eddies lead to quasi-steady Lagrangian transport patterns in a weak western boundary current

The Brazil Current (BC) is a weak western boundary current flowing along the Southwestern Atlantic Ocean. It is frequently described as a flow with intense mesoscale activity and relatively low volume transport between 5.0 and 10.0 Sv. We use a 13-year eddy-resolving primitive-equation simulation to show that the presence of persistent meanders and eddies leads to characteristic quasi-steady Lagrangian transport patterns, aptly extracted through climatological Lagrangian Coherent Structures (cLCS). The cLCS position the surface expression of the BC core along the 2000-m isobath, in excellent visual agreement with high resolution satellite sea-surface temperature and the model Eulerian mean velocity. The cLCS deformation pattern also responds to zonally persistent cross-shelf SSH transition from positive (high) values near coastline to low (negative) values between 200- and 2000-m and back to positive (high) offshore from the 2000-m isobath. Zonally-paired cyclonic and anticyclonic structures are embedded in this transition, also causing the cLCS to deform into chevrons. An efficient transport barrier is identified close to the 200-m isobath confirmed by limited inshore movement of drogued drifters and accurately indicated by an along slope maxima of climatological strength of attraction. We also show that the persistent cyclonic and anticyclonic structures may induce localized cross-shelf transport. Regions of low climatological strength of attraction coincide with large shelves and with stagnant synthetic trajectories. We also show that cLCS accurately depict trajectories initiated at the location of Chevron’s spill (November 2011) as compared to synthetic and satellite-tracked trajectories, and the outline of the oil from that accident. There is also an agreement between the large-scale oil slicks reaching the Brazilian beaches (from August 2019 to February 2020) and the strength of climatological attraction at the coast. Our work also clarifies the influence of persistent mesoscale structures on the regional circulation. The identification and quantitative description of climatological Lagrangian coherent structures is expected to improve the effectiveness of future emergency response to oil spills, contingency planning, rescue operations, larval and fish connectivity assessment, drifter launch strategies, waste pollutant and marine debris dispersion and destination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Combined Effect of the Wegener–Bergeron–Findeisen Mechanism and Large Eddies on Microphysics of Mixed-Phase Stratiform Clouds

Abstract The process of glaciation in mixed-phase stratiform clouds was investigated by a novel Lagrangian–Eulerian model (LEM) in which thousands of adjoining Lagrangian parcels moved within a turbulent-like velocity field with statistical parameters typical of the Arctic boundary layer. We used detailed bin microphysics to describe the condensation/evaporation processes in each parcel, in which droplets, aerosols, and ice particles were described using size distributions of 500 mass bins. The model also calculated aerosol mass inside droplets and ice particles. Gravitational sedimentation of droplets and ice particles was also accounted for. Assuming that droplet freezing is the primary source of ice particles, the Arctic clouds observed in Indirect and Semi-Direct Aerosol Campaign (ISDAC) were successfully simulated. The model showed that at a low ice particle concentration typical of ISDAC, large vortices (eddies) led to a quasi-stationary regime, in which mixed-phase St existed for a long time. The large eddies controlled the water partitioning in the mixed-phase clouds. Droplets formed and grew in updrafts, typically reaching the cloud top, and evaporated in downdrafts. Ice particles grew in updrafts and downdrafts. The Wegener–Bergeron–Findeisen (WBF) mechanism was efficient in downdrafts and some parts of updrafts, depending on ice concentration and vertical velocity. At low ice concentrations, the effect of ice on the phase partitioning was negligible. In this regime, liquid droplets were found near the cloud top, whereas ice particles precipitated through the cloud base. When ice concentration exceeded about 10 L −1 , the WBF mechanism led to glaciation of almost the entire cloud, with the exception of narrow cloud regions associated with strong updrafts. At ice particle concentrations of a few tens per liter, the oscillatory regime took place due to the ice–liquid interaction. The microphysical structure of mixed-phase St forms as a combined effect of cloud dynamics (large eddies) and the WBF mechanism.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Numerical experiments on the role of radiative processes in the development and maintenance of upper level clouds

A time-dependent, two-dimensional Eulerian model is presented whose purpose is to obtain more realistic parameterizations of extended high level cloudiness, and the results of a numerical experiment using the model are reported. The model is anelastic and the Bousinesque assumption is invoked. Unresolved subgrid scale processes are parameterized as eddy diffusion processes. Two phases of water are incorporated and equilibrium between them is assumed. The effects of infrared radiative processes are parametrically represented. Two simulations were conducted with identical initial conditions; in one of them, the radiation term was never turned on. The mean values of perturbation potential temperature at each level in the domain are plotted versus height after 15, 30, and 60 minutes of simulated time. The influence of the radiative term is seen to impose a cooling trend, leading to an increased generation of ice water and an increased generation of turbulent kinetic energy in the cloud layer.

Starr, D. O'C.↗
NASA NTRS1987

Evidence for methane and ammonia in the coma of comet P/Halley

Methane and ammonia abundances in the coma of Halley are derived from Giotto ion mass spectrometer data using an Eulerian model of chemical and physical processes inside the contact surface to simulate Giotto high-intensity spectrometer ion mass spectral data for mass-to-charge ratios (m/q) from 15 to 19. The ratio m/q = 19/18 as a function of distance from the nucleus is not reproduced by a model for a pure water coma. It is necessary to include the presence of NH3, and uniquely NH3, in coma gases in order to explain the data. A ratio of production rates Q(NH3)/Q(H2O) = 0.01 = 0.02 results in model values approximating the Giotto data. Methane is identified as the most probable source of the distinct peak at m/q = 15. The observations are fit best with Q(CH4)/Q(Q2O) = 0.02. The chemical composition of the comet nucleus implied by these production rate ratios is unlike that of the outer planets. On the other hand, there are also significant differences from observations of gas phase interstellar material.

Allen, M.↗
NASA NTRS1988

Predicting the velocity and azimuth of fragments generated by the range destruction or random failure of rocket casings and tankage

A calculational method is described which provides a powerful tool for predicting solid rocket motor (SRM) casing and liquid rocket tankage fragmentation response. The approach properly partitions the available impulse to each major system-mass component. It uses the Pisces code developed by Physics International to couple the forces generated by an Eulerian-modeled gas flow field to a Lagrangian-modeled fuel and casing system. The details of the predictive analytical modeling process and the development of normalized relations for momentum partition as a function of SRM burn time and initial geometry are discussed. Methods for applying similar modeling techniques to liquid-tankage-overpressure failures are also discussed. Good agreement between predictions and observations are obtained for five specific events.

Eck, Marshall↗
NASA NTRS2016

Using CATS Near-Real-time Lidar Observations to Monitor and Constrain Volcanic Sulfur Dioxide (SO2) Forecasts

An eruption of Italian volcano Mount Etna on 3 December 2015 produced fast-moving sulfur dioxide (SO2) and sulfate aerosol clouds that traveled across Asia and the Pacific Ocean, reaching North America in just 5 days. The Ozone Profiler and Mapping Suite's Nadir Mapping UV spectrometer aboard the U.S. National Polar-orbiting Partnership satellite observed the horizontal transport of the SO2 cloud. Vertical profiles of the colocated volcanic sulfate aerosols were observed between 11.5 and 13.5 km by the new Cloud Aerosol Transport System (CATS) space-based lidar aboard the International Space Station. Backward trajectory analysis estimates the SO2 cloud altitude at 7-12 km. Eulerian model simulations of the SO2 cloud constrained by CATS measurements produced more accurate dispersion patterns compared to those initialized with the back trajectory height estimate. The near-real-time data processing capabilities of CATS are unique, and this work demonstrates the use of these observations to monitor and model volcanic clouds.

space↗
DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Episodes of Cross-Polar Transport in the Arctic Troposphere During July 2008 as Seen from Models, Satellite, and Aircraft Observations

During the POLARCAT summer campaign in 2008, two episodes (2 5 July and 7 10 July 2008) occurred where low-pressure systems traveled from Siberia across the Arctic Ocean towards the North Pole. The two cyclones had extensive smoke plumes from Siberian forest fires and anthropogenic sources in East Asia embedded in their associated air masses, creating an excellent opportunity to use satellite and aircraft observations to validate the performance of atmospheric transport models in the Arctic, which is a challenging model domain due to numerical and other complications. Here we compare transport simulations of carbon monoxide (CO) from the Lagrangian transport model FLEXPART and the Eulerian chemical transport model TOMCAT with retrievals of total column CO from the IASI passive infrared sensor onboard the MetOp-A satellite. The main aspect of the comparison is how realistic horizontal and vertical structures are represented in the model simulations. Analysis of CALIPSO lidar curtains and in situ aircraft measurements provide further independent reference points to assess how reliable the model simulations are and what the main limitations are. The horizontal structure of mid-latitude pollution plumes agrees well between the IASI total column CO and the model simulations. However, finer-scale structures are too quickly diffused in the Eulerian model. Applying the IASI averaging kernels to the model data is essential for a meaningful comparison. Using aircraft data as a reference suggests that the satellite data are biased high, while TOMCAT is biased low. FLEXPART fits the aircraft data rather well, but due to added background concentrations the simulation is not independent from observations. The multi-data, multi-model approach allows separating the influences of meteorological fields, model realisation, and grid type on the plume structure. In addition to the very good agreement between simulated and observed total column CO fields, the results also highlight the difficulty to identify a data set that most realistically represents the actual pollution state of the Arctic atmosphere.

Sodemann, H.↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2023

Modeling a Miniature Voitenko-Like Explosively-Driven Shock Tube Using FLAG

Using the FLAG hydrocode, hypervelocity plasma jets observed in miniature explosively–driven shock tube experiments have been modeled with a relative error in velocity of less than 5% (a reduction from previous attempts which yielded errors of 37%). This suggests that FLAG may be well-suited to model HE systems involving the acceleration of gases through complex crack systems in HE devices: A phenomenon not accurately captured using current modeling capabilities. The Eulerian hydrocode model presented here leverages AMR to reduce computational cost and, importantly, it was found that mesh resolution and model geometry each significantly affected the average jet velocity. The model and the results are detailed here together with comments on challenges during its development and avenues for future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Comprehensive modeling of turbulent particulate flows using Eulerian and Lagrangian schemes

This paper addresses turbulent particle dispersion and modulation effects in dilute gas-particle turbulent flows using the Eulerian and Lagrangian modeling approaches. Gradient diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangian formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. For the particle size and loading considered, the turbulence transport equations must be modified to account for the modulation effects. Models are proposed for both Eulerian and Lagrangian schemes. Comparisons and predictions are made in fully developed gas-solid pipe flow and confined coaxial jets laden with particles. For the monodispersed system investigated, Eulerian approach is less expensive and gives more consistent results than the Lagrangian approach. The Lagrangian technique should be further developed to eliminate current inherent inconsistencies especially with regard to the symmetry boundary condition.

Adeniji-Fashola, A.↗
NASA NTRS2018

The NASA Eulerian Snow on Sea Ice Model (NESOSIM) v1.0: initial model development and analysis

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a new, open-source snow budget model that is currentlyconfigured to produce daily estimates of the depth and density of snow on sea ice across the Arctic Ocean throughthe accumulation season. NESOSIM has been developed in a three-dimensional Eulerian framework and includes two(vertical) snow layers and several simple parameterizations (accumulation, wind packing, advection_divergence, blowingsnow lost to leads) to represent key sources and sinks of snow on sea ice. The model is forced with daily inputs of snowfall and near-surface winds (from reanalyses), sea ice concentration (from satellite passive microwave data) and sea ice drift (from satellite feature tracking) during the accumulation season (August through April). In this study, we present the NESOSIM formulation, calibration efforts, sensitivity studies and validation efforts across an Arctic Ocean domain (100 km horizontal resolution). The simulated snow depth and density are calibrated with in situ data collected on drifting ice stations during the 1980s. NESOSIM shows strong agreement with the in situ seasonal cycles of snow depth and density, and shows good (moderate) agreement with the regional snow depth (density) distributions. NESOSIM is run for a contemporary period (2000 to 2015), with the results showing strong sensitivity to the reanalysis derived snowfall forcing data, with the Modern-Era Retrospective analysis for Research and Applications (MERRA) and the Japanese Meteorological Agency 55-year reanalysis (JRA-55) forced snow depths generally higher than ERA-Interim,and the Arctic System Reanalysis (ASR) generally lower. We also generate and force NESOSIM with a consensus "median" daily snowfall dataset from these reanalyses.The results are compared against snow depth estimates derived from NASA's Operation IceBridge (OIB) snow radar data from 2009 to 2015, showing moderate_strong correlations and root mean squared errors of 10 cm depending on the OIB snow depth product analyzed, similar to the comparisons between OIB snow depths and the commonly used modified Warren snow depth climatology. Potential improvements to this initial NESOSIM formulation are discussed in the hopes of improving the accuracy and reliability of these simulated snow depths and densities.

Petty, Alek A.↗
DOE OSTI2024

Application of a Comprehensive Lagrangian–Eulerian Spark-Ignition Model to Different Operating Conditions

Increasing engine efficiency is essential to reducing emissions, which is a priority for automakers. Unconventional modes such as boosted and highly dilute operation have the potential to increase engine efficiency but suffer from stability concerns and cyclic variability. To aid engineers in designing ignition systems that reduce cyclic variability in such engine operation modes, reliable and accurate spark-ignition models are necessary. Here, in this article, a Lagrangian–Eulerian spark-ignition (LESI) model is used to simulate electrical discharge, spark channel elongation, and ignition in inert or reacting crossflow within a combustion vessel, at different pressures, flow speeds, and dilution rates. First the model formulation is briefly revisited. Then, the experimental and simulations setups are presented. The results showcase the model’s ability to predict the secondary circuit voltage, current, and power signals, in addition to the spark channel elongation, for the inert cases, or flame front growth, for the reacting cases. The results also compare simulation spark channel and flame growth plots to experimental Schlieren images at different instants in time. This work serves to highlight LESI’s ability to predict the characteristics of discharge and ignition across a variety of operating conditions.

42 ENGINEERING↗
NASA NTRS2011

Assessment of Turbulence-Chemistry Interaction Models in the National Combustion Code (NCC) - Part I

This paper describes the implementations of the linear-eddy model (LEM) and an Eulerian FDF/PDF model in the National Combustion Code (NCC) for the simulation of turbulent combustion. The impacts of these two models, along with the so called laminar chemistry model, are then illustrated via the preliminary results from two combustion systems: a nine-element gas fueled combustor and a single-element liquid fueled combustor.

Wey, Thomas Changju↗