Search NASASearch

SEARCH · Search NASA

Results for “Numerical modelling”

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 91 records · Page 5

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham

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

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Solar Decarbonization of Paraffin Dehydrogenation Through Particle Heat Carriers (Final Technical Report)

This project focuses on solutions to decarbonize high-temperature catalytic processes using solar thermal heat. The primary project goal is to show the validity of a moving packed bed reactor for propane dehydrogenation using catalyst particles as the heat carrier for the reaction, which can be heated by concentrated solar energy in a particle receiver. This concept, if further developed, may provide a cost-effective pathway for converting lower value gases to important chemical precursors for industrial materials using only renewable energy. The project was divided into six tasks. In Task 1, DFT calculations were performed to understand the role of Pt and Sn in the catalytic dehydrogenation reaction. In Task 2 chemical kinetics measurements were made for several catalyst formulations at high temperatures. In Task 3, the solar absorptance of catalyst particles was compared to the absorptance of commonly used materials in particle receivers. In Task 4, numerical models were developed which could predict performance of the complete system and predict specific temperatures in the system. In Task 5, a prototype system was designed, fabricated, and tested to show the validity of the concept. Task 6 concerned project management activities. Experiments with the prototype showed repeatable thermal performance at temperatures targeted for the reaction. A limited set of tests were done with active catalyst and propane dehydrogenation, showing conversion of propane to propylene with a range of conversions and selectivities. The results are promising, and the prototype designed was reliable during testing, and the team expects that further development of the prototype would yield improved results. A numerical model framework based on coupled fluid and particle mechanics was developed with high computational efficiency using GPU calculations. The model may prove highly useful for evaluating other high-temperature particle systems. However, it was determined that simpler porous media models were good fits for the needs of the current moving packed bed concept. Data showing strong solar absorption of the particles validates the plan of using existing solar particle receivers with only a change in the particle type. Catalyst investigation showed that Pt 1 Sn 3 is the most viable candidate for developing PtSn catalysts for high temperature propane dehydrogenation, considering the balance of activity, selectivity, and deactivation. This project completed an initial study of various factors needed to incorporate a moving bed catalytic reactor for propane dehydrogenation into a concentrated solar thermal particle system. Future developments may allow this technology to be scaled up and help to use solar thermal energy to decarbonize not only the propane dehydrogenation reaction, but other gas-solid catalytic reactions at similar temperatures.

14 SOLAR ENERGY

TEAMER: Crossflow Turbine Fairing Geometry Optimization - Report and CFD Modeling Files

The dataset includes computational fluid dynamics (CFD) models and simulation files for crossflow turbines as well as a detailed project report. The report documents the project undertaken by the Ocean Renewable Power Company (ORPC) to design and optimize a modular fairing for the Modular RivGen Marine Hydrokinetic (MHK) turbine, which enhances the efficient deployment and operation of turbine arrays. The project focused on optimizing the hydrodynamic performance of the fairing using CFD, with an emphasis on two key geometric parameters: the fairing's cross-sectional shape and the spacing between the rotor and the fairing. The analysis aimed to maximize net power output while also assessing discretized loading to evaluate ultimate and fatigue loads on the turbine components. The numerical modeling was conducted using both the commercial CFD software Star-CCM+ and the open-source code openFOAM, with the latter utilizing the actuator line library, turbinesFOAM. This dual-code approach was intended to increase confidence in the results and demonstrate the viability of using open-source tools for high-fidelity marine energy modeling. This dataset includes all necessary files for actuator line simulations in openFOAM, as well as 2D blade-resolved CFD results, along with Python and Java scripts for setting up and post-processing simulations.

16 TIDAL AND WAVE POWER

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES

Design and Performance Evaluation of a Resistive Control Using a Hydraulic PTO System for the TALOS Wave Energy Converter

