Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface estimation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Field‐representative evaluation of PID‐polarization in TOPCon PV modules by accelerated stress testing

Abstract Potential‐induced degradation‐polarization (PID‐p) can reduce module power, but how to project the extent to which PID‐p may occur in field conditions considering the factors of system voltage, condensed moisture, temperature, and illumination has not been clarified. Using tunnel oxide passivated contact (TOPCon) modules, this work demonstrates a method to test full‐size crystalline silicon PV modules for PID‐p to provide field‐representative results. In initial screening tests with positive or negative 1000 V electrical bias applied at 60°C for 96 h using Al foil electrodes on the glass surfaces, the module type exhibited reversible PID‐p only on the front face when the cell circuit was in negative voltage potential. No PID was detected on the rear after testing in either polarity. We then evaluated the PID‐p sensitivity on the front side under different UV irradiances while maintaining the glass surface wet to estimate real‐world susceptibility to PID‐p. The magnitude of the observed behavior was fit using a previously developed charge transfer and depletion by light model. Whereas power loss with −1000 V applied to the cell circuit at 60°C for 96 h in the dark was about 30%, testing the module front under 0.051 W·m −2 nm −1 at 340 nm UVA irradiation using fluorescent tubes, the mean degradation was only 3%. When the modules were tested in the dark for PID‐p with in situ dark current–voltage (I‐V) characterization, the thermal activation energy for degradation was 0.71 eV; for recovery in the dark, it was 0.58 eV. Whereas recovery from the degraded state at 60°C in the dark without voltage bias was 5% absolute in 38 h, rapid recovery of about 5% absolute was observed with 1000 W·s/m 2 exposure at 25°C using a flash tester.

14 SOLAR ENERGY↗
DOE OSTI2023

Evolution of titanium particle combustion in potassium perchlorate and air

Understanding titanium particle combustion processes is critical not only for characterizing existing pyrotechnic systems but also for creating new igniter designs. In order to characterize titanium particle combustion processes, morphologies, and temperatures, simultaneous spatially-resolved electric field holography and imaging pyrometry techniques were used to capture post-ignition data at up to 7 kHz. Due to the phase and thermal distortions present in the combustion cloud, traditional digital in-line holography techniques fail to capture accurate data. Here in this work, electric field holography techniques are used in order to cancel distortions and capture the three-dimensional spatial locations and diameters of the particles. In order to estimate the projected surface temperatures of the titanium particles, an imaging pyrometry method that ratios emission at 750 and 850 nm is utilized. Using these diagnostics, joint statistics are collected for particle size, morphology, velocity, and temperature. Results show that, early in the combustion process, the titanium particles are primarily oxidized by potassium perchlorate inside the igniter cup, resulting in projected surface temperatures near 3000 K. Later in the process, the particles interact with ambient air, resulting in lower surface temperatures around 2400 K and the formation of flame zones. These results are consistent with adiabatic flame temperature predictions as well as particle morphology observations of a titanium core with a TiO 2 surface. Late stage particle expansion, star fragmentation, and molten droplet breakup events are also observed using the time-resolved morphology and temperature diagnostics. These results illustrate the different stages of titanium particle combustion in pyrotechnic environments, which can be used to inform improvements in next-generation igniters.

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

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A critical look at the prediction of the temperature field around a laser-induced melt pool on metallic substrates

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

42 ENGINEERING↗
DOE OSTI2020

H i imaging of dwarf star-forming galaxies: masses, morphologies, and gas deficiencies

