Search NASA⌕ Search

SEARCH · Search NASA

Results for “Long Time Integration”

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

Detection and Diagnosis of Data Integrity Attacks in Solar Farms Based on Multilayer Long Short-Term Memory Network

Photovoltaic (PV) systems are becoming more vulnerable to cyber threats. In response to this emerging concern, developing cyber-secure power electronics converters has received increased attention from the IEEE Power Electronics Society that recently launched a cyber-physical-security initiative. Here this letter proposes a deep sequence learning based diagnosis solution for data integrity attacks on PV systems in smart grids, including dc–dc and dc–ac converters. The multilayer long short-term memory networks are used to leverage time-series electric waveform data from current and voltage sensors in PV systems. The proposed method has been evaluated in a PV smart grid benchmark model with extensive quantitative analysis. As a comparison, we have evaluated classic data-driven methods, including K-nearest neighbor, decision tree, support vector machine, artificial neural network, and convolutional neural network. Comparison results verify performances of the proposed method for detection and diagnosis of various data integrity attacks on PV systems.

42 ENGINEERING↗
DOE OSTI2026

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Spatio-Temporal Denoising Graph Autoencoders with Data Augmentation for Missing Photovoltaic Data Imputation

The integration of the global Photovoltaic (PV) market with real time data-loggers has enabled large scale PV data analytical pipelines for power forecasting and long-term reliability assessment of PV fleets. Nevertheless, the performance of PV data analysis heavily depends on the quality of PV timeseries data. This paper proposes a novel Spatio-Temporal Denoising Graph Autoencoder (STD-GAE) framework to impute missing PV Power Data. STDGAE exploits temporal correlation, spatial coherence, and value dependencies from domain knowledge to recover missing data. It is empowered by two modules. (1) To cope with sparse yet various scenarios of missing data, STD-GAE incorporates a domain-knowledge aware data augmentation module that creates plausible variations of missing data patterns. This generalizes STD-GAE to robust imputation over different seasons and environment. (2) STD-GAE nontrivially integrates spatiotemporal graph convolution layers (to recover local missing data by observed “neighboring” PV plants) and denoising autoencoder (to recover corrupted data from augmented counterpart) to improve the accuracy of imputation accuracy at PV fleet level. We have evaluated our proposed model on two realworld PV datasets. Experimental results show that STD-GAE can achieve a gain of 43.14% in imputation accuracy and remains less sensitive to missing rate, different seasons, and missing scenarios, compared with state-of-the-art data imputation methods such as MIDA and LRTC-TNN.

Fan, Yangxin↗
DOE OSTI2026

Comparability of Liquid Chromatography Tandem Mass Spectrometry Analysis of Dissolved Organic Matter across Laboratories

Non-targeted liquid chromatography tandem highresolution mass spectrometry (LC−MS/MS) is increasingly applied for the structure-resolved chemical analysis of dissolved organic matter (DOM). With new developments in MS instrumentation and analysis software, the approach has gained substantial momentum over the past decade. However, achieving high-quality analytical data that is reproducible and comparable across laboratories can be a bottleneck in non-targeted metabolomics and organic matter chemical analysis, especially for data reuse in repository-scale analyses. Understanding the capabilities as well as challenges of comparing LC−MS/MS data from different laboratories is necessary for inferring global trends from public data sets. To illuminate instrumentation factors that drive differences and variability, we used a standardized data analysis pipeline, including classical (CMN) and featurebased molecular networking (FBMN), to analyze data from a ring trial by 24 laboratories on identical sample sets of algal and DOM extracts that were mixed in predefined concentrations and spiked with standards. Our results showed that data sets from similar mass spectrometer types with unified instrument parameters were qualitatively comparable, resolving the same general trends and shared mass spectral features. Interlaboratory comparability was best for high-intensity features, while low-intensity features showed greater detection variability. Our analysis also highlights challenges when comparing data from instruments with different acquisition rates or operating with less standardized methods. Lastly, we provide recommendations for data integration, public data sharing, standardization, and best practices for standardized LC−MS/MS data acquisition, which will be critical for long-term time series and intercomparability of DOM chemical analyses.

DOM↗
DOE OSTI2020

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗
DOE OSTI2022

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗
DOE OSTI2024

Replica symmetry breaking in spin glasses in the replica-free Keldysh formalism

