Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

iharm3D: Vectorized General Relativistic Magnetohydrodynamics

iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM (“High Accuracy Relativistic Magnetohydrodynamics”) algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Prediction of Alpha-Particle-Immune Gate-All-Around Field-Effect Transistors (GAA-FET) Based SRAM Design

Alpha particles are known to be a major source of particles creating soft errors in semiconductor devices, such as content flipping in Static Random-Access Memory (SRAM). Recent advancements in transistor nodes have led to the introduction of Gate-All-Around Field Effect Transistors (GAA-FETs), which have better gate control, thus better electrostatics. Moreover, the introduction of bottom dielectric isolation (BDI) eliminates substrate leakage and thus is expected to enhance its radiation hardness. It is thus important to explore if one can design an SRAM that is completely radiation-hard to alpha particles. In this paper, using 3D Technology Computer-Aided-Design (TCAD) simulations, we show that it is possible to design an SRAM using GAA-FET technology so that it is immune to single alpha particle radiation error. In other words, with the design, there will be no single-event upset (SEU) due to alpha particles. We first use ab initio calculations in PHITS to show that there is a maximum linear energy transfer (LET), LET max , for the alpha particle in Si and Si x Ge 1-x . Based on that, by de signing a sub-7nm GAA-FET-based SRAM with BDI, we show that the SRAM does not flip even if the particle strike is in the worst-case scenario for LET > LET max .

42 ENGINEERING↗
DOE OSTI2021

Fast digital lossy compression for X-ray ptychographic data

Increases in X-ray brightness from synchrotron light sources lead to a requirement for higher frame rates from hybrid pixel array detectors (HPADs), while also favoring charge integration over photon counting. However, transfer of the full uncompressed data will begin to constrain detector design, as well as limit the achievable continuous frame rate. Here a data compression scheme that is easy to implement in a HPAD's application-specific integrated circuit (ASIC) is described, and how different degrees of compression affect image quality in ptychography, a commonly employed coherent imaging method, is examined. Using adaptive encoding quantization, it is shown in simulations that one can digitize signals up to 16383 photons per pixel (corresponding to 14 bits of information) using only 8 or 9 bits for data transfer, with negligible effect on the reconstructed image.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Laser-direct-drive fusion target design with a high- Z gradient-density pusher shell

Laser-direct-drive fusion target designs with solid deuterium-tritium (DT) fuel, a high-Z gradient-density pusher shell (GDPS), and a Au-coated foam layer have been investigated through both 1D and 2D radiation-hydrodynamic simulations. Compared with conventional low-Z ablators and DT-push-on-DT targets, these GDPS targets possess certain advantages of being instability-resistant implosions that can be high adiabat (α ≥ 8) and low hot-spot and pusher-shell convergence (CR hs ≈ 22 and CR PS ≈ 17), and have a low implosion velocity (v imp < 3 × 10 7 cm/s). Using symmetric drive with laser energies of 1.9 to 2.5 MJ, 1D LILAC simulations of these GDPS implosions can result in neutron yields corresponding to >~50–MJ energy, even with reduced laser absorption due to the cross-beam energy transfer (CBET) effect. Two-dimensional DRACO simulations show that these GDPS targets can still ignite and deliver neutron yields from 4 to ~10 MJ even if CBET is present, while traditional DT-push-on-DT targets normally fail due to the CBET-induced reduction of ablation pressure. If CBET is mitigated, these GDPS targets are expected to produce neutron yields of >20 MJ at a driven laser energy of ~2 MJ. Here are the key factors behind the robust ignition and moderate energy gain of such GDPS implosions are as follows: (1) The high initial density of the high-Z pusher shell can be placed at a very high adiabat while the DT fuel is maintained at a relatively low-entropy state; therefore, such implosions can still provide enough compression ρR >1 g/cm 2 for sufficient confinement; (2) the high-Z layer significantly reduces heat-conduction loss from the hot spot since thermal conductivity scales as ~1/Z; and (3) possible radiation trapping may offer an additional advantage for reducing energy loss from such high-Z targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Model for Radiolysis in a Flowing-Water Target during High-Intensity Proton Irradiation

