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 433 records · Page 24

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

NUMERICAL PREDICTIONS OF MEAN PERFORMANCE AND DYNAMIC BEHAVIOR OF A 10 MWe SCO2 COMPRESSOR WITH TEST DATA VALIDATION

High-fidelity aerodynamic analysis has been demonstrated for mean performance and unsteady dynamics in a sCO2 compressor designed by Hanwha Power Systems Americas for a 10 MWe Concentrating Solar Power plant. Simulations were performed with CRUNCH CFD® software tool that was matured to accurately model near critical real fluid effects in sCO2. Pre-test predictions for mean performance were validated with test data collected later. Performance predictions were accurate and captured sensitivity of the efficiency to inlet temperature of CO2 as well as steep drop-off at high flow rates due to condensation in the inlet throat. Detailed analysis was performed to understand the source of these performance losses at near critical conditions. Unsteady dynamic effects in the compressor at off-design conditions were also identified and quantified. In particular a system wide “condensation surge” condition was detected at high flow coefficients that results in large amplitude pulsations with accompanying mass flow fluctuations at low frequencies and has potential to cause damage in closed loop systems.

hosangadi, ashvin↗

Flow Induced Corrosion in 6061-T6 Aluminum Pipes in De-Ionized Water Environments

Flow induced corrosion/erosion of 6061-T6 aluminum in de-ionized (DI) water environments has not been studied widely. Especially, the long-term effects of corrosion/erosion in seemingly benign flow velocity, temperature, and resistivity ranges of 8 ft/s, 85 oF, 3-5 MOhm-cm, respectively. This study concludes that the flow induced corrosion/erosion in the above parameter ranges is minimal. This is detailed by presenting a literature survey, measuring pipe wall samples from a system that has operated in the above parameter range, predict the loss of material in mm/year at the above velocities and at higher temperatures using electric potential values from other experimental studies coupled with wall shear stress simulated using Computational Fluid Dynamic (CFD) analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preliminary Modeling of Triply Periodic Minimal Surface (TPMS) Structures with RELAP5-3D

This poster goes over preliminary work performed on modeling triply periodic minimal surfaces (TPMS) with RELAP5-3D. The poster starts off with the need for advanced nuclear fuel designs for future reactors, and how TPMS are being researched for this purpose. There is then a focus on what a TPMS is and then moves onto the motivation for why this geometry is being modeled and looked into. The methodology is then shown for how data from different experiments will be used for pre-test predictions and eventual creation of empirical coefficients and multipliers to be used in models. Then there is a show of how data would be compared between the experimental tests, CFD, and RELAP5-3D values. The results from the simulations and comparisons are then presented to show how the CFD and RELAP5-3D compare to each other currently, with no experimental data to compare to for one test model. Results from the second test model are also presented to show that there is stability in the solving of that model.

42 - ENGINEERING↗

Effect of the J-Shaped Wind Turbine Airfoil Opening Ratio and Thickness on the Performance of Symmetrical Airfoils

Abstract In this work, the effect of the inner opening ratio on the J-shaped airfoils aerodynamic performance was studied and documented for symmetrical airfoils. Three different airfoil thicknesses were investigated: small (NACA0008), medium (NACA0015), and large (NACA0024). For each airfoil thickness, effects of three inner opening ratios were analyzed: one-third, one-half, and two-thirds. The performance of each opening ratio was compared with the performance of the solid airfoil “zero opening ratio” for different angles of attack between 5 deg and 20 deg. All designs were simulated using the computational fluid dynamics (CFD) technology against experimental results for solid NACA4412 airfoil in the University of Wisconsin-Milwaukee (UWM) wind tunnel facility and other published experimental data. It was found that large eddy simulation yields accurate solutions with a smaller number of mesh cells compared to the k–ω turbulence model but with much longer computational time. The lift-to-drag ratio for all studied airfoils has a maximum value for solid airfoils compared to those equipped with openings. For airfoils equipped with 00.00% opening ratio “solid,” NACA0015 airfoil has the maximum lift-to-drag ratio. Furthermore, it was found that NACA0008 equipped with a 33.33% opening ratio has the best performance of all studied J-shaped airfoils.

