Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cavity 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 109 records · Page 6

Design Evaluation of a Next-Generation High-Temperature Particle Receiver for Concentrating Solar Thermal Applications

High-temperature particle receivers are being developed to achieve temperatures in excess of 700 °C for advanced power cycles and solar thermochemical processes. This paper describes designs and features of a falling particle receiver system that has been evaluated and tested at the National Solar Thermal Test Facility at Sandia National Laboratories. These advanced designs are intended to reduce heat losses and increase the thermal efficiency. Novel features include aperture covers, active air flow, particle flow obstructions, and optimized receiver shapes that minimize advective heat losses, increase particle curtain opacity and uniformity, and reduce cavity wall temperatures. Control systems are implemented in recent on-sun tests to maintain a desired particle outlet temperature using an automated closed-loop proportional–integral–derivative controller. These tests demonstrate the ability to achieve and maintain particle outlet temperatures approaching 800 °C with efficiencies between 60 and 90%, depending on incident power, mass flow, and environmental conditions. Lessons learned regarding the testing of design features and overall receiver operation are also presented.

Mills, Brantley H.↗

Analysis of Flow Area in the Axial Cross-Vessel Annular Duct Following a Breach of the Primary Helium Boundary of an HTGR

The lock-exchange flow phenomenon has been heavily studied for decades, including the lock-exchange phenomenon that takes place during an air ingress accident in a High-Temperature Gas-cooled Reactor (HTGR) during a break on the axial cross-vessel duct. The present study aims to characterize the gas dynamics in the cavity region of an HTGR, particularly the axial cross-vessel outer duct region. This study computes the velocity and flow area of the heavy and light currents when a break on the primary helium boundary occurs. This analysis is supported using experimental data from experiments conducted in a 1/20th length-scale HTGR. The following study summarizes the conditions of the experiment, the approach adopted to predict the flow area, and future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗

DEVELOPMENT OF PLASMA SIMULATION TOOL FOR THE ACCEL-ERATING CAVITIES

Plasma processing of superconducting radio frequency (SRF) cavities has shown an improvement in accelerating gradient by reducing the radiation due to field emission and multipacting. Plasma processing is a common tech-nique where the free oxygen produced by the plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the second-ary emission coefficient. The hydrocarbon fragments of H2, CO, CO2, and H2O are removed from the system with the process gas which is flowing through the system. Here, we present COMSOL for the first time to simulate the plasma processing of an SRF cavity. In this work, we use Jefferson Lab's C75 SRF cavities design as our case study. Using simulation, we predict the condition of plasma ignition inside the SRF cavity. The simulation provides information about the optimal rf coupling to the cavity, mode for plasma ignition, choice of gas concen-tration, power, and pressure.

Senevirathne, I.↗

Development of a Reference Model for Molten-Salt-Cooled Pebble-Bed Reactor Using SAM

To support the development and utilization of the SAM code for fluoride-salt-cooled high-temperature pebble-bed reactor (PB-FHR) safety analysis and licensing, an effort was devoted to developing reference models for PB-FHR primary loop and the reactor cavity cooling system (RCCS). A reference standard problem of a prototypical reactor design is foundational to NRC and industry to verify the adequacy of computer codes and evaluation models for a specific reactor type. The SAM code was enhanced for the reference PB-FHR model development, including a 1D-3D flow coupling scheme and conjugate heat transfer between porous media and solid structures. The reference FHR primary loop model was developed based on the UC Berkeley Mk1 FHR design with additional design features from the Kairos Power’s KP-FHR core. The PB-FHR core is modeled by a two-dimensional porous medium while the rest of the primary loop is represented by a SAM one-dimensional model. A SAM one-dimensional primary loop model was also developed for comparison using a single channel approach for the reactor core. Furthermore, the water-based NSTF cavity was modeled using SAM multi-dimensional flow module, to assist in evaluating SAM capabilities of modeling the emergency heat removal systems relying on RCCS concepts.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Measurement of atomic oxygen densities using TALIF on a dielectric barrier discharge: insights into the volume above a micro cavity plasma array

