Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat flow”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2022

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions. Part I: Correlation Development

Abstract High-temperature supercritical CO2 Brayton cycles are promising possibilities for future stationary power generation and hybrid electric propulsion applications. Heat exchangers are critical components in supercritical CO2 thermal cycles and require accurate correlations and comprehensive performance modeling under extreme temperatures and pressures. In this paper (Part I), new Colburn and friction factor correlations are developed to quantify shell-side heat transfer and friction characteristics of flow within heat exchangers in the shell-and-tube configuration. Using experimental and computational fluid dynamics (CFD) data sets from existing literature, multivariate regression analysis is conducted to achieve correlations that capture the effect of multiple critical geometric parameters. These correlations offer superior accuracy and versatility as compared to previous studies and predict the thermohydraulic performance of about 90% of the existing experimental and CFD data within ±15%. Supplementary thermohydraulic performance data are acquired from CFD simulations with supercritical CO2 as working fluid to validate the developed correlations and demonstrate its capability to be applied to supercrtical CO2 heat exchangers.

Engineering↗
DOE OSTI2020

EGS Collab Modeling and Simulation Working Group Teleconference Series (1-98)

This submission contains the presentation slides and recordings from the first 98 EGS Collab Modeling and Simulation Working Group teleconferences. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Dynamic Behavior of Oval-Twisted Helical Tube Heat Exchanger: Numerical Study with RELAP5-3D

Convective heat transfer characteristics and theoretical thermal stress behaviors are numerically calculated using RELAP5-3D for the helical-coiled once-through steam generator (H-OTSG) and the novel heat exchanger design known as the oval-twisted helically coiled heat exchanger (OTHCHX) under (1) fluctuating wall temperature conditions, (2) square-wave pulsating flow conditions, and (3) the combined effects of fluctuating wall temperature and square-wave pulsating flow conditions. Heat transfer coefficient models for the H-OTSG and OTHCHX were developed based on existing data and implemented into RELAP5-3D, successfully capturing the N⁢uavg behavior within 8% to 10% of the reported data. Under fluctuating wall temperature conditions, the OTHCHX displayed higher N⁢u avg behavior than the H-OTSG. As 𝑓 increased, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ decreased. The $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ was higher for the OTHCHX than for the H-OTSG under fluctuating wall temperature conditions. Under pulsating flow conditions, the H-OTSG and OTHCHX displayed much higher 𝑁⁢𝑢 𝑎𝑣𝑔 than under constant flow conditions. The H-OTSG displayed a higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ over the OTHCHX. Under combined fluctuating wall temperature and pulsating flow conditions, the augmented heat transfer behavior from the pulsating flow was counteracted by the wall temperature fluctuations, producing slightly higher 𝑁⁢𝑢 𝑎𝑣𝑔 over constant wall temperature, constant flow conditions, but much lower than only constant pulsating flow under constant wall temperature conditions. The effects of simultaneous wall temperature fluctuations and square-wave pulsating flow caused higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ than that of only wall temperature fluctuations or pulsating flow. As the Reynolds number (Re) increased, $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ increased. However, when 𝑓=𝑓$_{\dot{m}}$, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ showed decreasing values as Re increased. In conclusion, the results indicate that thermal-fluid resonance can help mitigate thermal stresses.

Thermal stress↗
DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗
DOE OSTI2024

Experimental analysis of convective drying of paper and board

Conventional multi-cylinder drying of paper and board involves a mixture of conductive drying from steam-heated dryer cylinders and convective drying by the flow of heated air over the surface of the paper web in the pockets. Pocket ventilation is a critical component in assisting heat and mass transfer during the drying process but is the primary contributor towards removing evaporated water from the web. Air temperature, velocity, and humidity are critical parameters involved in the convective drying process. This paper covers an experimental study involving the design and development of a small lab-scale setup for convective drying of various grades of paper and board, monitoring multiple parameters like paper temperature, moisture content, air humidity, temperature, and velocity measured in situ as the drying proceeds with continuous and accurate sampling capabilities for all parameters in the sample and the system. Instantaneous drying rates, heat, and mass transfer coefficients were also deduced for every time step till the paper completely dried. Furthermore, the coefficients obtained were also reported in the form of dimensionless correlations, and the results were compared against traditional correlations used in the modelling of paper drying. Furthermore, this data will be useful in process development, modelling, design, and the paper drying process simulation.

42 ENGINEERING↗
DOE OSTI2021

Systems and methods for expanding flow in a waste heat recovery system

A waste heat recovery system includes a first heat exchanger, a second heat exchanger, and an expander. The first heat exchanger receives working fluid from a first portion of a first loop and provides the working fluid to a second portion of the first loop. The second heat exchanger receives the working fluid from a first portion of a second loop and provides the working fluid to a second portion of the second loop. The expander provides the working fluid to a first portion of a common line. The expander includes a stator. The stator includes a first inlet and a second inlet. The common line provides the working fluid to both the first loop and the second loop upstream of the first portion of the first loop and upstream of the first portion of the second loop.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

An Experimental Investigation of Sintered Particle Effect on Heat Transfer Performance in an “Annular Flow” Evaporation Tube

