SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A hybrid porous model for full reactor core scale CFD investigation of a prismatic HTGR

Three-dimensional (3-D) Computational Fluid Dynamics (CFD) analysis of a whole nuclear reactor core is a tremendous challenge due to the large geometric volume and complex structures. Here, this research presents a hybrid porous (HP) model to simplify a prismatic High Temperature Gas-cooled Reactor (HTGR) core, so 3-D CFD investigation can be performed on a full reactor core scale. In the HP model, the prototypic small coolant channels in the nuclear fuel blocks are lumped together to form multiple equivalent large coolant channels, and then the porous medium flow model is applied to each of them. Therefore, heat transfer in fuel blocks is computed by a hybrid combination of solid energy and porous flow energy equations. The similarity between the HP model and prototypic model is achieved by deriving the porous flow permeability, inertial resistance factor, and artificial thermophysical properties. Compared with the widely used whole porous (WP) flow model, the HP model preserves more realistic geometric structures, and therefore more accurate physical processes. The General Atomics' Modular High Temperature Gas-cooled Reactor (MHTGR) design was chosen as a prototype to demonstrate the methodology. Simulations were performed using the prototypic CFD model and HP model at steady-state forced circulation, steady-state natural circulation, and transient conditions that correspond to normal operation, extended period of pressurized cool down, and short-term transients after reactor shutdown, respectively. The comparison shows good agreement between the HP model and prototypic model in the maximum fuel temperature, average solid temperature, and helium flow rate, which demonstrates the potential applicability of the HP model for a full reactor core scale simulation in the future. As a benefit, the HP model reduces the mesh quantity by a factor of 50 from a prototypic model. Correspondingly, the computation time was reduced by a factor of at least 30.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗

DOE OSTI2024

Lie–Poisson Neural Networks (LPNets): Data-based computing of Hamiltonian systems with symmetries

An accurate data-based prediction of the long-term evolution of Hamiltonian systems requires a network that preserves the appropriate structure under each time step. Every Hamiltonian system contains two essential ingredients: the Poisson bracket and the Hamiltonian. Hamiltonian systems with symmetries, whose paradigm examples are the Lie–Poisson systems, have been shown to describe a broad category of physical phenomena, from satellite motion to underwater vehicles, fluids, geophysical applications, complex fluids, and plasma physics. The Poisson bracket in these systems comes from the symmetries, while the Hamiltonian comes from the underlying physics. We view the symmetry of the system as primary, hence the Lie–Poisson bracket is known exactly, whereas the Hamiltonian is regarded as coming from physics and is considered not known, or known approximately. Using this approach, we develop a network based on transformations that exactly preserve the Poisson bracket and the special functions of the Lie–Poisson systems (Casimirs) to machine precision. We present two flavors of such systems: one, where the parameters of transformations are computed from data using a dense neural network (LPNets), and another, where the composition of transformations is used as building blocks (G-LPNets). We also show how to adapt these methods to a larger class of Poisson brackets. We apply the resulting methods to several examples, such as rigid body (satellite) motion, underwater vehicles, a particle in a magnetic field, and others. The methods developed in this paper are important for the construction of accurate data-based methods for simulating the long-term dynamics of physical systems.

Data-based modeling↗

DOE OSTI2023

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗

DOE OSTI2023

Modeling Heat Pipes with Non-Condensable Gases

The report showcases advancements in modeling heat pipes, considering both scenarios with and without non-condensable gases. To achieve this, two distinct modeling approaches are juxtaposed. Firstly, the effective conduction model, implemented in the MOOSE-based code Sockeye, is employed. Secondly, a first-of-a-kind two-phase Euler-Euler Computational Fluid Dynamics (CFD) model is developed using the STAR-CCM+ code. Both models undergo validation against experimental data, acknowledging the inherent uncertainties associated with each modeling assumption. Interestingly, the non-tuned CFD model surpasses the performance of the calibrated conduction model for heat pipes operating under both conditions: with and without non-condensable gases. It's worth noting, however, that the CFD models entail significantly longer runtimes compared to the conduction models. Nevertheless, the insights garnered from the CFD model shed invaluable light on the intricate operational dynamics of heat pipes. Future work involves broadening the validation scope of these models and continuing their development to enhance their utility as robust tools for heat pipe design and operational support.

42 ENGINEERING↗

DOE OSTI2022

Defect-driven anomalous transport in fast-ion conducting solid electrolytes