Dielectric barrier discharges, particularly micro cavity plasma arrays, offer significant potential for plasma-catalytic research due to their ability to ignite plasma in direct contact with a catalytic surface, enabling the observation of plasma-surface interactions. A key factor in their application is the generation of reactive species, such as atomic oxygen, within the cavities. These species can interact with both the surface (e.g. for activation or cleaning) and the gas being treated (e.g. for oxidation). Given the central role of oxygen atoms in plasma catalysis and their use as a model for more complex species, this work investigates the transport of these atoms out of the cavities. Two-photon absorption laser-induced fluorescence spectroscopy with picosecond laser excitation is performed in the volume above the cavities. The results are compared with a basic diffusion model. The reactor operates with a He/O 2 mixture at a flow rate of 1 slm and atmospheric pressure. Densities of up to 10 16 cm -3 are measured near the surface. Time-dependent measurements show that, at a distance of 350 µm from the surface, a density equilibrium is reached within less than 3 ms of reactor operation. Decay times due to ozone formation after the reactor is turned off are on a similar scale. Spatially resolved measurements show that the oxygen density decreases exponentially from the surface but remains detectable up to approximately 1 mm above the surface, indicating significant application potential. Variations in the O 2 admixture show a density maximum at 0.4%, confirming previous helium state enhanced actinometry measurements within the cavities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Flux Spreading in a Light-Trapping, Planar-Cavity Receiver for Enclosed Solar Particle Heating

Concentrating solar thermal power (CSP) technology development has recently focused on increasing the operating temperatures to accommodate high efficiency power cycles and thermochemical processes. Inert solid particles as heat transfer media enable solar receivers to operate above 700 degrees Celsius resulting in increased system thermal efficiency compared to the conventional molten salt based CSP system. An open-cavity falling-particle solar receiver that can efficiently heat particles by direct heating from concentrated solar radiation faces challenges with large particle losses from wind and unable to support thermochemical reactions. A light-trapping, planar cavity reiver (LTPCR) where particles are indirectly heated can significantly minimize the particle losses during the operation, support thermochemical reactions, and offer scalability potential. The LTPCR features an array of vertical planar receiver/absorber panels arranged within a cavity configuration. Concentrated solar radiation from heliostats is focused onto the receiver walls, where heat is indirectly transferred to solid particles flowing inside the receiver channels. Heat transfer occurs through direct contact between the receiver panel walls and particles, and can be enhanced by fluidizing particles with air. This fluidization increases particle-wall contact and extends particle residence time, maximizing heat transfer efficiency. The unique vertical planar receiver structure originated from a near-blackbody tubular light absorber, effectively distributing the incoming solar beam spread across the panel walls and trapping light. This flux spreading effect, driven by cosine projection, converts high incident solar flux into a lower, more uniform heat flux on the panel walls. This redistribution enhances heat transfer efficiency between particle-wall or reaction gases-wall, while preventing localized overheating of the receiver panel. Indirect planar cavity solar receivers completely separate solid particles from the ambient environment that can greatly reduce the thermal losses in heated particles resulting in high efficiency at high temperatures above 700 degrees Celsius. This design ensures no particle losses to the environment during the operation while open-cavity designs can experience significant particle losses from wind. An experimental investigation was conducted to observe flux spreading on the receiver panel wall. A lab-scale prototype planar receiver, fabricated using Haynes 230 alloy, was tested under direct concentrated solar radiation using the high-flux solar furnace (HFSF) facility at NREL. The experiment was performed under normal peak radiative heat fluxes ranging from 800 to 1900 kW/m2. A temperature distribution on the panel wall was measured using a thermal imaging camera (FLIR A 6600). To prevent overheating at the receiver front tip, prism-shaped heat shields (Zircar UNIFROM C1) were placed in front of the receiver, and their influence on flux spreading was also studied. Absorbed flux distribution on the panel wall was modeled using SolTrace. The total solar power and flux distributions delivered from HFSF were determined based on the heliostat mirror optical properties, direct normal irradiance (DNI) on the on-sun testing days, peak flux measurement during the on-sun testing, and shutter/attenuator settings Due to the large incident angles of the solar beam on the panel wall, the angular optical properties of Haynes 230 alloy and Zircar heat shields were incorporated into the model. This flux distribution model was then integrated into a computational fluid dynamics (CFD) simulation to predict the receiver panel wall temperature, which was compared with the experimental measurements. Both prediction and measurements identified a temperature hotspot at the backside of the panel, indicating that the incident solar beam can fully reach to the rear of the receiver. The heat shields positioned at the front of the receiver effectively reduced the excessive temperature rise at the receiver front tip. Overall, the temperature was well distributed over the panel wall, with a minor hotspot at the back of the receiver. The model slightly overpredicted the temperature, possibly due to discrepancies in optical properties of the panel and an underprediction of thermal loss in the receiver. The advancement of the particle LTPCR offers a viable alternative to open-cavity receivers by addressing particle loss issues. Additionally, it presents a pathway for enabling solar thermochemical processes, extending CSP technology beyond power generation to fuel and chemical production.

