Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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 1,225 records · Page 68

Computational identification of HCV neutralizing antibodies with a common HCDR3 disulfide bond motif in the antibody repertoires of infected individuals

Despite recent success in hepatitis C virus (HCV) treatment using antivirals, an HCV vaccine is still needed to prevent reinfections in treated patients, to avert the emergence of drug-resistant strains, and to provide protection for people with no access to the antiviral therapeutics. The early production of broadly neutralizing antibodies (bNAbs) associates with HCV clearance. Several potent bNAbs bind a conserved HCV glycoprotein E2 epitope using an unusual heavy chain complementarity determining region 3 (HCDR3) containing an intra-loop disulfide bond. Isolation of additional structurally-homologous bNAbs would facilitate the recognition of key determinants of such bNAbs and guide rational vaccine design. Here we report the identification of new antibodies containing an HCDR3 disulfide bond motif using computational screening with the Rosetta software. Using the newly-discovered and already-known members of this antibody family, we review the required HCDR3 amino acid composition and propose determinants for the bent versus straight HCDR3 loop conformation observed in these antibodies.

59 BASIC BIOLOGICAL SCIENCES↗

Quantum fluctuations lead to glassy electron dynamics in the good metal regime of electron doped KTaO 3

One of the central challenges in condensed matter physics is to comprehend systems that have strong disorder and strong interactions. In the strongly localized regime, their subtle competition leads to glassy electron dynamics which ceases to exist well before the insulator-to-metal transition is approached as a function of doping. Here, we report on the discovery of glassy electron dynamics deep inside the good metal regime of an electron-doped quantum paraelectric system: KTaO 3 . We reveal that upon excitation of electrons from defect states to the conduction band, the excess injected carriers in the conduction band relax in a stretched exponential manner with a large relaxation time, and the system evinces simple aging phenomena—a telltale sign of glassy dynamics. Most significantly, we observe a critical slowing down of carrier dynamics below 35 K, concomitant with the onset of quantum paraelectricity in the undoped KTaO 3 . Our combined investigation using second harmonic generation technique, density functional theory and phenomenological modeling demonstrates quantum fluctuation-stabilized soft polar modes as the impetus for the glassy behavior. This study addresses one of the most fundamental questions regarding the potential promotion of glassiness by quantum fluctuations and opens a route for exploring glassy dynamics of electrons in a well-delocalized regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetic field mapping of inaccessible regions using physics-informed neural networks

A difficult problem concerns the determination of magnetic field components within an experimentally inaccessible region when direct field measurements are not feasible. In this paper, we propose a new method of accessing magnetic field components using non-disruptive magnetic field measurements on a surface enclosing the experimental region. Magnetic field components in the experimental region are predicted by solving a set of partial differential equations (Ampere’s law and Gauss’ law for magnetism) numerically with the aid of physics-informed neural networks (PINNs). Prediction errors due to noisy magnetic field measurements and small number of magnetic field measurements are regularized by the physics information term in the loss function. We benchmark our model by comparing it with an older method. The new method we present will be of broad interest to experiments requiring precise determination of magnetic field components, such as searches for the neutron electric dipole moment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Observing flow of He II with unsupervised machine learning

Abstract Time dependent observations of point-to-point correlations of the velocity vector field (structure functions) are necessary to model and understand fluid flow around complex objects. Using thermal gradients, we observed fluid flow by recording fluorescence of $${\text{He}}_{2}^{*}$$ He 2 ∗ excimers produced by neutron capture throughout a ~ cm 3 volume. Because the photon emitted by an excited excimer is unlikely to be recorded by the camera, the techniques of particle tracking (PTV) and particle imaging (PIV) velocimetry cannot be applied to extract information from the fluorescence of individual excimers. Therefore, we applied an unsupervised machine learning algorithm to identify light from ensembles of excimers (clusters) and then tracked the centroids of the clusters using a particle displacement determination algorithm developed for PTV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compositional analysis of organosolv poplar lignin by using high-performance liquid chromatography/high-resolution multi-stage tandem mass spectrometry

