Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded systems”

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 487 records · Page 27

BETO 2021 Peer Review - Inverse Bioproduct Design Through Machine Learning and Molecular Simulation

This work aims to identify performance advantaged bioproducts (PABPs) through property prediction, which will guide experimental synthesis. The impact of this work will be faster market adoption of bioproducts with greater performance relative to incumbent products. We have identified >106 bioproduct candidates, but only some will have superior performance to create a market pull. High-throughput property prediction, enabled by machine learning, and elucidation of structure-function relationships, enabled by molecular simulation, provide a hypothesis driven approach for down selection of candidate biomolecules to pursue experimentally. To enable machine learning and molecular simulation for bioproduct discovery, automated structure generation and embedding must capture relevant features for prediction, databases must cover domains applicable to biobased products, and best practices for simulation of polymer systems must be developed. To address these challenges, we have established bioproduct relevant datasets, developed high-throughput polymer structure generation, and built end-to-end neural networks that have predicted 8 properties for >1.4 x 106 biopolymers. A molecular simulation pipeline for building, running, and analyzing polymers and polymer additives is being used to predict performance and develop design principles of biobased products. In collaboration with the PABP synthesis project, these computational tools are guiding synthesis and informing design of PABPs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A Chemist’s Guide to Solid-State Qubits: Diamond Color Centers as Embedded Vacancy-Centered Molecules

Diamond point defects or color centers are at the crux of solid-state quantum technology applications, from computing to sensing. Modeling molecular qubits has advanced by leaps and bounds mainly due to the advanced quantum mechanical (QM) tools in the arsenal of quantum chemists for modeling (isolated) molecular systems, which should be easily translatable for analogous systems in the solid state. Here, we simulate a series of diamond defects: a carbon vacancy (V C ) and second row substitutions that are adjacent to a V C or D C V C , with D = B, N, and O, to evaluate their magnetic and optical properties. We applied molecular orbital principles, an embedded cluster method, and a multiconfigurational QM theory to computationally characterize the above-mentioned defects for a multitude of spin, charge, and excited states. Our computational approach delivered high-quality QM insights that compare well with experiments. Furthermore, with this framework, we confirm the spin-triplets B C V C – and O C V C 0 to be analogous to N C V C – , with calculated bright transitions that respectively correspond to UV (3.80 eV) and blue (2.65 eV) light absorption and emission lifetimes (τ) of 2.3 and 5.6 ns, which are higher in energy and more short-lived (and thus brighter) than that for N C V C – : green (2.26 eV), τ = 9.5 ns.

Martirez, John Mark P. [Princeton Plasma Physics L↗

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupling SCALE with DAKOTA for Axial Burnup Profiles Assessment in Burnup Credit

This paper presents a computational study that demonstrates the application of the SCALE code system in conjunction with the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) for the analysis of key factors influencing the evaluation of burnup credit (BUC) in pressurized water reactors (PWRs). The primary objective of this analysis is to characterize the model by utilizing parameterization, uncertainty quantification, and optimization studies. Using this approach, we can comprehensively assess the system and conduct informed predictive studies. This study highlights the effectiveness of the SCALE code system integrated within the DAKOTA framework in terms of efficiency and capability. With the coupling of the burnup code ORIGAMI with the CSAS or TSUNAMI-3D sequence embedded in a DAKOTA analysis, we can characterize the factors that influence the k eff of PWR 17x17 spent nuclear fuel (SNF) in the GBC-32 computational benchmark cask for the assessment of BUC in criticality safety analysis. The coupling methodology used in this study is not exclusive to BUC analysis. However, the choice to apply this methodology to the BUC problem is particularly significant because of the diverse range of aspects it encompasses in nuclear criticality safety analyses. This problem presents a unique opportunity to explore and address multiple facets of such analyses related to BUC and illustrates the capability of the SCALE code system with DAKOTA. This analysis makes use of historical reference data for the axial burnup profile, where the entire space within the bounds is considered. Both SCALE and DAKOTA are currently integrated in the Nuclear Energy Advanced Modeling Simulation (NEAMS) Workbench code system, which has a user-friendly graphical interface that simplifies the setup of simulations and configuration of input parameters as well as the visualization of simulation results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗

Double photoionization of atomic oxygen: Feshbach resonances in the two-electron continuum

