Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 91 records · Page 5

Tool Chain for Simulations of Bi-Filar Coil Winding for Fast Quench Protection

The advancement of high-field magnets utilizing high-temperature superconductors (HTS) brings about complex challenges, especially in quench detection and protection. Traditional methods often fall short due to the inherently slow quench propagation in HTS materials. One promising approach to overcome this involves using a bifilar winding configuration, where two conductors are placed side by side. Under normal operation, they function in series, but during a quench event, they switch to an anti-parallel mode. This shift reduces the differential inductance of the coil to near zero, enabling rapid current oscillations through a capacitor discharge. The resulting high-frequency current flow leads to swift, uniform heating, triggering a full-coil quench within microseconds. Moreover, the strong mutual coupling between the two windings significantly reduces electrical noise in voltage measurements. In this work, we explore the viability of this concept by designing, constructing, and testing a REBCO bifilar racetrack coil in liquid nitrogen. We also present a validated simulation model that closely mirrors the coil's dynamic behavior under these conditions, aligning well with experimental observations.

Jayathilaka, Rehan [Northern Illinois U.]↗

Global geodesic acoustic mode in an ideal magnetohydrodynamic tokamak plasma

A concise and transparent second order ordinary differential equation (ODE) describing the radial structure of the global geodesic acoustic mode (GAM) is analytically presented in a low-β tokamak plasma. The large-aspect-ratio and circular cross section are assumed to linearize the ideal magnetohydrodynamic equations. We show clearly how finite β-dependent terms affect the global GAM frequency and radial mode structure. Here, a typical Wentzel–Kramers–Brillouin form of solution is found for some reversed shear equilibria. For some other equilibria with lower β, even also in a reversed shear tokamak, the GAM continuum is upraised by the high order β-dependent terms so that its maximum is beyond ω G , where ω G is the classical local frequency of GAM. As a result, no self-consistent solution to the ODE can be found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparing Transcriptomic Points of Departure to Apical Effect Concentrations For Larval Fathead Minnow Exposed to Chemicals with Four Different Modes Of Action

It is postulated that below a transcriptomic-based point of departure, adverse effects are unlikely to occur, thereby providing a chemical concentration to use in screening level hazard assessment. The present study extends previous work describing a high-throughput fathead minnow assay that can provide full transcriptomic data after exposure to a test chemical. One-day post-hatch fathead minnows were exposed to ten concentrations of three representatives of four chemical modes of action: organophosphates, ecdysone receptor agonists, plant photosystem II inhibitors, and estrogen receptor agonists for 24 h. Concentration response modeling was performed on whole body gene expression data from each exposure, using measured chemical concentrations when available. Transcriptomic points of departure in larval fathead minnow were lower than apical effect concentrations across fish species but not always lower than toxic effect concentrations in other aquatic taxa like crustaceans and insects. The point of departure was highly dependent on measured chemical concentration which were often lower than the nominal concentration. Differentially expressed genes between chemicals within modes of action were compared and often showed statistically significant overlap. In addition, reproducibility between identical exposures using a positive control chemical (CuSO 4 ) and variability associated with the transcriptomic point of departure using in silico sampling were considered. Results extend a transcriptomic-compatible fathead minnow high-throughput assay for possible use in ecological hazard screening.

59 BASIC BIOLOGICAL SCIENCES↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Test of lepton flavor universality and search for lepton flavor violation in B → Kℓℓ decays