ABSTRACT The Giant Meter-wave Radio Telescope observations of the H i 21 cm-line emission from 13 nearby dwarf star-forming galaxies are presented. These galaxies are selected from the catalogues of Wolf−Rayet galaxies having very young (≤10 Myr) star formation. The ranges of star formation rates and stellar masses of the sample galaxies are 0.03–1.7 M⊙ yr−1 and 0.04–22.3 × 108 M⊙, respectively. The H i line emission is detected from 12 galaxies with peak column density >1 × 1021 cm−2. The 3σ H i column density sensitivities per channel width of 7 km s−1 for low (60 arcsec × 60 arcsec) resolution images are in the range 0.8–1.9 × 1019 cm−2. The H i channel images, moment images, global profiles, and mass surface density profiles are presented here. The average value of the peak H i mass surface density is estimated to be ∼2.5 M⊙ pc−2, which is significantly less compared to that in massive spiral galaxies. The scaling relations of $(M_{stars} + M_{\rm H\, I} + M_{\rm He})$versus Mdyn, gas fraction versus MB, $M_{\rm H\, I}$versus Mstars, H i-to-stellar mass ratio versus Mstars, and $M_{\rm H\, I}$versus $D_{\rm H\, I}$for the sample galaxies are estimated. These scaling relations can be used to constraint the key parameters in the galaxy evolution models. These galaxies are residing in group environment with galaxy density up to eight galaxy Mpc−3. An H i mass deficiency (with DEFH i > 0.3) is noticed in majority of galaxies for their optical diameters as compared to galaxies in field environments. Clear signatures of tidal interactions in these galaxies could be inferred using the H i images. Isolated H i clouds without known optical counterparts are seen in the vicinity of several galaxies. H i emission envelope is found to be having an offset from the optical envelope in several galaxies. Consistent with the previous studies on galaxy evolution in group environments, tidal interactions seem to play an important role in triggering recent star formation.

Jaiswal, S.↗
DOE OSTI2020

Temperature Measurement of Laser-Irradiated Metals Using Hyperspectral Imaging

Accurate noncontact surface-temperature measurements during laser-based materials processing remain challenging due to the difficulty of establishing reliable emissivity values as a function of temperature and wavelength. Direct measurement of emissivity is difficult, as the emissivity may be changing constantly in the laser-material interaction region, where the temperature gradients are extreme and surface displacement can complicate the measurement. Here, we present a hyperspectral imaging method using a multiwavelength camera to capture the spectral radiance in eight different bands and fit the wavelength-dependent radiance to Planck’s law. Time-resolved temperature measurements during microsecond pulsed-laser irradiation of a metal plate made of the titanium alloy Ti -6 Al -4 V provide temperature information about the irradiated surface with an estimated accuracy of ± 10 % . The extracted wavelength-dependent emissivity slope compares well with reported results for macroscale titanium melts obtained using IR spectroscopy, which measure temperature without any built-in assumptions about the emissivity. The results are directly applicable to temperature monitoring during welding, additive manufacturing, and semiconductor growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: ​​Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗
DOE OSTI2022

On-Demand Designing of Cathode Internal Surface Architecture for Dramatic Enhancement of SOFC Performance and Durability