At the Facility for Rare Isotope Beams (FRIB), interactions between heavy-ion beams and beam-dump water will create a wide variety of radionuclides which can be accessed by a technique known as “isotope harvesting”. However, irradiation of water is always accompanied by the creation of numerous radical, ionic, and molecular radiolysis products. Some of the radiolysis products have sufficiently long lifetimes to accumulate in the irradiated water and affect the harvesting chemistry. Here we investigate the formation of hydrogen peroxide, molecular hydrogen, and molecular oxygen during a high-intensity proton irradiation of a flowing-water isotope-harvesting target and compare the experimental results to simulations. The simulations kinetically model the chemical reactions occurring in the homogeneous phase of radiolysis in flowing water and establish an “effective yield”. In both the experiment and simulations, the bulk quantities of H 2 , H 2 O 2 , and O 2 are considerably lower than predicted by primary radiolysis yields (escape yields), meaning that in the high beam intensity regime the homogeneous phase reactions have a considerable impact on the overall chemical composition of the water. Further, it could be shown that for radiation which is characterized by a limited linear energy transfer, such as the here applied protons, the bulk outcome of the microscopic kinetic modeling could be estimated by a simplified steady-state model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Data/code for Mlawer et al. 2024 (A more transparent infrared window)

The infrared window region (780–1,250 cm –1 , 12.8 to 8.0 μm) is of great importance to Earth's climate due to its high transparency and thermal energy. We present here a new investigation of the transparency of this spectral region based on observations by interferometers of downwelling surface radiance at two DOE Atmospheric Radiation Measurement program sites. We focus on the dominant source of absorption in this region, the water vapor continuum, and derive updated values of spectral absorption coefficients for both the self and foreign continua. Our results show that the self continuum is too strong in the previous version of Mlawer-Tobin_Clough-Kneizys-Davies (MT_CKD) water vapor continuum model, a result that is consistent with other recent analyses, while the foreign continuum is too weak in MT_CKD. In general, the weaker self continuum derived in this study results in an overall increase in atmospheric transparency in the window, although in atmospheres with low amounts of water vapor the transparency may slightly decrease due to the increase in foreign continuum absorption. These continuum changes lead to a significant decrease in downwelling longwave flux at the surface for moist atmospheres and a modest increase in outgoing longwave radiation. The increased fraction of surface-leaving radiation that escapes to space leads to a notable increase (~5–10%) in climate feedback, implying that climate simulations that use the new infrared window continuum will show somewhat less warming than before. This study also points out the possibly important role that aerosol absorption may play in the longwave radiative budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impacts of active satellite sensors' low-level cloud detection limitations on cloud radiative forcing in the Arctic

Previous studies revealed that satellites sensors with the best detection capability identify 25 %–40 % and 0 %–25 % fewer clouds below 0.5 and between 0.5–1.0 km, respectively, over the Arctic. Quantifying the impacts of cloud detection limitations on the radiation flux are critical especially over the Arctic Ocean considering the dramatic changes in Arctic sea ice. In this study, the proxies of the space-based radar, CloudSat, and lidar, CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations), cloud masks are derived based on simulated radar reflectivity with QuickBeam and cloud optical thickness using retrieved cloud properties from surface-based radar and lidar during the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment. Limitations in low-level cloud detection by the space-based active sensors, and the impact of these limitations on the radiation fluxes at the surface and the top of the atmosphere (TOA), are estimated with radiative transfer model Streamer. The results show that the combined CloudSat and CALIPSO product generally detects all clouds above 1 km, while detecting 25 % (9 %) fewer in absolute values below 600 m (600 m to 1 km) than surface observations. These detection limitations lead to uncertainties in the monthly mean cloud radiative forcing (CRF), with maximum absolute monthly mean values of 2.5 and 3.4 Wm –2 at the surface and TOA, respectively. Cloud information from only CALIPSO or CloudSat lead to larger cloud detection differences compared to the surface observations and larger CRF uncertainties with absolute monthly means larger than 10.0 Wm –2 at the surface and TOA. The uncertainties for individual cases are larger – up to 30 Wm –2 . These uncertainties need to be considered when radiation flux products from CloudSat and CALIPSO are used in climate and weather studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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

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

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Robust direct laser acceleration of electrons with flying-focus laser pulses