We report solid-state ionic conduction is a key enabler of electrochemical energy storage and conversion. The mechanistic connections between material processing, defect chemistry, transport dynamics, and practical performance are of considerable importance, but remain incomplete. Here, inspired by studies of fluids and biophysical systems, we re-examine anomalous diffusion in the iconic two-dimensional fast-ion conductors, the β- and β"-aluminas. Using large-scale simulations, we reproduce the frequency dependence of alternating-current ionic conductivity data. We show how the distribution of charge-compensating defects, modulated by processing, drives static and dynamic disorder, which lead to persistent sub-diffusive ion transport at macroscopic timescales. We deconvolute the effects of repulsions between mobile ions, the attraction between the mobile ions and charge-compensating defects, and geometric crowding on ionic conductivity. Our quantitative framework based on these model solid electrolytes connects their atomistic defect chemistry to macroscopic performance with minimal assumptions and enables mechanism-driven 'atoms-to-device' optimization of fast-ion conductors.

25 ENERGY STORAGE↗

DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Non-Nuclear Test Facility Development to Support Verification and Validation Needs

This study focuses on developing a non-nuclear test facility to help meet the verification and validation (V&V) requirements of a water-cooled small modular reactor (SMR) system. The scaling analysis necessary for reactor system experimentation and assessment model V&V considers various reactor system accident scenarios such as main steam line break (MSLB), steam generator tube rupture (SGTR), and containment condensation and heat transfer. The goal is to identify key physics based phenomena of interest (POI), specific figures of merit (FOMs), and the state of knowledge (SOK) by employing phenomena identification and ranking table (PIRT) studies. Important physics phenomena such as heat transfer, fluid dynamics, structural integrity, and material selection must be addressed when scaling a thermal-hydraulics system. Alongside these considerations, modifications to the facility may be needed to align the scaled system with the available footprint (height and space) of the facility. These modifications may include infrastructure upgrades, changes to the piping layout, implementation of safety measures, and installation of new instrumentation and control systems. The test facility aims to provide a scaled-down version of prototypic facilities for comprehensive V&V activities by overcoming these challenges and making the necessary modifications. This approach minimizes scaling distortion, simulates prototypic plant conditions, and accurately models accident event progressions. Thus, the facility will be capable of adequately modeling a scaled SMR system to support the design, development and licensing needs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Investigation of Farm-to-Farm Interactions and Blockage Effects from AWAKEN using Large-Scale Numerical Simulations