This project is aimed to design and modify the internal surfaces of porous composite cathode from currently commercially viable Solid Oxide Fuel Cells (SOFCs), using additive manufacturing process of Atomic Layer Deposition (ALD). The material systems being investigated are commercial composite electrodes complex three-dimensional topographies. In term of the chemistry of the ALD layer applied on the internal surface of the porous cathode, this project has employed commercially relevant electrolyte, electrocatalyst and noble metal materials set. Such materials are fully compatible with the commercial fuel cells, and this project has developed special nanostructure on the surface of the commercial composite cathodes. The formation of the designed nanoarchitecture on the surface of SOFC cathode has been achieved through precise control of ALD parameters and their effect on overall cell performance and resultant electrochemical reaction mechanism of cathodes has been investigated through full cell electrochemical performance testing and nanostructure characterization by transmission electron microscopy (TEM). Under the support of this award, following has been achieved: (1). For cathode materials in solid oxide fuel cells (SOFCs), such as perovskite mixed conductor La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-x (LSCF), cation surface segregation and consequently losing conductivity and active sites for the oxygen reduction reaction (ORR) are problematic. To mitigate the cation segregation and enhance SOFC durability, further decorating the internal backbone surface using the desired electrocatalysts could be a promoting approach. Commonly, the cation segregation such as Sr is very volatile, so the effective surface decoration is ideally conformal. Nevertheless, the conformal surface coating would inevitably alter the ORR pathways that initially take place on the surface of the backbone. To reveal the impact of the conformal coating on both the catalytic activity and the conductivity of the cathode, the unary electrocatalyst of Pt or CoO x , was applied to the LSCF/SDC composite electrode of inherently functional SOFCs, respectively. Both ALD coating layers evolve strong interaction with the LSCF composite cathode. Upon operations, the Pt coating layer remains conformal on LSCF grain surfaces but turns into discrete particles on SDC grain surfaces. Meanwhile, CoO x conformal coating grows to be the discrete nanograins on both the LSCF and SDC grains. ALD coating of the cathode alone reduces the ohmic resistance up to 28 % for the entire cells. The increased conductivity induced by the ALD coating of Pt or CoO x is ascribed to different mechanisms. For the inherent functional SOFCs, the present study presents a novel and feasible approach to apply a conformal, dense coating layer on the surface of a mixed conductor, simultaneously increasing the conductivity and durability of the SOFC cathode. (2). High resistance of the oxygen electrode still significantly hinders the state-of-the-art Solid Oxide Fuel Cells (SOFCs). In particular, for an oxygen electrode consisting of mixed electronic and ionic conductors, such as perovskite lanthanum strontium cobalt ferrite (LSCF), it deteriorates due to its low chemical stability of the grain surface. Such degradation is often associated with the segregation of cations. To prevent the cation surface segregation and its resultant perovskite phase decomposition, we demonstrate a conformal ultra-thin (7-10 nm) surface heterogeneous coating layer consisting of subjacent discrete Pt nanoparticles capped with a superjacent fully dense conformal CoO x layer. The performance studies indicate the ALD coating reduces the cell series resistance by up to 40 %. The conformal CoO x layer consists of randomly orientated but single-layered nanograins, with high-density intergranular and surface grain boundaries serving as the electrochemical reaction sites and facilitating mass transport. The conformal coating layer appears to have successfully suppressed the Sr outward diffusion and confined the Sr enriched layer to a ~ 2 nm interface perovskite phase between the coating layer and the LSCF grain surface. Moreover, this ultra-thin Sr enriched perovskite layer presumably possesses high oxygen vacancy and high ionic conductivity and further imposes tensile strain to the LSCF grain surfaces. With the combination of a conformal CoO x nanoionics, Sr enriched layer, and its strained interface, the ALD coating induced surface layer is estimated to have a conductivity of ~ 1.27x10 4 S/cm, which is over two orders magnitude of that from LSCF at 750 ºC.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Methodological guidelines on concrete degradation based on predictive models and the release of MOSAIC for industry use

In light-water reactors (LWRs), the concrete biological shield (CBS) designates the concrete structure directly facing the reactor pressure vessel (RPV). The primary function of the CBS is to protect equipment and personnel from the neutron and gamma radiation exiting the RPV. Most CBSs in operation in the United States also provide structural support to the RPV. Based on a long-term structural analysis, the structural integrity of the CBS to transfer the in-service passive load to the foundation, dynamic loading during an earthquake, and thermal loading during a loss-of-coolant accident (LOCA) must be maintained. Based on test reactor data for accelerated conditions, the mechanical properties of irradiated concrete are affected when exposed to neutron fluence higher than ≈ 10 19 n.cm -2 (E > 0.1 MeV). At 80 years of operation, the fluence estimates at the surface of the CBS range between 1×10 19 and 7×10 19 n.cm -2 (E > 0.1 MeV): the fluence value depends mostly on the design. Although fluence attenuates through the CBS by virtue of the shielding properties of concrete, about 5% to 10% of the CBS wall depth is affected by high fluence. The main degradation mechanism is the so-called radiation-induced volumetric expansion (RIVE) caused by the amorphization of aggregate-forming minerals. RIVE causes large deformation of up to 18% in pure quartz, and it also reduces the structural properties of concrete. The structural effects of long-term exposure to fast neutron irradiation in LWRs was studied in the report entitled “Assessment of the Effect of the Irradiation-Induced Degradation on the In-Service Structural Performance of the Concrete Biological Shields” (M3LW-23OR0403013). This report illustrates that (1) irradiation-induced damage expands further than the region subject to a cumulative fluence of 10 19 n.cm -2 (E > 0.1 MeV) to reach the steel reinforcement location, and (2) the irradiation-induced damage region shows no substantial residual bearing capacity and should not be accounted for in the integrity assessment of the CBS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

