Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD results”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗
DOE OSTI2023

Heat transfer enhancement due to cold cap motion from bubbling in a waste glass melter

In this study, a computational fluid dynamics (CFD) model was developed to model the motion of a solid cold cap in a waste glass melter. Forced convection bubblers at the base of the melter release air into the molten glass, which forms large bubbles that travel upwards to the cold cap and augment heat transfer from the glass to the cold cap. The CFD model employs the Navier- Stokes equations to solve for the fluctuating flowfield using a rigid body motion dynamic fluid body interaction module. This allows for movement of the floating body in response to the bubbling forces calculated at each time step. Further, the heat flux delivered to the cold cap by the convective bubbling is studied as a function of the bubbling rate. Results for the moving cold cap are compared with the computed heat flux trends for a stationary cold cap. The heat flux delivered to the cold cap from the glass is 25% higher for the case with the moving cold cap. The heat flux was found to be proportional to v0.6 as opposed to v0.9 (where v is the normalized bubbling rate) for the stationary cold cap.

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

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Row spacing as a controller of solar module temperature and power output in solar farms

We report that when the temperature of solar photovoltaic modules rises, efficiency drops and module degradation accelerates. The spatial arrangement of solar modules can affect convective cooling and, consequently, module temperatures. However, the impact of row spacing on convective cooling in realistic solar farms has not yet been studied. Here, we develop six solar farm arrangements consisting of a fixed number of rows with varying streamwise row spacing. We model the flow and heat transfer of each solar farm using high-resolution large-eddy simulations. Results indicate that increasing row spacing can enhance convective cooling by 14.8%, which reduces module temperature by 6.6 °C and increases power output by 4.0% on average.

14 SOLAR ENERGY↗
DOE OSTI2023

Reducing the Flow Maldistribution in Heat Exchangers through a Novel Polymer Manifold: Numerical Evaluation

The maldistribution of working fluid is one of the issues in heat exchangers that causes a reduction in performance of not only the heat exchanger but also the entire HVAC system. One of the methods to reduce such maldistribution is to improve manifold design to evenly distribute the flow. In the present work, an advanced maldistribution reduction manifold, which was based on a preliminary maldistribution reduction manifold, was designed to further improve the flow distribution in the heat exchanger. In the design, spiral baffles are used to create vortices in the tubes to regulate the flow in each tube. The design also keeps the tubes away from the manifold inlet to avoid direct flow from the inlet. Due to the complexity, the design of the advanced maldistribution reduction manifold is for AM only, which cannot be fabricated by traditional manufacturing. To evaluate the design, a computational fluid dynamic model is developed to study flow distribution in heat exchanger manifolds. The simulation results reveal that the relative standard deviation of the tubes in the advanced maldistribution reduction design is half of the preliminary maldistribution reduction design and about 1/20 of the reference design.

42 ENGINEERING↗
DOE OSTI2025

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗
DOE OSTI2025

Glovebox Fires Phase 1: Simulations of Open Burner Experiments

Understanding the duration of a fire required to breach a glovebox via glove ports and/or windows within a facility is important in developing mitigation strategies for safety. To facilitate this understanding, both simulation and experiments are utilized. The approach is to first validate the Computation Fluid Dynamics (CFD) code, Fire Dynamics Simulator (FDS), with experimental data collected by New Mexico Tech (NMT) and then perform simulations of full-scale rooms containing gloveboxes to assess numerous scenarios which would otherwise be cost-prohibitive experimentally. This report provides comparison to first-phase experiments involving a fire without a glovebox. A subsequent report will provide comparison to experiments involving a glovebox. The results indicate good agreement with FDS tending to over predict the pre-mixed and diffusion-mode tests by 6% and 10%, respectively.

42 ENGINEERING↗
DOE OSTI2021

CFD-Guided Evaluation of Spark-Assisted Gasoline Compression Ignition for Cold Idle Operation