A large-scale numerical computation of five wind farms was performed as a part of the American Wake Experiment (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100km × 100km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions with wind farms with large scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

Cheung, Lawrence C.↗

DOE OSTI2024

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗

DOE OSTI2020

PyDA: A Hands-On Introduction to Dynamical Data Assimilation with Python

Dynamic data assimilation offers a suite of algorithms that merge measurement data with numerical simulations to predict accurate state trajectories. Meteorological centers rely heavily on data assimilation to achieve trustworthy weather forecast. With the advance in measurement systems, as well as the reduction in sensor prices, data assimilation (DA) techniques are applicable to various fields, other than meteorology. However, beginners usually face hardships digesting the core ideas from the available sophisticated resources requiring a steep learning curve. In this tutorial, we lay out the mathematical principles behind DA with easy-to-follow Python module implementations so that this group of newcomers can quickly feel the essence of DA algorithms. We explore a series of common variational, and sequential techniques, and highlight major differences and potential extensions. We demonstrate the presented approaches using an array of fluid flow applications with varying levels of complexity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗

DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Computationally-Aided Design of a Small-Scale Radioactive Waste Glass Melter

To provide mission support to the Hanford Waste Immobilization and Treatment Plant for the vitrification of legacy nuclear tank waste, a reduced-scale vitrification pilot system is being designed to process simulated and actual radioactive tank wastes. The tank waste will be separated into high-level waste (HLW) and low-activity waste (LAW) fractions, where the majority by mass (~90%) is LAW and by activity (>95%) is HLW. The first tank waste to be processed at the WTP will be LAW. Since Hanford tank waste and the resultant melter feed compositions are known to vary widely, pilot-scale operations are essential to identify potential problems, provide needed data, confirm assumptions, and determine impacts to full-scale melters and off-gas systems from these different feeds. In addition, the reduced-scale melter system must be capable of quickly providing results to operations without the typically high costs of radioactive operations and provide an engineering platform in close proximity to local operations staff. The objective is to produce similar process conditions to those encountered in the full-scale LAW melters and minimize the volume of radioactive waste necessary for the evaluations. To this end, melter design features are being explored to increase production without significantly increasing surface area, melter glass volumes or compromising data quality. In support of these objectives, a set of computational fluid dynamic models have been developed to provide insight into the flow patterns within these small melters and evaluate design features to increase production without significantly increasing glass inventory. Computed velocities were evaluated at the interface between the cold cap and the molten glass pool and within the bulk glass pool. A method was developed to correlate results to an estimated melt rate and down-select design features that produce the highest gains.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Computationally-Aided Design of a Small-Scale Radioactive Waste Glass Melter

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Numerical Investigation of Butterfly Valve Performance in Variable Valve Sizes, Positions and Flow Regimes

Reliability and efficiency of valves are necessary for precise control and sufficient heat-flow to heat application plants for the integrated energy systems of nuclear power plants (NPPs). Strategic Management Analysis Requirement and Technology (SMART) valves’ ability to control flow and assess environmental parameters stands out for these requirements. Their ability to sustain the downstream flow rate, prevent reverse flow, and maintain pressure in the heat transport loop is much more efficient with the integration of sensors and intelligent algorithms. For assessing valve performance and monitoring, mechanical design and operating conditions are two important parameters. In this study, the butterfly valves of three different sizes are simulated with water and steam using STAR-CCM+ in various flow regimes and positions to analyze performance parameters to strategize an automated control system for efficiently balancing the heat–transport network. Also, flow behavior is studied using velocity and pressure fields for valve–body geometry optimization. It can be observed, through performance parameters, that the valves are suitable for operation between 30° and 90° positions with significantly low loss coefficients and high flow coefficients, and the performance parameters follow a certain pattern in both water and steam flow in each scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Barriers and variable spacing enhance convective cooling and increase power output in solar PV plants

When the temperature of solar photovoltaic (PV) modules rises, efficiency drops and module degradation accelerates. Thus, it is beneficial to reduce module operating temperatures. Previous studies of solar power plants have illustrated that incoming flow characteristics, turbulent mixing, and array geometry can strongly impact convective cooling, as measured by the convective heat transfer coefficient h. In the fields of heat transfer and plant canopy flow, previous work has shown that system-scale arrangement modifications—e.g., variable spacing, barriers, or windbreaks—can passively alter the flow, enhance turbulent mixing, and influence convection. However, researchers have not yet explored how variable spacing or barriers might enhance convective cooling in solar power plants. Here, high-resolution large-eddy simulations model the air flow and heat transfer through solar power plant arrangements modified with missing modules and barrier walls. We then perform a control volume analysis to evaluate the net heat flux and compute h, which quantifies the influence of these spatial modifications on convective cooling and, thus, module temperature and power output. Installing barrier walls yields the greatest improvements, increasing h by 3.4%, reducing module temperature by an estimated 2.5 °C, and boosting power output by an estimated 1.4% on average. These findings indicate that incorporating variable spacing or barrier-type elements into PV plant designs can reduce module temperature and, thus, improve PV performance and service life.

14 SOLAR ENERGY↗

DOE OSTI2020

kynema-sgf [SWR-20-85]

Kynema-SGF (formerly AMR-Wind), wherein SGF stands for structured-grid fluid dynamics, is a massively parallel, block-structured adaptive-mesh, incompressible flow solver. The codebase was initiated in 2019 from incflo. The solver is built on top of the AMReX library. AMReX library provides the mesh data structures, mesh adaptivity, as well as the linear solvers used for solving the governing equations. Kynema-SGF is actively developed and maintained by a dedicated multi-institutional team from Lawrence Berkeley National Laboratory, National Laboratory of the Rockies, and Sandia National Laboratories. The primary applications for Kynema-SGF are: performing large-eddy simulations (LES) of atmospheric boundary layer (ABL) flows, simulating wind farm turbine-wake interactions using actuator disk or actuator line models for turbines, and as a background solver when coupled with a near-body solver (e.g., Kynema-UGF) with overset methodology to perform blade-resolved simulations of multiple wind turbines within a wind farm. For offshore applications, the ability to model the air-sea interaction effects and its impact on the ABL characteristics is another focus for the code development effort. As with other codes in the Kynema ecosystem, Kynema-SGF shares the following objectives: *an open, well-documented implementation of the state-of-the-art computational models for modeling wind farm flow physics at various fidelities that are backed by a comprehensive verification and validation (V&V) process; *be capable of performing the highest-fidelity simulations of flow fields within wind farms; and *be able to leverage the high-performance leadership class computing facilities available at DOE national laboratories.

Ananthan, Shreyas↗