Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 703 records · Page 39

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Fine Scale Stress Test Modelling

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

15 GEOTHERMAL ENERGY↗

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding flow dynamics in membrane distillation: Effects of reactor design on polarization

Optimization and design of full-scale membrane distillation (MD) systems usually require Sherwood and Nusselt correlations that are developed from lab-scale systems. However, entrance effects in lab-scale systems can significantly impact heat, mass and momentum transfer in the reactor, therefore affect the accuracy of the developed experimental Sherwood and Nusselt correlations. Here, Computational Fluid Dynamics (CFD) simulations using OpenFOAM are performed to understand the effects of right-angled bends and inlet design on flow dynamics, temperature and concentration polarization in MD systems. Simulation results show that the presence of right-angled bends and inlets with sudden expansions lead to the formation of Dean vortices. Dean vortices enhance perpendicular mixing in MD systems and reduce both temperature and concentration polarization. Temperature and concentration polarization coefficients in MD systems with right-angled bends and inlets with sudden expansions vary significantly for the same volumetric flow rate. Our studies show that lab-scale systems with the same volumetric flow rate but different designs lead to significantly different Nusselt and Sherwood correlations. Here, this study demonstrates the importance of CFD-informed design of lab-scale systems to minimize entrance effects and suppress Dean vortices for consistent model development and calibration across multiple scales.

42 ENGINEERING↗

Ocean Energy Sandia TEAMER WEC Simulation Results

Computational fluid simulations for wave energy converters and supporting materials from Ocean Energy's WEC Buoy TEAMER project in collaboration with Sandia National Laboratories. Each file includes images and video of simulation along with the simulation files for the case. The videos and images include air and water simulations of the buoy in 2D and 3D views. Data spreadsheets of the simulation outputs are also included in the files. The data are represented as different cases in which the wave height(Hs) and peak period(Tp) were different.

16 TIDAL AND WAVE POWER↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer

A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.

09 BIOMASS FUELS↗

Design of high deflection foils for MHK applications - CFD files

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. OpenFOAM V1912 files for straight foil model scale turbines in the University of New Hampshire tow tank. Strut Locations = (0.13, 0.225, 0.450, 0.675, 0.900) [m] Tip speed ratio = 2.40

16 TIDAL AND WAVE POWER↗

A reacting multiphase computational flow model for 2,3-butanediol synthesis in industrial-scale bioreactors

The microbe, Zymomonas mobilis, can efficiently convert sugars to 2,3 butanediol (BDO), an important intermediate for downstream chemical products, only in a well-controlled microaerated environment. However, controlling oxygen distribution in industrial-scale bioreactors is challenging, and at-scale BDO production is hard to optimize using traditional engineering methods. This study takes a step towards addressing this problem through a computational model for reacting multiphase flows in large-scale bioreactors. A phenomenological metabolic model was first developed and validated against experiments, and then it was coupled to a multiphase computational fluid dynamics (CFD) solver using a subcycling algorithm for simulating long conversion times (~30 hours). Large-scale (~500 m 3 ) bubble column simulations using this coupled approach demonstrated a 25% improvement in BDO yield with low sparging rates (0.002 m/s) and low height-to-diameter ratio (0.875) compared to a baseline bubble column (35 m height, 5 m diameter) with a superficial gas velocity of 0.06 m/s.

09 BIOMASS FUELS↗

A simulation-based parametric study of CLOU chemical looping reactor performance

Chemical looping with oxygen uncoupling (CLOU) is a variant on chemical looping combustion in which the oxygen carrier releases gaseous O 2 in the fuel reactor, making it well-suited for solid fuels, since the released gaseous O 2 readily reacts with solid char. This work presents several computational fluid dynamic (CFD) simulations of copper-based CLOU in a dual fluidized bed system, each with different operating conditions. The modeling predicted that coal particle sizes as large as 1000 μm did not significantly affect performance. Increased oxygen carrier copper loading resulted in an excess of gaseous oxygen in the product gas stream. Decreasing the oxygen carrier bed mass as well as reducing the air reactor fluidizing velocity did not supply enough oxygen to the fuel reactor to complete combustion of the coal. This generated a failure state in which the temperature continued to decrease in the fuel reactor from the lack of combustion, which in turn reduced the O 2 equilibrium partial pressure, further lowering the amount of combustion possible. Sufficient O 2 can be maintained in the fuel reactor by ensuring a high enough air reactor velocity and a large enough supply of oxygen carrier inventory to handle the chosen coal feed rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat Exchanger Intensification Through Powder Processing and Enhanced Design (HIPPED)

