Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System Simulation”

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

Plastic homogeneity in nanoscale heterostructured binary and multicomponent metallic eutectics: An overview

Heterostructured materials comprised of relatively soft/hard disparate phases typically exhibit composite strengthening but lack plastic deformability at ambient temperatures. However, heterostructured systems comprised of nanoscale phases can simultaneously enhance yield strength and strain hardening, thereby promoting uniform distribution of plastic flow. In this review, the atomic-scale deformation mechanisms in model systems of eutectic alloys, Al-Al 2 Cu and Al-Si, refined to nanoscales via laser rapid solidification are discussed, and compared with literature on multi-component (high entropy) eutectics such as Ni-Al-Fe-based with Cr and/or Co additions. The nano-lamellar Al-Al 2 Cu structures exhibit unit defect mechanisms not reported in monolithic Al 2 Cu intermetallic: localized shear on {0 1 1} and shear-induced faults on {1 2 1} planes, constrained by closely-spaced dislocation arrays in Al confined by Al/Al 2 Cu interfaces. The unexpected plasticity mechanisms are enabled by slip continuity in nanoscale Al-Al 2 Cu eutectics associated with the orientation relationship and interface habit planes. In nano-fibrous Al-Si eutectic, tensile ductility at strength approaching 600 MPa is observed resulting from dislocation plasticity in the nano-Al channels and cracking in Si nanofibers. Molecular dynamics simulations show that Al dislocations easily cross-slip (screw) or climb (edge) along Al-Si interfaces, making slip transmission difficult. The propagation of nano-cracks is suppressed by surrounding strain hardening Al, retaining good ductility of the sample, in spite of lack of direct slip transmission. Finally, the critical unit mechanisms of slip transmission and interface-enabled plasticity observed in nanoscale eutectic binary systems can also explain the strength-ductility relationship in multi-component eutectics and homogeneously distributed plastic flow with increasing microstructural heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterizing the Variation and Covariation of Cloud Microphysical Properties and Implications for Simulation of Subgrid-scale Warm-Rain Processes in Earth System Models (Final DOE-ASR Report)

Warm marine boundary layer (MBL) clouds constitute an important component in the global climate system, and precipitation plays a central role in controlling the water budget, radiative effects, and lifetime of these MBL clouds. Unfortunately, because of the relatively coarse effective grid resolution of the current generation of Earth system models (ESMs), the variety of cloud microphysical processes occurring inside an ESM grid cell are often oversimplified or unconstrained by observations. For example, the warm rain processes (e.g., autoconversion and accretion) are usually parameterized as nonlinear functions of grid-mean cloud properties. Because of the nonlinear nature of these functions, neglecting variability within the ESM grid volume can lead to substantial biases in precipitation production, cloud cover, and surface radiative fluxes. In state-of-the art ESMs, the influence of subgrid-scale variability is represented as an enhancement factor (EF) coefficient to the autoconversion, and accretion rates calculated from the model variables. However, EF is typically taken to be a constant or even used as a knob to tune model cloud properties to match observations, an ad hoc approach that may yield a desired cloud outcome yet introduce compensating errors. In this project, we used the combination of in situ cloud microphysics measurements from the ACE-ENA field campaign and large-eddy simulations (LES) to characterize and understand subgrid-scale variations and co-variations of cloud microphysical properties and use the results to evaluate and improve the representation of subgrid warm-rain processes in ESMs, in particular the EF used to tune the autoconversion and accretion processes. In this final report, we summarize our research activities and main findings in Section 2, provide a list of publications (Section 3) and presentations (Section 4) resulted from our research, and briefly discuss the student activities supported by this project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Using High-Resolution NSRDB Data to Evaluate Cloud Mask Forecast from WRF-Solar EPS

