Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 217 records · Page 12

Numerical modeling based machine learning approach for the optimization of falling - film evaporator in thermal desalination application

Scale formation that drastically increases thermal resistance and reduces freshwater production remains a critical challenge in thermal desalination. Novel designs of falling film evaporator and optimal operating condition hold great promise to mitigate scale formation, and increase heat transfer performance and fresh water production. In this work, CFD simulation based machine learning and multi-objective optimization are performed to identify optimal conditions and tube arrangement for evaporator. Non-dominated sorting genetic algorithm is adopted to determine and analyze the optimal pareto front for multiple objectives in desalination criteria. The errors of training, validation, and testing set are computed to identify an optimal hyperparameter set. For performance ratio, fouling resistance, and water production rate, the average relative error is 2.26%, 3.67%, and 3.24%. At pareto front, both performance ratio and water production rate increase at high temperature with fouling resistance (thermal resistance of the fouling layer) increasing as well. Tradeoffs between mitigating scale formation and enhancing desalination performance are evaluated in optimizations for different objectives. Finally, potential optima are identified and can be applied as guidelines to determine evaporator design and system operating conditions.

42 ENGINEERING↗

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗

Effect of coaxial nozzle wear on catchment efficiency in direct energy deposition built components

Laser based Direct Energy Deposition (DED) systems using metallic powder feedstock are recognized as a promising manufacturing method for their ability to shorten production cycles and create complex part geometries. Components are built by generating a melt pool with a high-power laser beam while material is coaxially injected and left to solidify. An impediment to large scale use of DED lies in poor powder catchment efficiency, the condition in which a portion of injected powder escapes the melt pool resulting in a ratio of decreased printed material mass to mass of supplied feedstock. The wear state of a coaxial nozzle on a DED system within a hybrid manufacturing machine tool has been observed to decrease catchment efficiency over time. Herein this study investigates this effect by adapting flow visualization techniques to an in-situ process monitoring format, the implementation of a Computational Fluid Dynamics (CFD) simulation, and deposition testing. Nozzle geometric defects due to wear are identified and categorized, and the impact of nozzle tip wear, resulting in axial tip reduction, on powder catchment efficiency is proven by multiple calculation methods. A linear correlation between catchment efficiency and powder stream diameter is identified, causing a 15–20% loss in efficiency sustained over incremental nozzle tip reduction up to − 1 mm. These results provide a foundation for further study of wear effects and Zero defect manufacturing solutions for powder fed DED systems.

42 ENGINEERING↗

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗

Status and performance of the GlueX DIRC

In support of the future physics program at Ce+BAF, we propose a prototype positron target. The prototype target consists of a tungsten disk mounted on a water-cooled copper support structure to enhance heat removal. The primary goals of this work are to benchmark CFD simulations, ensure vacuum integrity at the 10^-6 Pa level, achieve stable target rotation up to 10 Hz, and implement efficient water-based cooling for sustained high-power operation. To manage the thermal load, the target is designed to rotate, allowing the deposited heat to spread over a larger region, and thereby reducing localized thermal stress while maintaining the maximum temperature below 1000 K. To experimentally evaluate the thermal performance and structural stability of the target, we plan to conduct high-power laser heating tests at Laser Lab 5, LERF Building, JLab. These tests will serve as a stepping-stone for the positron production target at CEBAF.

