Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependence analysis”

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

Empirical probability and machine learning analysis of m , n = 2, 1 tearing mode onset parameter dependence in DIII-D H-mode scenarios

m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRS1979

A Simplified Shuttle Payload Thermal Analyzer /SSPTA/ program

A simple thermal analysis program for Space Shuttle payloads has been developed to accommodate the user who requires an easily understood but dependable analytical tool. The thermal analysis program includes several thermal subprograms traditionally employed in spacecraft thermal studies, a data management system for data generated by the subprograms, and a master program to coordinate the data files and thermal subprograms. The language and logic used to run the thermal analysis program are designed for the small user. In addition, analytical and storage techniques which conserve computer time and minimize core requirements are incorporated into the program.

Bartoszek, J. T.↗
DOE OSTI2021

Orientational Ordering within Semiconducting Polymer Fibrils

Abstract Due to a general paucity of suitable characterization methods, the internal orientational ordering of polymer fibrils has rarely been measured despite its importance particularly for semi‐conducting polymers. An emerging tool with sensitivity to bond orientation is polarized resonant soft X‐ray scattering (P‐RSoXS). Here, P‐RSoXS reveals the molecular arrangement within fibrils (if type I or type II fibrils), the extent of orientation in the fibril crystal, and an explicit crystal‐amorphous interphase. Neat films as well as binary blends with a fullerene derivative are characterized for three different polymers, that are prototypical materials widely used in organic electronics applications. Anisotropic P‐RSoXS patterns reveal two different fibril types. Analysis of the q ‐dependence of the anisotropy from simulated and experimental scattering patterns reveal that neat polymer fibrillar systems likely comprise more than two phases, with the third phase in addition to crystal and amorphous likely being an interphase with distinct density and orientation. Intriguingly, the fibril type correlates to the H‐ or J‐aggregation signature in ultraviolet‐visible (UV–vis) spectroscopy, revealing insight into the fibril formation. Together, the results will open the door to develop more sophisticated structure‐function relationships between chemical design, fibril type, formation pathways and kinetics, interfacial ordering, and eventually device functions.

Mukherjee, Subhrangsu↗
DOE OSTI2025

Contracted, Not Converted: Photoluminescence Thermochromism in Ultrabright Pure Blue Emissive Hybrid Copper(I) Halide

Thermochromic luminescent materials are promising for applications such as sensors or display technologies. Understanding the mechanism of thermochromic luminescence in these materials plays an important role for their targeted applications. Here, a 0D copper(I) halide (TMP)2CuBr3 and 1D (TMP)Ag2Br3 (TMP = tetramethylphosphonium) are reported. While (TMP)Ag2Br3 is a weak light emitter, (TMP)2CuBr3 demonstrates ultrabright pure blue emission with photoluminescence (PL) quantum yield exceeding 90% within the temperature range 300–80 K. The emission in (TMP)2CuBr3 is attributed to the radiative recombination of self-trapped excitons (STEs) that arise localized on [CuBr3]2− anions. (TMP)2CuBr3 demonstrates excitation-selective PL thermochromism at low temperatures. In the literature, the origin of multiple STEs is often unclear in luminescent metal halides. Here, a detailed temperature-dependent spectroscopic and structural analysis suggests that the PL thermochromism in (TMP)2CuBr3 is associated with the appearance of a second STE state. The appearance of second STE is expected to have a structural origin, and it may be associated with the observed anomaly in thermal contraction near 200 K. The distinct PL features of melt-processible (TMP)2CuBr3 coupled with its remarkable photo- and thermal stability allow its consideration for practical optical applications, including temperature sensing.

Popy, Dilruba [University of Oklahoma, Norman]↗
DOE OSTI2022

Generalized parton distributions through universal moment parameterization: zero skewness case