A closed-cycle, three-dimensional (3D) computational fluid dynamics (CFD) analysis campaign was conducted to evaluate the performance of using spark plugs to assist gasoline compression ignition (GCI) combustion during cold idle operations. A conventional spark plug using single-sided J-strap design was put at a location on the cylinder head to facilitate spray-guided spark assistance. Ignition was modeled with an L-type energy distribution to depict the breakdown and the arc-to-glow phases during the energy discharge process. Several key design parameters were investigated, including injector clocking, number of nozzle holes, spray inclusion angle, number of fuel injections, fuel split ratio, and fuel injection timings. The study emphasized the region around the spark gap, focusing on flame kernel formation and development and local equivalence ratio distribution. Flame kernel development and the ignition process were found to correlate strongly with the fuel stratification and the flow velocity near the spark gap. The analysis results showed that the flame kernel development followed the direction of the local flow field. In addition, the local fuel stratification notably influenced early-stage flame kernel development due to varying injection spray patterns and the fuel injection strategies. Among these design parameters, the number of nozzle holes and fuel injection timing had the most significant effects on the engine combustion performance.

42 ENGINEERING↗
DOE OSTI2024

Sensitivity Analyses of Natural Convection in the AP1000 Passive Containment Cooling System Following LBLOCA Using CFD

The AP1000 power plant has implemented multiple layers of passive safety systems to enhance reactor safety and ensure system integrity during postulated scenarios, such as loss of coolant accident (LOCA) and main steam line break (MSLB). Among these, the passive containment cooling system (PCCS) is a safety-related system designed to prevent the containment from exceeding the design limits of both pressure and temperature following a postulated design basis accident, by transferring heat from the steel containment vessel to the atmosphere. During LOCA, natural air convection plays a vital role in removing heat from the steel containment vessel to the atmosphere. This is integrated with the condensate film from the Passive Containment Cooling Water Storage Tank (PCCWST). Air natural convection is also the only heat removal mechanism after the dry-out of the PCCWST and during the loss of shutdown decay heat removal (LOSDHR) event. Therefore, it is essential to investigate the thermal hydraulics behavior of the containment under transient conditions to ensure its safety and integrity. This paper presents a simplified CFD/ANSYS FLUENT model developed to analyze the impact of different parameters, such as air relative humidity, air temperature, and steel containment temperature, on the natural convection capability to remove the decay heat following LOCA. The model aims to examine the thermal behavior of the containment and provides a comprehensive understanding of the system's natural convection capability during postulated accidents. The results obtained from the model can be used to improve the safety and integrity of the containment system and provide valuable insights for future research in the field.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

OVERVIEW OF PARTICLE DEPOSITION MODELS FOR SPENT NUCLEAR FUEL STORAGE SYSTEMS

Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.

Suffield, Sarah R.↗
DOE OSTI2024

Going beyond BEM with BEM: an insight into dynamic inflow effects on floating wind turbines

Blade element momentum (BEM) theory is the backbone of many industry-standard wind turbine aerodynamic models. To be applied to a broader set of engineering problems, BEM models have been extended since their inception and now include several empirical corrections. These models have benefitted from decades of development and refinement and have been extensively used and validated, proving their adequacy in predicting aerodynamic forces of horizontal-axis wind turbine rotors in most scenarios. However, the analysis of floating offshore wind turbines (FOWTs) introduces new sets of challenges, especially if new-generation large and flexible machines are considered. In fact, due to the combined action of wind and waves and their interaction with the turbine structure and control system, these machines are subject to unsteady motion and thus unsteady inflow on the wind turbine's blades, which could put BEM models to the test. Consensus has not been reached on the accuracy limits of BEM in these conditions. This study contributes to the ongoing research on the topic by systematically comparing four different aerodynamic models, ranging from BEM to computational fluid dynamics, in an attempt to shed light on the unsteady aerodynamic phenomena that are at stake in FOWTs and whether BEM is able to model them appropriately. Simulations are performed on the UNAFLOW 1:75 scale rotor during imposed harmonic surge and pitch motion. Experimental results are available for these conditions and are used for baseline validation. The rotor is analyzed in both rated operating conditions and low wind speeds, where unsteady aerodynamic effects are expected to be more pronounced. Results show that BEM, despite its simplicity, can adequately model the aerodynamics of FOWTs in most conditions if augmented with a dynamic inflow model.

