Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamic (CFD)”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗
DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2022

Pressure Drop in Seven-Pin Wire-Wrapped Rod Bundle for the Sodium Cartridge Loop in Versatile Test Reactor

This work studies the hydrodynamics of the seven pin wire-wrapped rod bundle in the sodium cartridge loop for the Versatile Test Reactor (VTR) through scaled water experiments and computational fluid dynamics (CFD) simulations. Here, the scaling analysis is first performed to demonstrate the hydrodynamic similarity between water and sodium flows at the same Reynolds number Re. A separate-effects test facility is designed and constructed based on the scaling analysis. Detailed experimental data on the pressure drop covering a wide range of Re values (1165 to 27 689) are obtained, which are used to evaluate existing correlations for friction factor and to benchmark CFD simulations. The experimentally determined friction factors f Exp agree well with the Upgraded Cheng and Todreas Detailed Correlation and Pacio-Chen-Todreas Detailed Model within ±7% but are significantly underpredicted by Rehme's correlation by 25%. Various CFD near-wall treatment methods are tested using ANSYS Fluent and evaluated by experimental data. It is found that when the recommended wall y + values are met, most of the near-wall treatment methods can give accurate friction factor predictions. The resolved near-wall method (y + ~ 1) with the Shear Stress Transport k - ω turbulence model and the scalable wall functions (y + > 12) with the realizable k - ε turbulence model can predict f Exp within ±5% The standard wall functions (y + > 30) and nonequilibrium wall functions (y + > 30) with the realizable k - ε model can predict f Exp within ±10%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Modelling Airborne Transmission and Ventilation Impacts of a COVID-19 Outbreak in a Restaurant in Guangzhou, China

Computational fluid dynamics (CFD) modelling was performed to simulate spatial and temporal airborne pathogen concentrations during an observed COVID-19 outbreak in a restaurant in Guangzhou, China. The reported seating configuration, overlap durations, room ventilation, layout, and dimensions were modelled in the CFD simulations to determine relative exposures and probabilities of infection. Results showed that the trends in the simulated probabilities of infection were consistent with the observed rates of infection at each of the tables surrounding the index patient. Alternative configurations that investigated different boundary conditions and ventilation conditions were also simulated. Increasing the fresh-air percentage to 10%, 50%, and 100% of the supply air reduced the accumulated pathogen mass in the room by an average of ~30%, ~70%, and ~80%, respectively, over 73 min. Overall, the probability of infection was reduced by ~10%, 40%, and 50%, respectively.

60 APPLIED LIFE SCIENCES↗
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↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Meta-uncertainty for particle image velocimetry

Uncertainty quantification for Particle Image Velocimetry (PIV) is critical for comparing experimentally measured flow fields with Computational Fluid Dynamics (CFD) results, and model design and validation. However, PIV features a complex measurement chain with coupled, non-linear error sources, and quantifying the uncertainty is challenging. Multiple assessments show that none of the current methods can reliably measure the actual uncertainty across a wide range of experiments, and estimates can vary. Because the current methods differ in assumptions regarding the measurement process and calculation procedures, it is not clear which method is best to use for an experiment where the error distribution is unknown. To address this issue, we propose a method to estimate an uncertainty method's sensitivity and reliability, termed the Meta-Uncertainty. The novel approach is automated, local, and instantaneous, and based on perturbation of the recorded particle images. We developed an image perturbation scheme based on adding random unmatched particles to the interrogation window pair considering the signal-to-noise (SNR) of the correlation plane. Each uncertainty scheme's response to several trials of random particle addition is used to estimate a reliability metric, defined as the rate of change of the inter-quartile range (IQR) of the uncertainties with increasing levels of particle addition. We also propose applying the meta-uncertainty as a weighting metric to combine uncertainty estimates from individual schemes, based on ideas from the consensus forecasting literature. We use planar and stereo PIV measurements across a range of canonical flows to assess the performance of the uncertainty schemes. Further, a novel method is introduced to assess an uncertainty scheme's performance based on a quantile comparison of the error and uncertainty distributions, generalizing the current method of comparing the RMS of the two distributions. Here, the results show that the combined uncertainty method outperforms the individual methods, and this work establishes the meta-uncertainty as a useful reliability assessment tool for PIV uncertainty quantification.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Measurement of a radial flow profile with eddy current flow meters and deep neural networks

Eddy current flow meters (ECFMs) measure flows of conductive fluids. Recent interest in ECFMs has increased due to applications in advanced nuclear reactors. ECFMs are well suited for such applications, as they can provide non-invasive measurements of flow in fluids that are often difficult to measure. Traditionally, ECFMs are operated using an alternating current at a single frequency, limiting ECFMs to measure average fluid velocities, blockages, or voids. Here, we expand the capabilities of ECFMs by measuring the fluid radial velocity profile of liquid mercury. To accomplish this, we made several ECFM sensitivity measurements at a range of frequencies. Different frequencies vary the electromagnetic skin depth of the device. By adjusting frequencies, we probed the fluid velocity at various radial locations and constructed a flow-velocity profile. The relationship between the ECFM measurements and velocity profile is nonlinear and requires solving an inverse problem. Using electromagnetic finite-element simulations to train a deep neural network (DNN), we created a model that provides a stable general relationship between the sensitivity measurements of an ECFM and the fluid velocity profile. Using ECFM measurements of liquid mercury, our DNN model calculates a flow profile that agrees well with computational fluid dynamics (CFD) simulations. This technique has potential to improve flow monitoring for optimization, safe operation of conductive fluid loops, and/or validating complex CFD models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Impact of high-fidelity temperature feedback modeling on fusion blanket neutronics

Accurate modeling and simulation of neutron and gamma transport and interactions within fusion tokamaks represents a critical area of research for fusion energy development. This paper presents a method to couple computational fluid dynamics (CFD) results to high-fidelity neutronics analyses, to explore the fusion neutronics’ feedback to blanket temperature and density. This method is demonstrated on a sector model of a tokamak with immersion molten salt blanket, and the relative differences in tritium breeding ratio (TBR) with spatially-resolved temperature modeling are explored. Additionally, the influence of higher fidelity neutronics results on the subsequent CFD analysis is explored. Results show that global TBR changes $<$0.1% due to temperature spatial distribution modeling. However, local TBR in areas with high temperature variability can have up to ∼4% difference.

blanket fusion↗
DOE OSTI2020

Enhancement of Unsteady and 3D Aerodynamics Models using Machine Learning

Unsteady aerodynamics will be an important part of the floating wind turbines of the future operating under high shear across the rotor disk coupled with platform motion and atmospheric turbulence. We develop an unsteady aerodynamics and dynamic stall model using a long short-term memory variant of recurrent neural networks. The neural network model is trained using the oscillating airfoil data set from Ohio State University. The predictions from our machine learning (ML)-based model show good agreement with the experimental data and other state-of-the-art dynamic stall models for a wide range of airfoils, Reynolds numbers and reduced frequencies. In some cases the predictions are better than the Beddoes-Leishman model implementation in OpenFAST, when using the default coefficients. The ML-based model is also able to capture the key physics associated with dynamic stall, such as the precedence of moment stall before lift stall and cycle-to-cycle variations in the aerodynamic response. The new unsteady aerodynamics model is expected to improve prediction of fatigue loads for yaw-based wake-steering control scenarios in actuator-line and actuator-disc simulations of wind farms. Our methodology for training the ML-model provides a pathway for improving design level tools using high-fidelity computational fluid dynamics (CFD) simulations in the future.

17 WIND ENERGY↗
DOE OSTI2025

Cryogenic Design and Thermal Analysis of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

The Electron Ion Collider (EIC) physics program utilizes a 2.0 T superconducting magnet at the heart of its ePIC detector system. This approximately 3.5 m long and 2.84 m diameter warm bore magnet has a 20 tons cold mass which is conduction-cooled using liquid helium at 4.5 K. A closed loop active thermosiphon system is chosen to facilitate the cooling and to maintain a minimum of 2 K temperature margin on the peak operating temperature (4.7 K) of the superconductor. Here, this paper presents the cryogenic design of the cooling system and the thermal analysis of the solenoid. A Computational fluid dynamics (CFD) model was developed to quantify the performance of the two–phase closed thermosiphon system and predict the temperature gradient on the cold mass.

