Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static analysis”

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

Heterogeneous Polymer Dynamics Explored Using Static 1H NMR Spectra

NMR spectroscopy continues to provide important molecular level details of dynamics in different polymer materials, ranging from rubbers to highly crosslinked composites. It has been argued that thermoset polymers containing dynamic and chemical heterogeneities can be fully cured at temperatures well below the final glass transition temperature (Tg). In this paper, we described the use of static solid-state 1H NMR spectroscopy to measure the activation of different chain dynamics as a function of temperature. Near Tg, increasing polymer segmental chain fluctuations lead to dynamic averaging of the local homonuclear proton-proton (1H-1H) dipolar couplings, as reflected in the reduction of the NMR line shape second moment (M2) when motions are faster than the magnitude of the dipolar coupling. In general, for polymer systems, distributions in the dynamic correlation times are commonly expected. To help identify the limitations and pitfalls of M2 analyses, the impact of activation energy or, equivalently, correlation time distributions, on the analysis of 1H NMR M2 temperature variations is explored. It is shown by using normalized reference curves that the distributions in dynamic activation energies can be measured from the M2 temperature behavior. An example of the M2 analysis for a series of thermosetting polymers with systematically varied dynamic heterogeneity is presented and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Implementation and Verification of Cable Bending Stiffness in MoorDyn: Preprint

Cable bending stiffness capabilities have been added to MoorDyn to enable modeling of dynamic power cables. The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. MoorDyn is a lumped-mass mooring dynamics model that is used for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim). Its existing line elasticity formulation considers axial stiffness only. To properly capture the dynamics of power cables, a cable bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent line elements. The resulting bending moment is realized by applying forces on adjacent nodes, leaving the underlying lumped-mass formulation unchanged. The bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2025

Cost-optimized energy storage operation for a grid-connected solar PV system at community and individual scales

This study provides a comparative analysis of grid-connected PV-integrated battery storage at individual and community scales. The paper addresses the challenge of managing energy demand-generation mismatch by using a battery energy storage optimization algorithm, which minimizes operational costs while accounting for battery degradation. Also, this work introduces a broader evaluation basis that includes seasonal variability, grid exchange smoothness, and scalability across different battery capacities. Results show that community-scale storage more effectively dampens grid exchange power fluctuations and reduces system costs, particularly with moderate price differences between electricity buying and selling prices and low battery capacities. The paper also analyzes the impacts of static control versus cost-optimized battery system management. Here, it is shown that the gap in system costs between the cost-optimized and static control scenarios widens as the price difference increases.

25 ENERGY STORAGE↗
DOE OSTI2025

Multiscale Physics of Atomic Nuclei from First Principles

Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding are approached near the neutron and proton drip lines, traditional shell structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe this multiscale physics starting from nuclear forces with roots in quantum chromodynamics. Here, we achieve this within a unified and nonperturbative quantum many-body framework that captures both short- and long-range correlations starting from modern nucleon-nucleon and three-nucleon forces from chiral effective field theory. The short-range (dynamic) correlations which account for the bulk of the binding energy are included within a symmetry-breaking framework, while long-range (static) correlations (and fine details about the collective structure) are included by employing symmetry projection techniques. Our calculations accurately reproduce—within theoretical error bars—available experimental data for low-lying collective states and the electromagnetic quadrupole transitions in 20−30 Ne. In addition, we reveal coexisting spherical and deformed shapes in 30 Ne, which indicates the breakdown of the magic neutron number 𝑁 = 20 as the key nucleus 28 O is approached, and we predict that the drip line nuclei 32,34 Ne are strongly deformed and collective. By developing reduced-order models for symmetry-projected states, we perform a global sensitivity analysis and find that the subleading singlet 𝑆-wave contact and a pion-nucleon coupling strongly impact nuclear deformation in chiral effective field theory. The techniques developed in this work clarify how microscopic nuclear forces generate the multiscale physics of nuclei spanning collective phenomena as well as short-range correlations and allow one to capture emergent and dynamical phenomena in finite fermion systems such as atom clusters, molecules, and atomic nuclei.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Design Report on the MiniFuel Instrumented Test Apparatus for Understanding Radiation Effects