This study is focused on developing a numerical model to evaluate the performance of a hydraulic PTO system for the TALOS Wave Energy Converter. The WEC device is described and the architecture of the hydraulic PTO system is presented with detail. The WEC is modeled using WEC-Sim, and the PTO is modeled using the Simscape Fluids library from Simulink. The hydraulic PTO is based on a constant pressure configuration that is suitable for WEC passive control. The hydraulic system is composed by a set of rectifying valves and two hydraulic accumulators that reduce the stiffness of the system and also serve as energy storage devices. One of the advantages of this hydraulic PTO architecture is the possibility of controlling the electric generator to operate around the optimal efficiency operating point. The main components of the hydraulic PTO are off-the-shelf devices that are commercially available, which will facility a future deployment of the designed system. The design variables used for this study are the accumulator size, the maximum pressure in the accumulators, the hydraulic motor maximum displacement, and the shaft speed in the electric generator. The performance of the system is evaluated individually, using sinusoidal inputs that replicates regular wave conditions. In addition to this, the numerical model of the PTO is coupled to a WEC-Sim simulation of the TALOS Wave Energy Converter with six PTOs to generate a wave-to-wire model. The main objective of this work is to present a comprehensive design methodology that could serve as a guideline for future research efforts focused on implementing control algorithms on multi degree of freedom WECs.

hydraulic systems

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences

naturf: a package for generating urban parameters for numerical weather modeling

The Neighborhood Adaptive Tissues for Urban Resilience Futures tool (NATURF) is a Python workflow that generates files readable by the Weather Research and Forecasting (WRF) model. NATURF uses geopandas and hamilton to calculate 132 building parameters from shapefiles with building footprint and height information. These parameters can be collected and used in many formats, and the primary output is a binary file configured for input to WRF. This workflow is a flexible adaptation of the National/World Urban Database and Access Portal Tool (NUDAPT/WUDAPT) that can be used with any study area at any spatial resolution. The climate modeling community and urban planners can identify the effects of building/neighborhood morphology on the microclimate using the urban parameters and WRF-readable files produced by NATURF. More information on the urban parameters calculated can be found in the documentation.

54 ENVIRONMENTAL SCIENCES

Techno-Economic Optimization of Fixed Bed and Moving Bed Contactors for CO2 Capture from NGCC Plants Using a Functionalized Metal-Organic Framework

Amine-appended metal–organic frameworks (MOFs) have strong potential for post-combustion CO2 capture from NGCC plants. In this work, a specific tetramine-appended MOF, Mg2(dobpdc)(3-4-3) (dobpdc4− = 4,4′-dioxidobiphenyl-3,3′dicarboxylate; 3-4-3 = N,N′bis(3-aminopropyl)-1,4-diaminobutane), is evaluated due to its superior performance compared to other candidate diamine-appended and tetraamine-appended MOFs. Techno-economic optimization is performed to minimize the cost of capture by using the Framework for Optimization, Quantification of Uncertainty, and Surrogates (FOQUS), which is capable of optimizing systems by using derivative-free optimization solvers where process models are built using numerous modeling platforms. The technical and economic performances of both fixed bed and moving bed processes are compared with state-of-the-art amine-based solvent systems.

Hughes, Ryan

Air Force Institute of Technology Dust Cloud Model Vertical Dispersion

The Air Force Institute of Technology (AFIT) has developed and refined a nuclear dust cloud model for modeling nuclear fallout, and the dust mass and radiation encountered by aircraft flying through a nuclear dust cloud. Recently, a question was raised requesting an explanation of why the vertical mass and activity distributions had an anomalous shape for yields above 1 Mt. This report explains that the anomalous shape occurs because the dust cloud model is a discrete numeric model, and not a continuous model.

54 ENVIRONMENTAL SCIENCES

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING

Analysis of heat and mass transfer potential of a dew-point cooling tower in different climatic conditions

In this study, the performance of the Dew-Point Cooling Tower (DPCT) was analyzed for different factors in a variety of climate conditions. For this purpose, a dedicated numerical model describing heat and mass transfer processes was developed and validated. The results of the numerical simulations allowed to analyze the potential of utilizing the heat and mass transfer process with the dew-point phenomenon for water cooling. It was established that the operational parameters that have a high impact on the performance of the DPCT are: inlet water temperature and inlet air humidity ratio. It was also established that DPCT achieves the highest COP and Specific Cooling Capacity for cold subtropical highland climates and that it achieves highest Wet-bulb Effectiveness for monsoon-influenced humid subtropical climates. In conclusion, the regions where all the efficiency factors achieved above-average values included warm, arid, and desert climates.

42 ENGINEERING

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY