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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development and validation of a mathematical model of heart rate response to fluid perturbation

Physiological closed-loop controlled (PCLC) medical devices monitor and automatically adjust the patient’s condition by using physiological variables as feedback, ideally with minimal human intervention to achieve the target levels set by a clinician. PCLC devices present a challenge when it comes to evaluating their performance, where conducting large clinical trials can be expensive. Virtual physiological patients simulated by validated mathematical models can be utilized to obtain preclinical evidence of safety and assess the performance of the PCLC medical device during normal and worst-case conditions that are unlikely to happen in a limited clinical trial. A physiological variable that plays a major role during fluid resuscitation is heart rate (HR). For in silico assessment of PCLC medical devices regarding fluid perturbation, there is currently no mathematical model of HR validated in terms of its predictive capability performance. This paper develops and validates a mathematical model of HR response using data collected from sheep subjects undergoing hemorrhage and fluid infusion. The model proved to be accurate in estimating the HR response to fluid perturbation, where averaged between 21 calibration datasets, the fitting performance showed a normalized root mean square error (NRMSE) of 7.41 ± 2.8%. The model was also evaluated in terms of model predictive capability performance via a leave-one-out procedure (21 subjects) and an independent validation dataset (6 subjects). Two different virtual cohort generation tools were used in each validation analysis. The generated envelope of virtual subjects robustly met the defined acceptance criteria, in which 95% of the testing datasets presented simulated HR patterns that were within a deviation of 50% from the observed data. In addition, out of 16000 and 18522 simulated subjects for the leave one-out and independent datasets, the model was able to generate at least one virtual subject that was close to the real subject within an error margin of 9.56 ± 3.15% and 11.1 ± 1.22% NRMSE, respectively. In conclusion, the model can generate valid virtual HR physiological responses to fluid perturbation and be incorporated into future non-clinical simulated testing setups for assessing PCLC devices intended for fluid resuscitation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Constraints on quasinormal modes and bounds for critical points from pole-skipping

We consider a holographic thermal state and perturb it by a scalar operator whose associated real-time Green’s function has only gapped poles. These gapped poles correspond to the non-hydrodynamic quasinormal modes of a massive scalar perturbation around a Schwarzschild black brane. Relations between pole-skipping points, critical points and quasinormal modes in general emerge when the mass of the scalar and hence the dual operator dimension is varied. First, this novel analysis reveals a relation between the location of a mode in the infinite tower of quasinormal modes and the number of pole-skipping points constraining its dispersion relation at imaginary momenta. Second, for the first time, we consider the radii of convergence of the derivative expansions about the gapped quasinormal modes. These convergence radii turn out to be bounded from above by the set of all pole-skipping points. Furthermore, a transition between two distinct classes of critical points occurs at a particular value for the conformal dimension, implying close relations between critical points and pole-skipping points in one of those two classes. We show numerically that all of our results are also true for gapped modes of vector and tensor operators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Instability seeding mechanisms due to internal defects in inertial confinement fusion targets

Performance degradation in laser-driven inertial confinement fusion (ICF) implosions is caused by several effects, one of which is Rayleigh–Taylor instability growth. Defects in ICF targets, such as internal voids and surface roughness, create instability seeds in the shell as shocks propagate through the target. A comprehensive understanding of seeding mechanisms is essential to characterize the impact of target defects on inflight shell integrity and mass injection into the central, lower-density vapor region. An analysis of early-time behavior of both single-mode shell mass modulations and isolated voids is performed by examining the evolution of the acoustic waves launched by these target imperfections. A systematic study of localized perturbation growth as a function of defect placement and size is presented. Furthermore, the use of low-density ablator materials (such as foams) is suggested as a potential mitigation strategy to improve target robustness against the impact of defect-initiated growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inverse design for waveguide dispersion with a differentiable mode solver

Inverse design of optical components based on adjoint sensitivity analysis has the potential to address the most challenging photonic engineering problems. However, existing inverse design tools based on finite-difference-time-domain (FDTD) models are poorly suited for optimizing waveguide modes for adiabatic transformation or perturbative coupling, which lies at the heart of many important photonic devices. Among these, dispersion engineering of optical waveguides is especially challenging in ultrafast and nonlinear optical applications involving broad optical bandwidths and frequency-dependent anisotropic dielectric material response. In this work, we develop gradient back-propagation through a general-purpose electromagnetic eigenmode solver and use it to demonstrate waveguide dispersion optimization for second harmonic generation with maximized phase-matching bandwidth. This optimization of three design parameters converges in eight steps, reducing the computational cost of optimization by ∼100x compared to exhaustive search and identifying new designs for broadband optical frequency doubling of laser sources in the 1.3–1.4 µm wavelength range. Furthermore, we demonstrate that the computational cost of gradient back-propagation is independent of the number of parameters, as required for optimization of complex geometries. This technique enables practical inverse design for a broad range of previously intractable photonic devices.

