Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD Modeling”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2024

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Modeling Induction Stirring and Particle Tracking in Molten Uranium

Two independent numerical models have been developed to simulate the behavior of carbon impurities in molten uranium metal. Informed by experimental parameters, one model was created using the commercial software Star-CCM+ and compared with another developed using open-source codes, including OpenFOAM, Finite Element Method Magnetics (FEMM) and a First Passage Kinetic Monte Carlo (FPKMC) approach. The target experimental system features a 404g uranium metal charge containing an average carbon concentration of 139 ppm which was melted in a vacuum induction furnace at 1400° C then resolidified. The microstructures of the uranium and its impurities before and after melting have been characterized and reported separately. Prior to simulating the uranium-carbon system described, the numerical models were validated using a previously published nonradioactive experimental system to ensure agreement with expected output values. Focus has been placed on modeling velocity fields under induction stirring and impurity particle trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Effect of fuel composition and EGR on spark-ignited engine combustion with LPG fueling: Experimental and numerical investigation

This paper presents an experimental and numerical investigation of a spark-ignited (SI) cooperative fuel research (CFR) engine fueled with different liquefied petroleum gas (LPG) fuels and exhaust gas recirculation (EGR). Here, the effects of LPG fuel composition on engine combustion characteristics are initially evaluated at two different compression ratios (CR). Results show normal combustion at CR 7 and heavy knocking combustion at CR 10 for all the tested fuels, with a more substantial impact for the LPG fuel with high proportions of n-butane species. The Livengood-Wu (LW) integral method is then used to analyze the knock occurrence risk of individual fuel based on the reactivity of the tested fuels. The introduction of EGR then demonstrates the potential of knock intensity reduction below the borderline knock limit. A zonal-based kinetic interactions study is also performed to understand the knock mitigation effectiveness of EGR over the pressure–temperature domain relevant to SI engine operation. Finally, a multidimensional, computational fluid dynamics (CFD) simulation model is shown to predict the LPG combustion characteristics and presents the evolution of in-cylinder temperature and chemical species to demonstrate the development of end-gas autoignition events without and with EGR.

02 PETROLEUM↗
DOE OSTI2023

Evaluation of Atmospheric Conditions Leading to a Fumigation Event

During the afternoon of 30 January 2022, the Savannah River Site experienced unusual temperature conditions leading to a fumigation event that triggered safety alarms and caused considerable confusion about the cause of the event. Normally, it is assumed that fumigation events occur early in the day once surface heating has begun. While most fumigation events are related to the breakup of a nocturnal inversion, this event was related to synoptic atmospheric conditions, which provided a more unique scenario that led to the fumigation event. The unusual synoptic atmospheric conditions led to the downwash and fumigation of the elevated plume, causing a pollutant to mix rapidly to the surface. These conditions could have potentially harmed workers within the facility, since the plume was directed toward a building air intake system. Here, we seek to outline the conditions that led to this unusual fumigation event and provide results of two-dimensional (2D) wind modeling of the event to characterize these conditions for future operational guidance of the facility air intake systems. This work sets the groundwork for future high-resolution modeling to explore the mechanisms and thresholds affecting fumigation on the facility-specific short distance scales and improve forecasting of non-standard fumigation events to protect human health.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗
DOE OSTI2021

A study on computational fluid dynamics modeling of a refrigerated container for COVID-19 vaccine distribution with experimental validation

A key issue with the distribution of vaccines to prevent COVID-19 is the temperature level required during transport, storage, and distribution. Typical refrigerated transport containers can provide a temperature-controlled environment down to -30 °C. However, the Pfizer vaccine must be carefully transported and stored under a lower temperature between -80 °C and -60 °C. One way to provide the required temperature is to pack the vaccine vials into small packages containing dry ice. Dry ice sublimates from a solid to a gas, which limits the allowable transport duration. This can be mitigated by transporting in a -30 °C refrigerated container. Moreover, because the dry ice will sublimate and thereby release CO 2 gas into the transport container, monitoring the CO 2 concentration within the refrigerated container is also essential. In the present work, a 3D computational fluid dynamics model was developed based on a commercially available refrigerated container and validated with experimental data. The airflow, temperature distribution, and CO 2 concentration within the container were obtained from the simulations. The modeling results can provide guidance on preparing experimental setups, thus saving time and lowering cost, and also provide insight into safety precautions needed to avoid hazardous conditions associated with the release of CO 2 during vaccine distribution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Partially cracked ammonia mixture combustion in gas turbine combustors