In this work, we describe a joint experimental and theoretical investigation on oxygen double photoionization—the emission of two electrons from atomic oxygen following single photon absorption. High-resolution experimental measurements were performed at the Advanced Light Source, revealing sharp resonance structure superimposed on the more familiar Wannier-like, nearly-linear background. These resonance features are attributed to ionization-plus-excitation Feshbach resonances embedded in the double ionization continuum, doubly-excited states that lie above the double-ionization threshold. Such features are absent in the double photoionization cross section of He, or other quasi-two-electron systems, for which the doubly-ionized atomic core remains inert. For a corresponding theoretical analysis, the R-matrix with pseudostates (RMPS) method was invoked by calculating final-state, two-electron resonances-plus-continua wavefunctions and corresponding single-photon absorption cross sections. Overall agreement is found in the direct, background double photoionization cross section. However, the RMPS method, using a small basis due to practical computational limitations, was unable to reproduce quantitatively the smooth background or the sharper resonance features observed in the measurements, showing instead large-scale oscillations about the experimental background, and characteristic pseudoresonance jitter, associated with an insufficient convergence of the pseudostate representation to the true two-electron infinite series of Feshbach resonances embedded in the two-electron continuum. The prominent resonance structure observed highlights the need to consider multiple excitation processes in atoms more complex than He or quasi-two-electron systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Covalent Drug Binding in Live Cells Monitored by Mid-Infrared Quantum Cascade Laser Spectroscopy: Photoactive Yellow Protein as a Model System

The detection of drug-target interactions in live cells enables analysis of therapeutic compounds in a native cellular environment. Recent advances in spectroscopy and molecular biology have facilitated the development of genetically encoded vibrational probes like nitriles that can sensitively report on molecular interactions. Nitriles are powerful tools for measuring electrostatic environments within condensed media like proteins, but such measurements in live cells have been hindered by low signal-to-noise ratios. In this study, we design a spectrometer based on a double-beam quantum cascade laser (QCL)-based transmission infrared (IR) source with balanced detection that can significantly enhance sensitivity to nitrile vibrational probes embedded in proteins within cells compared to a conventional FTIR spectrometer. Here, using this approach, we detect small-molecule binding in Escherichia coli, with particular focus on the interaction between para-Coumaric acid (pCA) and nitrile-incorporated photoactive yellow protein (PYP). This system effectively serves as a model for investigating covalent drug binding in a cellular environment. Notably, we observe large spectral shifts of up to 15 cm –1 for nitriles embedded in PYP between the unbound and drug-bound states directly within bacteria, in agreement with observations for purified proteins. Such large spectral shifts are ascribed to the changes in the hydrogen-bonding environment around the local environment of nitriles, accurately modeled through high-level molecular dynamics simulations using the AMOEBA force field. Our findings underscore the QCL spectrometer’s ability to enhance sensitivity for monitoring drug–protein interactions, offering new opportunities for advanced methodologies in drug development and biochemical research.

chromophores↗

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

Granular metals with SiN x dielectrics

Understanding and controlling nanoscale interface phenomena, such as band bending and secondary phase formation, is crucial for electronic device optimization. In granular metal (GM) studies, where metal nanoparticles are embedded in an insulating matrix, the importance of interface phenomena is frequently neglected. Here, we demonstrate that GMs can serve as an exemplar system for evaluating the role of secondary phases at interfaces through a combination of x-ray photoemission spectroscopy (XPS) and electrical transport studies. We investigated SiN x as an alternative to more commonly used oxide-insulators, as SiN x -based GMs may enable high temperature applications when paired with refractory metals. Comparing Co-SiN x and Mo-SiN x GMs, we found that, in the tunneling-dominated insulating regime, Mo-SiN x had reduced metal-silicide formation and orders-of-magnitude lower conductivity. XPS measurements indicate that metal-silicide and metal-nitride formation are mitigatable concerns in Mo-SiN x . Given the metal-oxide formation seen in other GMs, SiN x is an appealing alternative for metals that readily oxidize. Furthermore, SiN x provides a path to metal-nitride nanostructures, potentially useful for various applications in plasmonics, optics, and sensing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Light detection with power and signal transmission over fiber

The Deep Underground Neutrino Experiment is a next generation long baseline (1300 km) neutrino oscillation experiment. The neutrino beam measurements will be performed by a near detector and far detector. The far detector will consist of four modules, installed 1.5 km deep underground, based on Liquid Argon Time Projection Chamber (LArTPC) technology to detect particles. The Vertical Drift (VD) LArTPC is a recent technology proposed by the DUNE collaboration for the second FD module. In VD, light collection will be optimized by embedding photon detectors within the LArTPC cathode, which is biased at −300 kV. As a result, power must arrive to the Photon Detection System and signal must be transmitted via non-conductive material. The proposed solution is to use Power-over-Fiber and Signal-over-Fiber. An intense validation of the system is being performed by the collaboration at the CERN Neutrino Platform. The design of the system and the results of the validation collected over the first half of 2022 are presented here.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dynamically controlled random lasing with colloidal titanium carbide MXene