Gray, Dodd (ORCID:000000030469599X)↗
DOE OSTI2022

Influence of Polar Mutations on the Electronic and Structural Properties of Q A –· in Bacterial Reaction Centers

Reaction centers from Rb. sphaeroides with residue M265 mutated from isoleucine to threonine, serine, and asparagine (M265IT, M265IS, and M265IN, respectively) in the Q A –· state are studied by high-resolution ESEEM and ENDOR spectroscopies to investigate the structural characteristics of these mutants influencing the redox properties of the Q A site. All three mutants decrease the redox midpoint potential (E m ) of Q A by ~0.1 V, yet the mechanism for this drop in E m is unclear. In this work, we examine (i) the hydrogen bonding interactions between Q A –· and residues histidine M219 and alanine M260, (ii) the electron spin density distribution of the semiquinone, and (iii) the orientations of the ubiquinone methoxy substituents. 13 C measurements show no significant contribution of methoxy dihedral angles to the observed decrease in E m for the Q A mutants. Instead, 14 N three-pulse ESEEM data suggest electrostatic or hydrogen bond formation between the mutated M265 side chain and His-M219 N δ may be involved in the observed lowering of the Q A midpoint potential. For mutant M265IN, analysis of the proton hyperfine couplings reveals a weakened hydrogen bond network, resulting in an altered Q A –· spin density distribution. Here, the magnetic resonance study presented here is therefore most consistent with an electrostatic or structural perturbation of the His-M219 N δ hydrogen bond in these mutants as a mechanism for the ~0.1 V decrease in Q A E m .

14 SOLAR ENERGY↗
DOE OSTI2022

A deep redox proteome profiling workflow and its application to skeletal muscle of a Duchenne Muscular Dystrophy model

Perturbation to the redox state accompanies many diseases and its effects are viewed through oxidation of biomolecules, including proteins, lipids, and nucleic acids. The thiol groups of protein cysteine residues undergo an array of redox post-translational modifications (PTMs) that are important for regulation of protein and pathway function. To better understand what proteins are redox regulated following a perturbation, it is important to be able to comprehensively profile protein thiol oxidation at the proteome level. Herein, we report a deep redox proteome profiling workflow and demonstrate its application in measuring the changes in thiol oxidation along with global protein expression in skeletal muscle from mdx mice, a model of Duchenne Muscular Dystrophy (DMD). In depth coverage of the thiol proteome was achieved with >18,000 Cys sites from 5608 proteins in muscle being quantified. Compared to the control group, mdx mice exhibit markedly increased thiol oxidation, where ~2% shift in the median oxidation occupancy was observed. Further, pathway analysis for the redox data revealed that coagulation system and immune-related pathways were among the most susceptible to increased thiol oxidation in mdx mice, whereas protein abundance changes were more enriched in pathways associated with bioenergetics. This study illustrates the importance of deep redox profiling in gaining greater insight into oxidative stress regulation and pathways/processes that are perturbed in an oxidizing environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Grassmannian Shape Representations for Aerodynamic Applications

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗
DOE OSTI2022

Skewing the CMB×LSS: a fast method for bispectrum analysis

Upcoming cosmic microwave background (CMB) lensing measurements and tomographic galaxy surveys are expected to provide us with high-precision data sets in the coming years, thus paving the way for fruitful cross-correlation analyses. Here, in this paper we study the information content of the weighted skew-spectrum, a nearly-optimal estimator of the angular bispectrum amplitude, as a means to extract non-Gaussian information on both bias and cosmological parameters from the bispectra of galaxies cross-correlated with CMB lensing, while gaining significantly on speed. Our results show that for the combination of the Planck satellite and the Dark Energy Spectroscopic Instrument (DESI), the difference in the constraints on bias and cosmological parameters from the skew-spectrum and the bispectrum is at most 17%. We further compare and find agreement between our theoretical skew-spectra and those estimated from N-body simulations, for which it is important to include gravitational non-linearities beyond perturbation theory and the post-Born effect for CMB lensing. We define an algorithm to apply the skew-spectrum estimator to the data and, as a preliminary step, we use the skew-spectra to constrain bias parameters and the amplitude of shot noise from the simulations through a Markov chain Monte Carlo likelihood analysis, finding that it may be possible to reach percent-level estimates for the linear bias parameter b 1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Resonant Raman in armchair graphene nanoribbons from first-principles

