Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONVECTIVE 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

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamics of Downdrafts Around a Growing Convective Cloud: A Numerical Study

We examine the dynamics of cloud-edge downdrafts over the growth phase of isolated cumuli, combining Eulerian and Lagrangian analyses. As in previous studies, our results show that growing cumuli are surrounded by downdrafts linked to cloud-scale quasi-toroidal circulations at all times at middle and upper cloud levels consistent with the thermal chain description of convective clouds. These toroidal circulations are responsible for the most intense cloud-edge downdrafts in our simulations. In the upper cloud half, roughly 30%–50% of the upward mass flux is typically compensated within a radius of about twice the updraft radius in quasi-laminar simulations forced by a warm bubble in an initially quiescent flow. In a turbulent cloud forced by surface fluxes, this compensation fraction is around 10%–30% over the same region. In contrast to the buoyancy-centered view of subsiding shells, Eulerian and Lagrangian vertical momentum budget analyses show that the most intense cloud-edge downdrafts in the turbulent setup, and after spin-up of the toroidal circulation in the quasi-laminar experiments, are predominantly mechanically forced (i.e., driven by dynamic pressure accelerations). This is consistent throughout the entire growth phase of the cumulus clouds and across tests with varying assumptions, including drier and moister environments. Despite dynamic pressure perturbations being the main driver of toroidal downdrafts, the downdraft speed (relative to the corresponding updraft velocity) exceeds the prediction of the non-buoyant Hill's spherical vortex—a simple model frequently used for cloud-scale circulations—by more than 30%.

