Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extreme value”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Leveraging Extremal Dependence to Better Characterize the 2021 Pacific Northwest Heatwave

Abstract In late June, 2021, a devastating heatwave affected the US Pacific Northwest and western Canada, breaking numerous all-time temperature records by large margins and directly causing hundreds of fatalities. The observed 2021 daily maximum temperature across much of the U.S. Pacific Northwest exceeded upper bound estimates obtained from single-station temperature records even after accounting for anthropogenic climate change, meaning that the event could not have been predicted under standard univariate extreme value analysis assumptions. In this work, we utilize a flexible spatial extremes model that considers all stations across the Pacific Northwest domain and accounts for the fact that many stations simultaneously experience extreme temperatures. Our analysis incorporates the effects of anthropogenic forcing and natural climate variability in order to better characterize time-varying changes in the distribution of daily temperature extremes. We show that greenhouse gas forcing, drought conditions and large-scale atmospheric modes of variability all have significant impact on summertime maximum temperatures in this region. Our model represents a significant improvement over corresponding single-station analysis, and our posterior medians of the upper bounds are able to anticipate more than 96% of the observed 2021 high station temperatures after properly accounting for extremal dependence. Supplementary materials accompanying this paper appear online.

Zhang, Likun (ORCID:0000000154584556)↗
DOE OSTI2026

Machine learning method for enforcing variable independence in background estimation with LHC data: ABCDisCoTEC

A novel solution is presented for the problem of estimating the backgrounds of a signal search using observed data while simultaneously maximizing the sensitivity of the search to the signal. The 'ABCD method' provides a reliable framework for background estimation by partitioning events into one signal-enhanced region (A) and three background-enhanced control regions (B, C, and D) via two smoothly varying, statistically independent variables. In practice, even slight correlations between the two variables can significantly undermine the method's performance. Thus, choosing appropriate variables by hand can present a formidable challenge, especially when background and signal differ only subtly. To address this issue, the ABCD with distance correlation (ABCDisCo) method was developed to construct two learned variables via a neural network trained to provide strong signal-background discrimination with small values of the distance correlation (DisCo) measure between the two learned variables. However, relying solely on minimizing the DisCo can result in learned variables that may not have distributions of background events that are smoothly varying and localized at extreme values, as necessary for the validity of the background estimation. The ABCDisCo training enhanced with closure (ABCDisCoTEC) method is introduced to solve this issue by directly minimizing the nonclosure, expressed as a dedicated differentiable loss term. This extended method is applied to a data set of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS detector at the CERN Large Hadron Collider. Additionally, given the complexity of the minimization problem with constraints on multiple loss terms, the modified differential method of multipliers is applied and shown to greatly improve the stability and robustness of the ABCDisCoTEC method, compared to grid search hyperparameter optimization procedures.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗
DOE OSTI2025

Mechanisms of Projected Changes in Thunderstorm Downburst Environments Across the United States

Responses of downdraft convective available potential energy (DCAPE) to global warming were investigated using the Community Earth System Model (CESM2) under a high‐emission scenario through the year 2100. DCAPE is projected to increase by 5%–12% on average in most areas, independently of wind shear. A diagnostic of downdraft buoyancy is introduced to understand the mechanisms of DCAPE responses. Much of the increase in mean DCAPE is temperature‐driven, with additional contributions from changes in relative humidity and downdraft origin heights. However, extreme values increase at much faster rates than can be explained by local warming. In winter, the latitude of significant DCAPE and CAPE shifts poleward by more than 5° due to larger changes in downburst environments within midlatitude cyclones. The projected increase in cold‐season extremes indicates an interaction between weather events and warming trends that increases the potential for downbursts and straight‐line winds in winter.

