Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ATOMIC AND MICROSTRUCTURAL ORIGINS OF STORED ENERGY RELEASE IN NEUTRON-IRRADIATED SILICON CARBIDE

Here, we employ a combination of advanced synchrotron-based scattering characterization techniques to understand and unravel the atomic origins of the colossal stored energy release in neutron irradiated silicon carbide. The quantification of the neutron irradiation-induced defects, and their impact on the structureproperty relationship is important for the design and safe operation of advanced fission and fusion reactors. Our experimental results show that the atomic structure in the as irradiated samples is significantly perturbed due to a large fraction of vacancy and interstitial type defects that lead to complex microstructures and additional components in the x-ray diffraction and pair distribution function results. We directly correlate the stored energy release to the recovery of the sub-lattices, with pair distribution function analysis highlighting that the carbon interstitial and vacancy-type defects contribute to stored energy more than those of silicon. We find these results to be striking and believe our discoveries to be timely and noteworthy given the technological importance of silicon carbide to the nuclear fission and fusion communities.

Sprouster, D.J.↗
DOE OSTI2024

Machine Learning Driven Sensitivity Analysis of E3SM Land Model Parameters for Wetland Methane Emissions

Methane (CH 4 ) is globally the second most critical greenhouse gas after carbon dioxide, contributing to 16%–25% of the observed atmospheric warming. Wetlands are the primary natural source of methane emissions globally. However, wetland methane emission estimates from biogeochemistry models contain considerable uncertainty. One of the main sources of this uncertainty arises from the numerous uncertain model parameters within various physical, biological, and chemical processes that influence methane production, oxidation, and transport. Sensitivity Analysis (SA) can help identify critical parameters for methane emission and achieve reduced biases and uncertainties in future projections. This study performs SA for 19 selected parameters responsible for critical biogeochemical processes in the methane module of the Energy Exascale Earth System Model (E3SM) land model (ELM). The impact of these parameters on various CH 4 fluxes is examined at 14 FLUXNET- CH 4 sites with diverse vegetation types. Given the extensive number of model simulations needed for global variance-based SA, we employ a machine learning (ML) algorithm to emulate the complex behavior of ELM methane biogeochemistry. We found that parameters linked to CH 4 production and diffusion generally present the highest sensitivities despite apparent seasonal variation. Comparing simulated emissions from perturbed parameter sets against FLUXNET-CH 4 observations revealed that better performances can be achieved at each site compared to the default parameter values. This presents a scope for further improving simulated emissions using parameter calibration with advanced optimization techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Systematic planning of moving target defence for maximising detection effectiveness against false data injection attacks in smart grid

Abstract Moving target defence (MTD) has been gaining traction to thwart false data injection attacks against state estimation (SE) in the power grid. MTD actively perturbs the reactance of transmission lines equipped with distributed flexible AC transmission system (D‐FACTS) devices to falsify the attacker's knowledge about the system configuration. However, the existing literature has not systematically studied what influences the detection effectiveness of MTD and how it can be improved based on the topology analysis. These problems are tackled here from the perspective of an MTD plan in which the D‐FACTS placement is determined. We first exploit the relation between the rank of the composite matrix and the detecting effectiveness. Then, we rigorously derive upper and lower bounds on the attack detecting probability of MTDs with a given rank of the composite matrix. Furthermore, we analyse existing planning methods and highlight the importance of bus coverage by D‐FACTS devices. To improve the detection effectiveness, we propose a novel graph theory–based planning algorithm to retain the maximum rank of the composite matrix while covering all necessary buses. Comparative results on multiple systems show the high detecting effectiveness of the proposed algorithm in both DC‐ and AC‐SE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Not even 6 dB: Gaussian quantum illumination in thermal background

Abstract In analyses of target detection with Gaussian state transmitters in a thermal background, the thermal occupation is taken to depend on the target reflectivity in a way which simplifies the analysis of the symmetric quantum hypothesis testing problem. However, this assumption precludes comparison of target detection performance between an arbitrary transmitter and a vacuum state transmitter, i.e. ‘detection without illumination’, which is relevant in a bright thermal background because a target can be detected by its optical shadow or some other perturbation of the background. Using a target-agnostic thermal environment leads to the result that the oft-claimed 6 dB possible reduction in the quantum Chernoff exponent for a two-mode squeezed vacuum transmitter over a coherent state transmitter in high-occupation thermal background is an unachievable limiting value, only occurring in a limit in which the target detection problem is ill-posed. Further analyzing quantum illumination in a target-agnostic thermal environment shows that a weak single-mode squeezed transmitter performs worse than ‘no illumination’, which is explained by the noise-increasing property of reflected low-intensity squeezed light.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Sign-OPT: A Query-Efficient Hard-label Adversarial Attack

