Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistics of extremes”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

The Seyfert II nature of the IRAS source FSC 10214+4724

We have observed the rest-frame optical and UV spectra of the luminous, high redshift Infrared Astronomy Satellite (IRAS) source FSC 10214+4724. We find the optical emission lines to be characterized by ratios similar to those found in Seyfert II galaxies. We support the conclusion of previous work that the UV emission lines are similar to those attributed to Active Galactic Nuclei (AGN) rather than to star formation. The ratio H alpha/H beta greater than or equal to 20 (2 sigma lower limit) implies substantial reddening of the narrow line region with A(sub V) greater than 5.5, sufficient to hide a broad line region in our H alpha observations. Given this large inferred reddening and the strength of the UV continuum and emission lines, we conclude (as have others) that simple screen models of reddening are not appropriate for this object. These properties are very similar to those of the infrared luminous galaxies at lower redshift, suggesting that FSC 10244+4724 is the luminous extreme of the same population. We also present H band (1.6 micrometer) imaging polarimetry observations and find that the rest-frame optical emission is unpolarized (P = 3.2% +/- 22.0%). This deep image of the field shows FSC 10214+4724 to possess an unresolved core, with several companions located within 10 sec of the point source. We find it unlikely that this group of objects is physically associated with FSC 10214+4724 at z = 2.3, and we argue that their magnitudes and colors are more consistant with those expected with those expected for galaxies in a foreground group. While galaxy number counts would suggest that such a projection has a low probability of being observed randomly, a foreground group might gravitational lens the z = 2.3 source, making such random statistics inappropriate, and contribute to the large observed luminosity of FSC 10214+4724. Comparison of H band images taken on two occasions one year apart show that FSC 10214+4724 had varied by 0.16 +/- 0.03 mag relative to a nearby star during that time. The UV-optical emission line spectra, the bright dereddened continuum magnitude (K less than 12) of the unresolved core, its extremely large bolometric luminosity, strong UV polarization, and possible variability all suggest that FSC 10214+4724 harbors an obscured AGN. Given that we observe no spectral features attributable to star formation, it appears either that star formation in FSC 10214+4724 is more highly obscured at optical and UV wavelengths than the narrow line region of the AGN, or that the AGN (rather than star formation) dominates the large luminosity of this object.

Elston, Richard↗
DOE OSTI2022

Simulating Impacts of Extreme Events on Grids with High Penetrations of Wind Power Resources

As extreme weather events become more frequent and intense, the demand for connecting grid operation and infrastructure planning with extreme event models will increase as well. We present a methodology for creating damage contingencies and scenarios for electric transmission grids during a hurricane strike. Included is an example case study: Hurricane Dolly damaging a synthetic 2000 bus test system during its landing in Southern Texas. Using WIND Toolkit meteorological data in conjunction with fragility curves for various electric grid elements, we generate stochastic damage scenarios that can be used for short- and long-term planning problems including emergency asset management. We perform statistical analysis of damages and quantify topological effects on example synthetic grid. Also, we investigate loss-of-load events during two-day economic dispatch experiment. Finally, we point out various shortcomings of our method and suggest how it can be improved.

data-driven forecasting↗
DOE OSTI2021

Considering uncertainties expands the lower tail of maize yield projections

Crop yields are sensitive to extreme weather events. Improving the understanding of the mechanisms and the drivers of the projection uncertainties can help to improve decisions. Previous studies have provided important insights, but often sample only a small subset of potentially important uncertainties. Here we expand on a previous statistical modeling approach by refining the analyses of two uncertainty sources. Specifically, we assess the effects of uncertainties surrounding crop-yield model parameters and climate forcings on projected crop yield. We focus on maize yield projections in the eastern U.S.in this century. We quantify how considering more uncertainties expands the lower tail of yield projections. We characterized the relative importance of each uncertainty source and show that the uncertainty surrounding yield model parameters is the main driver of yield projection uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2005