Organosolv treatment is an efficient and environmentally friendly process to degrade lignin into small compounds. The capability of characterizing the individual compounds in the complex mixtures formed upon organosolv treatment is essential for the optimization of the further lignin conversion processes and for the rational genetic engineering of plants used to produce lignin in order to improve lignin properties. In this study, an organosolv poplar lignin sample was initially analyzed by high-resolution mass spectrometry coupled with negative-ion mode electrospray ionization ((-)ESI HRMS). Lignin monomers and dimers were found to constitute the majority of the compounds in the organosolv lignin sample. Larger lignin oligomers, such as trimers and tetramers, and some not lignin-related compounds, were also detected. A high-performance liquid chromatograph/linear quadrupole ion trap/orbitrap mass spectrometer capable of multi-stage high-resolution tandem mass spectrometry experiments (HRMS n ), equipped with an (-)ESI source (HPLC/(-)ESI HRMS n ), was employed to separate the unknown compounds in the organosolv mixture and to obtain structural information for the deprotonated compounds via collision-activated dissociation (CAD) HRMS n experiments. To improve the understanding of the CAD behavior of deprotonated lignin-related compounds, 16 deprotonated model compounds with different functionalities and linkage types were examined. This approach enabled the assignment of likely structures for several lignin monomers, dimers, trimers, and tetramers, and some not lignin-related compounds, most likely fatty acids. Based on the proposed structures, compounds in the organosolv lignin sample contain β-O-4, 5–5, β-5, and possibly also 4-O-5 linkages. Most compounds contain G- and S-monomeric units although a small amount of H-units were also detected.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling polymer crystallisation induced by a moving heat sink

Recent experimental work has shown that polymer crystallisation can be used to “move” and organize nanoparticles (NP). As a first effort at modeling this situation, we consider the classical Stefan problem but with the modification that polymer crystallisation does not occur at a single temperature. Rather, the rate of crystallisation is proportional to its subcooling, and here we employ a form inspired by the classical Avrami model to describe this functional form. Our results for the movement of the polymer crystallisation front, as defined as the point where the crystallinity is 50%, closely track the results of the classical Stefan problem. Thus, at this level of approximation, the crystallisation kinetics of the polymer do not cause qualitative changes to the physics of this situation. Inspired by this fact we study the more interesting situation where the directional recrystallisation of a polymer melt is considered, e.g. , through the application of a moving heat sink over an initially molten polymer, reminiscent of a processing technique termed zone annealing. The polymer crystallisation shows that a steady state exists for a range of sink velocities. The solid–melt interface moves slightly ahead of the sink but at the same velocity. The steady-state distance between the sink and the interface decreases with increasing sink velocity – this is a consequence of the excess cooling provided by the sink over what is required to crystallise the melt. The most interesting new result is that the temperature of the crystal–melt interface decreases with increasing sink velocity. This is in line with the ansatz of Lovinger and Gryte who suggested that larger zone annealing velocities correspond to progressively larger effective undercoolings at which polymer crystallisation occurs.

36 MATERIALS SCIENCE↗

Mechanistic insights into metal, nitrogen doped carbon catalysts for oxygen reduction: progress in computational modeling

We report metal and nitrogen doped carbon materials (denoted as M–N–C) synthesized through high-temperature pyrolysis have been found to exhibit activity for oxygen reduction reaction (ORR) approaching that of Pt and electrochemical stability higher than previous MN 4 -containing macrocyclic molecular catalysts. Tremendous efforts have thus been devoted to the advancement of M–N–C catalysts as an economical alternative to Pt-based catalysts for proton exchange membrane fuel cell cathodes with a focus on simultaneously improving activity and stability. To this end, novel computational modeling techniques have been developed and applied to acquire knowledge crucial for accelerating the pace of M–N–C catalyst development. In this review, recent progress in computational method development, as well as the predictions of chemical structure of active sites, reaction pathways, ORR kinetics, and catalyst stability in electrochemical environments, are critically surveyed. Moreover, the crucial role of computational modeling to elucidate the functional mechanism of M–N–C catalysts for ORR in acid media and enable rational design of M–N–C catalysts is discussed with a visionary outlook for the field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Exploiting electricity market dynamics using flexible electrolysis units for retrofitting methanol synthesis

Here we investigate the economic viability of integrating flexible electrolysis units to produce hydrogen in methanol synthesis processes. Specifically, we investigate whether this approach can help reduce methanol production costs by strategically exploiting dynamics of electricity markets. Our study integrates high-fidelity process simulations, optimization tools, and microkinetic modeling (informed by density functional theory) to conduct detailed techno-economic analyses and to compare performance against traditional processes that use hydrogen produced via steam-methane reforming (SMR). We also use this approach to estimate the levelized cost of hydrogen (LCOH) as a function of time-varying electricity prices (from day-ahead and real-time prices) and of key techno-economic parameters. Our results show that the proposed electrification framework is cost-competitive under certain electricity market conditions. Specifically, we find that, when the electrolysis system is operated in flexible mode (and can respond to dynamics of electricity markets), the associated electricity cost nearly collapses to zero. Conversely, when the unit is not flexible (and cannot respond to markets), the electricity cost comprises 60% of the total cost. Our results also reveal that the LCOH of the flexible electrolysis system participating in real-time electricity markets is 31% lower than the LCOH obtained from SMR. Overall, this indicates that exploiting the dynamics of electricity markets can make hydrogen production cost-competitive and this can lead to viable alternatives to electrify methanol production and other hydrogen-based processes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Co-orchestration of multiple instruments to uncover structure–property relationships in combinatorial libraries

