Researcher profile

Elia Merzari

Elia Merzari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

A Study on Convolution Neural Network for Reconstructing the Temperature Field of Wall-Bounded Flows

In the present study, the capabilities of a new Convolutional Neural Network (CNN) model are explored with the paramount objective of reconstructing the temperature field of wall-bounded flows based on a limited set of measurement points taken at the boundaries of the fluid domain. For that, we employ an algorithm that leverages the CNN capabilities provided with additional information about the governing equations of the physical problem. Great progress has been made in recent years towards reconstructing and characterizing the spatial distribution of physical variables of interest using CNNs. In principle, CNNs can represent any continuous mathematical function with a relatively reduced number of parameters. However, depending on the complexity imposed by the physical problem, this technique becomes unfeasible. The present study employs a Physics Informed Neuron Network technique featuring a data-efficient function approximator. As a proof of concept, the CNN is trained to retrieve the temperature of a heated channel based on a limited number of sensors placed at the boundaries of the domain. In this context, the training data are the temperature fields solutions considering various flows conditions at steady state, e.g varying the Reynolds and the Prandtl numbers. Additionally, a demonstration case considering the more complex geometry of a MSR is also provided. Assessment on the performance of the CNN is done by the mean L2 and the maximum Linf Euclidean norms stemmed from the difference between the actual solutions and the predictions made by the CNN. Finally, a sensitivity analysis is carried out such that the robustness of the CNN is tested considering a potential real application scenario where noise is inevitable. For that, the original test inputs are overlaid with a normal distribution of random numbers targeting to mimic different levels of noise in the measurement points.

preprint2022arXiv

Comparison of Pebble Bed Velocity Profiles Between High-Fidelity and Intermediate-Fidelity Codes

Recent interest for the development of high-temperature gas reactors has increased the need for more advanced understanding of flow characteristics in randomly packed pebble beds. A proper understanding of these flow characteristics can provide a better idea of the cooling capabilities of the system in both normal operation and accident scenarios. In order to enhance the accuracy of computationally efficient, intermediate fidelity modeling, high-fidelity simulation may be used to generate correlative data. For this research, NekRS, a GPU-enabled spectral-element computational fluid dynamics code, was used in order to produce the high-fidelity flow data for beds of 1,568 and 45,000 pebbles. Idaho National Lab's Pronghorn porous media code was used as the intermediate fidelity code. The results of the high-fidelity model were separated into multiple concentric regions in order to extract porosity and velocity averages in each region. The porosity values were input into the Pronghorn model and the resulting velocity profile was compared with that from NekRS. Both cases were run with a Reynolds number of 20,000 based on pebble diameter. The Pronghorn results were found to significantly overestimate the velocity in the outermost region indicating that changes in the porosity alone do not cause the difference in fluid velocity. We conclude that further work is necessary to develop a more effective drag coefficient correlation for the near-wall region and improve predictive capabilities of intermediate fidelity models.

preprint2022arXiv

Direct Numerical Simulation of High Prandtl Number Fluid Flow in the Downcomer of an Advanced Reactor

The passive safety is a crucial feature of advanced nuclear reactor (Gen IV) design. During loss of power scenarios, the downcomer plays a crucial role. The fluid-flow behavior in the downcomer can involve forced to mixed to natural convection and characterizing the heat transfer for these changing regimes is a daunting challenge. High-resolution heat transfer numerical database can potentially support the development of precise and affordable reduced resolution heat transfer models. These models can be designed based on a multiscale hierarchy developed as part of the recently DOE-funded center of excellence for thermal-fluids applications in nuclear energy. In this paper, the downcomer is simplified to heated parallel plates and High Prandtl number fluid (FLiBe) is considered for all simulations. The calculations are performed for a wide range of Richardson number from 0 to 400 at 2 different FLiBe Prandtl numbers (12 and 24) which result in 40 simulated Cases in total. Time averaged and time series statistics as well as Nusselt number correlations are investigated to illuminate the mixed convection behaviors. The calculated database will be instrumental in understanding the flow behaviors in the downcomer. Ultimately, we aimed to evaluate existing heat transfer correlations and some modifications will be proposed for Cases where no satisfactory choice is available.

preprint2022arXiv

Direct Numerical Simulation of high Prandtl number fluids and supercritical carbon dioxide canonical flows using the spectral element method

