Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Electric vehicle fast charging infrastructure planning in urban networks considering daily travel and charging behavior

Electric vehicles are a sustainable substitution to conventional vehicles. This work introduces an integrated framework for urban fast charging infrastructure to address the range anxiety issue. A mesoscopic simulation tool is developed to generate trip trajectories, and simulate charging behavior based on various trip attributes. The resulting charging demand is the key input to a mixed-integer nonlinear program that seeks charging station configuration. The model minimizes the total system cost including charging station and charger installation costs, and charging, queuing, and detouring delays. The problem is solved using a decomposition technique incorporating a commercial solver for small networks, and a heuristic algorithm for large-scale networks, in addition to the Golden Section method. The solution quality and significant superiority in the computational efficiency of the decomposition approach are confirmed in comparison with the implicit enumeration approach. Furthermore, the required infrastructure to support urban trips is explored for future market shares and technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Dynamic Nanocrystal Superlattices with Thermally Triggerable Lubricating Ligands

The size-dependent and collective physical properties of nanocrystals (NCs) and their self-assembled superlattices (SLs) enable the study of mesoscale phenomena and the design of metamaterials for a broad range of applications. However, the limited mobility of NC building blocks in dried NCSLs often hampers the potential for employing postdeposition methods to produce high-quality NCSLs. In this study, we present tailored promesogenic ligands that exhibit a lubricating property akin to thermotropic liquid crystals. The lubricating ability of ligands is thermally triggerable, allowing the dry solid NC aggregates deposited on the substrates with poor ordering to be transformed into NCSLs with high crystallinity and preferred orientations. The interplay between the dynamic behavior of NCSLs and the molecular structure of the ligands is elucidated through a comprehensive analysis of their lubricating efficacy using both experimental and simulation approaches. Coarse-grained molecular dynamic modeling suggests that a shielding layer from mesogens prevents the interdigitation of ligand tails, facilitating the sliding between outer shells and consequently enhancing the mobility of NC building blocks. Further, the dynamic organization of NCSLs can also be triggered with high spatial resolution by laser illumination. The principles, kinetics, and utility of lubricating ligands could be generalized to unlock stimuli-responsive metamaterials from NCSLs and contribute to the fabrication of NCSLs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Matched Guiding and Controlled Injection in Dark-Current-Free, 10-GeV-Class, Channel-Guided Laser-Plasma Accelerators

Here we measure the high-intensity laser propagation throughout meter-scale, channel-guided laser-plasma accelerators by adjusting the length of the plasma channel on a shot-by-shot basis, showing high-quality guiding of 500 TW laser pulses over 30 cm in a hydrogen plasma of density n 0 ≈ 1 X 10 17 cm -3 . We observed transverse energy transport of higher-order modes in the first ≈ 12 cm of the plasma channel, followed by quasimatched propagation, and the gradual, dark-current-free depletion of laser energy to the wake. We quantify the laser-to-wake transfer efficiency limitations of currently available petawatt-class lasers and demonstrate via simulation how control over the laser mode can significantly improve beam parameters. Using 21.3 J of laser energy, and triggering localized electron injection, we observed electron bunches with single, quasimonoenergetic peaks up to 9.2 GeV with charge extending beyond 10 GeV.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗
DOE OSTI2022

Interpreting machine learning prediction of fire emissions and comparison with FireMIP process-based models