We show that the algebra of Parisi ultrametric matrices is recovered by the real-time, replica-free, Dyson-Keldysh equations of infinite-range quantum spin glasses in the late time glassy limit. This connects to earlier results on classical and quantum systems showing how ultrametricity emerges from the persistent slow aging dynamics of the glass phase. The stationary spin glass state thereby spontaneously breaks thermal symmetry, or the Kubo-Martin-Schwinger relation of a state in global thermal equilibrium. We describe the Keldysh path integral of the infinite-range Ising model in transverse and longitudinal fields, and in the context of the Landau expansion of the action functional, show how the long-time limit connects to the full replica symmetry breaking obtained in the equilibrium formalism. We also illustrate our formalism by applying it to the spherical quantum p p -spin model, which only exhibits one-step replica symmetry breaking.

Lang, Johannes (ORCID:0000000283533392)↗
DOE OSTI2025

Key Performance Indicators for Vehicle Grid Integration

Electric vehicle (EV) sales account for a rapidly growing portion of the light-duty vehicle market and a portion of medium and heavy-duty fleet vehicles. However, in many locations, charging stations will require costly utility grid upgrades with long lead times. Today, there are a few methods of smart charge management (SCM) which can reduce the costs and wait times for electric vehicle supply equipment (EVSE) interconnection approvals, as well as reduce impacts of the charging stations on the grid and on EV driver transportation costs. It is crucial for the EV charging industry to understand the vehicle grid integration (VGI) requirements for EVSE to prevent adverse impacts from EVSE interconnections and assure that charging loads and interconnections are affordable.

33 - ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Implementation of Distributed Memory Computing in MOSAIC to Enable Large 3D Simulations of Irradiated Concrete

The concrete biological shield (CBS) of light-water reactors protects workers and the surrounding environment by absorbing neutron and gamma irradiation emitted from the reactor core. The radiation dose increases with the CBS’s operational time and, in the long term, becomes significant enough to raise the question of irradiation effects on concrete—and particularly on the structural integrity of the CBS. Irradiation-induced damage has been identified as one of the main degradation mechanisms in the CBS. Neutron radiation causes the swelling of aggregate-forming minerals at different rates and amplitudes depending on the mineral’s nature. Silicate-bearing minerals such as quartz are particularly sensitive to neutron radiation and experience up to 17.8% volumetric expansion. Aggregates comprise several minerals with different orientations and are, therefore, subject to cracking as a result of mismatch strains. Additionally, the swelling of aggregates creates significant stresses in the surrounding cement paste matrix, which also results in crack formation. In parallel with the collection of characterization and irradiation test data, development of modeling and simulation tools for irradiated concrete is ongoing with the support of the US Department of Energy Office of Nuclear Energy’s Light Water Reactor Sustainability (LWRS) program. This effort resulted in the development and application of the fast-Fourier transform (FFT)–based code Microstructure-Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) at Oak Ridge National Laboratory.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Evaluation of Techniques for Vegetation Removal from UAV-Based Photogrammetric DSM near WIPP Land Withdrawal Act Boundary, NM

The Waste Isolation Pilot Plant is located within New Mexico's karst landscape and is the United States' only deep geologic radioactive waste repository which isolates transuranic waste from defense activities underground in a bedded salt formation. Characteristic karst topography in this region is made up of material (such as limestone) that can be dissolved by water over a long period of time, and includes features such as springs, surface streams, sinking streams, caves, and sinkholes. Karst features can potentially impact the integrity of the waste repository in the future and compromise WIPP's performance, particularly when coupled with incompatible land-use activities within the Land Withdrawal Act (LWA) boundary as well as increased water withdrawals outside the LWA. In order to better capture the ground surface topography to support hydrology modeling efforts, the U.S. DOE has enlisted FIU to develop a high-resolution digital elevation model (DEM). In a previous field study, a high-resolution digital surface model (DSM) of an adjacent representative site, Basin 6, was developed by photogrammetry using aerial images captured with unmanned aerial vehicles. The DSM generated was based on point clouds representing the entire landscapes, including points of terrain, vegetation and infrastructure. These above-ground features need to be removed to create a DSM/DEM with good representation of the bare ground. In recent years, many techniques including both pre- and post-processing of DSMs, were developed for vegetation removal with varied success. This research focuses on applying and comparing various vegetation removal methods using point clouds generated of the Basin 6 pilot study area. The 3D RGB-based point cloud classification using Python scripts was a modification of the procedure by Themistocleous (2019). The method using Liblas classifies the point cloud three-dimensionally, whereas the method by Themistocleous (2019) is two-dimensional. The preliminary result of the method using Liblas, however, still had a large amount of vegetation remaining compared to the Pix4D Machine Learning method. The vegetation layer remains highly visible in the DEMs processed using the vegetation indices as opposed to the Pix4D Machine Learning method, which produces a DEM where mostly bare ground is seen. Further modifications to the input parameters are required to improve the vegetation vs. bare ground classification using this approach. The next step will be to improve the method using Liblas and to also test other vegetation removal methods found in the literature review to compare them against the two described here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

