Engineering topics

Busco, Giacomo

Publications and source records attributed to Busco, Giacomo.

4 records

DOE OSTI2023

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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 OSTI2020

Sneezing and asymptomatic virus transmission

The novel coronavirus disease (COVID-19) spread pattern continues to show that geographical barriers alone cannot contain a virus. Asymptomatic carriers play a critical role in the nature of this virus quickly escalating into a global pandemic. Asymptomatic carriers may transmit the virus unintentionally through sporadic sneezing. A novel Computational Fluid Dynamics (CFD) approach has been proposed with a realistic modeling of a human sneeze achieved by the combination of state-of-the-art experimental and numerical methods. This modeling approach may be suitable for future engineering analyses aimed at reshaping public spaces and common areas, with the main objective to accurately predict the spread of aerosol and droplets that may contain pathogens. This study shows that the biomechanics of a human sneeze, including complex muscle contractions and relaxations, can be accurately modeled by the angular head motion and the dynamic pressure response during sneezing. These have been considered as the human factors and were implemented in the CFD simulation by imposing a momentum source term to the coupled Eulerian–Lagrangian momentum equations. The momentum source was modeled by the measured dynamic pressure response in conjunction with the angular head motion. This approach eliminated the need to create an ad hoc set of inlet boundary conditions. With this proposed technique, it is easier to add multiple fixed and/or moving sources of sneezes in complex computational domains. Additionally, extensive sensitivity analyses based on different environmental conditions were performed, and their impact was described in terms of potential virus spread.

42 ENGINEERING↗

DOE OSTI2020

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