Energy & Fuels↗

Development of a Filtered CFD-DEM Drag Model with Multiscale Markers Using an Artificial Neural Network and Nonlinear Regression

Here, the accuracy of coarse-grained Euler-Lagrangian simulations of fluidized beds heavily depends on the mesoscale drag models to account for the influences of the unresolved sub-grid structures. Traditional filtered drag models are regressed with mesoscale markers such as voidage and slip velocities. In this research, a filtered drag was regressed with both mesoscale and macro-scale markers using fine grid Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) simulations. The traditional non-linear regression method was compared with machine learning regression using an Artificial Neural Network (ANN) implemented in PyTorch and coupled with MFiX. The new drag showed higher accuracy than the Wen-Yu drag and another filtered drag derived from the two-fluid model. The nonlinear regression shows slightly better results than ANN regression in cases with similar R 2 values. The utilization of the gas inlet velocity as an additional macro-scale marker reduced the errors by up to 55.3% in the tested cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

ExaWind: A multifidelity modeling and simulation environment for wind energy

We introduce the open-source ExaWind modeling and simulation environment for wind energy. The primary physics codes of ExaWind are Nalu-Wind and OpenFAST. Nalu-Wind is a wind-focused computational fluid dynamics (CFD) code that is coupled to the whole-turbine simulation code OpenFAST. The ExaWind environment was created under U.S. Department of Energy funding to achieve the highest-fidelity simulations of wind turbines and wind farms to date, with the goal of enabling disruptive changes to turbine and plant design and operation. Innovation will be gleaned through better understanding of the complex flow dynamics in wind farms, including wake evolution and the impact of wakes on downstream turbines and turbulent flow from complex terrain. High-fidelity predictive simulations employ hybrid turbulence models, geometry/boundary-layer-resolving CFD meshes, atmospheric turbulence, nonlinear structural dynamics, and fluid-structure interaction. While there is an emphasis on very high-fidelity simulations (e.g., blade resolved with full fluid-structure coupling), the ExaWind environment supports lower-fidelity modeling capabilities including actuator-line and -disk methods. Important in the development of ExaWind codes is that the codes scale well on today's largest petascale supercomputers and on the next-generation platforms that will enable exascale computing.

17 WIND ENERGY↗

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗

M-Star ® Software Test and Verification

Savannah River Mission Completion (SRMC) currently manages the risk for retained hydrogen in the Defense Waste Processing Facility (DWPF) vessels by implementing a Retained Hydrogen Program. The current program relies on Sludge Batch (SB) 8 Gas Chromatograph data and on conservative assumptions concerning gas release and retention to determine allowable vessel Quiescent time (Q-time). The authors identified M-Star ® CFD as a software that could simulate processes such as fluid flow, heat transfer, species transport, chemical reactions, particle transport, and retained hydrogen gas release. Preliminary simulation results suggest that more realistic assumptions on gas retention and release may be feasible for the DWPF retained hydrogen program. Because of the desire to use the M-Star ® software to perform analyses that support nuclear safety, SRMC has requested Savannah River National Laboratory (SRNL) to upgrade the software classification level of M-Star ® CFD from level D to level A to perform analyses that support nuclear safety.

08 HYDROGEN↗

Evaluation of adsorption and mechanical strength of 13X zeolite mixtures with phyllosilicate binders using molecular dynamics simulation and positron annihilation spectroscopy