17 WIND ENERGY↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2026

A Scaling Study for Incompressible Multispecies Solver in Vertex-CFD

Multispecies incompressible flows occur widely in engineering and environmental applications, such as chemical reactors, fuel cells, ocean mixing, and biomedical systems. However, accurately resolving the complex transport and mixing phenomena associated with multiple interacting species remains computationally challenging, especially for large-scale problems. In this study, we present a robust, high-performance computing--enabled multispecies incompressible Navier–Stokes solver integrated within the Vertex-CFD framework. Our solver employs a fully coupled, implicit, finite element--based formulation that accurately captures the advection, diffusion, and interaction of multiple species in incompressible flows by leveraging the Kokkos library for parallel computing to achieve high computational efficiency. For pressure coupling, the entropically damped artificial compressibility method is utilized. We validated the solver against canonical test cases, including multispecies advection, diffusion, and Bateman systems; the results demonstrate second- and third-order spatial accuracy and consistent convergence. Additionally, we demonstrated the strong and weak scaling study results obtained on the leadership-class high-performance computing system, Frontier at Oak Ridge National Laboratory.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
DOE OSTI2021

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗
DOE OSTI2020

Validation of blade-resolved computational fluid dynamics for a MW-scale turbine rotor in atmospheric flow

In this study, simulation results of two different computational fluid dynamics codes, Nalu-Wind and EllipSys3D, are presented for a wind turbine rotor in complex yawed and sheared inflow. The results are compared to measurements from the DanAero experiments, to validate computed pressures and azimuthal trends. Despite different code methodologies and grid setups, the codes agree well in computed pressures and integrated forces along the blade for all blade azimuthal positions, however with some discrepancy in the very yawed case. Additionally, both codes capture well the azimuthal trends and force levels seen in measurements. Investigation into discrepancies shows that expanding grids before the rotor, lead to smearing of the wind profiles, which is likely the main cause of the differences in the results between the codes. Additionally, omission of the ground constraint cause discrepancies in relative velocity seen by the passing blade, due to an over speeding beneath the rotor.

17 WIND ENERGY↗
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 OSTI2024

Simulation of Microcapsule Transport in Fractured Media Using Coupled CFD‐DEM

Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Piston geometry and stroke optimization for high efficiency propane spark ignition engines

Propane has unique properties and offers interesting characteristics for high-efficiency spark ignition engines. Its high volatility reduces or completely eliminates fuel-wall wetting and facilitates fuel air mixing. Furthermore, propane has a research octane number of 112 and a high octane sensitivity of 15. Finally, its laminar flame speed is on the same order as that of conventional gasoline, and it exhibits high dilution tolerance. Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, high stroke-to-bore ratios and high geometric compression ratios have been used in new engine designs. Therefore, propane’s relatively high laminar flame speeds, high knock resistance, and dilution tolerance make it an excellent candidate fuel for modern spark ignition engines. The objective of this work is to co-optimize the piston geometry and the engine stroke to maximize the efficiency of a spark-ignition engine fueled with propane. 3D computational fluid dynamics (CFD) simulations employing the extended coherent flamelet model were used to study the parametric effects of piston shape and stroke length. A piston geometry based on high performing pistons was parameterized using four controlling parameters. The piston geometry and engine stroke design space was explored using deterministic and quasi-random sampling techniques. In conclusion, a Gaussian process regression model was built using the simulation data to explain the results observed.

33 ADVANCED PROPULSION SYSTEMS↗