Most nuclear fuels irradiations at Oak Ridge National Laboratory (ORNL) over the past decade have been conducted using MiniFuel—a static capsule design employing subscale fuel specimens to collect separate-effects irradiated fuel performance data. Irradiation conditions for MiniFuel experiments are predicted pre-test using reactor physics, and thermal models are verified post-irradiation via SiC dilatometry and various spectrometry methods. Relevant fuel performance parameters are also observed post-irradiation in a hot cell, thereby providing a single data point for each parameter representing the cumulative effects of the irradiation conditions. Substantially more data can be harvested from a single test and within a shorter duration by instrumenting irradiation vehicles and measuring desired quantities in situ. This report presents the design and analysis of the MiniFuel INstrumented Irradiation Test Apparatus for Understanding Radiation Effects (MINITAURE)—an instrumented test rig based on the separate-effects MiniFuel concept that aims to capture fission gas release (FGR) and thermal conductivity degradation of fuel specimens during irradiation in the High Flux Isotope Reactor (HFIR). MINITAURE will be integrated with the Materials Irradiation Facility (MIF) located in the HFIR building outside the reactor containment. The MIF will act as the instrumentation and control center for the experiment, enabling real-time feedback from in situ sensors and control of irradiation temperatures via a gas delivery system. Two unique capsule designs were developed to capture each phenomenon: the thermal conductivity capsule, which uses a thermopile method to estimate fuel specimen thermal conductivity, and the fission gas release capsule, which will have continuous flowing gas communication to high-purity germanium detectors that are housed in the MIF for monitoring FGR. This report details the reactor physics and heat transfer modeling activities that were used to inform the experiment design and predict capsule performance. It also describes out-of-pile activities conducted to stand up this new capability and verify the measurement techniques. Modeling efforts to date have demonstrated the feasibility of the in situ measurement techniques and supported the development of the MINITAURE assembly configuration. Out-of-pile testing of the thermal conductivity measurement shows promise in capturing relative changes in thermal conductivity. However, significant errors exist in the measured absolute value, posing a need for further refinement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

ParFlow-CLM model simulation files for Maina et al., Journal of Hydrology, 2022

This dataset contains the files to run a ParFlow-CLM integrated hydrologic model simulation for Maina et al., HESS, 2022. It also contains the associated daily pressure output of those simulations. Simulations are for a hillslope-similarity approach based on seasonal groundwater changes in the East River, Colorado. Comparisons are made to different clustering, or functional zonation approaches. We assess the ability of these clustering approaches to identify and categorize hillslopes with similar static characteristics, hydroclimate, land surface processes, and subsurface dynamics in a mountainous watershed. ParFlow binary files (.pfb) are model specific, and can be modified/viewed with conversion to .sa, .silo, or .vtk file formats. ParFlow input scripts are in the .tcl file format. Associated post-processing scripts are included as .f90 files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Spin wave Hamiltonian and anomalous scattering in NiPS 3

Here, we report a comprehensive spin wave analysis of the semiconducting honeycomb van der Waal antiferromagnet NiPS 3 . Using single-crystal inelastic neutron scattering, we map out the full Brillouin zone and fit the observed modes to a spin wave model with rigorously defined uncertainty. We find that the third-nearest-neighbor exchange J 3 dominates the Hamiltonian, a feature which we fully account for by ab initio density functional theory calculations. We also quantify the degree to which the threefold rotation symmetry is broken and account for the Q = 0 excitations observed in other measurements, yielding a spin exchange model which is consistent across multiple experimental probes. We also identify a strongly reduced static ordered moment and reduced low-energy intensity relative to the linear spin wave calculations, signaling unexplained features in the magnetism which requires going beyond the linear spin wave approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Analysis and Characterization of Four-quadrant Switches based Commutation Cell