14 SOLAR ENERGY↗

Segmented annular combustion system

An annular combustion system includes a fuel nozzle, a panel fuel injector including a fuel plenum and at least one premixing channel therein, an inner liner including an inner band radially spaced from an outer band and an inner flow annulus formed therebetween, and an outer liner radially spaced from the inner liner and comprising an inner band radially spaced from an outer band and an outer flow annulus formed therebetween. The fuel nozzle is positioned radially between the inner liner and the outer liner, and the panel fuel injector extends radially between the inner liner and the outer liner. At least one of the inner flow passage and the outer flow passage is in fluid communication with at least one of a premix air plenum, a cooling air cavity and a cooling air plenum defined within the panel fuel injector.

42 ENGINEERING↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Laser-driven, ion-scale magnetospheres in laboratory plasmas. II. Particle-in-cell simulations

Ion-scale magnetospheres have been observed around comets, weakly magnetized asteroids, and localized regions on the Moon and provide a unique environment to study kinetic-scale plasma physics, in particular in the collision-less regime. In this work, we present the results of particle-in-cell simulations that replicate recent experiments on the large plasma device at the University of California, Los Angeles. Using high-repetition rate lasers, ion-scale magnetospheres were created to drive a plasma flow into a dipolar magnetic field embedded in a uniform background magnetic field. The simulations are employed to evolve idealized 2D configurations of the experiments, study highly resolved, volumetric datasets, and determine the magnetospheric structure, magnetopause location, and kinetic-scale structures of the plasma current distribution. We show the formation of a magnetic cavity and a magnetic compression in the magnetospheric region, and two main current structures in the dayside of the magnetic obstacle: the diamagnetic current, supported by the driver plasma flow, and the current associated with the magnetopause, supported by both the background and driver plasmas with some time-dependence. From multiple parameter scans, we show a reflection of the magnetic compression, bounded by the length of the driver plasma, and a higher separation of the main current structures for lower dipolar magnetic moments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗

Kinematic Evidence for Bipolar Ejecta Flows in the Galactic Supernova Remnant W49B

W49B is a unique Galactic supernova remnant with centrally peaked, “bar”-like ejecta distribution, which was once considered evidence for a hypernova origin that resulted in a bipolar ejection of the stellar core. However, chemical abundance measurements contradict this interpretation. Closely connected to the morphology of the ejecta is its velocity distribution, which provides critical details for understanding the explosion mechanism. We report the first ever observational constraint on the kinematics of the ejecta in W49B using the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM). Using XRISM/Resolve, we measured the line-of-sight velocity traced by the Fe Heα emission, which is the brightest feature in the Resolve spectrum, to vary by ±300 km s −1 with a smooth east-to-west gradient of a few tens of kilometers per second per parsec along the major axis. Similar trends in the line-of-sight velocity structure were found for other Fe-group elements Cr and Mn, traced by the Heα emission, and also for intermediate-mass elements Si, S, Ar, and Ca, traced by the Lyα emission. The discovery of the east–west gradient in the line-of-sight velocity, together with the absence of a twin-peaked line profile or enhanced broadening in the central region, clearly rejects the equatorially expanding disk model. In contrast, the observed velocity structure suggests bipolar flows reminiscent of a bipolar explosion scenario. An alternative scenario would be a collimation of the ejecta by an elongated cavity sculpted by bipolar stellar winds.