High intensity, supercritical CO 2 heat exchangers capable of operating at temperatures between 800 and 1100 ˚C and pressures greater than 80 bar are considered a critical component for high efficiency power generation and a range of next-generation industrial processes. A novel heat exchanger concept was proposed in this project. The proposed concept involved a plate-type architecture with three-dimensional fins, called twisted s-shaped fins, that capitalize on additive manufacturing. Four manufacturing processes were studied to realize the complex envisioned design: press and sinter of metallic powders, directed energy deposition, laser powder bed fusion, and a novel powder bed process involving a photoreactive polymer called SEAM (scalable expeditious additive manufacturing), that later was developed at MSU. The team was able to manufacture subscale heat exchangers using laser powder bed fusion (LPBF) and SEAM. The team settled on using SEAM as a final manufacturing approach due to its very low cost. The SEAM-manufactured heat exchanger, however, when tested under pressure, exhibited leaks and the process still requires further development beyond the timeline of the project. To support the heat exchanger design and testing, the team developed a new high temperature flow loop allowing testing under high temperature and pressure. Computational modeling of the heat exchanger was also performed for the twisted s-shaped fins and the computations were validated with experimental data. The new twisted fins resulted in enhanced heat transfer and reduced pressure drop. To further minimize pressure drop in the plate heat exchanger, the team also developed a simplified model for flow maldistribution in the manifolds, providing a new tool for designers that eliminates the need for complex CFD calculations. A technoeconomic analysis was also developed to contrast SEAM with LPBF, this tool indicated that SEAM is several times less costly and offer significant promises as a rapid low cost additive manufacturing process for metallic systems.

30 DIRECT ENERGY CONVERSION↗

An indirect approach to optimize the reaction rates of thermal NO formation for diesel engines

With stringent emission regulations, it has become more important for modern diesel engine manufacturers to accurately predict engine-out nitrogen oxide (NO x ) emissions across a wide range of operating conditions. Thermal NO is the major source of engine-out NO x in modern diesel engines. For thermal NO formation, several earlier studies have recommended the forward and reverse reaction rate coefficients of the rate-limiting reaction (O + N 2 ⇌ NO + N). However, due to deficiencies in sub-models and inadequacies of reduced chemical mechanisms to represent diesel combustion, these recommended values more often than not need to be adjusted in reduced order combustion models to accurately predict engine-out NO x . Hence, in this work a systematic and computationally efficient approach has been proposed to streamline the process of determining the optimum reaction rate coefficients. Here, to develop the optimization approach, four different production diesel engines with different operating conditions in terms of speed, load, and exhaust-gas recirculation have been considered. Numerical simulations have been performed using a detailed zero-dimensional velocity-composition-frequency transported probability density function (0D-VCF-tPDF) model that uses hundreds of notional particles to capture in-cylinder stratification. Four different combinations of hydrocarbon and NO x chemical mechanisms were used to represent chemistry. It was found that for the rate-limiting reaction, the pre-exponent factors (A f1 , A r1 ) and activation energies (E A,f1 , E A,r1 ) of the forward and reverse reaction rates follow a linear band in A f1 - E A,f1 and A r1 - E A,r1 space where predicted engine-out NO x match the measured values closely. By encompassing such bands from different engines and considering constraints on activation energies, a reduced search domain of pre-exponent factors and activation energies was constructed that is expected to be applicable to any diesel engine. Eventually, computationally efficient three-line and one-line search approaches were proposed to determine the optimum values of the pre-exponent factors and activation energies that led to a minimum error between measured and predicted engine-out NO x . Finally, these three-line and one-line NO x optimization approaches were applied to a fifth production diesel engine for which the 0D-VCF-tPDF model showed a very good predictive performance in terms of predicting peak pressure, 50% burn rate, and engine-out NO x when compared to measured and 3D-CFD values.

33 ADVANCED PROPULSION SYSTEMS↗

Comparison of URANS and LES predictions for the open phase of the OECD NEA CSNI fluid structure interaction CFD benchmark

The OECD NEA CSNI WGAMA CFD Task Group ran a benchmark in 2020 and 2021 to assess the predictive capabilities of coupled fluid structure interaction (FSI) CFD analysis methods. This paper presents the predictions made for the open phase of the benchmark using URANS and LES turbulence modelling approaches, and a comparison of the results to the experimental data. The benchmark comprised a channel containing two inline cylinders in cross-flow. The cylinders were fixed at one end, free at the other, and had measured resonant frequencies and damping properties. The URANS modelling used ANSYS Fluent 2-way coupled to ANSYS Mechanical. The LES modelling used Nek5000, 1-way coupled to Diablo. Comparisons with cross-channel velocity profiles are presented, both for the mean flow and its RMS. Comparisons are also made to the frequency spectra for point measurements of fluid velocity and pressure, and for the accelerations of the free end of each cylinder. URANS predicts the average velocity profiles relatively well, and is able to predict the velocity and acceleration spectra at the shedding frequency. However, the frequency content at the 4th harmonic of the shedding frequency is low in the URANS flow fields, and so does not excite accelerations at the resonant frequency of the cylinders. LES makes better predictions of the average profiles, and the velocity spectra agree well at both the shedding frequency and at higher frequencies. In conclusion, the 1-way coupled LES results show good agreement for acceleration spectra.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