There is growing interest in developing zeolites with novel internal structures that have optimal adsorptive capacity and high mechanical strength, while offering advantages, such as being light weight. We integrate computational and experimental methods to explore the effect of binder/zeolite types, and weight percentages on the mechanical strength of 13X zeolite and adsorption capacities of N 2 , H 2 O, and CO 2 for additive manufacturing (AM) applications with the goal of maximizing both adsorption and strength. Zeolite 13X mixtures and phyllosilicate binders (either bentonite or kaolin) are combined using molecular dynamics (MD) simulations to create structures with various binder/zeolite weight percentages. Further, adsorption capabilities and mechanical strength are assessed using the grand canonical Monte Carlo (GCMC) and ReaxFF modules, respectively. Our modeling shows that an optimized zeolite/binder ratio for N 2 adsorption is around 15 wt% for kaolin and roughly 10 wt% for bentonite. The resulting parameters can be applied to facilitate macro-scale computational fluid dynamics (CFD) and finite element method (FEM) simulations of an AM zeolite structure. We also performed Positron Annihilation Lifetime Spectroscopy (PALS) measurements on zeolite samples to explore the effect of changes in the internal volume. The results show an inverse relationship between the free volume and the solid loading and adding a binder changes the morphology of the zeolite-binder compound and decreases open-volume area significantly.

36 MATERIALS SCIENCE↗

A hybrid porous model for full reactor core scale CFD investigation of a prismatic HTGR

Three-dimensional (3-D) Computational Fluid Dynamics (CFD) analysis of a whole nuclear reactor core is a tremendous challenge due to the large geometric volume and complex structures. Here, this research presents a hybrid porous (HP) model to simplify a prismatic High Temperature Gas-cooled Reactor (HTGR) core, so 3-D CFD investigation can be performed on a full reactor core scale. In the HP model, the prototypic small coolant channels in the nuclear fuel blocks are lumped together to form multiple equivalent large coolant channels, and then the porous medium flow model is applied to each of them. Therefore, heat transfer in fuel blocks is computed by a hybrid combination of solid energy and porous flow energy equations. The similarity between the HP model and prototypic model is achieved by deriving the porous flow permeability, inertial resistance factor, and artificial thermophysical properties. Compared with the widely used whole porous (WP) flow model, the HP model preserves more realistic geometric structures, and therefore more accurate physical processes. The General Atomics' Modular High Temperature Gas-cooled Reactor (MHTGR) design was chosen as a prototype to demonstrate the methodology. Simulations were performed using the prototypic CFD model and HP model at steady-state forced circulation, steady-state natural circulation, and transient conditions that correspond to normal operation, extended period of pressurized cool down, and short-term transients after reactor shutdown, respectively. The comparison shows good agreement between the HP model and prototypic model in the maximum fuel temperature, average solid temperature, and helium flow rate, which demonstrates the potential applicability of the HP model for a full reactor core scale simulation in the future. As a benefit, the HP model reduces the mesh quantity by a factor of 50 from a prototypic model. Correspondingly, the computation time was reduced by a factor of at least 30.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tank 27 Solids Disturbance Testing Using M-Star for the Impact of CSMP Usage on Blend Tank Mixing Effectiveness

Savannah River Mission Completion (SRMC) plans to operate Tank 27as a blend tank for future waste disposal operations. It is desirable to avoid additional worker exposure by using the planned commercial submersible mixing pump (CSMP) installations in Risers B2 and B4 to function as the blend pumps for prolonged blend tank operations. The purpose of this task is to evaluate the impacts to mixing across a range of pump installation heights from 1 to 201 inches above the tank floor. The CSMPs shall maintain adequate blending within a 24 hour maximum timeframe and an adequate nozzle discharge velocity as outlined in U-ESR-G-00030, using the specifications of the planned CSMPs listed in PO SRRA149971 and M-DS-G-00081. Specifications include: maximum pump diameter of D=22.5 in., dual nozzles placed tangentially opposed with nominal inside diameter of 2 in., total flow range (both nozzles combined) of 500 to 1950 GPM and maximum RPM of 1800. Tank 27 is expected to have a fill height maximum of 360". This analysis began by using results from the 2010 testing and applying them to the Tank 27 design to calculate a mixing time as a function of pump flow rate. These results were modified to account for the higher liquid level in Tank 27. Once the liquid level was increased to 360 inches, the results were modified to account for a higher viscosity in Tank 27 than in the testing. Once the viscosity was increased, the effect of changing pump elevation was applied to calculations of the mixing time in Tank 27. Finally, the uncertainty from the testing was included in the analysis to recommend a blend time as a function of pump flow rate and elevation. The conclusion from this study is that the particles in tank 27 will be suspended at both the 600 and 900 gpm flow rates. Only particles above 90 micron will not be suspended in significant amounts when the flow rate is 600 gpm. At 900 gpm even 100 micron particles were thoroughly suspended in the tank. From the series of simulations it has been concluded that it will not be possible to mix tank 27 at the proposed fluid flow rates and not suspend the solid particles at the bottom of the tank.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Assessment of Thermal Stratification in Versatile Test Reactor Transients