The rapid growth of automated and autonomous instrumentation brings forth opportunities for the co-orchestration of multimodal tools that are equipped with multiple sequential detection methods or several characterization techniques to explore identical samples. This is exemplified by combinatorial libraries that can be explored in multiple locations via multiple tools simultaneously or downstream characterization in automated synthesis systems. In co-orchestration approaches, information gained in one modality should accelerate the discovery of other modalities. Correspondingly, an orchestrating agent should select the measurement modality based on the anticipated knowledge gain and measurement cost. Herein, we propose and implement a co-orchestration approach for conducting measurements with complex observables, such as spectra or images. The method relies on combining dimensionality reduction by variational autoencoders with representation learning for control over the latent space structure and integration into an iterative workflow via multi-task Gaussian Processes (GPs). This approach further allows for the native incorporation of the system's physics via a probabilistic model as a mean function of the GPs. We illustrate this method for different modes of piezoresponse force microscopy and micro-Raman spectroscopy on a combinatorial Sm-BiFeO3 library. However, the proposed framework is general and can be extended to multiple measurement modalities and arbitrary dimensionality of the measured signals.

47 OTHER INSTRUMENTATION↗

SLAB: simultaneous labeling and binding affinity prediction for protein–ligand structures

Machine learning models are often used as scoring functions to predict the binding affinity of a protein–ligand complex. These models are trained with limited amounts of data with experimentally measured binding affinity values. A large number of compounds are labeled inactive through single-concentration screens without measuring binding affinities. These inactive compounds, along with the active ones, can be used to train binary classification models, while regression models are trained using compounds with binding affinities only. However, the classification and regression tasks are often handled separately, without sharing the learned feature representations. In this paper, we propose a novel model architecture that jointly performs regression and classification objectives, aiming to maximize data utilization and improve predictive performance by leveraging two complementary tasks. In our setup, the regression yields the binding affinity, whereas the classification task yields the label as active or inactive. We demonstrate our method using PDBbind, the standard 3D structure database, as well as a dataset of flavivirus protease compounds with binding affinity data. Our experiments show that the new joint training strategy improves the accuracy of the model, increasing applicability in various practical drug screening scenarios.

Biological and medical sciences↗

MIRACL Co‐Simulation platform for control and operation of distributed wind in microgrid

Abstract This paper presents a co‐simulation platform (CSP) developed as a part of the “Microgrids, Infrastructure Resilience and Advanced Controls Launchpad (MIRACL)” project, hereafter called MIRACL‐CSP, to allow simulation‐based testing, demonstration, and evaluations of distributed wind under various grid operating conditions. MIRACL‐CSP provides modular interactions among the power distribution system, distributed wind, and utility decision‐making framework. A hierarchical engine for large‐scale infrastructure co‐simulations (HELICS) is used as the core engine of MIRACL‐CSP to establish time and information coordination among the MIRACL‐CSP modules. In this study, MIRACL‐CSP performance is demonstrated using the IEEE 123‐node test distribution grid modelled in GridLAB‐D and the utility decision support application modelled in Python. The functionality of MIRACL‐CSP is demonstrated through various grid operational scenarios in a microgrid and networked microgrid environment. A co‐simulation study is described that demonstrated MIRACL‐CSP capabilities for the microgrid operational scenario in the IEEE 123‐node test system. MIRACL‐CSP is a generic platform that facilitates distributed wind research for different test systems, applications, and valuation of distributed wind.

17 WIND ENERGY↗

Influence of globally spin-aligned vector mesons to the measurements of the chiral magnetic effect in heavy-ion collisions

The chiral magnetic effect (CME) in high-energy heavy-ion collisions arises from the interplay between the chirality imbalance and the intense magnetic field and will cause a charge separation along the magnetic field direction. While the CME search is still ongoing in experiments, the non-CME contributions need to be excluded from the CME observables. In this work, we examine the influence of globally spin-aligned ρ mesons on the γ 112 correlator, the R Ψ 2 (ΔS) correlator, and the signed balance functions, via a toy model and a multiphase transport model (AMPT). We find that the CME observables are sensitive to the 00-component of the spin density matrix, ρ 00 : they receive positive (negative) contributions when ρ 00 is larger (smaller) than 1/3.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigating the Impact of Temporal and Directional Traffic Distribution on Crash Frequencies