Ammonia is gaining popularity for single fuel use in existing gas turbines, as it is carbonless, has higher hydrogen content per unit volume than liquid hydrogen, and facilitates easy storage and handling. However, the low flame speed of ammonia makes it prone to flame blow off and remains the chief concern for its use in power generation. Following previous studies that have shown improved combustion performance by partially cracking the fuel, in this computational study, efficacy of partially cracked ammonia as a replacement fuel for natural gas is evaluated. A high-pressure optically accessible combustor operated at ~10 bar pressure with non-vitiated heated air was used as the test platform. This combustor was equipped with a novel multifuel multi-tube micromixing injector (M3 injector). High resolution simulations of the flow fields both in the injector as well as the combustor were performed to understand the fuel mixing in the injector and in the inlet region of the combustor. These modeling results showed that the injector was highly effective for achieving mixture homogeneity at the entrance of the combustor for ammonia and natural gas injection, whereas, for hydrogen an earlier injection facilitating longer residence time would be beneficial. Also, these simulations show that while harnessing waste heat the ideal fuel dissociation fraction is between 0.4 and 0.55 for optimal operation.

Ammonia↗
DOE OSTI2021

End-gas autoignition fraction and flame propagation rate in laser-ignited primary reference fuel mixtures at elevated temperature and pressure

Knock in spark-ignited (SI) engines is initiated by autoignition of the unburned gasses upstream of spark-ignited, propagating, turbulent premixed flames. Knock propensity of fuel/air mixtures is typically quantified using research octane number (RON), motor octane number (MON), or methane number (MN; for gaseous fuels), which are measured using single-cylinder, variable compression ratio engines. In this study, knock propensity of SI fuels was quantified via observations of end-gas autoignition (EGAI) in unburned gasses upstream of laser-ignited, premixed flames at elevated pressures and temperatures in a rapid compression machine. Stoichiometric primary reference fuel (PRF; n-heptane/isooctane) blends of varying reactivity (50 ≤ PRF ≤ 100) were ignited using an Nd:YAG laser over a range of temperatures and pressures, all in excess of 545 K and 16.1 bar. Laser ignition produced outwardly-propagating premixed flames. High-speed pressure measurements and schlieren images indicated the presence of EGAI. The fraction of the total heat release attributed to EGAI (i.e., EGAI fraction) varied with fuel reactivity (i.e., octane number) and the time-integrated temperature of the end-gas prior to ignition. Flame propagation rates, which were measured using schlieren images, were only weakly correlated with octane number but were affected by turbulence caused by variation in piston timing. Under conditions of low turbulence, measured flame propagation rates approached one-dimensional premixed laminar flame speed computations performed at the same conditions. Experiments were simulated with a three-dimensional CONVERGE™ model using reduced chemical kinetics (121 species, 538 reactions). The simulations accurately captured the measured flame propagation rates, as well as the variation in EGAI fraction with fuel reactivity and time-integrated end-gas temperature. The simulations also revealed low-temperature heat release as well as formaldehyde and hydrogen peroxide formation in the end-gas upstream of the propagating flame, which increased the temperature and degree of chain branching in the end-gas, ultimately leading to EGAI.

42 ENGINEERING↗
DOE OSTI2025

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Developing Science-based fueling protocols for 250-bar hydrogen tanks onboard hydrogen ferries: Experiments and modeling

Combined modeling and experimental studies are reported of the fueling of a large (28 kg capacity) 250-bar Type IV hydrogen tank of the type being deployed on early hydrogen ferries, such as the MV Sea Change. The primary goal was to determine how such tanks can be successfully fueled with hydrogen (state of charge greater than 97%) within 45 minutes without exceeding the 82 °C temperature limit for such tanks. The modeling studies show that a gas injector is needed to avoid thermal stratification during hydrogen fueling which can result in potential hot spots. Empirically, precooling of the hydrogen to 0 °C was found to be needed in some of the cases examined, as ambient conditions greatly affected the need for a precooling to achieve the 45-minute fill time desired by end users. The experimental results afforded a calibration of the engineering model SOFIL for these large 250-bar tanks, which now enables using SOFIL to predict volume-averaged hydrogen fueling temperatures to an accuracy of ±2.7°C for these tanks. The model can therefore be used to evaluate potential scenarios for development of a standardized fueling methodology for ferries utilizing large Type-IV tanks.

08 HYDROGEN↗
DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Design of a gas-solid-solid separator to remove ash from circulating fluidized bed reactors