On Building Predictive Digital Twin Incorporating Wave Predicting Capabilities: Case Study on UMaine Experimental Campaign - FOCAL

The response of floating wind turbines (FWT) are susceptible to stochastic wave variations. For the optimal operation of FWT, a comprehensive understanding of the phaseresolved wave dynamics and the consequential system response is crucial for real-time monitoring and control. A multi-variate, multi-step, long short term memory (MLSTM), a type of recurrent neural network (RNN) is used to capture complex system dynamics for real-time application. Results indicate that the integration of a wave prediction-reconstruction (WRP) model substantially enhances prediction accuracy by 50% on average relative to the baseline model. The improvement is consistent across various wave extremity and prediction horizons, thereby significantly broadening the scope for timely and precise predictive capabilities.

17 WIND ENERGY↗
DOE OSTI2022

An extrapolation method for strain ranges and hold times in developing the EPP+SMT creep-fatigue design curves for Alloy 617

Experimental and numerical studies in developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY2022. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue (CF) design curves at long hold times and low strain ranges. In this study, the available CF failure data on Alloy 617 at 950°C were analyzed to determine a set of CF failure criteria. At very low strain ranges and long hold times, CF failure data are not accessible by experiments because of the extraordinarily long test durations and the inability of the test machines to accurately control these small strain ranges. A CF experimental approach with the concept of block-strain range CF testing protocol was developed. Tests using this protocol were conducted to generate the needed information for calibrating material parameters of the numerical material models. The Time Fraction based method and Dissipated Energy method were used to extrapolate the CF life curves to low strain ranges and long hold times. Based on the new experimental approach and CF life prediction methods, the CF life curves with various hold times were developed for Alloy 617 at 950°C. In addition, an experiment was designed and is being performed to verify the predicted CF curves at 950°C. The extrapolation procedure will be applied at lower temperatures to complete the development of the EPP+SMT CF design curves for Alloy 617 in F2023.

36 MATERIALS SCIENCE↗
DOE OSTI2022

36 Hz integral linewidth laser based on a photonic integrated 4.0 m coil resonator

Laser stabilization sits at the heart of many precision scientific experiments and applications, including quantum information science, metrology, and atomic timekeeping. Many of these systems narrow the laser linewidth and stabilize the carrier by use of Pound–Drever–Hall (PDH) locking to a table-scale, ultrahigh quality factor (Q), vacuum spaced Fabry–Perot reference cavity. Integrating these cavities to bring characteristics of PDH stabilization to the chip scale is critical to reducing their size, cost, and weight, and enabling a wide range of portable and system-on-chip applications. We report a significant advance in integrated laser linewidth narrowing, stabilization, and noise reduction by use of a photonic integrated 4.0 m long coil resonator to stabilize a semiconductor laser. We achieve a 36 Hz 1 / π <#comment/> -integral linewidth, Allan deviation of 1.8 × <#comment/> 10 − <#comment/> 13 at 10 ms measurement time, and a 2.3 kHz/s drift—to the best of our knowledge, the lowest integral linewidth and highest stability demonstrated for an integrated waveguide reference cavity. This performance represents over an order of magnitude improvement in integral linewidth and frequency noise over previous integrated waveguide PDH stabilized reference cavities and bulk-optic and integrated injection locked approaches, and over two orders of magnitude improvement in frequency and phase noise than integrated injection locked approaches. Two different wavelength coil designs are demonstrated, stabilizing lasers at 1550 nm and 1319 nm. The resonator is bus-coupled to a 4.0 m long coil, with a 49 MHz free spectral range, mode volume of 1.0 × <#comment/> 10 10 µ <#comment/> m 3 , and 142 million intrinsic Q , fabricated in a CMOS compatible, ultralow loss silicon nitride waveguide platform. Our measurements and simulations show that the thermorefractive noise floor for this particular cavity is reached for frequencies down to 20 Hz in an ambient environment with simple passive vibration isolation and without vacuum or thermal isolation. The thermorefractive noise limited performance is estimated to yield an 8 Hz 1 / π <#comment/> -integral linewidth and Allan deviation of 5 × <#comment/> 10 − <#comment/> 14 at 10 ms, opening a stability regime that heretofore has been available only in fundamentally non-integrated systems. These results demonstrate the potential to bring the characteristics of laboratory-scale stabilized lasers to the integrated, wafer-scale compatible chip scale, and are of interest for a number of applications in quantum technologies and atomic, molecular, and optical physics, and with further developments below 10 Hz linewidth, can be highly relevant to ultralow noise microwave generation.