Resonant Raman spectra of armchair graphene nanoribbons (AGNRs) are computed using Density Functional Theory (DFT) and third-order perturbation theory. Results are benchmarked against available experimental data and compared to previously used theoretical approaches based on the Placzek approximation. Comparable agreement with experiments is found for both previously and presently used methods. In addition, a numerical analysis is carried out to provide a justification for the resonant modeling method based on the use of the frequency-dependent dielectric tensor in the Placzek approximation. Finally, this work also provides additional predictions and references for wide AGNRs that might be investigated with Raman scattering experiments in the future.

42 ENGINEERING↗
DOE OSTI2020

Why matter effects matter for JUNO

In this paper we focus on the Earth matter effects for the solar parameter determination by a medium baseline reactor experiment such as JUNO. We derive perturbative expansions for the mixing angles θ 12 and θ 13 as well as the Δm$^{2}_{21}$ and Δm$^{2}_{31}$ in terms of the matter potential relevant for JUNO. These expansions, up to second order in the matter potential, while simple, allow one to calculate the electron antineutrino survival probability to a precision much better than needed for the JUNO experiment. We use these perturbative expansions to semi-analytically explain and confirm the shift caused by the matter effects on the solar neutrino mixing parameters θ 12 and Δm$^{2}_{21}$ which were previously obtained by a purely numerical χ 2 analysis. Since these shifts do not satisfy the naive expectations and are significant given the precision that can be achieved by the JUNO experiment, a totally independent cross check using a completely different method is of particular importance. We find that these matter effect shifts do not depend on any of the details of the detector characteristics apart from the baseline and earth mass density between reactor(s) and detector, but do depend on the normalized product of reactor neutrino spectrum times the inverse-beta decay cross-section. The results of this manuscript suggests an alternative analysis method for measuring sin 2 θ 12 and Δm$^{2}_{21}$ in JUNO which would be a useful cross check of the standard analysis and for the understanding of the Wolfenstein matter effect. The explanation of these shifts together with a quantitative understanding, using a semi-analytical method, is the principal purpose of this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Revisiting Interior Water Mass Responses to Surface Forcing Changes and the Subsequent Effects on Overturning in the Southern Ocean

Abstract Two coupled climate models, differing primarily in horizontal resolution and treatment of mesoscale eddies, were used to assess the impact of perturbations in wind stress and Antarctic ice sheet (AIS) melting on the Southern Ocean meridional overturning circulation (SO MOC), which plays an important role in global climate regulation. The largest impact is found in the SO MOC lower limb, associated with the formation of Antarctic Bottom Water (AABW), which in both models is enhanced by wind and weakened by AIS meltwater perturbations. Even though both models under the AIS melting perturbation show similar AABW transport reductions of 4–5 Sv (50%–60%), the volume deflation of AABW south of 30°S is four times greater in the higher resolution simulation (−20 vs. −5 Sv). Water mass transformation (WMT) analysis reveals that surface‐forced dense water formation on the Antarctic shelf is absent in the higher resolution and reduced by half in the lower resolution model in response to the increased AIS melting. However, the decline of the AABW volume (and its inter‐model difference) far exceeds the surface‐forced WMT changes alone, which indicates that the divergent model responses arise from interactions between changes in surface forcing and interior mixing processes. This model divergence demonstrates an important source of uncertainty in climate modeling, and indicates that accurate shelf processes together with scenarios accounting for AIS melting are necessary for robust projections of the deep ocean's response to anthropogenic forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Black hole perturbations and quasinormal modes in hybrid metric-Palatini gravity

The rapid advancement of gravitational wave astronomy in recent years has paved the way for the burgeoning development of black hole spectroscopy, which enhances the possibility of testing black holes by their quasinormal modes (QNMs). In this paper, the axial gravitational perturbations and the QNM frequencies of black holes in the hybrid metric-Palatini gravity (HMPG) are investigated. Additionally, the HMPG theory is characterized by a dynamical scalar degree of freedom and is able to explain the late-time accelerating expansion of the universe without introducing any ad hoc screening mechanism to preserve the dynamics at the Solar System scale. We obtain the master equation governing the axial gravitational perturbations of the HMPG black holes and calculate the QNM frequencies. Moreover, in the scrutiny of the black holes and their QNMs, we take into account the constraints on the model parameters based on the post-Newtonian analysis, and show how the QNM frequencies of the HMPG black holes would be altered in the observationally consistent range of parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Picometer sensitivity metrology for EUV absorber phase