The design of advanced nuclear reactors (Gen IV) involves an array of challenging fluid-flow issues that affect safety and performance. Currently, these problems are addressed in an ad-hoc manner at varying scales which are time-consuming and expensive. The creation of a high-resolution heat transfer numerical database has the potential to help develop to accurate and inexpensively reduced resolution heat transfer models. Such models can help address industrial-driven issues associated with the heat transfer behavior of advanced reactors. The models can be developed using the multiscale hierarchy developed as part of the recently DOE-funded center of excellence for thermal-fluids applications in nuclear energy. Ultimately this can lead to fast-running reliable models, thus accelerating the deployment of advanced reactors. In this paper, we performed a series of Direct Numerical Simulation using the spectral element codes Nek5000 and NekRS to investigate heat transfer in mixed convection conditions. First, we investigate the heat transfer of the flow in heated parallel plates for high Prandtl number fluids. The calculated database will eventually be used to evaluate existing heat transfer correlations and some modifications will be proposed for cases where no satisfactory choice is available. We have also investigated the heated transfer alteration phenomena in a straight heated tube for supercritical carbon dioxide (sCO2). The low-Mach-number approximation is used to decouple thermal and dynamic pressure, as pressure drop is negligible in this problem. The properties of sCO2 are calculated using multi-region polynomials. We observed that the heat transfer deterioration occurred in combination with the property changes of sCO2 and the depreciation of turbulence kinetic energy (TKE) for upward flow. Whereas, in downward flow, the heat transfer is enhanced thanks to the increase of TKE.

preprint2022arXiv

Direct Numerical Simulation of Low and Unitary Prandtl Number Fluids in Reactor Downcomer Geometry

Buoyancy effect on low-flow condition convective heat transfer of non-conventional coolants, such as liquid metal and molten salts, is a crucial safety factor to advanced reactors under transient or accidental scenarios. The distinct heat transfer characteristics of non-unitary Prandtl fluids and the inherent complexity of the low-flow mixed convection phenomena requires the development of novel turbulent and heat transfer models that are adaptive to different spatiotemporal scales involved in the mixed convection heat transfer. In this work, direct numerical simulation of low-flow mixed convection is carried out at low-to-unitary Prandtl numbers that are of industrial interest. Time-averaged statistics, turbulent Prandtl number, as well as time signals are analyzed to investigate mixed convection phenomenon. From the time-averaged statistics, buoyant plume altered velocity boundary layer as well as the intensity of the fluctuation near both walls and channel centerline. Buoyancy effect also rendered different degree of convective heat transfer enhancement and impairment depends on Prandtl and Richardson number. Analysis of time series was conducted on the sodium mixed convection case to emphasize on the low-Pr mixed convection behavior at transition region. Resulting power spectra density and wavelet spectrogram suggests possible large convective structure in transition region. Future work will focus on providing broader data coverage on Pr-Re-Ri parameter space to facilitate more comprehensive analysis of mixed convection.

preprint2022arXiv

Toward Development of an Improved Friction Correlation for the Near-Wall Region of Pebble Bed Systems

The development of nuclear reactors that utilize pebble fuel has drastically increased the demand for improving the capabilities to simulate the packed beds found in these reactors. The complex flow fields found in a pebble bed make computational fluid dynamics (CFD) simulations time consuming and costly. Intermediate fidelity porous media models, however, are capable of approximating these flow fields in a much more computationally efficient manner. These models require the use of closures to capture the effects of complex flow phenomena without modeling them explicitly. This research employs data obtained from high-fidelity CFD simulations of a pebble bed to improve the drag closures used in porous media models in the near-wall region of the bed. Specifically, NekRS, a GPU-enabled spectral element CFD code, was used to simulate a bed of 1,568 pebbles at multiple Reynolds numbers. The case was divided into five concentric subdomains to extract radial profiles of the average porosity, velocity, and wall shear in each subdomain. A model consistent with the high-fidelity model was created in Idaho National Laboratory's Pronghorn porous media code and the KTA correlation was chosen as the drag closure of comparison. It was found that the KTA correlation overestimates the velocity in the near-wall region. An investigation of the drag coefficients between the two codes revealed that the KTA correlation underestimated the form factor in the outermost region while overestimating it in the inner four regions. This analysis in this work has revealed the underlying inaccuracy in the near-wall region of the KTA correlation and has set up the process for using high-fidelity simulation to predict more accurate drag coefficients a priori, rather than with a manual velocity-matching approach. This process will allow for the development of an improved drag closure for use in porous media models.