Studies of Coronae and Large Volcanoes on Venus: Constraining the Diverse Outcomes of Small-Scale Mantle Upwellings on Venus

Proxemy Research had a grant from NASA to perform science research on upwelling and volcanism on Venus. This was a 3 year Planetary Geology and Geophysics grant to E. Stofan, entitled Coronae and Large volcanoes on Venus. This grant closes on 12/31/05. Here we summarize the scientific progress and accomplishments of this grant. Scientific publications and abstracts of presentations are indicated in the final section. This was a very productive grant and the progress that was made is summarized. Attention is drawn to the publications and abstracts published in each year. The proposal consisted of two tasks, one examining coronae and one studying large volcanoes. The corona task (Task 1) consisted of three parts: 1) a statistical study of the updated corona population, with Sue Smrekar, Lori Glaze, Paula Martin and Steve Baloga; 2) geologic analysis of several specific groups of coronae, with Sue Smrekar and others; and 3) determining the histories and significance of a number of coronae with extreme amounts of volcanism, with Sue Smrekar. Task 2, studies of large volcanoes, consisted of two subtasks. In the first, we studied the geologic history of several volcanoes, with John Guest, Peter Grindrod, Antony Brian and Steve Anderson. In the second subtask, I analyzed a number of Venusian volcanoes with evidence of summit diking along with Peter Grindrod and Francis Nimmo.

Stofan, Ellen R.↗
DOE OSTI2023

Data Mining and Machine Learning for Power System Monitoring, Understanding, and Impact Evaluation

This chapter presents results from the Big Data analysis framework to improve power system situational awareness and system reliability. For this purpose, a dataset with real-world phasor measurement unit data and historical transmission system outage data has been created and used to carry out the analysis. Several statistical analysis and machine learning methods have been developed and implemented for event and anomaly detection and modeling. Detection and analysis results for actual examples of power system events are presented. Finally, data-driven characterization and risk assessment methods for weather-related extremes in power systems are developed and demonstrated on the Bonneville Power Administration system. These applications demonstrate the capability of Machine Learning (ML) methods to monitor system abnormalities, to predict system events, and to characterize the impact of extreme events on power grid

data mining, power grid, machine learning, anomaly↗
DOE OSTI2022

Chemomechanical weakening of muscovite quantified with in situ liquid nanoindentation

Abstract Chemomechanical weakening of layered phyllosilicate muscovite mica was studied as a function of chemical environment via in situ liquid-nanoindentation under four conditions (dry, deionized water, and two NaCl solutions of different pH). While traditional indentation analyses of layered materials with extreme mechanical anisotropy have been limited due to displacement bursts (pop-ins), here the bursts were used as proxies for delamination, fracture, and spalling events. Since displacement bursts during an indent represent a stochastic process, 120 indents were conducted for each condition to generate statistically significant amounts of data. In total, over 9000 bursts were assessed using a load–displacement threshold criterion, classifying this as a high-throughput nanoscale fracture technique. For each burst, initiation load, initiation displacement, plastic zone volume at initiation, and energy dissipation were analyzed. A power-law relationship between the burst load and burst energy was noted which separated the bursts into two continuous distributions: (1) bursts due only to the mechanics of the indent and (2) bursts due to both the mechanics of the indent and the environment. By using a cumulative probability distribution, it was found that the NaCl solutions decreased the minimum plastic zone volume necessary to initiate a displacement burst by an order of magnitude relative to the dry condition. Finally, the underlying mechanisms explaining the trends in initiation volume as a function of environment were discussed, with a focus on the chemomechanical degradation processes via chemical attack and cation exchange.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
NASA NTRSDate not supplied

Formation and Propagation of Atmospheric River and Its Impact on Extreme Precipitation Events in the North Pacific and the Western North America