The Versatile Test Reactor (VTR) is a fast-spectrum test reactor currently being developed in the United States under the direction of the US Department of Energy. The conceptual design of the 300 MWth pool-type sodium-cooled fast reactor (SFR) has been led by the US National Laboratories in collaboration with General Electric-Hitachi and Bechtel National Inc. Safety performance analysis for the VTR conceptual design is being performed with the systems thermal-hydraulics (SYSTH) module of the SAS4A/SASSYS-1 liquid-metal reactor safety analysis code system. Since the current model of the VTR employs a simple perfect mixing model for large plena like the hot pool, it is not able to predict temperature variations that may develop during the transient. Prior work simulating the response of SFRs to postulated events like the Protected Station Blackout (PSBO) has shown that the phenomenon of thermal stratification, where stable thermal layers accumulate in the hot pool, may delay the transition to natural circulation and thus impact the predicted transient progression. Thus, an effort has begun to model this transient by integrating a Computational Fluid Dynamics (CFD) model of the hot pool into the SAS4A/SASSYS-1 model of the Primary Heat Transport System during the simulation of the PSBO event. A three-dimensional Volume-Of-Fluid CFD model of the VTR hot pool has been developed for the co-simulation of SAS4A/SASSYS-1 with CFD. In this work, the standalone SAS4A/SASSYS-1 calculation and the standalone CFD calculation based on the SAS4A/SASSYS-1 calculation result were produced on the Idaho National Laboratory High Performance Computing cluster, SAWTOOTH. At this time, only the standalone CFD and SAS4A/SASYS-1 simulation results are provided.

99 GENERAL AND MISCELLANEOUS↗

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗

Modeling and simulation of gas separations with spiral–wound membranes

Models for gas separations with spiral–wound membranes are developed and found to exhibit good agreement with experiments performed on N 2 /O 2 mixtures. The two–dimensional (2D) model can be accurately approximated by a one–dimensional (1D) surrogate model when the spacer widths are chosen to make the channel pressure drops small. Subsequently, the separation of propane/propylene mixtures from the recycle purge stream of a polypropylene reactor is investigated. Assuming ideal gas is found to lead to significant overestimations in membrane stage cuts (sometimes more than 10%), an extent comparable to that associated with extrapolating constant olefin permeance from a low–pressure condition. While olefin permeance can change significantly with pressure, using a constant–permeance formulation can result in a small (< 2.5%) underprediction in stage cut if the value for the permeance is taken from the feed condition. Finally, membrane properties and costs necessary for a viable separation process are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-resolution simulations of a turbulent boundary layer impacting two obstacles in tandem

High-fidelity large-eddy simulations of the flow around two rectangular obstacles are carried out at a Reynolds number of 10 000 based on the freestream velocity and the obstacle height. The incoming flow is a developed turbulent boundary layer. Mean-velocity components, turbulence fluctuations, and the terms of the turbulent-kinetic-energy budget are analyzed for three flow regimes: skimming flow, wake interference, and isolated roughness. Three regions are identified where the flow undergoes the most significant changes: the first obstacle’s wake, the region in front of the second obstacle, and the region around the second obstacle. In the skimming-flow case, turbulence activity in the cavity between the obstacles is limited and mainly occurs in a small region in front of the second obstacle. In the wake-interference case, there is a strong interaction between the freestream flow that penetrates the cavity and the wake of the first obstacle. Furthermore, this interaction results in more intense turbulent fluctuations between the obstacles. In the isolated-roughness case, the wake of the first obstacle is in good agreement with that of an isolated obstacle. Separation bubbles with strong turbulent fluctuations appear around the second obstacle.

97 MATHEMATICS AND COMPUTING↗