Wicking structures have been widely used within passive heat transfer devices with high heat fluxes, such as heat pipes, to enhance their thermal performance. While wicking structures promote capillary pumping of the working fluid and thin film evaporation, they can result in capillary evaporation and further enhance the evaporation heat transfer. In this study, a 0.5 mm thick layer of 105 µm sintered copper particles was added to the inner wall of a copper tube, aiming to form an “annular flow” and enhance the heat transfer characteristics by taking advantage of thin film and capillary evaporation. Acetone was chosen as the working fluid, and the performance of an evaporation tube was tested for power inputs of 10, 30, 50, and 70 W. For each power input, trials were run at inclination angles varying from -90 deg to 90 deg to investigate the capillary effects. The temperature measurements showed that the temperature distribution along the evaporation tube is always downward sloping, meaning the temperature at the fluid inlet is larger than the outlet. In conclusion, results show that an “annular flow” formed by a thin layer of sintered particles can promote thin film and capillary evaporation and, therefore, boost the evaporation heat transfer coefficient.

42 ENGINEERING↗
DOE OSTI2020

Elucidating the Impact of Flow on Material-Sensitive Critical Heat Flux and Boiling Heat Transfer Coefficients: An Experimental Study with Various Materials

In this work, steady-state internal flow boiling critical heat flux (CHF) experiments were carried out at 10°C inlet subcooling and atmospheric pressure (84 kPa at an altitude of about 5300 feet) with mass flux ranging from 200 to 2,000 kg/m 2 -s on seven tube materials, including the accident-tolerant fuel cladding candidate FeCrAl (C26M) alloy. At a relatively low flow rate, the sensitivity of CHF to different materials was appreciable (the average difference among the tested materials was 19% at 200 kg/m 2 -s). As the absolute CHF differences stayed constant for higher flow rates, the increased flow rate led to reduced relative CHF differences among the tested materials (average relative difference of 15% at 1,000 kg/m 2 -s). FeCrAl alloy claddings that were oxidized under simulated Light Water Reactor (LWR) conditions yielded CHF comparable to as-received FeCrAl alloy cladding, with less than 1.8% difference. Experimental CHF results imply marginal material sensitivity impacts on the departure from nucleate boiling ratio in steady-state nuclear reactor cores. Yet, for accident progression analysis, material sensitivity may hold importance, as boiling behavior variation at low mass flux is appreciable. Macroscopic material parameters (surface wettability, surface roughness, and thermal properties) do not provide a suitable explanation for observed material-sensitive CHFs. Yet, surface roughness and thermal effusivity revealed power and linear relations, respectively, to the nucleate boiling heat transfer coefficients at low mass flux. At high mass flux, material sensitivity on nucleate boiling heat transfer significantly decreased due to the increased shear stress on boiling surface overwhelming material-related boiling implications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Heat Transfer Characteristics of Particle Flow Through Additively Manufactured (SS 316L) Lattice Frame Material Based on Octet-Shape Topology

Abstract This paper presents experimentally obtained heat transfer characteristics of additively manufactured lattice frame material based on Octet-shaped unit cell. Binder jetting technology was used to print lattices in Stainless Steel 316L material. Lattice porosities ranging from 0.75 to 0.9 were investigated where thermal transport characteristics were obtained for void space occupied by air and particles. Particle diameters were varied from 266–966 microns. Effective thermal conductivity and averaged heat transfer coefficient was calculated through steady-state experiments. It was found that presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200–400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape based lattice frame material for particle flow through them makes it a very promising candidate for particle-to-supercritical carbon dioxide (sCO2) heat exchanger in concentrating solar power (CSP) application.

Aider, Youssef↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Process intensification approach to enhancing heat and mass transfer: Radio frequency (RF) assisted drying of paper and board

Conventional multi-cylinder drying of paper and board typically relies on both conductive drying from steam-heated cylinders and convective drying, where heated air flows over the paper web surface. Conduction primarily contributes to heat transfer, while convection is the main driver of mass transfer. However, conventional drying systems are heavily dependent on steam, usually powered by fossil fuels, and are often energy-inefficient with high levels of waste. Additionally, these surface-driven processes result in a lower percentage of energy absorption compared to the energy supplied, leading to significant energy losses. To improve this long-standing process, an experimental system was developed to investigate a process intensification approach involving the integration of Radio Frequency (RF) heating, a volumetric electromagnetic technology, alongside traditional conduction and convection drying methods. This study also emphasizes the use of sensors to continuously monitor key parameters such as moisture content, supply system temperatures, sample temperatures, air flows, and drying rates in real-time. The effect of RF as an auxiliary energy source in localized environments at varying moisture levels was explored to optimize industrial drying systems, quantify potential improvements, and provide insights for future studies. Experimental results from trials combining RF with convection and with the base case alternating conduction-convection drying processes are presented. It is shown that RF is a viable process intensification approach for paper drying improving the drying rate and energy intensity at higher moisture contents. These findings offer valuable insights for process intensification and contribute to the process development, modeling, and simulation of advanced paper drying techniques.

42 ENGINEERING↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