preprint2020arXiv

Direct Numerical Simulation of Fluid Flow in a 5x5 Square Rod Bundle Using Nek5000

Rod bundle flows are commonplace in nuclear engineering, and are present in light water reactors (LWRs) as well as other more advanced concepts. Inhomogeneities in the bundle cross section can lead to complex flow phenomena, including varying local conditions of turbulence. Despite the decades of numerical and experimental investigations regarding this topic, and the importance of elucidating the physics of the flow field, to date there are few publicly available direct numerical simulations (DNS) of the flow in multiple-pin rod bundles. Thus a multiple-pin DNS study can provide significant value toward reaching a deeper understanding of the flow physics, as well as a reference simulation for development of various reduced-resolution analysis techniques. To this end, DNS of the flow in a square 5x5 rod bundle at Reynolds number of 19,000 has been performed using the highly-parallel spectral element code Nek5000. The geometrical dimensions were representative of typical LWR fuel designs. The DNS was designed using microscales estimated from an advanced Reynolds-Averaged Navier-Stokes (RANS) model. Characteristics of the velocity field, Reynolds stresses, and anisotropy are presented in detail for various regions of the bundle. The turbulent kinetic energy budget is also presented and discussed

preprint2020arXiv

The effect of varying viscosity in turbulent channel flow

In this article we examine channel flow subject to spatially varying viscosity in the streamwise direction. The Reynolds number is imposed locally with three different ramps. The setup is reminiscent of transient channel flow, but with a space dependent viscosity rather than a time dependent viscosity. It is also relevant to various applications in nuclear engineering and in particular in test reactors, where the viscosity changes significantly in the streamwise direction, and there is a severe lack of Direct Numerical Simulation (DNS) data to benchmark turbulence models in these conditions. As part of this work we set up a novel benchmark case: the channel is extended in the stream-wise direction up to 20p. The viscosity is kept constant in the first 4p region. This inlet region is used as a cyclic region to obtain a fully developed flow profile at the beginning of the ramping region. In the ramping region the Reynolds number is linearly increased along the channel. The flow is homogenous in the spanwise direction, while it is nonhomogenous in the stream-wise and wall-normal direction. We perform here Direct Numerical Simulation (DNS) with Nek5000, a spectral-element computational fluid dynamics (CFD) code developed at Argonne National Laboratory. In this study, specific focus is given to the investigation of turbulence properties and structures in the near-wall region along the flow direction. Turbulent statistics are collected and investigated. Similarly to transient channel flow, the results show that a variation in the Reynolds across a channel does not cause an immediate change in the size of turbulent structures in the ramp region and a delay is in fact observed in both wall shear and friction Reynolds number. The results from the present study are compared with a correlation available in the literature for the friction velocity and as a function of the Reynolds number.

preprint2020arXiv

Towards Direct Numerical Simulation of a 5x5 Rod Bundle

Rod bundle flows are prevalent in nuclear engineering for both light water reactors (LWR) and advanced reactor concepts. Unlike canonical channel flow, the flow in rod bundles presents some unique characteristics, notably due to the inhomogeneous cross section which can present different local conditions of turbulence as well as localized effects characteristic of external flows. Despite the ubiquity of rod bundle flows and the decades of experimental and numerical knowledge acquired in this field, there are no publicly available direct numerical simulations (DNS) of the flow in multiple pin rod bundles with heat transfer. A multiple pin DNS study is of great value as it would allow for assessment of the reliability of various turbulence models in the presence of heat transfer, as well as allow for a deeper understanding of the flow physics. We present work towards DNS of the flow in a square 5x5 rod bundle representative of LWR fuel. We consider standard configurations as well as configurations where the central pin is replaced with a guide thimble. We perform simulations in STAR-CCM+ to design the numerical DNS, which is to be conducted using the open source spectral element code Nek5000. Large Eddy Simulations are also performed in Nek5000 to confirm that the resolution requirements are adequate. We compare results from STAR-CCM+ and Nek5000, which show very good agreement in the wide gaps with larger discrepancies in the narrow gaps. In particular, evidence of a gap vortex street is seen in the edge subchannels in LES but is not predicted by STAR-CCM+.