Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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↗
DOE OSTI2023

A PIC Bootstrapping Strategy for Exascale CFD-DEM Simulation

The exascale computing era comes with the release of MFIX-Exa, a new code for studying gas-particle fluidization using CFD-DEM and PIC models on massively parallel, heterogeneous high-performance computing architectures. However, as compute resources continue to grow in scale and complexity, so too does the impetus to use them efficiently. Here, we propose a novel bootstrapping method to minimize neglected simulation time and cast the approach more broadly among a growing field of multi-fidelity methods.

Porcu, Roberto↗
DOE OSTI2021

Benchmarking a Pronghorn Model of the Natural Convection Shutdown Heat Removal Test Facility Against STAR-CCM+ (Poster)

Summary Verification of Pronghorn’s Net Radiation Transfer via Code-to-Code comparison with STAR-CCM+. Comparison based on Natural Convection Shutdown Heat-Removal Test Facility (NSTF). Near-perfect agreement if conjugate heat transfer inside the riser channels is replaced by fixed temperature. Difference of up to 5K if conjugate heat transfer is modeled (correlation vs. full CFD).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Method for designing a combustion system with reduced environmentally-harmful emissions

A method for designing a combustion system which emits less of at least one environmentally-harmful emission is presented. In a describing step, an injector which introduces a fuel into a combustion chamber is described via a CFD code. In a modeling step, combustion kinetics of the fuel are modeled via a pre-processing code as the fuel mixes and reacts with an oxidizer. In a first selecting step, at least one primary scalar is derived during the modeling of the combustion kinetics. In a performing step, a table look-up is performed to obtain at least one data from a look-up database based on the primary scalar. In a second selecting step, at least one secondary scalar is selected in addition to the primary scalar(s). In a specifying step, at least one chemical pathway of formation or destruction for the secondary scalar is specified via a chemistry manager wherein the secondary scalar is representative of the environmentally-harmful emission(s) of the chemical pathway(s). In a utilizing step, the data is utilized to evaluate the chemical pathway(s) to quantify the environmentally-harmful emission(s). In an identifying step, an improvement to the combustion system is identified which reduces the environmentally-harmful emission(s).

Zambon, Andrea C.↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
DOE OSTI2024

Implementation of Detailed Polyethylene Pyrolysis Kinetics into CFD Simulations using Machine Learning

Municipal solid waste (MSW) and waste plastics have received significant attention due to the issues of waste generation and storage, as well as their potential as an energy resource. High-density polyethylene (HDPE) makes up a large portion of plastic waste and has been the subject of several conversion studies. However, the mechanisms associated with converting HDPE through pyrolysis and gasification are extensive and complex making them difficult to implement into high-fidelity computational fluid dynamic (CFD) simulations. For this project, a primary pyrolysis mechanism containing 42 unique species and 737 heterogeneous reactions was used to generate kinetic data over a range of operating conditions. A machine learning (ML) model was developed to replicate the results of the detailed pyrolysis mechanism while significantly increasing the computational efficiency. A deep operator network (DeepONet) architecture was adopted to train the model using time steps relevant to CFD simulations. The ML used physics-based loss functions to ensure mass conservation. The ML model has been deployed in simple MFiX CFD simulations, single particle, and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗
DOE OSTI2022

Validation of NEAMS Tools Using MSRE Data

This report documents the benchmarks developed and performed for the tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program utilizing the MSRE experiment data. There are three focus areas explored in this work: fuel evolution, CFD, and system transient analyses. Fuel depletion is important for MSR reactor physics and fuel cycle analysis. Fuel depletion capability is currently being developed in the Griffin reactor physics code. To support this, a benchmark problem is developed using the MSRE data and an existing depletion code with flowing fuel analysis functionality (ADDER) to validate this new capability in Griffin. The ADDER calculation is compared with the TRITON-MSR results as well as the MSRE experimental data. This benchmark can be directly used to verify and validate the Griffin fuel depletion capability in future years. For CFD simulation, the MSRE CAD model is obtained from the Copenhagen Atomics and a mesh is generated for the MSRE inlet plenum. The Nek5000 code will be employed to find the flow distribution in the MSRE inlet plenum based on the mesh generated. For system analysis, a simple benchmark problem for three types of MSR system transients is developed. The SAM and SPECTRA codes are used to solve the simple problem, and the results from the two codes agree well. The SAM code is also employed to obtain the frequency response of the MSRE system at the 5 and 8 MWt power levels. The results agree well with the experiment measurements. These verification and validation efforts gives SAM users more confidence in MSR system transient simulations using SAM. Lastly, the SAM-Griffin Multiphysics coupling work is initiated. A simple one-way coupling of evaluating the precursor drift effect on the MSRE k-eff is demonstrated. This preliminary coupling work serves as the basis for future more complex Multiphysics simulation and modelling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

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

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal Modelling of Advanced Test Reactor Fuel in a Generalized Dry Storage System