We study the most practical problem setup for evaluating adversarial robustness of a machine learning system with limited access: the hard-label black-box attack setting for generating adversarial examples, where limited model queries are allowed and only the decision is provided to a queried data input. Several algorithms have been proposed for this problem but they typically require huge amount (>20,000) of queries for attacking one example. Among them, one of the state-of-the-art approaches (Cheng et al., 2019) showed that hard-label attack can be modeled as an optimization problem where the objective function can be evaluated by binary search with additional model queries, thereby a zeroth order optimization algorithm can be applied. In this paper, we adopt the same optimization formulation but propose to directly estimate the sign of gradient at any direction instead of the gradient itself, which enjoys the benefit of single query. Using this single query oracle for retrieving sign of directional derivative, we develop a novel query-efficient Sign-OPT approach for hard-label black-box attack. We provide a convergence analysis of the new algorithm and conduct experiments on several models on MNIST, CIFAR-10 and ImageNet. We find that Sign-OPT attack consistently requires 5X to 10X fewer queries when compared to the current state-of-the-art approaches, and usually converges to an adversarial example with smaller perturbation.

Cheng, Minhao↗
DOE OSTI2025

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗
DOE OSTI2023

The persistent nonperturbative charm enigma

The question of the existence and possible magnitude of nonperturbative (often called “intrinsic”) charm in the proton has long confounded attempts to cleanly isolate such a contribution in global analyses of high-energy experiments. In this letter, we show that the available (non)perturbative QCD theory and hadronic data have still not developed to a sufficient level to clearly resolve this problem. We highlight a number of challenging aspects that must be confronted in extracting nonperturbative charm in PDF fits, and in so doing, present an updated next-to-next-to-leading order CT analysis of fitted charm, CT18 FC, which we also compare to recent studies. We outline the theory developments and future data needed to make progress on this subject.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Perturbative QCD contribution to transverse single spin asymmetries in the Drell-Yan process and SIDIS

In a previous publication [Benić , ], we have computed the perturbative QCD contribution to transverse single spin asymmetries (SSAs) in semi-inclusive deep inelastic scattering (SIDIS) involving the g T ( x ) distribution. In this paper, we first present a more efficient derivation of the asymmetries which is applicable to both transverse and longitudinal SSAs and correct some inconsistencies in our previous calculation. We then adapt the method to compute transverse SSAs in the Drell-Yan process proportional to g T ( x ) and its gluonic counterpart and discuss the crossing symmetry between the results for SIDIS and the Drell-Yan process. Finally, we present numerical results for various asymmetries measurable at the EIC, RHIC, COMPASS, and Fermilab (SpinQuest), including also part of the genuine twist-three corrections to g T ( x ) from a recent global analysis. We find that the asymmetries can reach percent-level magnitude, if the kinematics predominantly probes the large- x (valence) region of the polarized proton, but remain at subpercent levels otherwise. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Measurements of the differential cross sections of the production of Z + jets and $\gamma$ + jets and of Z boson emission collinear with a jet in pp collisions at $ \sqrt{s} $ = 13 TeV

Measurements of the differential cross sections of Z + jets and γ + jets production, and their ratio, are presented as a function of the boson transverse momentum. Measurements are also presented of the angular distribution between the Z boson and the closest jet. The analysis is based on pp collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 35.9 fb$^{−1}$ recorded by the CMS experiment at the LHC. The results, corrected for detector effects, are compared with various theoretical predictions. In general, the predictions at higher orders in perturbation theory show better agreement with the measurements. This work provides the first measurement of the ratio of the differential cross sections of Z + jets and γ + jets production at 13 TeV, as well as the first direct measurement of Z bosons emitted collinearly with a jet.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

SST-TG-P1F4R3200: Decaying Stably-Stratified Turbulence (SST), Initialized Using Taylor-Green Vortices (TG) at Prandtl Number Pr=1, Froude Number Fr=4, Reynolds Number Re=3200