Pardo, Lianet Hernández [Goethe Univ., Frankfurt (↗

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Demonstration of gold nanorod systems for enhanced total efficiency: Experimental and numerical analysis

This study investigates the photothermal performance of gold nanorods engineered to exhibit longitudinal plasmon resonances at 695 nm, 780 nm, and 970 nm. The work combines synthesis, structural characterization, extinction measurements, numerical modeling, and controlled temperature experiments to quantify how nanorod geometry, resonance tuning, concentration, and chamber shape jointly influence heat generation. Transmission electron microscopy confirms that increasing nanorod aspect ratio systematically shifts the longitudinal plasmon peak toward the near-infrared region. Extinction measurements show strong agreement with theoretical predictions, validating the numerical model across two independent datasets. Three chamber geometries were tested under laser excitation at 640 nm, 808 nm, and 980 nm: an ascending stepped base, a flat base, and a descending stepped base. Without nanorods, the ascending geometry produced the highest efficiency due to enhanced natural convection. After introducing gold nanorods, all geometries exhibited substantial thermal enhancement, with total efficiencies exceeding 20%. The strongest improvement was obtained for nanorods resonant at 780 nm with a mass concentration of 4.6 mg/mL implemented on the descending stepped-base geometry. This performance resulted from the combined effect of spectral overlapping with the 808 nm laser, the highest nanorod concentration, and localized heat accumulation that intensified buoyancy-driven flow. The findings demonstrate that total efficiency is governed by a synergistic interplay between optical resonance, nanoparticle concentration, and macroscopic chamber design, revealing the system-level coupling between nanoscale plasmonic absorption and macroscale heat-transfer phenomena. The results provide a validated framework for tuning nanoscale plasmonic absorbers and optimizing thermal systems for applications requiring efficient light-to-heat conversion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Computational Fluid Dynamics Analysis of the Molten Salt Tritium Transport Experiment Test Section

Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of Cathode Cooling in Three-Dimensional Simulations of an Atmospheric Pressure Glow Discharge

The Atmospheric Pressure Glow Discharge (APGD) is a relatively simple and versatile plasma source used in a wide range of applications. Active cooling of the cathode can effectively mitigate instabilities, leading to glow-to-arc transitions. This study investigates the effect of varying the degree of cathode cooling in APGD with a planar cathode in helium. The plasma flow model incorporates mass conservation, chemical species transport, momentum conservation, conservation of thermal energy of heavy species and of electrons, and electrostatics. The model is applied to time-dependent simulations through a three-dimensional computational domain describing the whole discharge, without geometric symmetry or steady-state assumptions. Simulations of an experimentally characterized APGD explore the effects of electric current and cathode cooling—ranging from thermally insulated to extreme convective cooling. Results show the formation of an annular region with high electric field over the cathode surface under conditions of high current and low cooling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

Multi-physics Modeling of Radiative Heat Transfer and Fluid Flow for the Reactor Cavity Cooling System

High-temperature gas-cooled reactors (HTGRs) are notable for their high thermal efficiency and potential for combined heat and power applications. These reactors are particularly appealing due to their advanced passive safety features. HTGRs utilize passive safety systems that function without requiring active components like pumps or compressors during emergencies. These reactor designs depend on a Reactor Cavity Cooling System (RCCS) to manage decay heat removal from the reactor pressure vessel (RPV) during accident conditions. The RCCS consists of vertical rectangular channels known as "risers" or riser ducts positioned around the RPV. These risers receive heat from the RPV through both convective and radiative heat transfer mechanisms. Understanding the interplay of multiple physical phenomena, such as fluid dynamics, heat transfer, and neutron interactions, is essential for the effective design and operation of nuclear reactors, particularly for systems like the RCCS. Multi-physics simulations provide a comprehensive approach to studying these interactions, offering detailed insights and enhancing accuracy. They are especially important in RCCS designs, where the interaction between radiative and convective heat transfer can significantly impact system performance. By leveraging multi-physics simulations, complex reactor behaviors can be modeled without compromising the fidelity of the underlying physical processes. This work aims to establish a robust methodology for coupling multiple physical processes in an air-cooled RCCS. By focusing on validating this multi-physics approach, the study involves designing test cases that simulate various conditions to verify the numerical models employed. The outcomes of this research will provide critical insights for accurately modeling and optimizing complex nuclear systems like the RCCS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal Performance of Triply Periodic Minimal Surface Lattice Structures in Single-Phase Dielectric Fluid Cooling of Power Electronics

Additive manufacturing has transformed thermal management by enabling the production of complex, optimized geometries that conventional manufacturing methods cannot achieve. This study investigates the single-phase convective heat transfer performance of gyroid triply periodic minimal surface (TPMS) lattice structures with functional porosity. TPMS structures provide high surface area to volume ratios and are amenable to 3D printing. A gyroid numerical model was created and validated against an existing experimental study with a similar feature size to the investigated geometries. The TPMS structure has a periodic width of 1.6 mm, a length of 10 mm, and a height of 4 mm, with a functional porosity ranging from 0.5 to 0.8, decreasing with distance from the heated surface. Three different flow configurations were examined for an inlet fluid temperature of 70 °C. The inlet velocities range from 0.01 to 1.2 m/s, corresponding to a Reynolds number range of 10–900 with a heat flux of 50 W/cm 2 applied at the base. AmpCool ® AC-110 dielectric fluid (Prandtl number 59.5) was used as the coolant. Thermal performance and friction characteristics were studied for the three flow orientations. The parallel flow configuration was identified as the most efficient for heat removal. A detailed analysis of the numerical results highlights the underlying physics behind the thermal performance differences among the flow configurations.

33 ADVANCED PROPULSION SYSTEMS↗

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗

Enhancing particle-to-sCO2 heat exchanger effectiveness through novel high-porosity metallic foams

In this project, a novel concept of Octet-shape based lattice frame material was developed for deployment in particle-to-sCO2 heat exchangers for CSP application. A comprehensive experimental and numerical program was performed to investigate and characterize the thermal and flow transport behavior in high porosity lattice frame material. The thermal properties include effective thermal conductivity and interstitial heat transfer coefficient and flow properties including permeability and inertial coefficient. Above quantities were determined for both air flow and heated particle flow through the porous channels which were representative of the hot channel in a typical counter flow heat exchanger. The cold side of the heat exchanger feature minichannels which are known to be very efficient in thermal transport. The goal of this project was to enhance the overall heat transfer coefficient (Up) on the hot side by introducing Octet-shape based lattice frame material through enhancement in effective thermal conductivity of the channel and through enhancement in the interstitial heat transfer coefficient.

14 SOLAR ENERGY↗

Ekman-driven salt transport as a key mechanism for open-ocean polynya formation at Maud Rise

Open-ocean polynyas formed over the Maud Rise, in the Weddell Sea, during the winters of 2016–2017. Such polynyas are rare events in the Southern Ocean and are associated with deep convection, affecting regional carbon and heat budgets. Using an ocean state estimate, we found that during 2017, early sea ice melting occurred in response to enhanced vertical mixing of heat, which was accompanied by mixing of salt. The melting sea ice compensated for the vertically mixed salt, resulting in a net buoyancy gain. An additional salt input was then necessary to destabilize the upper ocean. This came from a hitherto unexplored polynya-formation mechanism: an Ekman transport of salt across a jet girdling the northern flank of the Maud Rise. Such transport was driven by intensified eastward surface stresses during 2015–2018. Our results illustrate how highly localized interactions between wind, ocean flow and topography can trigger polynya formation in the open Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

Experimental investigation of Multi-Mode heat transfer to a Free-Falling dilute particle cloud in a heated vertical tube

The development of dilute particle heat exchangers and reactors for advanced energy systems requires an understanding of the multi-mode heat transfer from a heated wall to falling particles. This study presents experimental results of the overall heat transfer coefficient for a free-falling, dilute flow of particles with solid volume fraction from 0.0005 to 0.006 corresponding to feed rates from 3.7 kg s -1 m -2 to 44 kg s -1 m -2 in a vertical, heated tube containing quiescent air at atmospheric pressure. Tube wall temperatures are varied between 300°C to 900°C while keeping the particle inlet temperature constant. The experimental results show that the overall heat transfer coefficient is a strong function of particle feed rate and surface temperature. Good agreement was obtained with prior studies conducted at comparable temperatures but lower particle feed rates (< 4 kg m -2 s -1 ). The established correlations for particle-to-wall radiation and particle-to-gas convection were used to estimate the wall-to-gas convective contribution from the measured overall heat transfer coefficient. The experimental results indicated a 4 to 6 times improvement in the wall convection in the solid-gas mixture compared to that expected from natural convection in a single-phase gas. Furthermore, the data presented here are applicable to characterize heat transfer in dilute particle heat exchangers, furnaces, and solar receivers.

14 SOLAR ENERGY↗