We present a global analysis program for the generalized parton distributions (GPDs) based on conformal moment expansion. We apply the strategy of universal moment parameterization to fit both the collinear parton distribution functions (PDFs) from phenomenology and generalized form factors from lattice calculations, and show that the parameterization is flexible enough to accommodate these constraints. In addition, we can also fit direct lattice calculations of GPDs from large-momentum effective theory. In this work we focus on the analysis of t-dependent PDFs which correspond to GPDs in the ξ → 0 limit. The strategy also applies to the ξ ≠ 0 region with extra parameters, and therefore can be fitted to experimental observables in the future. With a demonstrative example of fitted GPDs, we exhibit the quark transverse angular momentum densities of the proton as well as the impact parameter space distributions of quarks in both unpolarized and transversely polarized protons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Impact of heating and cooling loads on battery energy storage system sizing in extreme cold climates

Efficient operation of battery energy storage systems requires that battery temperature remains within a specific range. Current techno-economic models neglect the parasitic loads heating and cooling operations have on these devices, assuming they operate at constant temperature. In this work, these effects are investigated considering the optimal sizing of battery energy storage systems when deployed in cold environments. Here, a peak shaving application is presented as a linear programming problem which is then formulated in the PYOMO optimization programming language. The building energy simulation software EnergyPlus is used to model the heating, ventilation, and air conditioning load of the battery energy storage system enclosure. Case studies are conducted for eight locations in the United States considering a nickel manganese cobalt oxide lithium ion battery type and whether the power conversion system is inside or outside the enclosure. The results show an increase of 42% to 300% in energy capacity size, 43% to 217% in power rating, and 43% to 296% increase in capital cost dependent on location. This analysis shows that the heating, ventilation, and air conditioning load can have a large impact on the optimal sizes and cost of a battery energy storage system and merit consideration in techno-economic studies.

25 ENERGY STORAGE↗
DOE OSTI2026

Optimizing Porous Transport Layer Porosity for Proton Exchange Membrane Water Electrolysis

An empirical model is presented that describes anode-side losses related to porous transport layer (PTL) morphology in proton exchange membrane water electrolysis (PEMWE). The model is based on an advanced voltage breakdown analysis that links various overpotentials to PTL morphology. Custom Ti PTLs, spanning uncommonly low porosities (22 - 31%), were fabricated and analyzed with X-ray CT to obtain pore and particle size distributions. Particle size distributions were consistent across samples with an average particle diameter of 12.0?..mu..m, whereas average pore diameters ranged from 6.0 to 7.0?..mu..m. The PTLs were tested in standard PEMWE cell assemblies with anode catalyst loadings of 0.1 mgIr cm-2 to obtain polarization curves, electrochemical impedance spectra, and augmented Tafel analysis. The PTL-dependent anode side losses were deconvoluted and assigned to excess utilization, concentration, ion transport resistance, and electrical contact resistance overpotentials. The data and model reveal an optimal 20 - 28% PTL porosity region where utilization and contact resistance overpotentials are minimized without triggering concentration and ion transport losses related to water deprivation. The optimal PTL porosity depends on the operating current density and is demonstrated at realistic PEMWE water flow rates to establish PTL design guidance for operation at scale.

08 HYDROGEN↗
DOE OSTI2022

Impact of acetic acid exposure on metal contact degradation of different crystalline silicon solar cell technologies

Degradation due to acetic acid in photovoltaic (PV) modules has been a commonly observed phenomenon for both damp-heat exposure and outdoor operations. Acetic acid is a degradation byproduct of ethylene-vinyl acetate (EVA), a common module encapsulant. To address this issue, robust metallization pastes and cell technologies are being developed. However, it is important to assess how these technologies perform in an acetic acid environment and withstand degradation before they are implemented in the solar market. In this work, we investigate the impact of acetic acid exposure on four different cell groups: monofacial passivated emitter and rear contact (PERC) cells with advanced telluride-based front contact pastes, bifacial PERC cells with novel aluminum rear contact pastes, bifacial tunnel oxide passivated contacts (TOPCon) cells, and silicon heterojunction (SHJ) cells. These cells were exposed to acetic acid for different time increments. The recombination losses were characterized by Suns-VOC, and multi-variate regression analysis of intensity-dependent photoluminescence (PL) images with Griddler AI. Resistive losses were tracked with the transmission line method (TLM). Samples showing severe performance degradation were selected for further materials characterization to understand the root cause. Top-down and cross-sectional scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), and X-ray photoelectron spectroscopy (XPS) were performed to investigate the change in materials properties. Our study shows that the front contacts of the bifacial TOPCon cells and monofacial PERC cells were significantly affected by acetic acid exposure. Here, the SHJ cells were found to be the most stable.