A four-quadrant switch (FQS) blocks either polarity voltage and controls current flow in both directions. Unlike voltage-source converters, in which two-quadrant switches operate over a narrow voltage range, four-quadrant switches are required to operate over a wide range of both voltage and current in applications such as matrix converters and current-source converters. Furthermore, matrix converters require multi-step commutation schemes compared to two-step schemes for current-bidirectional switch based voltage-source converters and voltage-bidirectional switch based current-source converters. This paper provides a generalized overview of commutation schemes used for two and four quadrant switches based two-level commutation cells, identifies comparison indices for FQS commutation schemes, and discusses the need for adaptive commutation-step times for wide voltage and current variation applications. Also, the static and dynamic characteristics of 1.2 kV rated FQS implementations utilizing commercial SiC MOSFETs from four different manufacturers and novel monolithic SiC BiDirectional Field Effect Transistor (BiDFET) have been reported.

Commutation cell, four quadrant switch, SiC Bidire↗
DOE OSTI2022

Modeling Heterogeneity in UO 2 Nanoparticles Using X‐ray Absorption Spectroscopy

Abstract EXAFS provides the capability to interrogate nanoparticle (NP) structure in atomistic detail without relying on long‐range crystallinity. There is a limitation in that EXAFS provides averaged structural information, making it difficult to separate a small amount of heterogeneous structure from bulk. In this work, models were developed to extract surface‐specific information from conventional EXAFS measurements collected on UO 2 NPs of varying size. Specifically, the surface terminating species of UO 2 NPs was determined from comparison of coordination numbers with geometric models while the origin of static disorder was interrogated from user‐defined simulations. Results show that the degree of oxygenation on the NP surface does not significantly deviate from bulk surface and that static disorder is highly enhanced in NP surface layers but cannot be attributed to surface relaxation effects alone. The approach described herein has the potential to be adapted to a range of inorganic NP systems to interrogate surface structure.

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

Mechanism and dynamics of fatty acid photodecarboxylase