High resolution spectroscopy↗

Bathymetric Influences on Antarctic Ice-Shelf Melt Rates

Ocean bathymetry exerts a strong control on ice sheet-ocean interactions within Antarctic ice-shelf cavities, where it can limit the access of warm, dense water at depth to the underside of floating ice shelves. However, ocean bathymetry is challenging to measure within or close to ice-shelf cavities. It remains unclear how uncertainty in existing bathymetry datasets affect simulated sub-ice-shelf melt rates. In this work, we infer linear sensitivities of ice-shelf melt rates to bathymetric shape with grid-scale detail by means of the adjoint of an ocean general circulation model. Both idealized and realistic-geometry experiments of sub-ice-shelf cavities in West Antarctica reveal that bathymetry has a strong impact on melt in localized regions such as topographic obstacles to flow. Moreover, response of melt to bathymetric perturbation is found to be non-monotonic, with deepening leading to either increased or decreased melt depending on location. Our computational approach provides a comprehensive way of identifying regions, where refined knowledge of bathymetry is most impactful, and also where bathymetric errors have relatively little effect on modeled ice sheet-ocean interactions.

54 ENVIRONMENTAL SCIENCES↗

High blowing ratio high effectiveness film cooling configurations

A film cooled component is provided that includes a body having a cavity defined therein, wherein a coolant is supplied to the cavity cooling apertures having at least one channel defined on an inside surface. The coolant exits the cavity through the plurality of cooling apertures thereby providing a cooling effect to the body and the channel creates an air flow such that the air coolant moves towards the body upon exiting the plurality of cooling apertures. The film cooled component may include cooling apertures configured to reverse blow coolant air.

42 ENGINEERING↗

Precision-controlled ultrafast electron microscope platforms. A case study: Multiple-order coherent phonon dynamics in 1T-TaSe2 probed at 50 fs–10 fm scales

We report on the first detailed beam tests attesting the fundamental principle behind the development of high-current-efficiency ultrafast electron microscope systems where a radio frequency (RF) cavity is incorporated as a condenser lens in the beam delivery system. To allow for the experiment to be carried out with a sufficient resolution to probe the performance at the emittance floor, a new cascade loop RF controller system is developed to reduce the RF noise floor. Temporal resolution at 50 fs in full-width-at-half-maximum and detection sensitivity better than 1% are demonstrated on exfoliated 1T-TaSe2 system under a moderate repetition rate. To benchmark the performance, multi-terahertz edge-mode coherent phonon excitation is employed as the standard candle. The high temporal resolution and the significant visibility to very low dynamical contrast in diffraction signals via high-precision phase-space manipulation give strong support to the working principle for the new high-brightness femtosecond electron microscope systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of self-magnetic pinch radiographic diode performance on RITS-6 at Sandia National Laboratories. II. Coupling between the inductive voltage adder and the SMP load

