Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMAL CONVECTION”

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 379 records · Page 21

Nonlinear proper orthogonal decomposition for convection-dominated flows

Autoencoder techniques find increasingly common use in reduced order modeling as a means to create a latent space. This reduced order representation offers a modular data-driven modeling approach for nonlinear dynamical systems when integrated with a time series predictive model. In this Letter, we put forth a nonlinear proper orthogonal decomposition (POD) framework, which is an end-to-end Galerkin-free model combining autoencoders with long short-term memory networks for dynamics. By eliminating the projection error due to the truncation of Galerkin models, a key enabler of the proposed nonintrusive approach is the kinematic construction of a nonlinear mapping between the full-rank expansion of the POD coefficients and the latent space where the dynamics evolve. We test our framework for model reduction of a convection-dominated system, which is generally challenging for reduced order models. Our approach not only improves the accuracy, but also significantly reduces the computational cost of training and testing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Science Plan for the Deployment of the Third ARM Mobile Facility to the Southeastern United States at the Bankhead National Forest, Alabama (AMF3 BNF)

In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land-surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land-atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow-to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land-atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.

54 ENVIRONMENTAL SCIENCES↗

High-fidelity Experiments of Turbulent Buoyant Jets from Simultaneous Particle Image Velocimetry and Laser Induced Fluorescence