14 SOLAR ENERGY↗
DOE OSTI2024

Using Data-Science Approaches to Unravel Insights for Enhanced Transport of Lithium Ions in Single-Ion Conducting Polymer Electrolytes

Solid polymer electrolytes have yet to achieve the desired ionic conductivity (>1 mS/cm) near room temperature required for many applications. This target implies the need to reduce the effective energy barriers for ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer–Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules 2023, 56, 15, 6051), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify anions leading to the energy barriers <30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers <20 kJ/mol, we have performed density functional theory (DFT) based calculations to connect the chemical structures of boron-based anions via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion hopping. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity >1 mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes.

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

Structurally Frustrated Remeika Phases: Tolerance Factor and Non-Kramers’ Ion Effects

Quantum materials that are found on the verge of structural, magnetic, and electronic instabilities are deep reservoirs for exotic phenomena. Our criteria to study the reported Ln 5 Ru 6 Sn 18 (Ln = Gd, Tb) are guided by the tolerance factor for Remeika phases─a class of intermetallic compounds with the general formula A 3 M 4 X 13 (where A = rare-earth element, M = transition metal, and X = tetrel), which can be viewed as structural analogues of pseudoperovskites (4 × ABX 3 ), with one of the tetrel atoms occupying the A-site. We have grown single crystals, up to 0.5 cm on the largest facet, of Gd 5 Ru 6 Sn 18 (a = 13.8052(12) Å) and Tb 5 Ru 6 Sn 18 (a = 13.7825(9) Å), which we expect to be on the verge of a formation instability. Bulk magnetic, heat capacity, and electrical transport measurements reveal unusual behavior that originates from the trivalent Gd and Tb ions, with deviations from typical metallic behavior and magnetic ordering that is short-range or disrupted in some other way. Furthermore, to better understand these phenomena, single-crystal neutron diffraction is investigated, which reveals short-range correlations among the partially occupied Tb and Sn sublattices, with 3D-ΔPDF analysis supporting direction-dependent local order. Taken together, these measurements show the utility of using tolerance factors in identifying materials that are likely to exhibit complex phenomena and uncover environments where a formation instability may strongly impact emergent bulk magnetic and electronic phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2020

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

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

Proton Transfer Dynamics in the Aprotic Proton Accepting Solvent 1-Methylimidazole

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) were investigated using fast fluorescence measurements. The closely related molecule, 8- methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as: protonated (P), associated I (A I ), associated II (A II ), and deprotonated (D). We find that the simple kinetic scheme of P → A I → A II → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into A I , which mainly goes on to deprotonation (D), but AI can also feed into A II . A II can return to A I or decay to the ground state, but does not deprotonate within experimental error. Finally, quantum chemistry and excitedstate QM/MM Born-Oppenheimer molecular dynamics simulations indicate that A I and A II are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial and equatorial, respectively.

14 SOLAR ENERGY↗
DOE OSTI2023

Plasmon-Induced Hot-Carrier Excited-State Dynamics in Plasmonic Semiconductor Nanocrystals