Safety Performance Functions (SPFs) are mathematical models that establish relationships between the frequency of various crash types and site-specific characteristics, serving as essential tools for traffic safety analysis and roadway design. Traditional SPFs, however, often overlook the temporal fluctuations in traffic flow (such as peak-hour surges) and directional imbalances between opposing traffic streams. These traffic patterns can exacerbate congestion, disrupt driver behavior, and create unexpected conflict points, potentially leading to increased crash frequencies and more severe accidents. In light of this gap, this study aims to explore the potential of incorporating K-factors (representing peak-hour traffic proportions) and D-factors (reflecting the imbalance of directional traffic) into the development of SPFs to assess whether these factors can effectively represent the impact of temporal and spatial traffic distribution on roadway safety. Using crash data from Pennsylvania urban-suburban collector roadways, it is found that the D-factor plays a significant role in predicting the frequency of total crashes, fatal + injury crashes, and angle crashes, with positive coefficient signs indicating that higher directional imbalances correspond to increased crash risks. Similarly, the K-factor emerges as a critical predictor for fatal + injury crashes and rear-end crashes, with negative coefficients suggesting that a more pronounced traffic peak is associated with a reduction in expected crash frequencies. These results highlight the importance of accounting for uneven traffic distribution in both time and direction when developing SPFs, offering deeper insights into crash patterns and supporting more effective safety interventions and roadway designs.

Xu, Guanhao [ORNL] (ORCID:0000000214326357)↗

Full waveform analysis for high pressure ultrasonic measurement

In this paper we report a data analysis protocol for ultrasonic velocity measurements carried out in a synchrotron x-ray facility based multi-anvil high-pressure apparatus. Synthetic ultrasonic signals for the time period between echoes from the two ends of the sample assembly are created with a few simple parameters as the echo waveform is modeled by the first echo from the anvil-assembly boundary. Each echo is modeled as a delta function with an amplitude and arrival time and convolved with this waveform. The final waveform, fit to the data, enables analysis of signals with overlapped echoes as they are common in small cell assemblies and high pressure.

47 OTHER INSTRUMENTATION↗

Solution phase high repetition rate laser pump x-ray probe picosecond hard x-ray spectroscopy at the Stanford Synchrotron Radiation Lightsource

We present a dedicated end-station for solution phase high repetition rate (MHz) picosecond hard x-ray spectroscopy at beamline 15-2 of the Stanford Synchrotron Radiation Lightsource. A high-power ultrafast ytterbium-doped fiber laser is used to photoexcite the samples at a repetition rate of 640 kHz, while the data acquisition operates at the 1.28 MHz repetition rate of the storage ring recording data in an alternating on-off mode. The time-resolved x-ray measurements are enabled via gating the x-ray detectors with the 20 mA/70 ps camshaft bunch of SPEAR3, a mode available during the routine operations of the Stanford Synchrotron Radiation Lightsource. As a benchmark study, aiming to demonstrate the advantageous capabilities of this end-station, we have conducted picosecond Fe K-edge x-ray absorption spectroscopy on aqueous [FeII(phen)3]2+, a prototypical spin crossover complex that undergoes light-induced excited spin state trapping forming an electronic excited state with a 0.6–0.7 ns lifetime. In addition, we report transient Fe Kβ main line and valence-to-core x-ray emission spectra, showing a unique detection sensitivity and an excellent agreement with model spectra and density functional theory calculations, respectively. Notably, the achieved signal-to-noise ratio, the overall performance, and the routine availability of the developed end-station have enabled a systematic time-resolved science program using the monochromatic beam at the Stanford Synchrotron Radiation Lightsource.

47 OTHER INSTRUMENTATION↗

Probing nonlinear velocity profiles of shear-thinning, nematic platelet dispersions in Couette flow using x-ray photon correlation spectroscopy

In this study, we report experiments employing x-ray photon correlation spectroscopy (XPCS) to characterize the velocity profiles of complex fluids in Couette flow. The approach involves modeling the XPCS correlation functions obtained with the incident x-ray beam passing tangentially through the Couette cell gap at various distances from the inner wall. We first demonstrate the technique with measurements on a dilute colloidal dispersion in the Newtonian liquid glycerol, where the expected linear velocity profiles are recovered. We then employ the technique to map the shear-rate-dependent velocity profiles of a shear-thinning dispersion of nematically ordered Gibbsite platelets. The nonlinear velocity profiles of the Gibbsite dispersion include a narrow slip region adjacent to the outer wall and a band with a small velocity gradient in the interior of the gap that evolves into a region increasingly resembling plug flow with increasing shear rate. Variations in the velocity profile along the vorticity direction indicate an instability in the interface between this region of small velocity gradient and a region of high velocity gradient near the inner wall. The analysis of the small-angle scattering patterns provides information about the spatial and temporal variations in the nematic order of the Gibbsite dispersion and their coupling to the velocity profile. Additional potential applications of this XPCS-based technique and comparisons with established methods for characterizing velocity profiles are discussed.

42 ENGINEERING↗