Validating spatiotemporal distributions of cloud forecasts using numerical weather prediction (NWP) models is difficult as this requires high-quality cloud-property at significantly high spatial and temporal resolution over extended periods of time. Observations of cloud properties, such as cloud mask, cloud optical thickness, and cloud type, are vital for assessing the capability of NWP models to forecast various types of clouds. Using the National Solar Radiation Database (NSRDB), this research evaluates ensemble cloud-mask predictions from the WRF-Solar ensemble prediction system (WRF-Solar EPS). From the WRF-Solar EPS, day-ahead solar forecasts for the contiguous United States (CONUS) for 2018 are simulated. Given the NSRDB data is accessible at a resolution of 2 km, we can calculate the cloud fraction across the 9-km grid of WRF-Solar EPS. This allows us to spatially assess the cloud-mask forecasts using two methods against the high-resolution NSRDB: (a) considering all 2-km NSRDB clouds in the forecast domain (EMAll), and (b) using a minimum cloud fraction threshold of 50% to designate a pixel as cloudy (EMP50). The low-resolution cloud masks from WRF-Solar EPS are evaluated directly against the cloud-resolving scale gridded observations from NSRDB using EMAll. With EMP50, we presume that scenes with less than 50% cloud cover from the 2-km NSRDB are clear. Thus, this assessment approach allows for a fair comparison with WRF-Solar EPS resolved to a 9-km grid. A method of point-by-point verification is used to evaluate dichotomous (yes/no) cloud mask predictions against the NSRDB. For each pixel of model extent, cloud frequency and traditional metrics (e.g., probability of detection, false alarm rate, and hit rate, etc.) are computed and compared with satellite-derived data sets. Mismatched cloud frequency (MCF) is computed to measure the present capability of WRF-Solar EPS in representing various types of clouds, which are categorized using three levels of cloud top height (CTH) and cloud optical depth (COD) across entire CONUS. Preliminary results show that the WRF-Solar EPS provides MCF values ranging from 9% to 46%, 16% to 33%, and 8% to 27% for low-level, middle-level, and high-level clouds, respectively, for three CTHs. The model produces MCFs ranging from 27% to 46%, 13% to 34%, and 8% to 19% for thin, medium-thickness, and thick clouds, respectively, for three CODs. The presentation will include a detailed description of the current outcomes as well as potential future extensions. The evaluation approach established in this study is readily extensible to the evaluation of cloud predictions from different ensemble NWP models. In addition, the findings of the suggested evaluation technique aid in identifying model weaknesses and will ultimately lead to advances in WRF-Solar EPS's skill in predicting clouds and solar irradiance.

cloud mask forecast↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unraveling electronic correlations in warm dense quantum plasmas

The study of matter at extreme densities and temperatures has emerged as a highly active frontier at the interface of plasma physics, material science and quantum chemistry with relevance for planetary modeling and inertial confinement fusion. A particular feature of such warm dense matter is the complex interplay of Coulomb interactions, quantum effects, and thermal excitations, making its rigorous theoretical description challenging. Here, we demonstrate how ab initio path integral Monte Carlo simulations allow us to unravel this intricate interplay for the example of strongly compressed beryllium, focusing on two X-ray Thomson scattering data sets obtained at the National Ignition Facility. We find excellent agreement between simulation and experiment with a very high level of consistency between independent observations without the need for any empirical input parameters. Our results call into question previously used chemical models, with important implications for the interpretation of scattering experiments and radiation hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Material Characterization for Large Scale Additive Manufacturing (AM)