Direct laser acceleration (DLA) offers a compact source of high-charge, energetic electrons for generating secondary radiation or neutrons. While DLA in high-density plasma optimizes the energy transfer from a laser pulse to electrons, it exacerbates nonlinear propagation effects, such as filamentation, that can disrupt the acceleration process. Here, we show that superluminal flying-focus pulses (FFPs) mitigate nonlinear propagation, thereby enhancing the number of high-energy electrons and resulting x-ray yield. Three-dimensional particle-in-cell simulations show that, compared to a Gaussian pulse of equal energy (1 J) and intensity (2 × 10 20 W/cm 2 ), an FFP produces 80 × more electrons above 100 MeV, increases the electron cutoff energy by 20%, triples the high-energy x-ray yield, and improves x-ray collimation. These results illustrate the ability of spatiotemporally structured laser pulses to provide additional control in the highly nonlinear, relativistic regime of laser-plasma interactions.

Laser-produced plasmas↗
DOE OSTI2024

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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↗
DOE OSTI2020

3D Printing Prototypes for a Reflector-Based Broadband Axion Search

This work sought to begin research and development on a horn antenna which can be used to transfer the energy of microwave photons produced by axions in a strong magnetic field to a coaxial cable for measurement. 3D prototypes of pyramidal horn antennas and a parabolic reflector were printed and then coated with conductive spray paint. The S 11 parameter was measured and simulated for both 3D printed horn antenna and the conical horn antenna both with and without the reflector present to determine how much of the radiated power was reflected back. Ultimately, for both the conical and pyramidal antennas, less than one percent of the energy emitted by the antennas was reflected back. When the pyramidal antenna was tilted by about 45 degrees, however, around 30-40 percent of the microwaves emitted by the antennas was reflected back.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Power handling in a highly-radiative negative triangularity pilot plant

Abstract This work explores detailed power handling solutions for a class of high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept (Rutherford et al 2024 Plasma Phys. Control. Fusion ). The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only ∼50 cm (∼40% of the minor radius). The UEDGE code package is used for the boundary plasma solution to track plasma temperatures and fluxes to the divertor targets. It is found that for P SOL = 25 MW and n sep = 0.96 × 10 20 m −3 , conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For a fixed impurity fraction of just 0.13% Ne in the plasma, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW m −2 , while the electron temperature T e remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, P SOL can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1% neon (Ne) impurity content, the divertor still experiences less than 10 MW m −2 at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550 ∘ C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01% Ne content.

Miller, M. A. (ORCID:0000000265406533)↗
DOE OSTI2024

Modeling In-Reactor Temperature Gradients in Lithium Aluminate Ceramic Pellets

In reactor thermal performance of lithium aluminate ceramic pellets are calculated to determine operating temperatures if internal components shift during irradiation. A 2D conduction model was setup to model this scenario with additional physics, including thermal expansion of components to reactor temperatures, creep down of the cladding due to differential pressure between the inside of the rod and reactor coolant system, and net radiation exchange between components. The non-linear coupling of the governing equations presented convergence challenges due to the very different effects from a small temperature change on the radiation vs radial and circumferential conduction equations. Initial temperature guesses for the multidimensional solution were provided from conservative 1D models simulating bounding scenarios. Convergence was greatly improved by taking small steps towards the root solutions and averaging with initial guesses. The final method robustly converged all 132 points at an axial plane after a few dozen iterations. The temperature results show that the additional physics have a secondary effect as the dominant heat transfer mechanism is radial conduction across the gas gap followed by circumferential conduction around components. Circumferential conduction in the lithium aluminate pellets is enhanced by the relatively high thermal conductivity of this ceramic.

Carstens, Nathan A.↗
DOE OSTI2021

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