Accurate models of turbulent buoyant flows are essential for the design of nuclear reactors thermal hydraulics and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We performed experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous Particle Image Velocimetry and Laser Induced Fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The Refractive Index Matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy driven flows (e.g., natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as basis for the development of a new one that captures turbulent mixing.

42 ENGINEERING↗

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗

High-Fidelity Velocity and Concentration Measurements of Turbulent Buoyant Jets

Accurate models of turbulent buoyant flows are essential for the design of the cooling circuit of nuclear reactors and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We analyze experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous particle image velocimetry and laser-induced fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The refractive index matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy-driven flows (e.g. natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in a uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as a basis for the development of a new model that captures turbulent mixing.

laser induced fluorescence↗

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗

Development and Validation of Resistance-Capacitance Model (RCM) for Phase Change Material (PCM) Embedded in 3D Periodic Structures

The low thermal conductivity of Phase Change Materials (PCM) can be improved with extended surfaces such as additively manufactured 3D periodic lattice structures. Three different aluminum alloy-based lattices (base sizes 10, 20, and 40 mm) with average porosity of 0.95 filled with paraffin wax, with a nominal phase change temperature of 55°C, were experimentally investigated. In this work, a computationally efficient 2D Resistance Capacitance-based model (RCM) was developed for predicting the thermal characteristics of these geometries. Non-uniform porosity in the PCM-metal domain was estimated using image processing and served as model input. The solver does not solve for higher-order physics as in CFD but still can provide a good prediction of thermal resistance and energy storage at a very low computational cost. The simulation-to-real-time factor for this geometry is of the order of 10-4, while CFD simulations typically have a real-time factor greater than 1. The model was validated against the experimental data for melting under three different heat fluxes (6250 W/m2, 12500 W/m2, and 18750 W/m2). The mean deviation of the predicted average PCM temperature was between 1.34 K-2.81 K for different cases. The maximum average temperature deviation of 5.45 K was observed for the 20 mm geometry at the highest heat flux test condition. The effects of natural convection were neglected in the model, but the predicted PCM temperature and energy storage still showed good agreement with the experimental data.

25 ENERGY STORAGE↗

Supersaturation Variability from Scalar Mixing: Evaluation of a New Subgrid-Scale Model Using Direct Numerical Simulations of Turbulent Rayleigh–Bénard Convection

We report that supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars - water vapor and temperature - leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

54 ENVIRONMENTAL SCIENCES↗

Nanoparticle-enhanced absorptivity of copper during laser powder bed fusion

We report laser powder bed fusion (LPBF) of pure copper for thermal and electrical applications is hampered by its low near-infrared absorptivity and high thermal diffusivity. These material properties make it very difficult to localize the thermal energy needed to produce high density 3D printed parts. Modification of metal powders via nanoparticle additives is a promising approach to increasing absorptivity, but the effect of nanoparticles on absorptivity and melting behavior during LPBF is not well understood. In this study, we developed an in situ calorimetry system to measure effective absorptivity during LPBF on copper substrates. We decorated copper substrates using three nanoparticle systems (CuS, TiB 2 , multilayer graphene flakes) and demonstrated an enhanced absorptivity of the decorated substrates relative to pure copper. Graphene nanoflakes resulted in the highest improved absorption relative to pure copper from 0.09 to 0.48, due to their stability at high laser scanning powers. A thermomechanical model with convective heat transfer provided confidence in the measurements by reproducing the experimental melt pool traces. Full 3D cylindrical prints demonstrated an improvement in relative density of the copper-graphene powder prints (in the range of 0.930–0.992), relative to that of as-purchased copper powder prints (in the range of 0.854–0.972). This work provides a fundamental study of nanoparticle-enabled LPBF of highly reflective metals and demonstrates a viable route for expanding the library of reliably printable metals.

36 MATERIALS SCIENCE↗

Analytical investigation of the dew point water recovery system with potential application for drinking water production

The article presents the analytical investigation of the Dew Point Water Recovery (DPWR) system with potential application for drinking water production. Four different configurations of the DPWR system were proposed, differing in the number and type of dew point cooling exchangers used, and their operation was described using a validated convective heat transfer coefficient mathematical model. It was found that the three stage configu-ration of the DPWR system with a regenerative exchanger has the best water recovery performance, while the three stage configuration with a counter-flow exchanger has the best ratio of energy supplied to the system in relation to the amount of water produced. It was also found that the DPWR system can operate in a wide range of operating conditions with a specific electric energy consumption reaching below 0.9 kWh/m 3 and GOR in the range of 0.36-2.16. What distinguishes the DPWR system from other methods of drinking water production is the ability to work without an additional source of thermal energy apart from the energy of the air supplied to the system. Thanks to this property, the system has the potential to fill the technological gap of a low-energy and easily accessible method of drinking water production.

42 ENGINEERING↗

RELAP5-3D Simulation of PG-27 Test at the HTTF Facility

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a quarter scaled integral-effect test facility designed to examine the transient phenomena occurring in a Modular Hight Temperature Gas-cooled Reactor (MHTGR). It is specifically designed to handle the pressurized and depressurized conduction cooldown (PCC and DCC) events. The PCC event is an accident scenario in which there is a loss of forced convection of the coolant through the system. However, the pressure boundary remains intact. This study focuses on the Reactor Excursions and Leak Analysis Program (RELAP)5-3D simulation of the PG-27 test (PCC phenomena) and the sensitivity of those simulations from the material properties and primary helium mass flow rate standpoint. The aim is to capture a reasonable one-dimensional picture of a three dimensional phenomena. The simulation will help in the code assessment and in understanding the various unknown parameters involved. The presented results will demonstrate the code is able to capture the detailed one-dimensional physics accurately. However, input parameters to the code (i.e., the primary mass flow rate, refractory material thermal conductivity, and heat capacity at elevated temperatures) impact the fuel temperature the most. This work will discuss the sensitivity of simulations to the identified input parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Code Validation of SAM Using Forced and Natural Circulation Data from NACIE-UP Benchmark

Heavy liquid metals (HLMs) are promising candidates as coolants of Generation IV fast reactors due to their thermophysical properties. In the last decade, experimental work has been proposed as part of research and development efforts to develop such systems. In this context, researchers from the Brasimone Research Center have conducted many experiments using the Natural Circulation Experiment Upgrade (NACIE-UP) facility to study the thermofluid dynamic behavior of HLMs in rod bundle configurations with or without wire wrappers. This facility consists of a rectangular loop operated with lead-bismuth eutectic. Sensors across the loop monitor relevant parameters, i.e. temperatures, heat transfer, and flow conditions. Here, in the present work, we carefully select published data from NACIE-UP to validate the System Analysis Module (SAM), a modern system analysis code developed at Argonne National Laboratory. We developed one SAM model using specifications of the facility geometry and materials existing in relevant papers and reports. On top of that, these references provided the boundary conditions for simulating natural circulation and forced convection experiments in either steady or transient conditions. The SAM model simulates five test cases with diverse operating conditions. Ultimately, the code is proven to predict temperatures and mass flow rates that closely match the experiments. The discrepancies between numerical predictions and diverse transients are limited to a few degrees Celsius, showcasing that SAM is well suited for analyzing nuclear systems relying on HLM coolants.

advanced reactors↗

A computational analysis of effective R-values of buried ducts – the dynamic performance of buried ducts

Here, this paper evaluates the thermal performance of ducts partially or fully buried in loose-fill attic insulation. The overall thermal resistance between the ducts and the attic is referred to as an effective R-value. This paper shows a strong dependency of assumed attic temperature on the effective R-value. Based on the results, the effective R-value can be about twice as much with an attic temperature of 130 °F [54.4 °C], compared to when the attic temperature is 80 °F [26.7 °C]. Thus, this paper provides a polynomial regression equation based on a large set of simulations to determine the effective R-value of buried ducts depending on attic temperature and whether the HVAC system runs in cooling or heating mode. Further, the work presented in this paper investigated the potential impact of convective airflow within the attic insulation, particularly around the ducts. The analysis was based on computational fluid dynamics (CFD) and indicated that convectional forces are presented around the exterior surface of the ducts, but with negligible impact on the overall heat balance between the duct and the attic space.

97 MATHEMATICS AND COMPUTING↗

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↗

THERMS Technology Commercialization Fund Final Report

The technology commercialization fund assisted with the design and development of an intermediate scale packed be thermal energy storage system. Presented in this final report is the design, testing, and modeling of a 100 kWh th radial packed-bed. Air is used as a heat transfer fluid and 3/8” pea gravel is used as the storage medium. Testing has demonstrated the formation of a thermocline within the packed-bed. An air charging temperature of 450 °C was achieved resulting in a charging power of 26.5 kWh th . The radial packed-bed was charged for 5 hours and discharged for 5 hours. Minimal temperature variation was observed at the bottom, middle, and top axial locations of the radial packed-bed which suggests that buoyancy and hydrodynamics in the injection well had a minimal effect on the uniformity of air flow through the bed during charging and discharging. This indicates that the forced convection in a radial direction is dominating the flow regime compared to the buoyancy. The findings of this study suggest an increased air duct diameter through the air heater will increase the charging temperature and increase air mass flow through the system. Increased mass flow through the packed-bed is hypothesized to shorten the thermocline length. Modeled and measured results are compared, and it is here determined that a more comprehensive CFD model will aid in the understanding of air flow and thermocline evolution within the packed-bed.

25 ENERGY STORAGE↗