Annual burned areas in the United States have increased 2-fold during the past decades. With more large fires resulting in more emissions of fine particulate matter, an accurate prediction of fire emissions is critical for quantifying the impacts of fires on air quality, human health, and climate. This study aims to construct a machine learning (ML) model with game-theory interpretation to predict monthly fire emissions over the contiguous US (CONUS) and to understand the controlling factors of fire emissions. The optimized ML model is used to diagnose the process-based models in the Fire Modeling Intercomparison Project (FireMIP) to inform future development. Results show promising performance for the ML model, Community Land Model (CLM), and Joint UK Land Environment Simulator-Interactive Fire And Emission Algorithm For Natural Environments (JULES-INFERNO) in reproducing the spatial distributions, seasonality, and interannual variability of fire emissions over the CONUS. Regional analysis shows that only the ML model and CLM simulate the realistic interannual variability of fire emissions for most of the subregions (r >0.95 for ML and r =0.14~0.70 for CLM), except for Mediterranean California, where all the models perform poorly (r =0.74 for ML and r <0.30 for the FireMIP models). Regarding seasonality, most models capture the peak emission in July over the western US. However, all models except for the ML model fail to reproduce the bimodal peaks in July and October over Mediterranean California, which may be explained by the smaller wind speeds of the atmospheric forcing data during Santa Ana wind events and limitations in model parameterizations for capturing the effects of Santa Ana winds on fire activity. Furthermore, most models struggle to capture the spring peak in emissions in the southeastern US, probably due to underrepresentation of human effects and the influences of winter dryness on fires in the models. As for extreme events, both the ML model and CLM successfully reproduce the frequency map of extreme emission occurrence but overestimate the number of months with extremely large fire emissions. Comparing the fire PM 2.5 emissions from the ML model with process-based fire models highlights their strengths and uncertainties for regional analysis and prediction and provides useful insights into future directions for model improvements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A review of future weather data for assessing climate change impacts on buildings and energy systems

The effectiveness of climate change impact assessments and the development of adaptation strategies depend on the availability of high-quality future weather data. However, significant gaps exist between the needs of the energy research community and the focus of the climate modeling community, primarily due to a historical lack of communication and collaboration between the two groups. Here, to address this issue, this work provides a comprehensive overview of the critical aspects involved in creating future weather data for building and energy system modeling, including emissions scenarios, general circulation models, downscaling methods, categories of future weather data, and uncertainties in climate simulations. Moreover, it critically evaluates the applicability and suitability of various types of future weather data in five key application scenarios: energy use analysis, resilience analysis, HVAC design, utility-scale analysis, and renewable energy analysis. Finally, this work presents recommendations for high-level actions and research directions to foster collaboration between the energy research and climate modeling communities and to promote the integration of future weather data into energy codes and the design practices of buildings and energy systems.

Climate change↗
DOE OSTI2025

Blinding scheme for the scale-dependence bias signature of local primordial non-Gaussianity for DESI 2024

The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics is important relative to the cosmological signal, and a detection of that signal would have significant implications. This is the case for local primordial non-gaussianity, as probed by the scale-dependent bias of the galaxy power spectrum at large scales that are heavily sensitive to the dependence of the target selection on the imaging quality, known as imaging systematics. We propose a blinding method for the scale-dependent bias signature of local primordial non-gaussianity at the density field level which consists in generating a set of weights for the data that replicate the scale-dependent bias. The applied blinding is predictable, and can be straightforwardly combined with other catalog-level blinding procedures that have been designed for the baryon acoustic oscillation and redshift space distortion signals. The procedure is validated through simulations that replicate data from the first year of observation of the Dark Energy Spectroscopic Instrument, but may find applications to other upcoming spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Penetration Through Slots in Cylindrical Cavities Operating at Fundamental Cavity Modes

In this article, we examine the coupling into an electrically short azimuthal slot on a cylindrical cavity operating at fundamental cavity modal frequencies. We first develop a matched bound formulation through which we can gather information for maximum achievable levels of interior cavity fields. Actual field levels are below this matched bound; therefore, we also develop an unmatched formulation for frequencies below the slot resonance to achieve a better insight on the physics of this coupling. Good agreement is observed between the unmatched formulation, full-wave simulations, and experimental data, providing a validation of our analytical models. We then extend the unmatched formulation to treat an array of slots, found again in good agreement with full-wave simulations. These analytical models can be used to investigate ways to mitigate electromagnetic interference and electromagnetic compatibility effects within cavities.