CERESMIP: a climate modeling protocol to investigate recent trends in the Earth's Energy Imbalance

The Clouds and the Earth's Radiant Energy System (CERES) project has now produced over two decades of observed data on the Earth's Energy Imbalance (EEI) and has revealed substantive trends in both the reflected shortwave and outgoing longwave top-of-atmosphere radiation components. Available climate model simulations suggest that these trends are incompatible with purely internal variability, but that the full magnitude and breakdown of the trends are outside of the model ranges. Unfortunately, the Coupled Model Intercomparison Project (Phase 6) (CMIP6) protocol only uses observed forcings to 2014 (and Shared Socioeconomic Pathways (SSP) projections thereafter), and furthermore, many of the ‘observed' drivers have been updated substantially since the CMIP6 inputs were defined. Most notably, the sea surface temperature (SST) estimates have been revised and now show up to 50% greater trends since 1979, particularly in the southern hemisphere. Additionally, estimates of short-lived aerosol and gas-phase emissions have been substantially updated. These revisions will likely have material impacts on the model-simulated EEI. We therefore propose a new, relatively low-cost, model intercomparison, CERESMIP, that would target the CERES period (2000-present), with updated forcings to at least the end of 2021. The focus will be on atmosphere-only simulations, using updated SST, forcings and emissions from 1990 to 2021. The key metrics of interest will be the EEI and atmospheric feedbacks, and so the analysis will benefit from output from satellite cloud observation simulators. The Tier 1 request would consist only of an ensemble of AMIP-style simulations, while the Tier 2 request would encompass uncertainties in the applied forcing, atmospheric composition, single and all-but-one forcing responses. We present some preliminary results and invite participation from a wide group of models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Marangoni Convection using a Velocity-Time Relationship

The Waste Treatment and Immobilization Plant (WTP) at the U.S. Department of Energy’s (DOE’s) Handford site is designed to treat 56 million gallons of radioactive waste through vitrification. The melting vessels physical and chemical integrity is essential to maintain, as failure has a high safety and financial penalty. Monofrax® K-3 refractory, the material widely used for the melting vessel, undergoes significant corrosion over time, concentrated at the triple-point junction between the gaseous atmosphere, refractory, and glass melt, producing a characteristic neck-like profile that significantly limits refractory service life. Accurate prediction of the corrosion within the melter can be used to prevent this and extend the lifetime of the melters. A computational fluid dynamics (CFD) model was previously developed, and predicted subsurface refractory corrosion, but failed to predict the neck-line profile observed experimentally. To address this limitation, the glass-melt meniscus geometry at the glass-refractory interface was estimated from experimentally measured neck profiles and wettability experiments. A curve was fitted to the neck profile to identify the triple point location, and the meniscus shape was characterized by a contact angle of 1°. A velocity profile accounting for Marangoni convection-driven corrosion was derived from the estimated meniscus geometry. Surface tension measurements from LORPM14R1 glasses of varying composition were coupled with K-3 coupon dissolution rates to establish a surface tension–time relationship. This analytical velocity-time model was integrated into the existing CFD framework to improve prediction of the neck corrosion profile. Work was done in a computational setting with national collaboration, resulting in professional development through presentational and technical writing growth.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Spectral solar irradiance on inclined surfaces: A fast Monte Carlo approach