Cyclones are one of the most common types of gas-solid separators used in circulating fluidized bed boilers. However, cyclones typically do not allow ash to leave the system through the cyclone exit, causing ash to build up in the fluidized bed and necessitating additional systems to remove ash that builds up in the bed. In this study, an alternative “disengager” gas-solid separator is proposed as a way of inherently separating small and large solids, resulting in a gas-solid-solid separation system where ash is allowed to leave the system along with gas while the desired fluidized bed material is retained. Unlike cyclones, which rely on centrifugal force to separate solids and gas, the disengager separates based on entrainment velocity of the particles. Upwards-flowing gas and particles strike a deflection plate and enter the disengaging chamber where particles with low terminal velocity such as ash fines flow with the gas, while larger particles such as sand fall to the bottom of the separator and are returned to the fluidized bed. In this study, several different proposed disengager configurations are simulated and compared to a typical cyclone using computational fluid dynamic (CFD) simulations. Furthermore, it was found that separation efficiency in the disengager is strongly influenced by the size of the deflection plate, rather than by the size of the unit itself. The predicted separation efficiency showed that compared to a cyclone, the disengager design allows significantly more ash to exit the system but retains a similar amount of desirable material. Additionally, the disengager was predicted to not suffer significantly more erosion that a cyclone.

42 ENGINEERING↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Comparison of Single-Phase Convection in Additive Manufactured Versus Traditional Metal Foams

Abstract Metal foams have been often used for thermal management due to their favorable characteristics including high specific surface area (SSA), high thermal conductivity, and low relative density. However, they are accompanied by shortcomings including the significant contact resistances due to attachment method, as well as the need for characterization of foam parameters such as pore diameter and SSA. Additive manufactured (AM) metal foams would eliminate the substrate/foam thermal resistance, decrease the need for pre-usage characterization, and allow for tailoring structures, while also taking advantage of the characteristics of traditionally manufactured foams. A commercial, aluminum foam (nominally 5 pores per inch (PPI), 86.5% porosity) was analyzed using X-ray microcomputed tomography, and a custom-designed metal foam based on the cell diameter and porosity of the commercial sample was subsequently manufactured. Reduced domain computational fluid dynamics/heat transfer (CFD-HT) models were compared against experimental data. Postvalidation, the flow behavior, effect of varying attachment thermal conductivities, and thermal performance were numerically investigated, demonstrating the usefulness of validated pore-scale models, as well as the potential for improved performance using AM metal foams over traditionally manufactured foams.

Engineering↗
DOE OSTI2023

Experimental and Numerical Investigation of Flow Boiling in Additive Manufactured Foam Structures With Vapor Pathways

Abstract The unique properties of metal foams make them potential candidates for a range of applications, including microsystem thermal management. Using additive manufacturing to create foam-type structures can improve upon prior thermal solutions by eliminating thermal interface materials and allowing for customization/local control of parameters. In the present investigation, flow boiling in additive-manufactured metal foams is investigated both experimentally and numerically. Two test samples, one with uniform structure and the other with pathways for vapor removal, are compared both experimentally and numerically. A conjugate computational fluid dynamics and heat transfer (CFD-HT) model utilizing a three-dimensional volume of fluid (VOF) model with accompanying evaporation/condensation model provided in-depth visualization of the boiling flow phenomena. The experiments generated the thermohydraulic performance over a range of heat fluxes, demonstrating that the sample incorporating dedicated vapor pathways performed better in both pressure and heat transfer performance metrics compared to the uniform foam. Additionally, negative system-level effects (i.e., hydraulic oscillations) were shown to be abated using the vapor removal structures. The numerical model yielded further insight into the factors contributing to the improved performance. Results indicated the pathways functioned as vapor removal channels, allowing the generated vapor to vent from the foam structure into the lanes. Further computational investigations demonstrated changes in flow regimes, where the addition of vapor channels caused the flow to change from churn to annular. Bubble behavior unique to the vapor pathway structure was studied, showing stagnant regions that eject vapor into the channel.

Engineering↗
DOE OSTI2022

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗
DOE OSTI2023

Evaluation of a CFD Commercial Package to the Modeling of Acoustic Wave Propagation in Bubbly-Liquid Column

This document reports on the considerations and development of an initial computational model based on the first principles for predicting acoustic propagation in a wellbore based on a logging-while-drilling (LWD) tool. This tool could be used to create a large database of simulated signals mimicking gas-influx-while-drilling events (so-called “kicks”). While other logging techniques, such as electrical resistivity and bulk density, are available in LWD tools, the acoustic logging method was identified in a complementary experimental endeavor as the most appropriate for early kick detection due to its high sensitivity to the presence of gas in liquids.

03 NATURAL GAS↗