Star-CCM+, a computational fluid dynamics (CFD) software was used to conduct modeling and simulation of the thermal performance of a dry storage configuration consisting of Department of Energy Standardized Canisters (DOESCs) loaded with aluminum-clad spent nuclear fuel (ASNF.) The configuration includes nine DOESCs loaded with Advanced Test Reactor ASNF contained within a stainless steel overcanister centered in a ventilated, concrete overpack. The simulations were used to estimate the maximum temperatures reached by backfill gases inside the overcanister and DOESCs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2020

Evaporation Submodel Development for Volume of Fluid (eVOF) Method Applicable to Spray-Wall Interaction Including Film Characteristics with Validation at High Pressure and Temperature Conditions

Internal combustion engines have seen a great evolution over the last several decades through application of high pressure direct injection, multiple injections, and other technologies to reduced fuel consumption, NOx, and PM. Although combustion systems with advanced injection strategies have been studied extensively, there exists a significant fundamental knowledge gap on the fuel-spray interactions with the piston surface and chamber walls. Advanced computational codes validated with experimental techniques have to be developed for accurate representation of the drop impingement, fuel film formation, and vaporization. Current engine CFD (Computational Fluid Dynamics) spray models utilize a Lagrangian framework for modeling which lacks critical considerations of the physics pertaining to these interactions and thus requiring extensive parameterization, tuning and validation. The team from Michigan Technological University, University of Massachusetts Dartmouth, and Argonne National Laboratory is composed of experts in sprays, combustion, engines and CFD with a wide spectrum of knowledge including specific expertise in the area under consideration. In the proposed work, a VOF (Volume of Fluid) modeling approach has been adopted for the spray-wall interaction, film formation and spreading, and vaporization. With the inclusion of a vaporization submodel, a more predictive and accurate simulation of the spray-film was performed without extensive need of parameterization and tuning. Extensive experimentation of the spray-wall interaction under the range of conditions matching the thermodynamic and surface temperatures that occur in diesel and gasoline engines were conducted to validate the SWI submodels and for development of the evaporation submodel, which has been implemented in the flow solver.

42 ENGINEERING↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Uncertainty Assessment of CFD Investigation of the Nonlinear Difference-Frequency Wave Loads on a Semisubmersible FOWT Platform

Current mid-fidelity modeling approaches for floating offshore wind turbines (FOWTs) have been found to underpredict the nonlinear, low-frequency wave excitation and the response of semisubmersible FOWTs. To examine the cause of this underprediction, the OC6 project is using computational fluid dynamics (CFD) tools to investigate the wave loads on the OC5-DeepCwind semisubmersible, with a focus on the nonlinear difference-frequency excitation. This paper focuses on assessing the uncertainty of the CFD predictions from simulations of the semisubmersible in a fixed condition under bichromatic wave loading and on establishing confidence in the results for use in improving mid-fidelity models. The uncertainty for the nonlinear wave excitation is found to be acceptable but larger than that for the wave-frequency excitation, with the spatial discretization error being the dominant contributor. Further, unwanted free waves at the difference frequency have been identified in the CFD solution. A wave-splitting and wave load-correction procedure are presented to remove the contamination from the free waves in the results. A preliminary comparison to second-order potential-flow theory shows that the CFD model predicted significantly higher difference-frequency wave excitations, especially in surge, suggesting that the CFD results can be used to better calibrate the mid-fidelity tools.

17 WIND ENERGY↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
DOE OSTI2025

Design Overview of a High-Pressure Helium Flow Visualization Apparatus for Blanket Cooling Studies

Cooling of the fusion blanket first wall remains a significant challenge given the adverse conditions of heat and particle flux encountered near the plasma. Helium emerges as an attractive cooling candidate because of its chemical and neutronic inertness and separability from hydrogenic species (e.g. tritium). Because of the low thermal mass of helium, optimization of these coolant channels is warranted to provide high heat transfer performance at low pumping costs. Increasingly, computational fluid dynamics (CFD) simulations are employed to model and optimize these flow channels, and accompanying experimental data are needed to validate the predictions of these models. To provide the aforementioned experimental data, a high-pressure helium flow visualization upgrade has been designed for the Helium Flow Loop Experiment facility. This apparatus was built to American Society of Mechanical Engineers boiler and pressure vessel standards to withstand operating pressure of 4 MPa and mated to high-pressure glass windows. Seedless flow visualization is performed via high-speed background oriented schlieren (BOS), with image correlation used for time-resolved two-dimensional velocimetry at frequencies in excess of 60 kHz. Rectangular flow channel test articles are additively manufactured via laser powder bed fusion and installed into this visualization apparatus, with one-sided heating supplied by resistive heaters. In conclusion, the chosen test geometries were informed by prior CFD simulations, and the helium flow structures observed via BOS (detachment, recirculation, etc.) will be used for the validation of these accompanying models, in support of the design and optimization of blanket cooling channel configurations.

Helium flow↗