42 ENGINEERING↗
DOE OSTI2022

Optimizing Error-Bounded Lossy Compression for Scientific Data With Diverse Constraints

Vast volumes of data are produced by today's scientific simulations and advanced instruments. These data cannot be stored and transferred efficiently because of limited I/O bandwidth, network speed, and storage capacity. Error-bounded lossy compression can be an effective method for addressing these issues: not only can it significantly reduce data size, but it can also control the data distortion based on user-defined error bounds. In practice, many scientific applications have specific requirements or constraints for lossy compression, in order to guarantee that the reconstructed data are valid for post hoc analysis. For example, some datasets contain irrelevant data that should be isolated in particular and users often have intuition regarding value ranges, geospatial regions, and other data subsets that are crucial for subsequent analysis. Existing state-of-the-art error-bounded lossy compressors, however, do not consider these constraints during compression, resulting in inferior compression ratios with respect to user's post hoc analysis, due to the fact that the data itself provides little or no value for post hoc analysis. In this work we address this issue by proposing an optimized framework that can preserve diverse constraints during the error-bounded lossy compression, e.g., cleaning the irrelevant data, efficiently preserving different precision for multiple value intervals, and allowing users to set diverse precision over both regular and irregular regions. We perform our evaluation on a supercomputer with up to 2,100 cores. Experiments with six real-world applications show that our proposed diverse constraints based error-bounded lossy compressor can obtain a higher visual quality or data fidelity on reconstructed data with the same or even higher compression ratios compared with the traditional state-of-the-art compressor SZ. Furthermore, our experiments also demonstrate very good scalability in compression performance compared with the I/O throughput of the parallel file system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

In-situ imaging of nucleation and growth of superlattices from nanoscale colloidal nanoparticles

Nanoparticle-based superlattices have attracted extensive research efforts due to their versatile composition and packing-dependent optical, magnetic, electronic, catalytic, mechanical properties. Yet the experimental understanding of the thermodynamics and pathways of the crystallization processes for these systems, which govern their ultimate morphology, size, and packing structure, has for long remained underexplored. Here, this minireview highlights recent integration of liquid-phase transmission electron microscopy, which allows for real-space, real-time direct imaging of the crystallization processes at the single-particle level, with simulation and the conceptual framework established by pioneers such as Dr. A. A. Chernov. Different aspects of superlattice formation were revealed by direct imaging for the first time, such as the existence of prenucleation precursor in nonclassical crystallization, surface energy-dependent growth, and coalescence, allowing for the charting of phase coordinates and thermodynamic quantities at the nanoscale. We discuss the similarities and differences of crystal growth behaviors in atomic and nanoparticle systems, as well as engineering opportunities to fashion, shape, and achieve quality control in superlattice formation from nanoparticles, for next-generation optical and mechanical metamaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward Mass-Production of Transition Metal Dichalcogenide Solar Cells: Scalable Growth of Photovoltaic-Grade Multilayer WSe2 by Tungsten Selenization

Semiconducting transition metal dichalcogenides (TMDs) are promising for high-specific-power photovoltaics due to their desirable band gaps, high absorption coefficients, and ideally dangling-bond-free surfaces. Despite their potential, the majority of TMD solar cells to date are fabricated in a nonscalable fashion, with exfoliated materials, due to the lack of high-quality, large-area, multilayer TMDs. Here, we present the scalable, thickness-tunable synthesis of multilayer WSe2 films by selenizing prepatterned tungsten with either solid-source selenium at 900 degrees C or H2Se precursors at 650 degrees C. Both methods yield photovoltaic-grade, wafer-scale WSe2 films with a layered van der Waals structure and superior characteristics, including charge carrier lifetimes up to 144 ns, over 14x higher than those of any other large-area TMD films previously demonstrated. Simulations show that such carrier lifetimes correspond to ~22% power conversion efficiency and ~64 W g-1 specific power in a packaged solar cell, or ~3 W g-1 in a fully packaged solar module. The results of this study could facilitate the mass production of high-efficiency multilayer WSe2 solar cells at low cost.