This dataset comprises direct numerical simulations (DNS) of decaying stably-stratified turbulence influenced by a linear background density gradient, initialized using an array of Taylor-Green vortices, as described in [Riley & de Bruyn Kops (2003)](https://doi.org/10.1063/1.1578077). The initial Prandtl, Froude, and Reynolds numbers are (Pr, Fr, Re) = (1, 4, 3200). A total of 15,000 snapshots are recorded at uniform time intervals, each with a spatial resolution of 512x512x256 grid points. Four flow variables are associated with each snapshot: the three velocity components (u,v,w) and the perturbed density field (rho) away from the background gradient. All fields are stored in binary format (32-bit little-endian), each with a size of 255 MB, yielding a total dataset size of 15.3 TB. Further details are referenced in the attached README file, and a current list of publications and associated analysis tools are provided at https://stratified-turbulence.github.io/web/.

42 ENGINEERING↗
DOE OSTI2025

SST-TG-P50F4R3200: Decaying Stably-Stratified Turbulence (SST), Initialized Using Taylor-Green Vortices (TG) at Prandtl Number Pr=50, Froude Number Fr=4, Reynolds Number Re=3200

This dataset comprises direct numerical simulations (DNS) of decaying stably-stratified turbulence influenced by a linear background density gradient, initialized using an array of Taylor-Green vortices, extending the Pr=1 simulations performed in [Riley and de Bruyn Kops (2003)](https://doi.org/10.1063/1.1578077). The initial Prandtl, Froude, and Reynolds numbers are (Pr, Fr, Re) = (50, 4, 3200). A total of 1,680 snapshots are recorded at uniform time intervals, each with a spatial resolution of 3584x3584x1792 grid points. Four flow variables are associated with each snapshot: the three velocity components (u,v,w) and the perturbed density field (rho) away from the background gradient. All fields are stored in binary format (32-bit little-endian), each with a size of 85.8 GB, yielding a total dataset size of 577 TB. Further details are referenced in the attached README file, and a current list of publications and associated analysis tools are provided at https://stratified-turbulence.github.io/web/.

42 ENGINEERING↗
DOE OSTI2025

SST-TG-P7F4R3200: Decaying Stably-Stratified Turbulence (SST), Initialized Using Taylor-Green Vortices (TG) at Prandtl Number Pr=7, Froude Number Fr=4, Reynolds Number Re=3200

This dataset comprises direct numerical simulations (DNS) of decaying stably-stratified turbulence influenced by a linear background density gradient, initialized using an array of Taylor-Green vortices, extending the Pr=1 simulations performed in [Riley and de Bruyn Kops (2003)](https://doi.org/10.1063/1.1578077). The initial Prandtl, Froude, and Reynolds numbers are (Pr, Fr, Re) = (7, 4, 3200). A total of 15,250 snapshots are recorded at uniform time intervals, each with a spatial resolution of 1280x1280x640 grid points. Four flow variables are associated with each snapshot: the three velocity components (u,v,w) and the perturbed density field (rho) away from the background gradient. All fields are stored in binary format (32-bit little-endian), each with a size of 4 GB, yielding a total dataset size of 244 TB. Further details are referenced in the attached README file, and a current list of publications and associated analysis tools are provided at https://stratified-turbulence.github.io/web/.

42 ENGINEERING↗
DOE OSTI2021

The multi-faceted inverted harmonic oscillator: Chaos and complexity

The harmonic oscillator is the paragon of physical models; conceptually and computationally simple, yet rich enough to teach us about physics on scales that span classical mechanics to quantum field theory. This multifaceted nature extends also to its inverted counterpart, in which the oscillator frequency is analytically continued to pure imaginary values. In this article we probe the inverted harmonic oscillator (IHO) with recently developed quantum chaos diagnostics such as the out-of-time-order correlator (OTOC) and the circuit complexity. In particular, we study the OTOC for the displacement operator of the IHO with and without a non-Gaussian cubic perturbation to explore genuine and quasi scrambling respectively. In addition, we compute the full quantum Lyapunov spectrum for the inverted oscillator, finding a paired structure among the Lyapunov exponents. We also use the Heisenberg group to compute the complexity for the time evolved displacement operator, which displays chaotic behaviour. Finally, we extended our analysis to N-inverted harmonic oscillators to study the behaviour of complexity at the different timescales encoded in dissipation, scrambling and asymptotic regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Contrasting Community Assembly Forces Drive Microbial Structural and Potential Functional Responses to Precipitation in an Incipient Soil System

Microbial communities in incipient soil systems serve as the only biotic force shaping landscape evolution. However, the underlying ecological forces shaping microbial community structure and function are inadequately understood. We used amplicon sequencing to determine microbial taxonomic assembly and metagenome sequencing to evaluate microbial functional assembly in incipient basaltic soil subjected to precipitation. Community composition was stratified with soil depth in the pre-precipitation samples, with surficial communities maintaining their distinct structure and diversity after precipitation, while the deeper soil samples appeared to become more uniform. The structural community assembly remained deterministic in pre- and post-precipitation periods, with homogenous selection being dominant. Metagenome analysis revealed that carbon and nitrogen functional potential was assembled stochastically. Sub-populations putatively involved in the nitrogen cycle and carbon fixation experienced counteracting assembly pressures at the deepest depths, suggesting the communities may functionally assemble to respond to short-term environmental fluctuations and impact the landscape-scale response to perturbations. We propose that contrasting assembly forces impact microbial structure and potential function in an incipient landscape; in situ landscape characteristics (here homogenous parent material) drive community structure assembly, while short-term environmental fluctuations (here precipitation) shape environmental variations that are random in the soil depth profile and drive stochastic sub-population functional dynamics.

16S amplicon sequencing↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
DOE OSTI2022

Rigorous justification for the space–split sensitivity algorithm to compute linear response in Anosov systems

Ruelle gave a formula for linear response of transitive Anosov diffeomorphisms. Recently, practically computable realizations of Ruelle's formula have emerged that potentially enable sensitivity analysis of certain high-dimensional chaotic numerical simulations encountered in the applied sciences. In this paper, we provide full mathematical justification for the convergence of one such efficient computation, the space–split sensitivity, or S3, algorithm. In S3, Ruelle's formula is computed as a sum of two terms obtained by decomposing the perturbation vector field into a coboundary and a remainder that is parallel to the unstable direction. Such a decomposition results in a splitting of Ruelle's formula that is amenable to efficient computation. Here, we prove the existence of the S3 decomposition and the convergence of the computations of both resulting components of Ruelle's formula.

97 MATHEMATICS AND COMPUTING↗