The variance of carrier relaxation pathways for WO 3–x plasmonic semiconductor nanocrystals (PSNCs) is monitored by transient absorption spectroscopy following excitation of the localized surface plasmon resonance (LSPR) versus the optical band gap (E g,opt ). Excitation of the LSPR leads to efficient hot carrier population above the Fermi level in WO 3–x via Landau damping, in analogy to noble metal LSPR relaxation mechanisms. Hot carrier depopulation occurs on the femtosecond timescale, observed as the concomitant recovery of an LSPR bleach with the appearance of discrete interband and intraband photoinduced absorption features. Further, by comparison, the direct excitation of E g,opt results in trion recombination at donor–acceptor sites within the WO 3–x NC, consistent with exciton decay dynamics observed for typical wide-band-gap semiconductor NCs. From the analysis of pump power dependency data, a hot-carrier electron–phonon coupling constant of 1.47 × 10 11 J K –1 s –1 cm –3 is extracted. The direct comparison of the decay dynamics following E g,opt versus LSPR excitation confirms that the observed plasmon in trioctylphosphine oxide passivated, spherical WO 3–x is a resonance state in which hot carriers are generated only from excitation on resonance with the LSPR frequency. This study on WO 3–x PSNCs provides a toolset that can be used to evaluate the role of hot carriers following LSPR excitation of n-type, plasmonic transparent conducting oxide NCs, where enhancement of photocatalysis, photovoltaic performance, and optical enhancement has been reported.

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

Investigating Cu-Site Doped Cu–Sb–S Nanoparticles Using Photoelectron and Electron Paramagnetic Resonance Spectroscopy

Tetrahedrite (Cu 12 Sb 4 S 13 ) and famatinite (Cu 3 SbS 4 ) are good candidates for green energy applications because they possess promising thermoelectric and photovoltaic properties as well as contain earth-abundant and nontoxic constituents. Herein, X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and electron paramagnetic resonance spectroscopy (EPR) methods examined inherent electronic properties and interatomic magnetic interactions of Cu-site doped tetrahedrite and famatinite nanomaterials. An energy-efficient modified polyol method was utilized for the synthesis of tetrahedrite and famatinite nanoparticles doped on the Cu-site with Zn, Fe, Ni, Mn, and Co. This is the first parallel study of tetrahedrite and famatinite nanomaterials with XPS, UPS, and EPR methods alongside a systematic analysis of dopant-dependent effects on the electronic structure and magnetic interactions for each material. XPS showed that the Cu and Sb species in tetrahedrite and famatinite possess different oxidation states, while UPS characterization reveals larger dopant-dependent shifts in the work function for tetrahedrite nanoparticles (4.21 to 4.79 eV) than for famatinite nanoparticles (4.57 to 4.77 eV). Finally, all famatinite nanoparticles display an EPR signal, indicating trace amounts of paramagnetic Cu(II) present below the detection limit of XPS. For tetrahedrite, EPR signatures were observed only for the Zn-doped and Mn-doped nanoparticles, suggesting signal broadening from Cu–Cu spin exchange or spin–lattice relaxation. This study demonstrates the complementary nature of XPS and EPR techniques for studying the oxidation states of metals in solid-state nanomaterials. Comparing the electronic and magnetic properties of tetrahedrite and famatinite while studying the impact of dopant incorporation will guide future endeavors in designing sustainable, high-performance materials for renewable energy applications.

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

Vibrational Dynamics and Phase Transitions of Hydrazine to 50 GPa

The high pressure behavior of hydrazine, N 2 H 4 , has been investigated to 50 GPa at room temperature using infrared and Raman spectroscopy to explore pressure induced phase transitions and changes in hydrogen bonding. Three solid–solid phase transitions were detected at 11, 21, and 32 GPa on room temperature compression through dramatic changes in the lattice vibration and N–H stretching regions with increasing pressure in both measurement techniques. The transition to phase IV, which appears at 32 GPa, exhibits increased hydrogen bonding with significant hysteresis, persisting to 9 GPa on decompression. This work presents a detailed analysis of the pressure dependence of mode shifts and calculations of mode Grüneisen parameters as well as a determination of an approximate thermodynamic Grüneisen parameter. We compare these results to the behavior of other small molecular materials such as ammonia and water and explore the evolution of hydrogen bonding in hydrazine toward the symmetrically hydrogen bonded state, which has previously been suggested by theoretical computations.

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

Importance of Spatially Continuous Urban Surface Properties in Urban‐Resolving Earth System Modeling

Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.

Cheng, Yifan [University of Illinois Urbana-Champa↗