With growing interest in EUV attenuated phase-shift masks due to their superior image quality for applications such as dense contact and pillar arrays, it is becoming critical to model, measure, and monitor the intensity and relative phase of multilayer and absorber reflections. We present a solution based on physical modeling of reflectometry data, which can achieve single picometer phase precision and sensitivity to changes in average film thickness below one atomic monolayer. We measure absorber and multilayer reflectivity to determine thin-film parameters with a multidimensional optimization and then acquire a new measurement of either multilayer or absorber to determine perturbations in surface contamination thickness. While it is difficult to assess the accuracy of the first step, the simplicity of the second step allows us to characterize our sensitivity to changes in contamination thickness. Thus, we apply this analysis using an initial set of measurements and repeated measurements after a period of storage. For the multilayer, the total contamination growth was 1068 pm, which occurred almost exclusively during storage (1085 pm) and decreased very slightly during repeated measurements (-17 pm). For the absorber, the behavior was quite different, with a total growth of 126 pm, which occurred much less during storage (28 pm) and primarily during repeated measurements (98 pm). Ultimately, the change in relative phase (absorber minus multilayer) was -0.86 deg for the multilayer and -1.12 deg for the absorber. We estimate the precision of the surface contamination measurement to be 3σ < 6 pm for measuring thickness and 3σ < 0.2 deg for measuring phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GradientGraph

Under this SBIR Phase II, Reservoir Labs has developed G2 Analytics, a new technology that allows network operators to analyze bottleneck and flow performance with high precision. G2 delivers a new analytical approach and framework to resolve a variety of key problems found in modern communication networks, including: traffic engineering, routing, flow scheduling, network design, capacity planning, resiliency analysis, network slicing, or service level agreement (SLA) management, among others. G2 leverages the bottleneck structure of congestion-controlled communication networks, a recent mathematical discovery by the Reservoir team [RL19b, RL20a, RL20b, RL21a]. Bottleneck structures reveal how perturbations on flows and links propagate through the network, providing an analytical framework to measure (qualitatively and quantitatively) the ripple effects induced as they traverse the network. Leveraging the mathematics of bottleneck structures, Reservoir Labs is developing the G2 technology to provide network operators with a framework to design, optimize and troubleshoot network performance. This delivery includes the G2 software stack.

Yellamraju, Sruthi↗
DOE OSTI2025

Using "AI Poincare" to analyze non-linear integrable optics

This study dives into the applicability of using automated discovery of conserved quantities in dynamical systems relevant to accelerator physics. Specifically, we explore the performance of AI Poincaré in analyzing numerical trajectory data obtained using the McMillan system of non-linear integrable optics. A comprehensive evaluation of the algorithm's performance is conducted through diverse methodologies. These include the analysis of the estimated number of conserved quantities embedded in a dataset and the deviation of interpolated points on the inferred manifold with respect to points in actually in the dataset. the investigation identifies an optimal range of perturbation distances where the underlying manifold extraction algorithm inside AI Poincaré exhibits optimal performance. Additionally, an improved neural network architecture is proposed based on the observed results. Finally, we apply the algorithm to preliminary experimental data from the Integrable Optics Test Accelerator at Fermilab to successfully infer the number of conserved quantities even in the presence of fast decoherence of the measured signal.

Osmanov, Lazare [Free U. Tbilisi]↗
DOE OSTI2022

A Strong-QCD Regime Measurement of the Proton’s Spin Structure

The theory of the strong force, Quantum Chromodynamics (QCD) remains one of the most important ways to understand the fundamental properties of ordinary matter. However, at low momentum transfer Q 2 , in the regime where the strong force becomes extremely strong, our understanding of QCD for ordinary nucleons becomes hazy. Several cutting edge theories such as Chiral Perturbation Theory (χPT) and Lattice QCD have provided valuable predictions in this regime, but Lattice QCD has not yet extended predictions of many important quantities to this kinematic region, and Chiral Perturbation Theory has faced several important disagreements with experimental data for the neutron over the last several decades. It is therefore of extreme importance to have a benchmark of experimental data in the low energy regime for the proton’s behavior, as a test of leading theories for the behavior of QCD in this regime. The E08-027 (g2p) experiment ran at Jefferson Lab in 2012 with the goal of collecting this valuable data, and though I was still completing my undergraduate studies at the time, I became involved in the analysis in 2015 and built on the previous work to complete it and analyze the exciting results. This experiment achieved a high precision measurement of the spin structure functions g1 and g2 for the proton, quantities which describe the internal spin structure of the proton. These measurements were taken in the valuable low Q 2 region described above, and used to extract several moments of these spin structure functions which can be directly compared to the cutting-edge predictions of Chiral Perturbation Theory. Though the experiment’s timeline was such that I didn’t have a chance to work directly on the experimental setup, I had the opportunity to acquire hands-on experience working on a polarized target at UNH which is very similar to the crucial polarized target used in the g2p experiment. Full details of the g2p experiment and my experimental work at the University of New Hampshire are presented in this thesis, as well as a detailed description of the analysis process and the exciting benchmark results, which serve as a direct test of all current and future theories of QCD in the low-Q 2 regime.

Ruth, David↗