Control of lasing properties through tailorable and dynamically tunable materials and reconfigurable compositions can augment the performance of random lasers for a wide range of applications. Here, a colloid of randomly dispersed weakly scattering single-layer titanium carbide (Ti 3 C 2 T x ) MXene flakes embedded within rhodamine 101 gain medium is experimentally shown to provide feedback for random lasing. Additionally, in contrast to previously reported random laser systems where the optical properties of scatterers are static, the relative permittivity of Ti 3 C 2 T x MXene flakes can be varied under optical pumping due to the saturable absorption properties. Numerical simulations indicate that the observed nonlinear response of Ti 3 C 2 T x MXene flakes enables dynamically tunable random lasing. Thus, pumping the Ti 3 C 2 T x MXene flakes with a second optical source decreases the gain threshold required to obtain random lasing. Also, using numerical simulations, it is shown that the control over the intensity of the second pump enables tuning the field distribution of the random lasing modes. Considering the diversity of the MXenes family, the proposed MXene colloidal metamaterial design opens up a new avenue to advanced control of lasing properties for photonic applications.

36 MATERIALS SCIENCE↗

Validation of a lightning protection system for a fusion-welded thermoplastic composite wind turbine blade tip

Fusion welding for thermoplastic composite wind turbine blades is being investigated as a replacement for adhesives in bond lines to reduce manufacturing time and cost and to promote ease of recyclability. Resistance welding uses an electrically conductive material, such as carbon fiber, embedded in the bond line to heat the joining surfaces and fuse them together under pressure. One main challenge with fusion joining for large wind blades is the ability to use this conductive heating element in the bond line while still protecting the blade from the catastrophic effects of a lightning strike. This paper outlines the design and manufacturing of a lightning protection system for a 5-m-long thermally welded blade tip. The outcomes of high voltage and high current lightning testing applied to the 5 m blade tip show that the proposed design was able to protect the fusion welded bond lines from significant lightning damage.

17 WIND ENERGY↗

Revisiting Power Systems Time-Domain Simulation Methods and Models

The changing nature of power systems dynamics is challenging present practices related to modeling and study of system-level dynamic behavior. While developing new techniques and models to handle the new modeling requirements, it is also critical to review some of the terminology used to describe existing simulation approaches and the embedded assumptions. This article provides a first-principles review of the simplifications and transformations commonly used in the formulation of time-domain simulation models. It introduces a taxonomy and classification of time-domain simulation models depending on their frequency bandwidth, network representation, and software availability. Furthermore, it focuses on the fundamental aspects of averaging techniques, and model reduction approaches that result in modeling choices, and discusses the associated challenges and opportunities of applying these methods in systems with large shares of Inverter Based Resources (IBRs). The article concludes with an illustrative simulation that compares the trajectories of an IBR-dominated system.

behavioral sciences↗

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗

In situ detection and characterization of alkali-silica reaction damage in concrete using contactless ultrasonic wavefield imaging

Here we present work to characterize cracking damage in concrete caused by alkali-silica reactivity distress. We apply contactless ultrasonic scan inspection, exploiting Rayleigh wave scattering from concrete inhomogeneities. The scanning system utilizes a multi-channel MEMS sensor array to collect ultrasonic wavefield data from concrete subjected to ASR-promoting environments. The location and extent of ASR damage are established by internal expansion measurements from embedded strain gauges, scanning electron microscope images of cored samples, and expansion measurements from companion samples. A wavefield data processing method extracts oscillatory fields from the ultrasonic wavefield data to detect distinct zones of ASR damage. A damage index is proposed to further characterize the extent of ASR damage. The results confirm the feasibility and accuracy of the approach to characterize ASR damage in concrete. The fully contactless ultrasonic scanning measurement system does not require separate material samples and enables in situ characterization of ASR damage within concrete structures.

36 MATERIALS SCIENCE↗

Netload Range Cost Curves for Coordinated Transmission-Distribution Planning Under DER Growth Uncertainty

The increasing penetration of distributed energy resources (DERs) requires better coordination between transmission and distribution (T&D) planning to ensure system security and cost efficiency. However, misaligned planning horizons, computational burdens, and privacy concerns hinder effective coordination, leading to either underutilized resources caused by overinvestments or reliability risks due to underinvestment. To address this challenge, we introduce netload range cost curves (NRCCs), a novel approach for managing long-term DER growth uncertainty through T&D coordination, while preserving existing data-sharing and regulatory structures. NRCCs provide pairs of (i) peak substation netload guarantees and (ii) corresponding distribution upgrade options and costs, enabling their seamless integration into transmission planning workflows. To compute NRCCs efficiently, we develop a transmission-aware distribution network planning (TADNP), which is subsequently integrated to an iterative computation procedure. These NRCCs are then embedded into an NRCC-informed transmission planning model to enable resource-efficient coordination. We illustrate our proposed approach with a case study based on realistic distribution and transmission systems in the San Francisco Bay Area, California. Our results indicate the possibility of dramatic savings in transmission investments by incorporating the proposed NRCC-integrated T&D coordination framework.

Li, Yujia↗