Liu, Kaikai (ORCID:0000000213609646)↗
DOE OSTI2023

High-density integrated delay line using extreme skin-depth subwavelength grating waveguides

Optical delay lines control the flow of light in time, introducing phase and group delays for engineering interferences and ultrashort pulses. Photonic integration of such optical delay lines is essential for chip-scale lightwave signal processing and pulse control. However, typical photonic delay lines based on long spiral waveguides require extensively large chip footprints, ranging from mm 2 to cm 2 scales. Here we present a scalable, high-density integrated delay line using a skin-depth engineered subwavelength grating waveguide, i.e., an extreme skin-depth (eskid) waveguide. The eskid waveguide suppresses the crosstalk between closely spaced waveguides, significantly saving the chip footprint area. Our eskid-based photonic delay line is easily scalable by increasing the number of turns and should improve the photonic chip integration density.

Ahmed, Ishtiaque↗
DOE OSTI2023

Looking at supersymmetric black holes for a very long time

We study correlation functions for extremal supersymmetric black holes. It is necessary to take into account the strongly coupled nature of the boundary supergraviton mode. We consider the case with {\cal N}=2 supercharges which is the minimal amount of supersymmetry needed to give a large ground state degeneracy, separated from the continuum. Using the exact solution for this theory we derive formulas for the two point function and we also give integral expressions for any n n -point correlator. These correlators are time independent at large times and approach constant values that depend on the masses and couplings of the bulk theory. We also explain that in the non-supersymmetric case, the correlators develop a universal time dependence at long times. This paper is the longer companion paper of arXiv:2207.00407.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Integrate FARM with PID controllers: IES Simulation Ecosystem Control System Development

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM was designed to support the HERON software module in the solution of the optimal dispatch problem for IES units. As the result of HERON-FARM dispatch simulation, the set-point trajectories are optimized to meet constraints on both the production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and the process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.) at a coarse time resolution (every 10 or 100 seconds) over long time horizons (several days or weeks). In case the operational constraints need to be met at finer time resolution, the computational burden of HERON-FARM would linearly increase with the sampling rate, and sub-optimal solutions might be obtained. System responses characterized by overshoots and damped oscillations temporarily violating the imposed constraints might occur during abrupt power transients. In this report, a hierarchical control system architecture for the operation of the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility constructed at INL was proposed. First, the preliminary studies on the proposed control strategy for operating the facility and the designed PI controllers were reviewed. In particular, the current approach for generating the set-point trajectories was studied, and its limits were identified. To this aim, the inclusion of a Supervisory Control layer embedding a modified version of the FARM algorithm for preserving the system safe operation over both long and real-time horizons was proposed. In this way, FARM would be applied twice, i.e., the original version (“FARM Validator”) aiding the solution of the power dispatch problem, and the modified version (“FARM Supervisory” coordinating the PI controllers to address the real-time control tasks. Despite the kernel of the two modules is the same algorithm, their roles, tasks, and capabilities are quite different. A detailed description of the role of FARM at addressing low-level control tasks is provided, along with tentative operational procedures for training the models embedded into the algorithm by using the collected experimental data.

42 ENGINEERING↗
DOE OSTI2022

Müllerian mimicry and the coloration patterns of sympatric coral snakes

Abstract Coral snakes in the genus Micrurus are venomous, aposematic organisms that signal danger to predators through vivid coloration. Previous studies found that they serve as models to several harmless species of Batesian mimics. However, the extent to which Micrurus species engage in Müllerian mimicry remains poorly understood. We integrate detailed morphological and geographical distribution data to investigate if coral snakes are Müllerian mimics. We found that coloration is spatially structured and that Micrurus species tend to be more similar where they co-occur. Though long supposed, we demonstrate for the first time that coral snakes might indeed be Müllerian mimics as they show some convergence in coloration patterns. Additionally, we found that the length of red-coloured rings in Micrurus is conserved, even at large geographic scales. This finding suggests that bright red rings may be under more substantial stabilizing selection than other aspects of coloration and probably function as a generalized signal for deterring predators.

Evolutionary Biology↗