Gopinath, Sandesh [Thomas Jefferson National Accel↗
DOE OSTI2022

MPEX High Heat Flux Plasma Dump Design, Manufacture, and Articles Test

Here, the Material Plasma Exposure eXperiment (MPEX) device is a steady-state linear plasma device currently in the final design phase at the Oak Ridge National Laboratory. This device will reach ion fluences up to 10 31 m -2 and will be used to expose neutron-irradiated materials to divertor-relevant plasmas and to study the effects of plasma-material interactions. These studies will elucidate the complex effects of plasmas with divertor candidate materials capable of withstanding high heat flux and high fluences for next-generation fusion devices. Bidirectional plasma will be generated using a high-power (200 kW) helicon source. Plasma will be confined by superconducting magnets. The last plasma-facing component on the upstream side of the MPEX device is the dump, which has been designed to intercept plasma and energetic particles. The dump will have a total heat load of 9.2 kW. A copper alloy (Glidcop AL-15) was selected for use in the water-cooled flange design because of its high thermal conductivity, its retaining strength at elevated temperatures, and its ability to be used in the high-temperature braze joints used in this application. Titanium-zirconium-molybdenum (TZM) tiles are brazed to the Glidcop AL-15 flange using high-temperature braze alloy. External water-cooling channels are used on the dump flange to prevent water leakage inside the vacuum space. This article discusses the details of the high heat flux dump design, including the computational fluid dynamics (CFD) and structural analyses performed to validate the design to meet the operational requirements of the MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Design and Analysis of Liquid Lithium Plasma Facing Components

In this study, liquid metal (LM) on plasma-facing components (PFCs) creates a renewable protective cover enhancing power exhaust and protecting the solid substrate. Additionally, energy confinement improvement through particle pumping can be achieved in the case of liquid lithium (LL). A new variant of PFCs was recently introduced at Princeton Plasma Physics Laboratory (PPPL) where a porous wall is used to stabilize the LM surface, while using magnetohydrodynamic (MHD) drive to push the LM flow inside the component. This arrangement allows efficient heat exhaust, and enhanced control of the LM surface temperature, leading to spatial control of evaporation and condensation of LL on the plasma interface. This feature is particularly attractive when vapor shielding is introduced to allow heat flux redistribution. This system has the advantage that as the heat flux increases, the evaporation rate will increase, while decreasing heat flux will decrease the evaporation rate of lithium, ideally creating a feedback effect which could self-regulate the amount of lithium evaporated. Analytical and numerical models for LL PFCs were developed in PPPL. To calculate the target temperature distribution for the case of evaporation from the divertor, we apply an iterative process allowing two-way coupling between the fluid-kinetic analysis of plasma using SOLPS-ITER code and the flow and heat transfer analysis of the PFC using an analytical model. In the last stage, the results are validated using computational fluid dynamics (CFD) analysis with customized version of the CFX code from ANSYS. CFX was modified at PPPL to allow MHD analysis at the high magnetic field typical for fusion applications. Results of the analysis for NSTX-U tokamak conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗
DOE OSTI2020

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

Engine Combustion System Optimization Using Computational Fluid Dynamics and Machine Learning: A Methodological Approach

Gasoline compression ignition (GCI) engines are considered an attractive alternative to traditional spark-ignition and diesel engines. Here, a Machine Learning-Grid Gradient Ascent (ML-GGA) approach was developed to optimize the performance of internal combustion engines. ML offers a pathway to transform complex physical processes that occur in a combustion engine into compact informational processes. The developed ML-GGA model was compared with a recently developed Machine Learning-Genetic Algorithm (ML-GA). Detailed investigations of optimization solver parameters and variable limit extension were performed in the present ML-GGA model to improve the accuracy and robustness of the optimization process. Detailed descriptions of the different procedures, optimization tools, and criteria that must be followed for a successful output are provided here. The developed ML-GGA approach was used to optimize the operating conditions (case 1) and the piston bowl design (case 2) of a heavy-duty diesel engine running on a gasoline fuel with a research octane number (RON) of 80. The ML-GGA approach yielded >2% improvements in the merit function, compared with the optimum obtained from a thorough computational fluid dynamics (CFD) guided system optimization. The predictions from the ML-GGA approach were validated with engine CFD simulations. This study demonstrates the potential of ML-GGA to significantly reduce the time needed for optimization problems, without loss in accuracy compared with traditional approaches.

33 ADVANCED PROPULSION SYSTEMS↗