Estimating spectral plane-of-array (POA) solar irradiance on inclined surfaces is an important step in the design and performance evaluation of both photovoltaic and concentrated solar plants. This work introduces a fast, line-by-line spectral, Monte Carlo (MC) radiative transfer model approach to simulate anisotropic distributions of shortwave radiation through the atmosphere as photon bundles impinge on inclined surfaces. This fast Monte Carlo approach reproduces the angular distribution of solar irradiance without the undesirable effects of spatial discretization and thus computes detailed POA irradiance values on surfaces at any orientation and also when surfaces are subjected to the anisotropic ground and atmospheric scattering. Polarization effects are also easily incorporated into this approach that can be considered as direct numerical simulation of the physics involved. Here, we compare our Monte Carlo radiative transfer model with the most widely used empirical transposition model, Perez4, under various conditions. The results show that the Perez4 model reproduces the more detailed Monte Carlo simulations with less than 10% deviation under clear skies for all relevant surface tilt and azimuth angles. When optically thin clouds are present, observed deviations are larger, especially when the receiving surface is strongly tilted. Deviations are also observed for large azimuth angle differences between the receiving surface and the solar position. When optically thick clouds are present, the two models agree within 15% deviation for nearly all surface orientation and tilt angles. The overall deviations are smaller when compared with cases for optically thin clouds. The Perez4 model performs very well (6.0% deviation) in comparison with the detailed MC simulations for all cases, thus validating its widespread use for practical solar applications. When detailed atmospheric profiles and cloud optical properties are available, the proposed fast Monte Carlo radiative model reproduces accurate spectral and angular POA irradiance levels for various atmospheric and cloud cover conditions, surface orientations, and different surface and ground properties.

42 ENGINEERING↗
DOE OSTI2020

Conserving Land–Atmosphere Synthesis Suite (CLASS)

Accurate estimates of terrestrial water and energy cycle components are needed to better understand climate processes and improve models’ ability to simulate future change. Various observational estimates are available for the individual budget terms; however, these typically show inconsistencies when combined in a budget. In this work, a Conserving Land–Atmosphere Synthesis Suite (CLASS) of estimates of simultaneously balanced surface water and energy budget components is developed. Individual CLASS variable datasets, where possible, 1) combine a range of existing variable product estimates, and hence overcome the limitations of estimates from a single source; 2) are observationally constrained with in situ measurements; 3) have uncertainty estimates that are consistent with their agreement with in situ observations; and 4) are consistent with each other by being able to solve the water and energy budgets simultaneously. First, available datasets of a budget variable are merged by implementing a weighting method that accounts both for the ability of datasets to match in situ measurements and the error covariance between datasets. Then, the budget terms are adjusted by applying an objective variational data assimilation technique (DAT) that enforces the simultaneous closure of the surface water and energy budgets linked through the equivalence of evapotranspiration and latent heat. Comparing component estimates before and after applying the DAT against in situ measurements of energy fluxes and streamflow showed that modified estimates agree better with in situ observations across various metrics, but also revealed some inconsistencies between water budget terms in June over the higher latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Identification of surface urban heat versus cool islands for arid cities depends on the choice of urban and rural definitions

The urban heat island (UHI) effect in arid cities can be small or even negative, the latter known as the urban cool island (UCI) effect. Differences in defining urban and rural areas can introduce uncertainties in detecting UHI or UCI, especially when the UHI signal is small. Here, we compared the surface UHI intensity (SUHII) estimated by a dozen different methods (with multiple urban and/or rural definitions) across 104 arid cities globally, providing a comprehensive evaluation of the uncertainty in SUHII estimates. Results show that the absolute difference in annual average SUHII (ΔSUHII) among methods exceeded 1°C in about half of the arid cities during both daytime and nighttime. Further, the overall annual mean ΔSUHII for all arid cities was 1.35°C during daytime and 1.03 °C at night. The uncertainty arising from simultaneous variations in urban and rural definitions was generally higher than that resulting from their individual changes. It was observed that, with varying definitions of urban and rural areas, nearly 50% of arid cities experienced a sign reversal in daytime SUHII estimates, while approximately 15% exhibited a sign reversal in nighttime SUHII. Variations in urban-rural differences in surface properties, such as vegetation index and albedo, due to differing urban and rural definitions, contributed strongly to the observed SUHII uncertainties. Overall, our results offer new insights into the ongoing debate on heat and cold islands in arid cities, emphasizing a critical need to standardize SUHII estimation frameworks.

54 ENVIRONMENTAL SCIENCES↗