Williams, Ian N. [Iowa State Univ., Ames, IA (Unit↗
DOE OSTI2024

Four Years of Atmospheric Boundary Layer Height Retrievals Using COSMIC-2 Satellite Data

This work aimed to study the atmospheric boundary layer height (ABLH) from COSMIC-2 refractivity data, endeavoring to refine existing ABLH detection algorithms and scrutinize the resulting spatial and seasonal distributions. Through validation analyses involving different ground-based methodologies (involving data from lidar, ceilometer, microwave radiometers, and radiosondes), the optimal ABLH determination relied on identifying the lowest refractivity gradient negative peak with a magnitude at least $τ$% times the minimum refractivity gradient magnitude, where $τ$ is a fitting parameter representing the minimum peak strength relative to the absolute minimum refractivity gradient. Different $τ$ values were derived accounting for the moment of the day (daytime, nighttime, or sunrise/sunset) and the underlying surface (land or sea). Results show discernible relations between ABLH and various features, notably, the land cover and latitude. On average, ABLH is higher over oceans (≈1.5 km), but extreme values (maximums > 2.5 km, and minimums < 1 km) are reached over intertropical lands. Variability is generally subtle over oceans, whereas seasonality and daily evolution are pronounced over continents, with higher ABLHs during daytime and local wintertime (summertime) in intertropical (middle) latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

How Frequent Will the Rarest Daily Rainfall Records of Hurricane Ida’s Remnants Be in the Future?

Abstract Gaining continued insights into the impact of global warming on the occurrence of hurricane-associated intense record downpours is essential for building climate resilient communities. This study investigates projected future changes in extreme rainfall over the Northeast United States, as represented by extreme daily amounts during Hurricane Ida in 2021. We used historical control simulations of Weather Research and Forecasting (WRF) Model generated from 40 years of weather events (1980–2014, 12 km) forced by the fifth generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis. These simulations are thermodynamically modified (2060–2100) via an imposed warming for the high-emission scenario of shared socioeconomic pathway (SSP585) from a range of general circulation models. Ground observations from the Global Historical Climatology Network (1950–2014) and WRF simulations (historical, 1980–2014, and future, 2060–2100) are integrated into a nonstationary generalized extreme value (GEV) framework to assess the frequency of Ida’s heaviest daily rain rates under the SSP585 scenario. Results show that Ida’s daily maximum rainfall recorded at different observation locations was higher than the single highest September daily maximum observed (1950–2014) for 5 out of 17 stations (∼30% of the stations). Ida-like extreme daily rain rates are projected to be, on average, more than 2 times more likely to occur at the end of the century in the simulations (with some regions as high as 5 times). This work demonstrates that integrating a high-resolution atmospheric model’s present-day and thermodynamically modified future simulations along with ground observations, within a nonstationary statistical framework, is crucial for understanding changing characteristics of extreme weather events. Significance Statement Daily scale extreme precipitation is expected to become more frequent and severe, as evidenced by observations and model simulations. While it is important to investigate how these intensifying heavy rainfall events affect current engineering standards, fewer studies have contextualized how warming impacts the most extreme rainfall from a single storm event relative to historical heavy downpours. In this study, we focused on the daily extreme rainfall associated with the extratropical transition of Hurricane Ida (2021), particularly over the northeastern United States—some of which exceeded the commonly used hydrologic design criteria for a 100-yr storm. Using a high-resolution atmospheric model simulation, we investigated how continued warming may influence the frequency of such daily rain rates. Under a high-emission scenario, these events are projected to become up to 5 times more likely at the end of the twenty-first century.

Dollan, Ishrat J↗
DOE OSTI2023

Projected changes in extreme streamflow and inland flooding in the mid-21st century over Northeastern United States using ensemble WRF-Hydro simulations

Study region: Northeastern United States (NEUS). Study focus: We investigate the potential impacts of climate change on precipitation, streamflow, and inland flooding in the NEUS during the mid-21st century. Dynamically downscaled climate projections from three global climate models for historical (1995-2004) and future (2045-2054) periods under business-as-usual scenarios were used to force the hydrologic model WRF-Hydro at 200-meter resolution and create ensemble hydrologic simulations. Additionally, an extreme value model was developed to project the risks associated with low-frequency hydrologic events. New hydrological insights for the region: Results from four major watersheds indicate a significantly wetter regime in winter months and potential drier conditions during late spring to early summer. Discharges in fall are projected to decrease in the northern watersheds and increase toward the south. Extreme flow and water depths resulting from extreme inland flooding are projected to increase by 5-20% and > 100%, respectively. The extent of the total flooded area is likely to be 20% greater by the mid-century. These increased risks can be attributed to (i) an approximate 25% increase in decadal mean and > 40% increase in decadal extreme precipitation intensity, (ii) up to 30% lower snow availability and 5-25% higher evapotranspiration throughout the year, and (iii) a projected 5% increase in soil moisture in all seasons except summer. Furthermore, rapid snow melting in winter will likely cause an earlier peak flow in the rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Is land use producing robust signals in future projections from Earth system models, all else being equal?

Abstract We use six Earth system models (ESMs) run under SSP3-7.0, a scenario characterized by a relatively large land use change (LUC) over the 21st century, and under a variant of the same scenario where a significantly different pattern of LUC, taken from SSP1-2.6, was used, all else being equal. Our goal is to identify changes in climate extremes between the two scenarios that are statistically significant and robust across the ESMs. The motivation for this study is to test a long-held assumption of the shared socio-economic pathway-representative concentration pathway (SSP-RCP) scenario framework: that the signal from LUC can be safely disregarded when pairing different SSPs to the compatible RCPs, where compatibility only considers global radiative forcing, predominantly determined by well-mixed greenhouse gasses emissions. We analyze extremes of daily minimum and maximum temperatures and precipitation, after fitting non-stationary generalized extreme value distributions in a way that borrows strength along the length of the simulation (2015–2100) and across initial condition ensembles. We consider changes in the 20 year return levels (RL20s) of these metrics by 2100, and focus on eight locations where LUC is large within each scenario, and strongly differs between scenarios, averaging the RL20s over a neighborhood characterized by the same LUC to enhance the signal to noise. We find that precipitation extremes do not show significant differences attributable to LUC differences. For temperature extremes (cold and hot) results are mixed, with some location-index combination showing significant results for some of the ESMs but not all, and not many coherent changes appearing for indices across regions, or regions across indices. These ESMs are representative of what is typically adopted as the source of climate information for impact studies, when the SSP-RCP framework is put to use. Overall, our analysis suggests that the hypothesis to pair SSPs to RCPs in a flexible fashion is overall defensible. However, the appearance of some coherence in a few locations and for some indices invites further investigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Constraining the Multiscale Structure of Geophysical Fields in Machine Learning: The Case of Precipitation

The use of deep-learning algorithms for estimating the value of geophysical variables from remotely sensed information is rapidly expanding. The typical objective function minimized in such algorithms is the mean square error (MSE), which is known to lead to smooth estimates with compressed dynamical range as compared to the true distribution of the variable of interest. Here, we introduce and evaluate alternative objective functions, focusing on the retrieval of precipitation rates from satellite passive microwave radiometric measurements using a deep convolutional neural network. For this testbed application, the results show that explicitly imposing the preservation of the statistical distribution and spatial wavelet power spectrum of the target variable allows to accurately reproduce extreme values and sharp gradients across multiple scales in this study.

58 GEOSCIENCES↗
DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2025

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗
DOE OSTI2023

Nanoscale fluctuation of stacking fault energy strengthens multi-principal element alloys

Chemical randomness and the associated energy fluctuation are essential features of multi-principal element alloys (MPEAs). Due to these features, nanoscale stacking fault energy (SFE) fluctuation is a natural and independent contribution to strengthening MPEAs. However, existing models for conventional alloys (i.e., alloys with one principal element) cannot be applied to MPEAs. The extreme values of SFEs required by such models are unknown for MPEAs, which need to calculate the nanoscale volume relevant to the SFE fluctuation. In the present work, we developed an analytic model to evaluate the strengthening effect through the SFE fluctuation, profuse in MPEAs. The model has no adjustable parameters, and all parameters can be determined from experiments and ab initio calculations. Further, this model explains available experimental observations and provides insightful guidance for designing new MPEAs based on the SFE fluctuation. It generally applies to MPEAs in random states and with chemical short-range order.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparing the cooling effectiveness of operationalisable urban surface combination scenarios for summer heat mitigation

Extreme summer heat in cities exacerbates the vulnerability of urban communities to heatwaves. Vegetative and reflective urban surfaces can help reduce urban heat. This study investigated the impacts of urban trees, green roofs and cool roofs on heat mitigation during average and extreme summer conditions in temperate oceanic Melbourne, Australia. We simulated the city climate using ‘The Air Pollution Model’ (TAPM) at a 1 km spatial resolution over 10 years, which according to our review of the literature, was the most prolonged period for simulation in Melbourne. During a widespread heatwave event, some of the tested scenarios with combined surface parameters could reduce the extreme values of the energy budget components- sensible heat, latent heat, and storage heat fluxes up to seasonal averages compared to the existing situation for Melbourne (control). The scenario with the highest (reasonable maximum) ground-level vegetation, green roofs, and cool roofs could reduce air temperatures up to 2.4 °C. The simulations suggest that a combined strategy with vegetative and high-albedo surfaces will deliver higher effectiveness with maximum cooling benefits and cost-effectiveness than individual strategies in cities. These results suggest the importance of collaborative strategic planning of urban surfaces to make cities healthier, sustainable, and liveable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2022

Anthropogenic Contributions to the 2021 Pacific Northwest Heatwave

Abstract Daily maximum temperatures during the 2021 heatwave in the Pacific Northwest United States and Canada shattered century old records. Multiple causal factors, including anthropogenic climate change, contributed to these high temperatures, challenging traditional methods of attributing human influence. We demonstrate that the observed 2021 daily maximum temperatures are far above the bounds of Generalized Extreme Value distributions fitted from historical data. Hence, confidence in Granger causal inference statements about the human influence on this heatwave is low. Alternatively, we present a more conditional hindcast attribution study using two regional models. We performed ensembles of simulations of the heatwave to investigate how the event would have changed if it had occurred without anthropogenic climate change and with future warming. We found that global warming caused a ∼0.8°C–1°C increase in heatwave temperatures. Future warming would lead to a ∼5°C increase in heatwave temperature by the end of the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Projecting Changes in the Frequency and Magnitude of Ozone Pollution Events Under Uncertain Climate Sensitivity

Abstract Climate change is projected to worsen ozone pollution over many populated regions, with larger impacts at higher concentrations. More intense and frequent ozone episodes risk setbacks to human health and environmental policy achievements. However, assessing these changes is complicated by uncertain climate sensitivity, closely related to climate model response, and internal variability in simulations projecting climate's influence on air quality. Here, leveraging a global modeling framework that one‐way couples a human activity model, an Earth system model of intermediate complexity, and an atmospheric chemistry model, we investigate the role of climate sensitivity in climate‐induced changes to high ozone pollution episodes in the United States using multiple greenhouse gas emissions scenarios, representations of climate sensitivity, and initial condition members. We bias correct and evaluate historical model simulations, identifying modeled and observed O 3 episodes using extreme value theory, and extend the approach to projections of mid‐ and end‐century climate impacts. Results show that the influence of climate sensitivity can be as significant as that of greenhouse gas emissions scenario absent precursor emissions changes. Climate change is projected to increase the magnitude of the highest annually occurring O 3 concentrations by over 2.3 ppb on average across the U.S. at mid‐century under a high climate sensitivity and moderate emissions scenario, but the increase is limited to less than 0.3 ppb under lower climate sensitivity. Further, we show that areas in the U.S. currently meeting air quality standards risk being pushed into non‐compliance due to a climate‐induced increase in frequency of high ozone days.

Environmental Sciences & Ecology↗
DOE OSTI2022

How galaxies populate halos in extreme density environments: An analysis of the halo occupation distribution in SDSS

Recent works have shown that the properties of galaxy populations in simulated dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments since the dynamical and astrophysical conditions prevailing in these regions may significantly affect the formation and evolution of their halos and residing galaxies, influencing the mean number of satellite galaxies. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the Sloan Digital Sky Server Data Release 12 (SDSS-DR12) and estimate the HOD within these super-structures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within cosmic voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these super-structures. We use a publicly available observational galaxy catalogue with information on redshifts, positions, magnitudes and other astrophysical features to build a volume complete galaxy sample and identify cosmic voids and FVS. Using a publicly available galaxy group catalogue as a proxy to dark matter halos, we compute the HOD within both types of super-structures for different absolute magnitude thresholds. We also study the dependence on the results on the main void and FVS properties, density and volume. We also analysed the main characteristics of the stellar content of galaxies inside these extreme-density regions such as the mean stellar age and the galaxy light concentration index. In all cases, we compare the results with those derived from the Field sample, defined by objects outside both types of environments. Inside cosmic voids, we find a strong decrease in HOD concerning the Field results. In the most extreme cases, the mean number of satellites fall to ~ 50%. Inside FVS, the HOD shows a significant increase to the Field, with a ~ 40% excess in the mean number of satellites. These results are present for the different galaxy luminosity ranges explored. In both environments, the differences with respect to the Field increases for the extreme values of the density environments. However, we obtain no signs of variations related to intrinsic characteristics of the super-structures, indicating that the effects depend mainly on the density of the large-scale environment. In addition, we find that the cumulative distribution of the mean age of stars of the central galaxy also varies in the different regions, this suggests that the history of the formation of the dark matter halos may be different. Finally, we explore the HOD for the 25% youngest (oldest) galaxies, based on the mean age of their stars. We find that for the low-mass groups the youngest galaxies are only present inside voids and are generally central galaxies. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field as previously mentioned. We find that cosmic voids lack a significant fraction of galaxies with the oldest stellar population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