We present measurements of the branching fractions for the decays B → Kμ + μ - and B → Ke + e - , and their ratio (R K ), using a data sample of 711 fb -1 that contains 772 × 10 6 B$\bar{B}$ events. The data were collected at the Υ(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e - collider. The ratio R K is measured in five bins of dilepton invariant-mass-squared (q 2 ): q 2 $\epsilon$ (0.1, 4.0), (4.00, 8.12), (1.0, 6.0), (10.2, 12.8) and (> 14.18) GeV 2 /c 4 , along with the whole q 2 region. The R K value for q 2 $\epsilon$ (1.0, 6.0) GeV 2 /c 4 is 1.03$^{+0.28}_{-0.24}$ ± 0.01. The first and second uncertainties listed are statistical and systematic, respectively. All results for RK are consistent with Standard Model predictions. We also measure CP-averaged isospin asymmetries in the same q 2 bins. The results are consistent with a null asymmetry, with the largest difference of 2.6 standard deviations occurring for the q 2 $\epsilon$ (1.0, 6.0) GeV 2 /c 4 bin in the mode with muon final states. The measured differential branching fractions, dB/dq 2 , are consistent with theoretical predictions for charged B decays, while the corresponding values are below the expectations for neutral B decays. We have also searched for lepton-flavor-violating B → Kμ ± e ∓ decays and set 90% confidence-level upper limits on the branching fraction in the range of 10 -8 for B + → K + μ ± e ∓ , and B 0 → K 0 μ ± e ∓ modes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ultrafast method for scanning tunneling spectroscopy

A scanning tunneling microscope (STM) combines unique capabilities in imaging and spectroscopy with atomic precision, and it can obtain energy-resolved spectroscopic data with atomic resolution. In this paper, we utilize a recently proposed modification to the STM feedback control loop to acquire high quality d 2 I/dV 2 images. We have developed a constant differential conductance imaging method by closing the STM feedback loop with a high precision dI/dV measurement. In this mode, the tip’s vertical position is adjusted so as to keep the differential conductance constant during raster scanning of the surface. Based on this imaging mode, we propose a new technique to acquire fast and reliable scanning tunneling spectroscopy (STS) data simultaneously with the imaging.

47 OTHER INSTRUMENTATION↗

Distinguishing fissile uranium isotopes using an active well neutron coincidence counter

Proposed thorium-based nuclear fuel cycles are likely to require quantification and verification of 233 U within nuclear material. Because of their similar fission cross sections, active neutron nondestructive assay (NDA) systems may respond similarly to 233 U and 235 U. Traditional safeguards equipment has been optimized for 235 U and 238 U quantification associated with conventional uranium/plutonium fuel cycles and may not be directly applicable to 233 U quantification when mixed with other actinides. This work used models of the large volume active well coincidence counter (LV-AWCC) at Oak Ridge National Laboratory to evaluate the performance of this neutron NDA system to differentiate fissile uranium isotopes. The models were developed to simulate NDA system performance in response to a number of triangular radiation signature training device sources within the central cavity or well. This work predicted that the LV-AWCC can effectively differentiate 233 U from 235 U in certain modes of operation. In active mode, the LV-AWCC with the cadmium liner results in different doubles count rates between the fissile isotopes for a given fissile uranium mass. Without the cadmium liner, the uranium isotopes provide a statistically indistinguishable doubles count rate response for the fissile masses considered in this work (up to approximately 150 g). The cadmium liner serves to harden the neutron interrogation spectrum, which better exploits the notable difference in the 233 U and 235 U fission cross sections at approximately 1 eV. In passive mode, the two fissile isotopes exhibit different doubles and singles count rates regardless of liner presence because the passive source strength of 233 U is approximately 2 orders of magnitude stronger than that of 235 U due to the shorter half-life and correspondingly higher (α, n) yield. We conclude that using neutron interrogation in the LV-AWCC, two measurements are needed to quantify 233 U content in mixed uranium items. The first measurement is used to determine the total fissile uranium mass using a mode that cannot distinguish fissile isotopes (i.e., where a similar response is observed for both fissile uranium isotopes such as active doubles without cadmium or using a thermal neutron interrogation source). In conclusion, the second measurement is used to determine the 233 U content by using a differentiating technique (e.g., passive doubles, passive doubles to singles ratio, active doubles with cadmium).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗

Nonlinear modeling and performance analysis of cracked beam microgyroscopes

Cracks constitute a common structural defect in microelectromechanical systems that may arise during manufacturing, mechanical fatigue, or shock loading. Areas weakened by cracks increase the flexibility of the microstructure. Furthermore, depending on the crack severity, the increased flexibility can cause dramatic changes in the static and dynamic behaviors of electrically actuated MEMS. In this study, a numerical investigation of the depth and location of multiple edge cracks on the performance of beam microgyroscopes is conducted. Applying the Griffith strain energy release theorem for an edge crack, the vibrational characteristics of microgyroscopes with multiple cracks including the static deflection, natural frequencies, and mode shapes are studied analytically. Then, the differential quadrature method is employed to model the nonlinear dynamic behaviors of the cracked gyroscope system. Cracks initiated along the driving or sensing directions of the microgyroscope are found to have differing effects on its nonlinear dynamic response. This numerical study reveals that the onset of severe single cracks or a series of multiple cracks can significantly affect the sensitivity of the microgyroscope to base rotations and therefore leads to a degradation in the performance of the damaged MEMS sensor. This study also shows that a potential design approach for broadband microgyroscope systems manufactured with cracks can be proposed.

42 ENGINEERING↗

Charge transport in topological graphene nanoribbons and nanoribbon heterostructures

Although it is generally accepted that structural parameters like width, shape, and edge structure crucially affect the electronic characteristics of graphene nanoribbons (GNRs), the exact relationship between geometry and charge transport remains largely unexplored. Here, in this paper, we present in situ through-transport measurements of various topological GNRs and GNR heterostructures by lifting the ribbon with the tip of a scanning tunneling microscope. At the same time, we develop a comprehensive transport model that enables us to understand various features, such as obscuring of localized states in through transport, the effect of topology on transport, as well as negative differential conductance in heterostructures with localized electronic modes. The combined experimental and theoretical efforts described in this paper serve to elucidate general charge transport phenomena in GNRs and GNR heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantification of the Radiative Effect of Aerosol–Cloud Interactions in Shallow Continental Cumulus Clouds

The indirect radiative effect of aerosol variability on shallow cumulus clouds is realized in nature with considerable concurrent meteorological variability. Large-eddy simulations constrained by observations at a continental site in Oklahoma are performed to represent the variability of different meteorological states on days with different aerosol conditions. The total radiative effect of this natural covariation between aerosol and other meteorological drivers of total cloud amount and albedo is quantified. The changes to these bulk quantities are used to understand the response of the cloud radiative effect to aerosol–cloud interactions (ACI) in the context of concurrent processes, as opposed to attempting to untangle the effect of individual processes on a case-by-case basis. Mutual information (MI) analysis suggests that meteorological variability masks the strength of the relationship between cloud drop number concentration and the cloud radiative effect. This is shown to be mostly due to variation in solar zenith angle and cloud field horizontal heterogeneity masking the relationship between cloud drop number and cloud albedo. By combining MI and more traditional differential analyses, a framework to identify important modes of covariation between aerosol, clouds, and meteorological conditions is developed. This shows that accounting for solar zenith angle variation and implementing an albedo bias correction increases the detectability of the radiative effects of ACI in simulations of shallow cumulus.

54 ENVIRONMENTAL SCIENCES↗

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗

Shot noise and differential conductance as signatures of putative topological superconductivity in FeSe 0.45 ⁢Te 0.55

We present a theory for the differential shot noise, dS/dV, as measured via shot-noise scanning tunneling spectroscopy, and the differential conductance, dI/dV, for tunneling into Majorana zero modes (MZMs) in the putative topological superconductor FeSe 0.45 Te 0.55 . We demonstrate that for tunneling into chiral Majorana edge modes near domain walls, as well as MZMs localized in vortex cores and at the end of defect lines, dS/dV vanishes whenever dI/dV reaches a quantized value proportional to the quantum of conductance. Furthermore, these results are independent of the particular orbital tunneling path, thus establishing a vanishing dS/dV concomitant with a quantized dI/dV, as universal signatures for Majorana modes in two-dimensional topological superconductors, irrespective of the material’s specific complex electronic band structure.

Iron-based superconductors↗

FLORA Equilibrium and Stability Code Archive

FLORA solves, in a 2-D domain (radial and axial dimensions) with a specified azimuthal Fourier mode, for the linearized stability of a long, thin, axisymmetric plasma equilibrium in an applied magnetic field. Before the stability equation is solved, FLORA solves a set of simple equations for pressure balance that specify the equilibrium magnetic field and plasma pressure in the long-thin limit given a simplified description of the magnetic coils. It uses an initial-value method for the linear stability problem in which an equilibrium is given an initial perturbation to its magnetic field, and the temporal behavior of the perturbation is followed. The perturbation has been Fourier expanded in the azimuthal coordinate; each azimuthal mode must be examined separately. The complex partial differential equation of motion for the perturbed radial displacement of the field lines is solved as a coupled system of two real p.d.e.'s and the solution consists of two parts, the real part and the imaginary part. The system is solved by bringing the coupling terms in each equation to the right side and using an iterative technique.

Cohen, Bruce [Lawrence Livermore National Laborato↗

Differential Privacy in Grid Kitchen: Implementation & Software Documentation

Sharing of power grid feeder models faces significant challenges due to the potential risk of exposing sensitive operational information. Traditional anonymization techniques have shown notable limitations in other sensitive domains, as evidenced by documented re-identification attacks that combine supposedly anonymized datasets with auxiliary information, raising concerns that similar vulnerabilities could affect power grid data. Consequently, there is a pressing need for a more rigorous privacy protection strategy that not only delivers formal mathematical guarantees but also preserves the analytical value of the shared models. To address this challenge, we have enhanced the Grid Kitchen framework by implementing differential privacy mechanisms within the distribution model dehydration pipeline. This implementation carefully calibrates and applies noise to sensitive attributes in feeder models according to configurable privacy levels—low, moderate, and high—each offering different balances between data utility and privacy protection. Our approach uses established noise functions (Gaussian for continuous data and Discrete Laplace for integer values) with parameters carefully calibrated so that the impact of individual data points is effectively masked in the final output. The integration leverages our Noise Catalog, which we developed to categorize feeder model properties by component type, data type, and sensitivity. This catalog guides the application of appropriate noise functions and privacy parameters ($\varepsilon$ and $\delta$) to each attribute, ensuring consistent privacy protection across the model while maintaining its structural integrity and analytical usefulness. This implementation also includes evaluation tools that allow model owners to assess the impact of privacy-preserving transformations before sharing data with external parties. This report provides documentation for the differential privacy capabilities added to the Grid Kitchen project. It includes a primer on differential privacy concepts and their importance in modern data sharing, details the architecture of our implementation, explains the privacy modes and parameter configurations, and offers practical guidance on using the code for applying differential privacy to grid feeder models. Through examples and code snippets, we demonstrate the effective application of these privacy-enhancing technologies, enabling utility operators and researchers to confidently share grid data while protecting sensitive information.

24 POWER TRANSMISSION AND DISTRIBUTION↗