Photoenzymes are rare biocatalysts driven by absorption of a photon at each catalytic cycle; they inspire development of artificial photoenzymes with valuable activities. Fatty acid photodecarboxylase (FAP) is a natural photoenzyme that has potential applications in the bio-based production of hydrocarbons, yet its mechanism is far from fully understood. RATIONALE To elucidate the mechanism of FAP, we studied the wild-type (WT) enzyme from Chlorella variabilis (CvFAP) and variants with altered active-site residues using a wealth of techniques, including static and time-resolved crystallography and spectroscopy, as well as biochemical and computational approaches. RESULTS A 1.8-Å-resolution CvFAP x-ray crystal structure revealed a dense hydrogen-bonding network positioning the fatty acid carboxyl group in the vicinity of the flavin adenine dinucleotide (FAD) cofactor. Structures solved from free electron laser and low-dose synchrotron x-ray crystal data further highlighted an unusual bent shape of the oxidized flavin chromophore, and showed that the bending angle (14°) did not change upon photon absorption (step 1) or throughout the photocycle. Calculations showed that bending substantially affected the energy levels of the flavin. Structural and spectroscopic analysis of WT and mutant proteins targeting two conserved active-site residues, R451 and C432, demonstrated that both residues were crucial for proper positioning of the substrate and water molecules and for oxidation of the fatty acid carboxylate by 1 FAD* (~300 ps in WT FAP) to form FAD ∙– (step 2). Time-resolved infrared spectroscopy demonstrated that decarboxylation occured quasi-instantaneously upon this forward electron transfer, consistent with barrierless bond cleavage predicted by quantum chemistry calculations and with snapshots obtained by time-resolved crystallography. Transient absorption spectroscopy in H 2 O and D 2 O buffers indicated that back electron transfer from FAD ∙– was coupled to and limited by transfer of an exchangeable proton or hydrogen atom (step 3). Unexpectedly, concomitant with FAD ∙– reoxidation (to a red-shifted form FAD RS ) in 100 ns, most of the CO 2 product was converted, most likely into bicarbonate (as inferred from FTIR spectra of the cryotrapped FAD RS intermediate). Calculations indicated that this catalytic transformation involved an active-site water molecule. Cryo-Fourier transform infrared spectroscopy studies suggested that bicarbonate formation (step 4) was preceded by deprotonation of an arginine residue (step 3). At room temperature, the remaining CO 2 left the protein in 1.5 μs (step 4'). The observation of residual electron density close to C432 in electron density maps derived from time-resolved and cryocrystallography data suggests that this residue may play a role in stabilizing CO 2 and/or bicarbonate. Three routes for alkane formation were identified by quantum chemistry calculations; the one shown in the figure is favored by the ensemble of experimental data. CONCLUSION Finally, we provide a detailed and comprehensive characterization of light-driven hydrocarbon formation by FAP, which uses a remarkably complex mechanism including unique catalytic steps. We anticipate that our results will help to expand the green chemistry toolkit.

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

Effect of Topology on Transient Dynamic and Shock Response of Polymeric Lattice Structures

Architected cellular materials, such as lattice structures, offer potential for tunable mechanical properties for dynamic applications of energy absorption and impact mitigation. In this work, the static and dynamic behavior of polymeric lattice structures was investigated through experiments on octet-truss, Kelvin, and cubic topologies with relative densities around 8%. Here, dynamic testing was conducted via direct impact experiments (25–70 m/s) with high-speed imaging coupled with digital image correlation and a polycarbonate Hopkinson pressure bar. Mechanical properties such as elastic wave speed, deformation modes, failure properties, particle velocities, and stress histories were extracted from experimental results. At low impact velocities, a transient dynamic response was observed which was composed of a compaction front initiating at the impact surface and additional deformation bands whose characteristics matched low strain-rate behavior. For higher impact velocities, shock analysis was carried out using compaction wave velocity and Eulerian Rankine–Hugoniot jump conditions with parameters determined from full-field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ConnectIt: a framework for static and incremental parallel graph connectivity algorithms

Connected components is a fundamental kernel in graph applications. The fastest existing multicore algorithms for solving graph connectivity are based on some form of edge sampling and/or linking and compressing trees. However, many combinations of these design choices have been left unexplored. In this paper, we design the ConnectIt framework, which provides different sampling strategies as well as various tree linking and compression schemes. ConnectIt enables us to obtain several hundred new variants of connectivity algorithms, most of which extend to computing spanning forest. In addition to static graphs, we also extend ConnectIt to support mixes of insertions and connectivity queries in the concurrent setting. We present an experimental evaluation of ConnectIt on a 72-core machine, which we believe is the most comprehensive evaluation of parallel connectivity algorithms to date. Compared to a collection of state-of-the-art static multicore algorithms, we obtain an average speedup of 12.4x (2.36x average speedup over the fastest existing implementation for each graph). Using ConnectIt, we are able to compute connectivity on the largest publicly-available graph (with over 3.5 billion vertices and 128 billion edges) in under 10 seconds using a 72-core machine, providing a 3.1x speedup over the fastest existing connectivity result for this graph, in any computational setting. For our incremental algorithms, we show that our algorithms can ingest graph updates at up to several billion edges per second. To guide the user in selecting the best variants in ConnectIt for different situations, we provide a detailed analysis of the different strategies. Finally, we show how the techniques in ConnectIt can be used to speed up two important graph applications: approximate minimum spanning forest and SCAN clustering.

Computer Science↗
DOE OSTI2025

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Dynamically tunable membrane metasurfaces for infrared spectroscopy and strong light-matter interactions

Mid-infrared spectroscopy enables biochemical sensing by identifying vibrational molecular fingerprints, but it faces limitations in instrumentation portability and analytical sensitivity. Optical metasurfaces with strong mid-infrared photonic resonances provide an attractive solution towards on-chip spectrometry and sensitive molecular detection, yet their static nature hinders their anticipated impact. Here, we introduce and demonstrate dynamically tunable silicon membrane metasurfaces exhibiting high-Q transmissive resonances in the fingerprint region. By harnessing silicon’s thermo-optical properties, we achieve continuous modulation of coupling-induced transparency (CIT) modes that emerge upon the interference of quasi-bound states in the continuum (q-BICs) and surface lattice modes (SLMs). We measure a spectral tuning rate of 0.06 cm −1 K −1 by continuously sweeping the sharp CIT resonances over a 23.5 cm −1 spectral range across a temperature range of 300–700 K. In the current proof‑of‑concept implementation, the dynamic transmission control enables non-contact chemical analysis of polymer films by detecting characteristic absorption bands of polystyrene (1450 and 1492 cm −1 ) and poly(methyl methacrylate) (1730 cm −1 ) without requiring conventional spectrometers. When analyte molecules fill the metasurface-generated photonic cavities, we demonstrate vibrational strong coupling between the poly(methyl methacrylate)’s carbonyl band and the CIT mode, manifested in a Rabi splitting of ~43 cm −1 . Our results establish a new photonic platform that unites spectral precision, strong field enhancement, and reconfigurability, offering diverse potential for compact mid-infrared spectroscopy, molecular sensing, and programmable polaritonic photonics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Semiglobal Safety-Filtered Extremum Seeking With Unknown CBFs

We introduce a safe extremum-seeking (Safe ES) algorithm which achieves the minimization of an unknown objective function while ensuring that an unknown, yet measured, control barrier function (CBF) remains above an arbitrarily small negative value for all time. In other words, “practical safety” is maintained during the entire period of convergence to the constrained extremum. Our design is based on quadratic program (QP) CBF style filters for safety, which is applied in an average and estimated sense. Using nonsmooth analysis tools, we guarantee semiglobal practical asymptotic (SPA) stability of the global constrained optimum, practical convergence to the safe set if starting in a condition violating the CBF, and practical safety for all time—semiglobally—if starting in safe set. The safety result of the paper is analogous with modern notions of SPA stability, guaranteeing that, for any small violation of safety, there exist design coefficients which guarantee that such a small violation is not exceeded. The paper outlines a set of sufficient conditions on the barrier and objective functions, and by way of a Lyapunov argument, we demonstrate that nonconvex constrained optimization problems can be solved. We present these results in the setting of a static map and a dynamical system. A simulation example illustrates the results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Thermally induced cracking on the massive concrete structure of the NSLS II synchrotron and its engineering remediation

Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.

42 ENGINEERING↗
DOE OSTI2021

Thermo-Oxidation Analysis of Structural Adhesives: An Experimental Study

Abstract Cross-linked polymers play an important role in structural applications as adhesives. In this article, we performed an experimental study to understand the effects of high temperature i.e. thermo-oxidative aging on the mechanical behavior of polymers. Subsequently, the constitutive behavior, i.e. stress-strain relation of these adhesives before and after the accelerated aging tests were studied on samples aged at different temperature and different time durations at zero humidity. The aging conditions and periods selected for the tests have been specifically designed to simulate the effects of aging under the working conditions that the materials will be subjected to in their standard applications. The focus of this study is on the investigation of the effects of aging on the polymeric adhesives used in automotive applications. Consequently, four different polymeric adhesives (one silicone-based, two polyurethane-based, and one acrylic type), which are being used extensively in these applications, have been chosen. The quasi-static tensile loading was used on these four adhesives and the changes in the constitutive behavior were measured. The failure tensile tests as well as the cyclic tensile tests were performed in our lab using a Test-Resources tensile testing machine. The temperatures that were used in this work are 60, 80 and 95°C for varying aging periods of 24hrs, 240hrs, and 720hrs. Then, all the data was collected and consequently analyzed, and the resulting behavior was compared to that of the unaged material, showing varying results depending on the material, temperature, and aging period. Different behaviors have been observed after aging. The polymers’ behavior was expanding from no effect to significant damage caused by aging conditions. This work provides the core understanding required for modeling the material behavior during aging and lays down the groundwork for further data-infused model-free approaches.

Alazhary, Sharif↗