SEARCH · Engineering Papers

Results for “openfoam”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

DOE OSTI2023

Entropy-driven Optimal Sub-sampling of Fluid Dynamics for Developing Machine-learned Surrogates

Optimal sub-sampling of large datasets from fluid dynamics simulations is essential for training reduced-order machine learned models. A method using Shannon entropy was developed to weight flow features according to their level of information content, such that the most informative features can be extracted and used for training a surrogate model. The method is demonstrated in the canonical flow over a cylinder problem simulated with OpenFOAM. Both time-independent predictions and temporal forecasting were investigated as well as two types of prediction targets: local per-grid-point predictions and global per-time-step predictions. When tested on training a surrogate model, results indicate that our entropy-based sampling method typically outperforms random sampling and yields more reproducible results in less iterations. Finally, the method was used to train a surrogate model for modeling turbulence in magnetohydrodynamic flows, which revealed various challenges and opportunities for future research.

Brewer, Wes↗

DOE OSTI2021

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗

DOE OSTI2021

densegranFoam

densegranFoam is a dense granular material flow solver developed using the OpenFoam Computational Fluid Dynamics (CFD) tool. It is developed using Euler-Euler framework. The transport model is implemented using the shear rate dependant constitutive model. The primary focus of the solver is to simulate shear rate dependency on the dense granular flow bulk behavior.

Ahsan, Syed↗

DOE OSTI2025

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗

DOE OSTI2025

Open-Datasets for WEC Simulation

SAND2025-11468O The Open-Datasets for WEC Simulation is a tool that uses datasets and simulation configuration files to conduct OpenFOAM wave energy converter (WEC) simulations. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Chartrand, Chris [Sandia National Lab. (SNL-CA), L↗

DOE OSTI2023

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗

DOE OSTI2021

Design of high-deflection foils MHK applications - CFD models - RivGen turbine

The Ocean Renewable Power Company's (ORPC's) goal is to design, develop, and test hydrofoils with large deflections. The effects of the deflections on cross-flow turbine performance would be evaluated in order to inform design considerations for full-scale water turbines and other marine hydrokinetic devices. CFD Models OpenFOAM V1912 RivGen turbine

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Design of high-deflection foils MHK applications - CFD models - Helical turbines

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Modeling Induction Stirring and Particle Tracking in Molten Uranium

Two independent numerical models have been developed to simulate the behavior of carbon impurities in molten uranium metal. Informed by experimental parameters, one model was created using the commercial software Star-CCM+ and compared with another developed using open-source codes, including OpenFOAM, Finite Element Method Magnetics (FEMM) and a First Passage Kinetic Monte Carlo (FPKMC) approach. The target experimental system features a 404g uranium metal charge containing an average carbon concentration of 139 ppm which was melted in a vacuum induction furnace at 1400° C then resolidified. The microstructures of the uranium and its impurities before and after melting have been characterized and reported separately. Prior to simulating the uranium-carbon system described, the numerical models were validated using a previously published nonradioactive experimental system to ensure agreement with expected output values. Focus has been placed on modeling velocity fields under induction stirring and impurity particle trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

MOSCATO Development and Integration in Fiscal Year 2022

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. We built upon the work for FY22 and improved the code from several different perspectives. First, we improved the user interface by adding a new component to the official Nek5000 input file (.par). This new part contains documents parameters like, salt properties (density, viscosity, Cp, thermal conductivity), diffusion coefficients, standard potential, etc. Second, we built a conversion script to extract salt properties from the MSTDB-TP salt database and write to MOSCATO input file. Third, we migrated the code to NekRS, which is the GPU branch of Nek5000 and suitable for next generation supercomputers. Verification and Validation (V&V) work was also continued in FY22. Two tasks were performed. The first V&V task involved the validation of the Poisson-Nernst-Planck equation solver and Butler-Volmer electrode kinetics, by comparing with numerical and experimental data about thermoelectric cells. The second task involved the comparisons to corrosion results from a thermal convection loop run during the MSRE era. Satisfactory agreement was obtained from both tasks.

Yuan, Haomin↗

DOE OSTI2023

Laser heating and evaporation of a single droplet

The laser technology is being abundantly studied for controlled energy deposition for a range of applications in aerodynamic flow control, material processing, ignition, and combustion. The absorption of laser radiation by liquid droplets affects further propagation of laser in the atmosphere and causes bleaching of suspended droplets while the ignition and combustion characteristics in combustors are influenced by the evaporation rate of the sprayed fuel. In this work, we present a multi-dimensional mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion dominated regime taking into account absorption of the laser radiation, evaporation process and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for the ethanol and water droplet heating and evaporation. For continuous heating the peak temperature is established by the balance of cooling, evaporation and heating and results in high temperature for larger droplets. It has been shown that for heating by a single laser pulse the maximum temperature of droplets depends only on the peak intensity of the laser radiation. Furthermore, for the peak irradiance close to the transition to the boiling regime, temporal dynamics of the droplet temperature is independent of the droplet size. With proper normalization of time, the dynamics of the droplet shrinkage and cooling is shown to be independent of droplet sizes and peak laser intensities. The influence of cooling and evaporation processes on droplet heating was found to be controlled by the pulse repetition rate for repeated pulse operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

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 OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Aerodynamic Characterization of 3D Scanned Wind Turbine Blades Using Experimental and Computational Methods

This study presents an aerodynamic characterization of 3D scanned wind turbine blades using both experimental and computational methods. The research was conducted by Gulf Wind Technology and Sandia National Laboratories. The primary objective was to investigate the aerodynamic impacts of leading-edge manufacturing defects on wind turbine blades. The study utilized the Stratasys NEO 800 3D Printer for high-precision manufacturing and the GWT Accelerator Wind Tunnel for experimental testing. Computational simulations were performed using COMSOL Multiphysics to model the wind tunnel and analyze flow characteristics and OpenFOAM to study the aerodynamic impacts of leading-edge defects. OpenFAST was used to estimate how these defects can lead to revenue losses for wind farm operators as high as 6%. The results demonstrated significant aerodynamic performance variations due to defects, with detailed analysis provided through wind tunnel and CFD data. The findings contribute to the understanding of defect impacts on wind turbine blade performance and offer insights for future design improvements.

17 WIND ENERGY↗