carrier lifetime↗
DOE OSTI2021

A Stochastic Multi-Criteria Decision-Making Algorithm for Dynamic Load Prioritization in Grid-Interactive Efficient Buildings

Increasing deployment of advanced sensing, controls, and communication infrastructure enables buildings to provide services to the power grid, leading to the concept of grid-interactive efficient buildings. Since occupant activities and preferences primarily drive the availability and operational flexibility of building devices, there is a critical need to develop occupant-centric approaches that prioritize devices for providing grid services, while maintaining the desired end-use quality of service. In this paper, we present a decision-making framework that facilitates a building owner/operator to effectively prioritize loads for curtailment service under uncertainties, while minimizing any adverse impact on the occupants. The proposed framework uses a stochastic (Markov) model to represent the probabilistic behavior of device usage from power consumption data, and a load prioritization algorithm that dynamically ranks building loads using a stochastic multi-criteria decision-making algorithm. The proposed load prioritization framework is illustrated via numerical simulations in a residential building use-case, including plug-loads, air-conditioners, and plug-in electric vehicle chargers, in the context of load curtailment as a grid service. Suitable metrics are proposed to evaluate the closed-loop performance of the proposed prioritization algorithm under various scenarios and design choices. Scalability of the proposed algorithm is established via computational analysis, while time-series plots are used for intuitive explanation of the ranking choices.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Evaluating the Effects of Nitrogen Doping and Oxygen Doping on SRF Cavity Performance

Superconducting radiofrequency (SRF) cavities are resonators with extremely low surface resistance that enable accelerating cavities to have extremely high quality factors (Q₀). High Q₀ decreases the capital required to keep the accelerators cold by reducing power loss. The performance of SRF cavities is largely governed by the surface composition of the first §I{100}{nm} of the cavity surface. Impurities such as oxygen and nitrogen have been observed to yield high Q₀, but their precise roles are still being studied. Here, we compare the performance of cavities doped with nitrogen and oxygen in terms of surface composition and heating behavior with field. A simulation of the diffusion of oxygen into the bulk of the cavity was built using COMSOL Multiphysics software. Simulated results were compared to the actual surface composition of the cavities as determined from secondary ion mass spectrometry analysis. Understanding how these impurities affects performance allows us to have further insight into the underlying mechanisms that enable these surface treatments to yield high Q₀.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Instabilities driven by the fundamental crabbing mode

The Electron Ion Collider (EIC) requires crab cavities to correct for crossing angle effects. These cavities produce a time varying horizontal kick that makes the collision head on in the center of mass frame of the bunches. This results in a luminosity that is very close to the the luminosity obtained without a crossing angle. Producing the kick requires a transverse RF field. This field requires a superconducting resonator and reasonable power levels require loaded quality factors of order a million. This results is a very large transverse impedance for the crabbing mode. The horizontal beta function is large and this, coupled with the large transverse impedance of the resonator can lead to strong instabilities. In the first section of the paper a simple model is used to describe the situation without feedback. There are large ranges of tune leading to strong instabilities in both the electron and hadron storage rings of the EIC. These estimates are confirmed with simulations. In the second section, RF feedback is used to reduce the apparent impedance of the crabbing mode. Simulating the actual feedback would require both a significant rewrite of the tracking code and a significant reduction in computational speed. Instead, the apparent impedance is fitted using one pole filters, which allows use of the existing simulation code. In the third section a transverse damping system is outlined.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Synthetic Scientific Image Generation with VAE, GAN, and Diffusion Model Architectures