The objective of this project was to enable prediction of residual stresses and warpage for extrusion deposition additive manufacturing (EDAM) fabricated carbon fiber thermoplastic matrix parts with three different materials systems produced by Techmer. The materials selected include Polyphenylene sulfide reinforced with 50% by weight of carbon fiber (CF-PPS), Polyethersulfone reinforced with 25% by weight of carbon fiber (CF-PESU), Polysulfone reinforced with 25% by weight of carbon fiber (CF-PSU). ADDITIVE3D © , a physics-based simulation workflow for EDAM, provided the simulation capabilities required for this project. Simulation predictions were validated against measurements carried out during and after the printing process carried out in the CAMRI and LSAM systems . Predictions for temperature, degree of crystallinity and deformation were carried out for the three material systems and two different geometries. Predictions for temperature were correlated very well with the experimental measurements for the two geometries printed using the three-material systems. The crystallinity level was verified for CF-PPS. The predictions for part deformation were in good agreement with the experimental measurements. In the best-case, predictions were within 8% of the maximum displacement observed in the 3-direction whereas predictions were within 14% for the worst-case. The adoption of this technology in commercial large-scale EDAM production processes would dramatically reduce the costs associated with producing articles by that process. Many thousands of dollars in materials, energy, and machine time could be saved by utilizing this simulation technology to develop articles rather than iterative printings to arrive at the optimal or correct design. The avoidance of a single failed print has the possibility of saving tens of thousands of dollars involved in the cost of material, machine and operators’ time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Short-range order and its impact on the CrCoNi medium-entropy alloy

Traditional metallic alloys are mixtures of elements in which the atoms of minority species tend to be distributed randomly if they are below their solubility limit, or to form secondary phases if they are above it. The concept of multiple-principal-element alloys has recently expanded this view, as these materials are single-phase solid solutions of generally equiatomic mixtures of metallic elements. This group of materials has received much interest owing to their enhanced mechanical properties 1,2,3,4,5 . They are usually called medium-entropy alloys in ternary systems and high-entropy alloys in quaternary or quinary systems, alluding to their high degree of configurational entropy. However, the question has remained as to how random these solid solutions actually are, with the influence of short-range order being suggested in computational simulations but not seen experimentally 6,7 . In this paper we report the observation, using energy-filtered transmission electron microscopy, of structural features attributable to short-range order in the CrCoNi medium-entropy alloy. Increasing amounts of such order give rise to both higher stacking-fault energy and hardness. These findings suggest that the degree of local ordering at the nanometre scale can be tailored through thermomechanical processing, providing a new avenue for tuning the mechanical properties of medium- and high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Brominated lipid probes expose structural asymmetries in constricted membranes

Lipids in biological membranes are thought to be functionally organized, but few experimental tools can probe nanoscale membrane structure. Using brominated lipids as contrast probes for cryo-EM and a model ESCRT-III membrane-remodeling system composed of human CHMP1B and IST1, we observed leaflet-level and protein-localized structural lipid patterns within highly constricted and thinned membrane nanotubes. These nanotubes differed markedly from protein-free, flat bilayers in leaflet thickness, lipid diffusion rates and lipid compositional and conformational asymmetries. Simulations and cryo-EM imaging of brominated stearoyl-docosahexanenoyl-phosphocholine showed how a pair of phenylalanine residues scored the outer leaflet with a helical hydrophobic defect where polyunsaturated docosahexaenoyl tails accumulated at the bilayer surface. Combining cryo-EM of halogenated lipids with molecular dynamics thus enables new characterizations of the composition and structure of membranes on molecular length scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗
DOE OSTI2020

Attosecond XUV probing of vibronic quantum superpositions in Br 2 +

Ultrafast laser excitation can create coherent superpositions of electronic states in molecules and trigger ultrafast flow of electron density on a few femtosecond timescale. While recent attosecond experiments have addressed real-time observation of these primary photochemical processes, the underlying roles of simultaneous nuclear motions and how they modify and disturb the valence electronic motion remain uncertain. In this work, we investigate coherent electronic-vibrational dynamics induced among multiple vibronic levels of ionic bromine (Br 2 + ), including both spin-orbit (X 2 Π 3/2,g – X 2 Π 1/2,g ) and valence (A 2 Π 3/2,u – 4 Σ 3/2,u ) electronic superpositions, using attosecond transient absorption spectroscopy. Decay, revival, and apparent frequency shifts of electronic coherences are measured via characteristic quantum beats on the Br-3d core-level absorption signals. Quantum-mechanical simulations attribute the observed electronic decoherence to broadened phase distributions of nuclear wave packets on anharmonic potentials. Molecular vibronic structure is further revealed to be imprinted as discrete progressions in electronic beat frequencies. These results provide a future basis to interpret complex charge-migration dynamics in polyatomic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamical Scaling as a Signature of Multiple Phase Competition in Yb 2⁢ Ti 2 ⁢O 7