Seasonal and interannual evolution patterns of atmospheric rivers (AR) in the North Pacific are examined as a function of the formation region where an AR is first detected using 43-year MERRA-2 reanalysis data. Integrated water vapor transport (IVT) is used to detect AR with latitude dependent thresholds of IVT to better detect AR-like features in the high latitudes. Based on 3-hourly AR statistics, three main AR genesis regions in the North Pacific (i.e., South China Sea (SCS), Western North Pacific (WNP) and Central Pacific Ocean and Hawaiian Islands (CPO) are identified. WNP is the main source of AR with 1475 ARs detected for 43 cold seasons (NDJFM). Over 70% of all AR formed in the WNP has longer than 2-day. On the other hand, AR from CPO tends to have shorter lifetime than those from WNP and SCS. While propagation patterns of AR from SCS and WNP are similar, AR from WNP tends to reach mature phase quicker and shows higher change of impacting west coast of North America. Longevity and strength of AR are also examined based on three large-scale circulation modes (e.g. ENSO, WP, and EAJS) over the East Asia identified from eigen analysis of upper-level zonal winds. During El Nino, the number of AR formed in the western Pacific increased by 20%. First two days, the average size and intensity shows little difference compared with ARs in La Nina years. AR appears to grow in size in El Nino vs. La Nina years. Positive phase of WP correlated with less, but larger and stronger AR formation over WNP and SCS regions. EAJS has little impact in the numbers of AR and its size, but ARs in a stronger EAJS tends to grow larger. Atmospheric circulation associated with the initial formation and propagation of AR from the different regions in the North Pacific as well as its impact on extreme precipitation events in the North Pacific and the west coast of North America will be also discussed.

Integrated water vapor transport↗
DOE OSTI2020

The ongoing need for high-resolution regional climate models: Process understanding and stakeholder information

Regional climate modeling addresses our need to understand and simulate climatic processes and phenomena unresolved in global models. High resolution models are generally more skillful in simulating extremes, such as heavy precipitation, strong winds, and severe storms. In addition, research has shown that fine-scale features such as mountains, coastlines, lakes, irrigation, land use, and urban heat islands can substantially influence a region’s climate and its response to changing forcings. Regional climate simulations explicitly simulating convection are now being performed, providing an opportunity to illuminate new physical behavior that previously was represented by parameterizations with large uncertainties. Regional and global models are both advancing toward higher resolution, as computational capacity increases. However, the resolution and ensemble size necessary to produce a sufficient statistical sample of these processes in global models has proven too costly for contemporary supercomputing systems. Regional climate models are thus indispensable tools that complement global models for understanding regional climate variability and change, and are critical for supporting societal responses to changing climate.

Gutowski, William↗
DOE OSTI2023

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗
DOE OSTI2020

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002

ABSTRACT RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3′ untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5′ and 3′ UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator. IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO 2 and sunlight into useful products.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The Multi-scenario Extreme Weather Simulator: Energy Resilience for Mission Assurance

The Multi-scenario extreme weather simulator (MEWS) is a stochastic weather generation tool. The MEWS algorithm uses 50 or more years of National Oceanic and Atmospheric Association (NOAA) daily summaries [1] for maximum and minimum temperature and NOAA climate norms [2] to calculate historical heat wave and cold snap statistics. The algorithm takes these statistics and shifts them according to multiplication factors provided in the Intergovernmental Panel on Climate Change (IPCC) physical basis technical summary [3] for heat waves.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

AMPTE/CCE CHEM observations of the energetic ion population at geosynchronous altitudes

The paper presents results of a statistical study of average characteristics of the energetic ion population at geosynchronous altitudes, using energetic-ion (1-300 keV/e) measurements from the CHEM spectrometer aboard the AMPTE Charge Composition Explorer between January 1985 and June 1987. Data were sorted into four MLT groups and two extreme geomagnetic activity levels ('very quiet' for AE less than 30 nT and 'very active' for AE greater than 700 nT). A clear quiet-time dayside feature found in the measurements was a dip in H(+) and He(2+) spectra, at 6.6 keV/e in the prenoon sector and at 3.5 keV/e in the postnoon sector. During active times, the ion fluxes increased (except for He(+)), and the O(+) contribution to the total energy density increased dramatically. The pitch angle distributions were normal during quiet times and isotropic or field-aligned during active times.

Daglis, Ioannis A.↗
NASA NTRS2011

Northern Hemisphere Atmospheric Influence of the Solar Proton Events and Ground Level Enhancement in January 2005

Solar eruptions in early 2005 led substantial barrage of charged particles on the Earth's atmosphere during the January 16-21 period. Proton fluxes were greatly increased during these several days and led to the production ofHO(x)(H, OH, BO2)and NO(x)(N, NO, NO2), which then caused the destruction of ozone. We focus on the Northern polar region, where satellite measurements and simulations with the Whole Atmosphere Community Climate Model (WACCM3) showed large enhancements in mesospheric HO(x) and NO(x) constituents, and associated ozone reductions, due 10 these solar proton events (SPEs). The WACCM3 simulations show enhanced short-lived OH throughout the mesosphere in the 60-82.5degN latitude band due to the SPEs for most days in the Jan.16-2l,2005 period, in reasonable agreement with the Aura Microwave Limb Sounder (MLS) measurements. Mesospheric HO2 is also predicted to be increased by the SPEs, however, the modeled HO2 results are somewhat larger than the MLS measurements. These HO(x) enhancements led to huge predicted and MLS-measured ozone decreases of greater than 40% throughout most of the Northern polar mesosphere during the SPE period. Envisat Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) measurements of hydrogen peroxide (H2O2) show increases throughout the stratosphere with highest enhancements of about 60 ppt y in the lowermost mesosphere over the Jan. 16-18, 2005 period due to the solar protons. WACCM3 predictions indicate H2O2 enhancements over the same time period of more than twice that amount. Measurements of nitric acid (HNO3) by both MLS and MIPAS show an increase of about 1 ppbv above background levels in the upper stratosphere during January 16-29, 2005. WACCM3 simulations show only minuscule HNO3 changes in the upper stratosphere during this time period. However due to the small loss rates during winter, polar mesospheric enhancements of NO(x) are computed to be greater than 50 ppbv during the SPE period. Computed NO(x)increases, which were statistically significant at the 95% level, lasted about a month past the SPEs. The SCISAT-I Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) NO(x) measurements and MIPAS NO, measurements for the polar Northern Hemisphere are in reasonable agreement with these predictions. An extremely large ground level enhancement (GLE) occurred during the SPE period on January 20, 2005. We find that protons of energies 300 to 20,000 MeV, not normally included in our computations, led to enhanced lower stratospheric odd nitrogen concentrations of less than 0.1% as a result of this GLE.

Jackman, C. H.↗
NASA NTRS2018

Determination of Uncertainties for Analytically Derived Material Properties to Be Used in Monte Carlo Based Orion Heatshield Sizing

Ablative materials are often used for spacecraft heatshields to protect underlying structures from the extreme environments associated with atmospheric reentry. NASA's Orion EM-1 capsule has been designed to use a molded Avcoat material system. In order to determine the required heatshield thickness, a Monte Carlo approach to the sizing process was proposed. To perform the Monte Carlo simulation, statistical uncertainties on all material property input parameters were required. Obtaining these values for measured properties is straightforward, however input parameters that are derived analytically have historically used uncertainties based on engineering judgment. A MATLAB program was created to use laboratory generated thermogravimetric analysis (TGA) data to calculate uncertainties on the Arrhenius parameters for molded Avcoat. Uncertainties associated with the normalized ablation rate and pyrolysis gas enthalpy were also generated using a wrapper script and the ACE code. These uncertainties could then be tied directly to measured values of individual elemental constituents. The resulting uncertainty values will allow for a probabilistic sizing approach on molded Avcoat with a higher level of confidence in the input parameters.

material property uncertainties↗
DOE OSTI2021

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