Generative AI (genAI) has emerged as a powerful tool for synthesizing diverse and complex image data, offering new possibilities for scientific imaging applications. This review presents a comprehensive comparative analysis of leading generative architectures, ranging from Variational Autoencoders (VAEs) to Generative Adversarial Networks (GANs) on through to Diffusion Models, in the context of scientific image synthesis. We examine each model's foundational principles, recent architectural advancements, and practical trade-offs. Our evaluation, conducted on domain-specific datasets including microCT scans of rocks and composite fibers, as well as high-resolution images of plant roots, integrates both quantitative metrics (SSIM, LPIPS, FID, CLIPScore) and expert-driven qualitative assessments. Results show that GANs, particularly StyleGAN, produce images with high perceptual quality and structural coherence. Diffusion-based models for inpainting and image variation, such as DALL-E 2, delivered high realism and semantic alignment but generally struggled in balancing visual fidelity with scientific accuracy. Importantly, our findings reveal limitations of standard quantitative metrics in capturing scientific relevance, underscoring the need for domain-expert validation. We conclude by discussing key challenges such as model interpretability, computational cost, and verification protocols, and discuss future directions where generative AI can drive innovation in data augmentation, simulation, and hypothesis generation in scientific research.

Generative Adversarial Networks↗
DOE OSTI2025

PtCoO 2 for Scaled Interconnects

Copper (Cu) interconnects are an increasingly important bottleneck in integrated circuits due to energy consumption and latency caused by the notable increase in Cu resistivity as dimensions decrease, primarily due to electron scattering at surfaces. Herein, the potential of a directional conductor, PtCoO 2 , which has a low bulk resistivity and a distinctive anisotropic structure that mitigates electron surface scattering is showcased. Thin films of PtCoO 2 of various thicknesses are synthesized by molecular beam epitaxy (MBE) coupled with a postdeposition annealing process and the superior quality of PtCoO 2 films is demonstrated by multiple characterization techniques. The thickness‐dependent resistivity curve illustrates that PtCoO 2 significantly outperforms effective Cu (Cu with TaN barriers) and Ru in resistivity below 20.0 nm with a more than 6x reduction compared to effective Cu below 6.0 nm, having a value of only 6.32 μΩ cm at 3.3 nm. It is determined that grain boundary scattering can still be improved for even lower resistivities in this material system through a combination of experiments and theoretical simulations. PtCoO 2 is therefore a highly promising alternative material for future interconnect technologies promising lower resistivities, better stability, and significant improvements in energy efficiency and latency for advanced integrated circuits.

Li, Yansong [Department of Electrical Engineering ↗
DOE OSTI2020

The Evolving Distribution of Relative Humidity Conditional Upon Daily Maximum Temperature in a Warming Climate

Abstract The impacts of heat waves in a warming climate depend not only on changing temperatures but also on changing humidity. Using 35 simulations from the Community Earth System Model Large Ensemble (CESM LENS), we investigate the long‐term evolution of the joint distribution of summer relative humidity (RH) and daily maximum temperature ( Tmax ) near four U.S. cities (New York City, Chicago, Phoenix, and New Orleans) under the high‐emissions Representative Concentration Pathway (RCP) 8.5. We estimate the conditional quantiles of RH given Tmax with quantile regression models, using functions of temperature for each city in July for three time periods (1990–2005, 2026–2035, and 2071–2080). Quality‐of‐fit diagnostics indicate that these models accurately estimate conditional quantiles for each city. As expected, each quantile of Tmax increases from 1990–2005 to 2071–2080, while mean RH decreases modestly. Conditional upon a fixed quantile of Tmax , the median and high quantiles of RH decrease, while those of the Heat Index (HI) and dew point both increase. This result suggests that, despite a modest decrease in median relative humidity, heat stress measured by metrics considering both humidity and temperature in a warming climate will increase faster than that measured by temperatures alone would indicate. For a fixed Tmax , the high quantiles of RH (and thus of HI and dew point) increase from 1990–2005 to 2071–2080 in all four cities. This result suggests that the heat stress of a day at a given Tmax will increase in a warming climate due to the increase of RH.

Yuan, Jiacan↗