Yb 2 ⁢Ti 2 ⁢O 7 is a celebrated example of a pyrochlore magnet with highly frustrated, anisotropic exchange interactions. To date, attention has largely focused on its unusual, static properties, many of which can be understood as coming from the competition between different types of magnetic order. Here, in this work, we use inelastic neutron scattering with exceptionally high energy resolution to explore the dynamical properties of Yb 2 ⁢Ti 2 ⁢O 7 . We find that spin correlations exhibit dynamical scaling, analogous to behavior found near to a quantum critical point. We show that the observed scaling collapse can be explained within a phenomenological theory of multiple-phase competition, and confirm that a scaling collapse is also seen in semiclassical simulations of a microscopic model of Yb 2 ⁢Ti 2 ⁢O 7 . These results suggest that dynamical scaling may be general to systems with competing ground states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Specific ion effects on ion transport in charged polymer membranes

Selective separation of like-charged ions is a central challenge in applications such as critical mineral recovery. Electrochemical membrane–based separations offer promising pathways to address this need, but a limited fundamental understanding of ion transport in charged polymer membranes hampers the development of highly selective materials. This study elucidates the role of specific ion effects (SIEs) in ion transport through charged polymer membranes by integrating experimental measurements of ion mobility and in situ ion/ion and ion/water interactions in model charged polymer membranes with results from molecular dynamics simulations on analogous systems. We demonstrate that solvent-mediated ion interactions drive pronounced SIEs within the studied membranes, with the ion softness, i.e., the malleability of ion hydration shells, from hard/soft acid/base (HSAB) theory emerging as a key predictor of transport properties. HSAB theory also explains the observed mechanism of solvent-mediated ion interactions. Our findings offer a mechanistic framework for designing membranes with tailored ion selectivity, potentially enabling efficient separations of chemically similar ions.

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

Model Coupling Through Learned Representations

Reliable climate predictions are important for making robust decisions in response to the changing climate. This project aims to reduce mis-modeling uncertainties arising from the representation of the land-atmosphere coupling in the Energy Exascale Earth System Model (E3SM) by using a machine learning approach. This approach will use an encoder-decoder architecture to represent the information that is developed in the land model and given to the atmosphere model. The simulated data will be taken from the E3SM simulation. However, the incorporation of observed data into the simulated dataset reduces mis-modeling uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pulse duration dependent effects of ultrafast laser induced damage on a 1030 nm multi-layer dielectric mirror for high repetition rate, high average power laser systems

High repetition rate, high peak, and average power laser systems are crucial for next-generation particle accelerators, inertial confinement fusion, and secondary particle sources. These applications demand durable laser optics, particularly interference coatings on optics lasting millions of shots at high fluence. This study focuses on designing, testing, and simulating multi-layer dielectric (MLD) mirrors for pulse durations of 260 fs, 77 fs, and 25 fs at 1030 nm wavelength and 45-degree incidence angle with p -polarization. S-on-1 laser-induced damage thresholds (LIDT) for varying pulse numbers were determined, with single-shot LIDT values of 0.98 Jcm -2 , 1.63 Jcm -2 , and 2.3 Jcm -2 for 25 fs, 77 fs, and 260 fs respectively. A strong correlation between blister shape and local fluence was observed, implying that the layer expansion in a blister depends on local fluence. We have also examined mechanisms responsible for laser-induced stress generation and energy release rates in blister formation. Damage mechanisms are further explored by finite-difference time-domain (FDTD) simulations, incorporating Keldysh strong field ionization, whose predictions were in excellent agreement with the onset of damage determined experimentally. These findings offer insights for enhancing MLD coating technology, promising more efficient and resilient laser systems for diverse scientific and industrial applications.

Noor, Mohamed Yaseen (ORCID:000000021036644X)↗