The self-magnetic pinch (SMP) diode is a type of radiographic diode used to generate an intense electron beam for radiographic applications. At Sandia National Laboratories, SMP was the diode load for the six-cavity radiographic integrated test stand inductive voltage adder (IVA) driver operated in a magnetically insulated transmission line (MITL). The MITL contributes a flow current in addition to the current generated within the diode itself. Extensive experiments with a MITL of 40Ω load impedance [T. J. Renk et al., Phys. Plasmas 29, 023105 (2022)] indicate that the additional flow current leads to results similar to what might be expected from a conventional high-voltage interface driver, where flow current is not present. However, when the MITL flow impedance was increased to 80 Ω, qualitatively different diode behavior was observed. This includes large retrapping waves suggestive of an initial coupling to low impedance as well as diode current decreasing with time even as the total current does not. A key observation is that the driver generates total current (flow + diode) consistent with the flow impedance of the MITL used. The case is made in this paper that the 80Ω MITL experiments detailed here can only be understood when the IVA-MITL-SMP diode is considered as a total system. Additionally, the constraint of fixed total current plus the relatively high flow impedance limits the ability of the diode (whether SMP or other type) to act as an independent load. An unexpected new result is that in tracking the behavior of the electron strike angle on the converter as a function of time, we observed that the conventional cIV x “Radiographic” radiation scaling (where x~2.2) begins to break down for voltages above 8 MV, and cubic scaling is required to recover accurate angle tracking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ocean Variability Beneath the Filchner-Ronne Ice Shelf Inferred From Basal Melt Rate Time Series

Fourteen phase-sensitive radars (ApRES) were deployed on the Filchner-Ronne Ice Shelf (FRIS) to measure variability in its basal melt rate. Melt rates from sites along the Ronne Depression vary seasonally, consistent with the dynamics of the propagation of seasonal dense water from the western ice front into the cavity. Several sites at the back of the FRIS cavity feature a signal with two seasonal maxima. Sub-ice shelf oceanographic data available from one of the sites indicate that this signal is caused by two different pathways followed by the same source water. Inter-annual variability is strongest along a direct flow pathway between western Ronne Ice Front and western Berkner Island. Highest melting occurred in 1999 and 2018, following anomalously low summer sea-ice concentrations in front of the ice shelf. Inter-annual melt rate variability at the back of the FRIS cavity is limited. If present, it is expressed as a suppression or delay in the arrival of the seasonal melt rate minimum, which can be understood in terms of inter-annual stratification changes and variable inflow pathways toward the sites. Long term mean ApRES melt rates agree with estimates from satellite data over eastern FRIS. However, the satellite estimates overstate the area of active basal freezing in the western part of the ice shelf. Furthermore, the temporal melt rate variability from the satellite estimates exaggerates the range of variability at both seasonal and inter-annual time scales with any correspondence between the in-situ and remotely derived inter-annual variability being limited to a single site.

54 ENVIRONMENTAL SCIENCES↗

A Simplified Numerical Approach to Characterize the Thermal Response of a Moving Bed Solar Reactor

Abstract Concentrated solar thermochemical storage in the form of a zero-emission fuel is a promising option to produce long-duration energy storage. Solar fuel is produced using a cavity reactor that captures concentrated solar radiation from a solar field of heliostats. In this paper, heat transfer model of a tubular plug-flow reactor designed and manufactured for a solar fuel production is presented. Experimental data collected from a fixed bed tubular reactor testing are used for model comparison. The system consists of an externally heated tube with counter-current flowing gas and moving solid particles as the heated media. The proposed model simulates the dynamic behavior of temperature profiles of the tube wall, gas, and particles under various gas flowrates and residence times. The heat transfer between gas–wall, solid particle–wall, and gas–solid particle is numerically studied. The model results are compared with the results of experiments done using a 4 kW furnace with a 150 mm heating zone surrounding a horizontal alumina tube (reactor) with 50.8 mm outer diameter and thickness of 3.175 mm. Solid fixed particles of magnesium manganese oxide (MgMn2O4) with the size of 1 mm are packed within the length of 250 mm at the center of the tube length. Simulation results are assessed with respect to fixed bed experimental data for four different gas flowrates, namely, 5, 10, 15, and 20 standard liters per minute of air, and furnace temperatures in the range of 200–1200 °C. The simulation results showed good agreement with maximum steady state error that is less than 6% of those obtained from the experiments for all runs. The proposed model can be implemented as a low-order physical model for the control of temperature inside plug-flow reactors for thermochemical energy storage applications.

Engineering↗