Investigation of countercurrent flow profile and liquid holdup in random packed column with local CFD data

Liquid holdup and mass transfer area are critical parameters used in packed column design and CO 2 capture efficiency prediction. Compared to experimental methods, 3D computational fluid dynamics (CFD) can provide high-resolution data for studying countercurrent flow hydrodynamics in a packed column. In this paper, a framework was established for modeling the liquid–gas countercurrent flow in a random packed column with pall rings. Besides the column-averaged information, the radial pall ring distribution, velocity, and liquid holdup profiles are obtained in the simulations. The radial two-phase velocity and liquid holdup profiles show significant fluctuations caused by the non-uniform pall ring distribution. The entrance effect and wall influence are also observed in the packed column. By using local countercurrent flow information, the packing specific area $a_p$ and liquid velocity $u_L$ range are significantly expanded for liquid holdup hL prediction. This enables us to develop the liquid holdup correlation with one packing size and six CFD runs only. The proposed power law of the liquid holdup correlation $h_L~u_L^{0.44}$ indicates the random packed column falls in a viscous to turbulent transition regime with a Reynolds Number range of [6.7–40.2]. The derived liquid holdup correlation is in good agreement with existing correlations that are developed using the column-averaged experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrodynamic Modelling for Tidal Energy Kites: Cooperative Research and Development Final Report, CRADA Number CRD-24-30123

The primary goal of the work is to improve accuracy of hydrodynamic loads acting on the lifting surface of the kite turbine due to complex air-water-lifting body interaction in a low-order numerical model using high-fidelity numerical input. The result of this study will help the Participant better design control algorithms for the kite turbine so that the turbine can harvest more energy while maintaining stability of the platform as well as its reliability. The use of hydrodynamic coefficients/load components from computational fluid dynamic simulation will inform better inputs for a medium-fidelity dynamic model (specifically an OrcaFlex model) which is intended for use in rapid design iterations in the design process.

16 TIDAL AND WAVE POWER↗

Editorial: Particulate Matter Emissions From Conventional and Reformulated Fuel Combustion: Advances in Experiments and Simulations

Particulate matter (PM or soot) emissions generated by combustion of conventional and reformulated fuels adversely impact the environment and human health (Bond et al., 2013; Landrigan et al., 2018; World Health Organization (WHO), 2016), generating public awareness and motivating efforts towards the mitigation of their harmful effects. Soot formation is one of the most complex phenomena in combustion, involving interactions between combustion chemistry, fluid mechanics, mass/heat transport, and particle dynamics, spanning different spatial and temporal scales. A fundamental understanding of soot formation process is thus necessary to achieve a strong reduction of PM emissions and design cleaner and more efficient combustion systems. This has motivated a longstanding and ongoing research activity aimed at improving our understanding of the physical and chemical processes involved in soot formation, well reviewed in these recent works (Wang and Chung, 2019; Michelsen et al., 2020; Martin et al., 2022). Despite the wide interest, the transition from gas-phase molecules to incipient soot particles is still elusive and the successive particle growth and oxidation processes are far from being fully understood, especially in conditions relevant to real-world applications. Recent advances in combustion PM emission diagnostic and computational capabilities helped in improving the predictability of fundamental chemical and aerosol models for practical applications, thus tackling some of the above-mentioned challenges. The aim of this research topic is to display the ongoing research efforts in addressing the existing gaps on particulate formation from various fuel sources and in conditions typical of practical combustion applications (e.g. flames, engines, pool fires), through both experimental and numerical approaches.

42 ENGINEERING↗

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Large Scale Grid

Mesh, properties, initial conditions, injection/withdrawal rates for modeling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, on large scale grid, with some results. First simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 1-5 pertain.

15 GEOTHERMAL ENERGY↗

An in-situ conductometric apparatus for physicochemical characterization of solutions and in-line monitoring of separation processes at elevated temperatures and pressures

Specific conductance and frequency-dependent resistance (impedance) data are widely utilized for understanding the physicochemical characteristics of aqueous and non-aqueous fluids and for evaluating the performance of chemical processes. However, the implementation of such an in-situ probe in high-temperature and high-pressure environments is not trivial. This work provides a description of both the hardware and software associated with implementing a parallel-type in-situ electrochemical sensor. The sensor can be used for in-line monitoring of thermal desalination processes and for impedance measurements in fluids at high temperature and pressure. Further, a comparison between the experimental measurements on the specific conductance in aqueous sodium chloride solutions and the conductance model demonstrate that the methodology yields reasonable agreement with both the model and literature data. A combination of hardware components, a software-based correction for experimental artifacts, and computational fluid dynamics (CFD) calculations used in this work provide a sound basis for implementing such in-situ electrochemical sensors to measure frequency-dependent resistance spectra.

47 OTHER INSTRUMENTATION↗