Stevens, J.R. [College of William and Mary, Willia↗

Experimental Insights into the Coupling of Methane Combustion and Steam Reforming in a Catalytic Plate Reactor in Transient Mode

The microstructured reactor concept is very promising technology to develop a compact reformer for distributed hydrogen generation. In this work, a catalytic plate reactor (CPR) is developed and investigated for the coupling of methane combustion (MC) and methane steam reforming (MSR) over Pt/Al 2 O 3 -coated microchannels in cocurrent and counter-current modes in transient experiments during start-up. A three-dimensional (3D) computational fluid dynamics (CFD) simulation shows uniform velocity and pressure distribution profiles in microchannels. For a channel velocity from 5.1 to 57.3 m/s in the combustor, the oxidation of methane is complete and self-sustainable without explosion, blow-off, or extinction; nevertheless, flashbacks are observed in counter-current mode. In the reformer, the maximum methane conversion is 84.9% in cocurrent mode, slightly higher than that of 80.2% in counter-current mode at a residence time of 33 ms, but at the cost of three times higher energy input in the combustor operating at ~1000 °C. Nitric oxide (NO) is not identified in combustion products, but nitrous oxide (N 2 O) is a function of coupling mode and forms significantly in cocurrent mode. Furthermore, this research would be helpful to establish the start-up strategy and environmental impact of compact reformers on a small scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature Homogeneity under Selective and Localized Microwave Heating in Structured Flow Reactors

Selective heating of different phases of multiphase systems via microwaves can result in energy savings and suppression of side reactions. However, materials properties and operating conditions that maximize temperature gradients are poorly understood. Here we utilize computational fluid dynamics (CFD) computations and temperature measurements in structured flow reactors (monoliths) in a monomodal microwave cavity to assess the temperature difference between the walls and the fluid and develop a simple lumped model to estimate when temperature gradients exist. We also explore the material’s thermal and electrical properties of structured reactors for isothermal catalyst conditions. We propose that CFD simulations can be used as a nonintrusive, predictive tool of temperature homogeneity. Importantly, we demonstrate that localized heating in the bed under several conditions rather than selective heating is responsible for the selectivity enhancement. Furthermore, our results indicate that structured beds made of high thermal conductivity materials avoid arcing and enable temperature homogeneity and low electrical conductivity materials allow microwaves to penetrate the domain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Flow Configurations on Response Time and Data Quality in Real-Time, In-Line Fourier Transform Infrared (FTIR) Monitoring of Viscous Flows

The real-time, in-line monitoring of continuous flow concentrations is widely conducted via infrared (IR) spectrometry by using a flow cell connected to a reactive flow stream. For protective purposes, the IR sensor tip is typically offset from the flow. This offset can cause the formation of a stagnant boundary layer above the sensor, especially when dealing with high-viscosity fluids. As a result, the IR signal response time is often controlled by the slow diffusional exchange of fluid in the boundary layer, as confirmed via 2D computational fluid dynamics (CFD) simulations. We evaluated several flow configuration modifications in a typical IR flow cell in order to identify the changes to the flow dynamics that enable improved response times with minimal changes to the cell configuration: the use of (i) vertical flow, where the standard horizontal flow over the sensor is redirected to contact vertically with the sensor, (ii) a static mixer to create radial flow momentum above the IR sensor, and (iii) horizontal or vertical nozzles to direct the flow toward the IR sensor. The vertical flow configuration did not show any significant improvement over the standard horizontal flow configuration. However, the static mixer, horizontal nozzle, and vertical nozzle configurations all resulted in markedly improved response times and signal quality, albeit at the expense of a higher pressure drop across the flow cell. These results point toward straightforward, user-accessible modifications of in-line IR flow cells that result in significant improvements in signal stability and acquisition times.

42 ENGINEERING↗

Chemical influence of ammonia on suppressing soot formation pathways

There is interest in mixing ammonia (NH3) with hydrocarbons, as this decreases overall CO2 emissions from the fuel mixture when burned. Co-firing NH3 with hydrocarbons like methane (CH4) also enables stable combustion over a larger set of operating conditions in engines and gas turbines than firing only NH3 permits. When NH3 is co-fired with hydrocarbons, the influence of NH3 on soot formation chemistry becomes relevant. In this work, soot concentrations in laminar nonpremixed co-flow flames of NH3/CH4 and N2/CH4 were analyzed using color-ratio pyrometry. The concentrations of combustion species (CH4, CO2, H2O, O2, C2H2, C6H6, HCN) across the different flames were also determined using on-line electron-impact mass spectrometry. The flames were also interrogated with 2D CFD simulations, which employ a detailed chemical mechanism tailored to capturing fuel-nitrogen chemistry. Addition of NH3 to CH4 flames was found to decrease soot concentrations more significantly than N2-addition (Fig. 1). The NH3/CH4 flames were found to have higher flame temperatures than the N2/CH4 flames, implying the suppression of soot afforded by NH3 is due to a chemical influence of NH3, as opposed to thermal or dilution effects. Based on experimental and simulated soot volume fractions and species concentrations, soot formation pathways affected by NH3 are postulated.

10 SYNTHETIC FUELS↗

Virtual sensing-enabled digital twin framework for real-time monitoring of nuclear systems leveraging deep neural operators

Abstract Real-time monitoring is a foundation of nuclear digital twin technology, crucial for detecting material degradation and maintaining nuclear system integrity. Traditional physical sensor systems face limitations, particularly in measuring critical parameters in hard-to-reach or harsh environments, often resulting in incomplete data coverage. Machine learning-driven virtual sensors offer a transformative solution by complementing physical sensors in monitoring critical degradation indicators. This paper introduces the use of Deep Operator Networks (DeepONet) to predict key thermal-hydraulic parameters in the hot leg of pressurized water reactor. DeepONet acts as a virtual sensor, mapping operational inputs to spatially distributed system behaviors without requiring frequent retraining. Our results show that DeepONet achieves low mean squared and Relative L2 error, making predictions 1400 times faster than traditional CFD simulations . These characteristics enable DeepONet to function as a real-time virtual sensor, synchronizing with the physical system to track degradation conditions and provide insights within the digital twin framework for nuclear systems.

Hossain, Raisa↗

UO 2 F 2 particulate formation in an impinging jet gas reactor

Computational fluid dynamics (CFD), chemical kinetics, and aerosol dynamics were combined to model the formation of uranyl fluoride particles from the hydrolysis of uranium hexafluoride gas. This chemical process has been studied for many decades, however, detailed measurements of aerosol formation have only become available in the past few years which provide a basis for model comparison. Here, CFD simulations predicted complicated flow patterns in the impinging jet gas reactor. Aerosol formation simulations also predicted higher mass concentrations than were observed experimentally. This suggested that experimental data provided a partial representation which was subsequently enhanced by modeling and simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unconventional gas-phase preparation of the prototype polycyclic aromatic hydrocarbon naphthalene (C 10 H 8 ) via the reaction of benzyl (C 7 H 7 ) and propargyl (C 3 H 3 ) radicals coupled with hydrogen-atom assisted isomerization

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the interstellar medium and in meteorites such as Murchison and Allende and signify the missing link between resonantly stabilized free radicals and carbonaceous nanoparticles (soot particles, interstellar grains). However, the predicted lifetime of interstellar PAHs of some 108 years imply that PAHs should not exist in extraterrestrial environments suggesting that key mechanisms of their formation are elusive. Exploiting a microchemical reactor and coupling these data with computational fluid dynamics (CFD) simulations and kinetic modeling, we reveal through an isomer selective product detection that the reaction of the resonantly stabilized benzyl (C 7 H 7 ) and the propargyl (C 3 H 3 ) synthesizes the simplest representative of PAHs – the 10π Hückel aromatic naphthalene (C 10 H 8 ) molecule – via the novel Propargyl Addition–BenzAnnulation (PABA) mechanism. The gas-phase preparation of naphthalene affords a versatile concept of the reaction of combustion and astronomically abundant propargyl radicals with aromatic radicals carrying the radical center at the methylene moiety (aromatic-CH 2 ) as a previously passed over source of aromatics in high temperature environments thus bringing us closer to an understanding of the aromatic universe we live in.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sneezing and asymptomatic virus transmission

The novel coronavirus disease (COVID-19) spread pattern continues to show that geographical barriers alone cannot contain a virus. Asymptomatic carriers play a critical role in the nature of this virus quickly escalating into a global pandemic. Asymptomatic carriers may transmit the virus unintentionally through sporadic sneezing. A novel Computational Fluid Dynamics (CFD) approach has been proposed with a realistic modeling of a human sneeze achieved by the combination of state-of-the-art experimental and numerical methods. This modeling approach may be suitable for future engineering analyses aimed at reshaping public spaces and common areas, with the main objective to accurately predict the spread of aerosol and droplets that may contain pathogens. This study shows that the biomechanics of a human sneeze, including complex muscle contractions and relaxations, can be accurately modeled by the angular head motion and the dynamic pressure response during sneezing. These have been considered as the human factors and were implemented in the CFD simulation by imposing a momentum source term to the coupled Eulerian–Lagrangian momentum equations. The momentum source was modeled by the measured dynamic pressure response in conjunction with the angular head motion. This approach eliminated the need to create an ad hoc set of inlet boundary conditions. With this proposed technique, it is easier to add multiple fixed and/or moving sources of sneezes in complex computational domains. Additionally, extensive sensitivity analyses based on different environmental conditions were performed, and their impact was described in terms of potential virus spread.

42 ENGINEERING↗

Effect of Nozzle Curvature on Supersonic Gas Jets Used in Laser-Plasma Acceleration

Supersonic gas jets produced by converging-diverging (C-D) nozzles are commonly used as targets for laser-plasma acceleration (LPA) experiments. A major point of interest for these targets is the gas density at the region of interaction where the laser ionizes the gas plume to create a plasma, providing the acceleration structure. Tuning the density profiles at this interaction region is crucial to LPA optimization. A "flat-top" density profile is desired at this line of interaction to control laser propagation and high energy electron acceleration, while a short high-density profile is often preferred for acceleration of lower-energy tightly-focused laser-plasma interactions. A particular design parameter of interest is the curvature of the nozzle's diverging section. We examine three nozzle designs with different curvatures: the concave "bell", straight conical and convex "trumpet" nozzles. We demonstrate that, at mm-scale distances from the nozzle exit, the trumpet and straight nozzles, if optimized, produce "flat-top" density profiles whereas the bell nozzle creates focused regions of gas with higher densities. An optimization procedure for the trumpet nozzle is derived and compared to the straight nozzle optimization process. We find that the trumpet nozzle, by providing an extra parameter of control through its curvature, is more versatile for creating flat-top profiles and its optimization procedure is more refined compared to the straight nozzle and the straight nozzle optimization process. Furthermore, we present results for different nozzle designs from computational fluid dynamics (CFD) simulations performed with the program ANSYS Fluent and verify them experimentally using neutral density interferometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resource Assessment for Distributed Wind Energy: An Evaluation of Best-Practice Methods in the Continental US

Current wind resources within the United States (US) indicate a potential to profitably install nearly 1,400 gigawatts of distributed wind (DW) capacity. This amount is equivalent to over half of the United States’ current energy demand from electricity, making it enough to power millions of homes and businesses and replace countless fossil fuel-based generating plants. Despite the potential growth of DW in the US, deployments are presently hindered by a lack of confidence in resource estimation methods. One potential challenge is that smaller-scale turbines, with hub heights of 40 meters or less, are disproportionately impacted by obstacles such as buildings and vegetation. These obstacles may produce complex wake effects, best modeled with high-fidelity complex fluid dynamics (CFD) models that are too computationally expensive to use for routine siting and resource assessment. Thus, installers today make use of heuristics and simple equations to approximate the impact of obstacles while also leveraging long-term resource data from commercial or publicly available atmospheric models. This study evaluates these historical and commonly used methods alongside new lower-order obstacle models produced from CFD simulations and measurement-based bias correction. The preliminary results from this study show the importance of taking care in the choice and application of mesoscale atmospheric models and the significant value of bias correction using measurements from nearby meteorological towers. Detailed obstacle modeling provides only modest additional gains in performance and, in some cases, can add error, especially at sites where turbines have already been located to avoid obvious impact from upwind obstacles. These findings reinforce the importance of collecting in situ measurements and suggest that obstacle models may be better applied in practice to automated or computer-aided siting, rather than in economic wind resource assessments.

17 WIND ENERGY↗

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗

Probing the Effects of NO x and SO x Impurities on Oxy-Fuel Combustion in Supercritical CO 2 : Shock Tube Experiments and Chemical Kinetic Modeling

The direct-fired supercritical carbon dioxide cycles are one of the most promising power generation methods in terms of their efficiency and environmental friendliness. Two important challenges in implementing these cycles are the high pressure (300 bar) and high CO 2 dilution (>80%) in the combustor. The design and development of supercritical oxy-combustors for natural gas require accurate reaction kinetic models to predict the combustion outcomes. The presence of a small amount of impurities in natural gas and other feed streams to oxy-combustors makes these predictions even more complex. During oxy-combustion, trace amounts of nitrogen present in the oxidizer is converted to NO x and gets into the combustion chamber along with the recirculated CO 2 . Similarly, natural gas can contain a trace amount of ammonia and sulfurous impurities that get converted to NO x and SO x and get back into the combustion chamber with recirculated CO 2 . In this work, a reaction model is developed for predicting the effect of impurities such as NO x and SO x on supercritical methane combustion. The base mechanism used in this work is GRI Mech 3.0. H 2 S combustion chemistry is obtained from Bongartz et al. while NOx chemistry is from Konnov. The reaction model is then optimized for a pressure range of 30–300 bar using high-pressure shock tube data from the literature. It is then validated with data obtained from the literature for methane combustion, H 2 S oxidation, and NO x effects on ignition delay. The effect of impurities on CH 4 combustion up to 16 atm is validated using NO x -doped methane studies obtained from the literature. In order to validate the model for high-pressure conditions, experiments are conducted at the UCF shock tube facility using natural gas identical mixtures with N 2 O as an impurity at ~100 bar. Current results show that there is a significant change in ignition delay with the presence of impurities. A comparison is made with experimental data using the developed model and predictions are found to be in good agreement. The model developed was used to study the effect of impurities on CO formation from sCO 2 combustors. It was found that NO x helps in reducing CO formation while the presence of H 2 S results in the formation of more CO. The reaction mechanism developed herein can also be used as a base mechanism to develop reduced mechanisms for use in CFD simulations.

Energy & Fuels↗

Evaluating the Effect of Vane Trailing Edge Flow on Turbine Rim Sealing

Abstract Modern gas turbine development continues to move toward increased overall efficiency, driven in part by higher firing temperatures that point to a need for more cooling air to prevent catastrophic component failure. However, using additional cooling flow bled from the upstream compressor causes a corresponding detriment to overall efficiency. A primary candidate for cooling flow optimization is purge flow, which contributes to sealing the stator–rotor cavity and prevents ingestion of hot main gas path (MGP) flow into the wheelspace. Previous research has identified that the external main gas path flow physics play a significant role in driving rim seal ingestion. However, the potential impact of other cooling flow features on ingestion behavior, such as vane trailing edge (VTE) flow, is absent in the open literature. This paper presents experimental measurements of rim cavity cooling effectiveness collected from a one-stage turbine operating at engine-representative Reynolds and Mach numbers. Carbon dioxide (CO2) was used as a tracer gas in both the purge flow and vane trailing edge flow to investigate flow migration into and out of the wheelspace. Results show that the vane trailing edge flow does in fact migrate into the rim seal and that there is a superposition relationship between individual cooling flow contributions. Computational fluid dynamics (CFD) simulations using unsteady Reynolds-averaged Navier–Stokes (URANS) were used to confirm VTE flow ingestion into the rim seal cavity. Radial and circumferential traverse surveys were performed to quantify cooling flow radial migration through the main gas path with and without vane trailing edge flow. The surveys confirmed that vane trailing edge flow is entrained into the wheelspace as purge flow is reduced. Local CO2 measurements also confirmed the presence of VTE flow deep in the wheelspace cavity.

Engineering↗