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

Decomposition of the photoabsorption Schumann-Runge continuum in O2

The photoabsorption Schumann-Runge continuum in O2 has been decomposed into the contribution from states of both valence and Rydberg character. The decomposition of the experimental photoabsorption (or high-energy, zero-angle, electron energy-loss) data contains contributions from three valence states, four Rydberg states, and one state of unknown character. Information on the perturbed valence-state potential-energy curves in the Franck-Condon region obtained by this analysis substantiate the recent theoretical conclusions of strong valence-Rydberg perturbations in the 8-12 eV excitation region. There is presently no unambiguous identification for the state at 8.93 plus or minus 0.04 eV, found in the optical and all the electron energy-loss data.

Cartwright, D. C.↗
NASA NTRS2003

Modeling beta-adrenergic control of cardiac myocyte contractility in silico

The beta-adrenergic signaling pathway regulates cardiac myocyte contractility through a combination of feedforward and feedback mechanisms. We used systems analysis to investigate how the components and topology of this signaling network permit neurohormonal control of excitation-contraction coupling in the rat ventricular myocyte. A kinetic model integrating beta-adrenergic signaling with excitation-contraction coupling was formulated, and each subsystem was validated with independent biochemical and physiological measurements. Model analysis was used to investigate quantitatively the effects of specific molecular perturbations. 3-Fold overexpression of adenylyl cyclase in the model allowed an 85% higher rate of cyclic AMP synthesis than an equivalent overexpression of beta 1-adrenergic receptor, and manipulating the affinity of Gs alpha for adenylyl cyclase was a more potent regulator of cyclic AMP production. The model predicted that less than 40% of adenylyl cyclase molecules may be stimulated under maximal receptor activation, and an experimental protocol is suggested for validating this prediction. The model also predicted that the endogenous heat-stable protein kinase inhibitor may enhance basal cyclic AMP buffering by 68% and increasing the apparent Hill coefficient of protein kinase A activation from 1.0 to 2.0. Finally, phosphorylation of the L-type calcium channel and phospholamban were found sufficient to predict the dominant changes in myocyte contractility, including a 2.6x increase in systolic calcium (inotropy) and a 28% decrease in calcium half-relaxation time (lusitropy). By performing systems analysis, the consequences of molecular perturbations in the beta-adrenergic signaling network may be understood within the context of integrative cellular physiology.

Non-NASA Center↗
NASA NTRS1983

Extrapolation of optimum design based on sensitivity derivatives

In many applications, it is useful to have information regarding the effect of the perturbations of the design problem constants (parameters) on the optimum solution of the problem. Sobieszczanski-Sobieski et al. (1981) have, therefore, proposed to utilize an optimum sensitivity analysis as a structural design tool to assess the parameter perturbation effects by extrapolations based on optimum sensitivity derivatives. The present investigation is concerned with an example of optimum sensitivity analysis involving a composite panel. The case was selected in connection with the nonlinearity of the constraints and, the tendency for a set of active constraints to change when the parameters are disturbed. Under these adverse conditions, sensitivity derivatives are used in linear extrapolations to extend optimization information over a limited range of variation for a typical parameter.

Barthelemy, J.-F. M.↗
DOE OSTI2022

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Enhanced simulation software for rocket turbopump, turbulent, annular liquid seals

One of the main objectives of this work is to develop a new dynamic analysis for liquid annular seals with arbitrary profile and to analyze a general distorted interstage seal of the space shuttle main engine high pressure oxygen turbopump (SSME-ATD-HPOTP). The dynamic analysis developed is based on a method originally proposed by Nelson and Nguyen. A simpler scheme based on cubic splines is found to be computationally more efficient and has better convergence properties at higher eccentricities. The first order solution of the original analysis is modified by including a more exact solution that takes into account the variation of perturbed variables along the circumference. A new set of equations for dynamic analysis are derived based on this more general model. A unified solution procedure that is valid for both Moody's and Hirs' friction models is presented. Dynamic analysis is developed for three different models: constant properties, variable properties, and thermal effects with variable properties. Arbitrarily varying seal profiles in both axial and circumferential directions are considered. An example case of an elliptical seal with varying degrees of axial curvature is analyzed in detail. A case study based on predicted clearances of an interstage seal of the SSME-ATD-HPOTP is presented. Dynamic coefficients based on external specified load are introduced to analyze seals that support a preload. The other objective of this work is to study the effect of large rotor displacements of SSME-ATD-HPOTP on the dynamics of the annular seal and the resulting transient motion. One task is to identify the magnitude of motion of the rotor about the centered position and establish limits of effectiveness of using current linear models. This task is accomplished by solving the bulk flow model seal governing equations directly for transient seal forces for any given type of motion, including motion with large eccentricities. Based on the above study, an equivalence is established between linearized coefficients based transient motion and the same motion as predicted by the original governing equations. An innovative method is developed to model nonlinearities in an annular seal based on dynamic coefficients computed at various static eccentricities. This method is thoroughly tested for various types of transient motion using bulk flow model results as a benchmark.

Padavala, Satya↗
NASA NTRS1999

Interaction Between Tropical Convection and its Embedding Environment: An Energetics Analysis of a 2-D Cloud Resolving Simulation

The phase relation between the perturbation kinetic energy (K') associated with the tropical convection and the horizontal-mean moist available potential energy (bar-P) associated with environmental conditions is investigated by an energetics analysis of a numerical experiment. This experiment is performed using a 2-D cloud resolving model forced by the TOGA-COARE derived vertical velocity. The imposed upward motion leads to a decrease of bar-P directly through the associated vertical advective cooling, and to an increase of K' directly through cloud related processes, feeding the convection. The maximum K' and its maximum growth rate lags and leads, respectively, the maximum imposed large-scale upward motion by about 1-2 hours, indicating that convection is phase locked with large-scale forcing. The dominant life cycle of the simulated convection is about 9 hours, whereas the time scales of the imposed large-scale forcing are longer than the diurnal cycle. In the convective events, maximum growth of K' leads maximum decay of the perturbation moist available potential energy (P') by about 3 hours through vertical heat transport by perturbation circulation, and perturbation cloud heating. Maximum decay of P' leads maximum decay of bar-P by about one hour through the perturbation radiative, processes, the horizontal-mean cloud heating, and the large-scale vertical advective cooling. Therefore, maximum gain of K' occurs about 4-5 hours before maximum decay of bar-P.

Li, Xiaofan↗
DOE OSTI2025

Post-inflationary enhancement of adiabatic perturbations in modular cosmology

We show that multi-field inflationary models with negligible turning in field space during inflation can lead to an effective sourcing of adiabatic from entropic perturbations afterthe end of inflation. We illustrate this general phenomenon with a detailed analysis of an inflationary model whose scalar potential is determined by modular invariance. Its entropic perturbations are frozen during inflation, but instead, they are converted into adiabatic perturbations in the first post-inflationary e-folds. The curvature power spectrum, giving rise to CMB fluctuations, reaches a novel and enhanced plateau in this process; we address the implications for the inflationary observables As , n s and r.

cosmological perturbation theory↗
DOE OSTI2025

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗
DOE OSTI2025

DECOVALEX-2023: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX initiative is an international research collaboration (www.decovalex.org), initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of the complex perturbations and coupled processes in geologic repositories and how these impact long-term performance predictions. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2023, started in early 2020 and ended in late 2023. More than fifty research teams associated with 17 international DECOVALEX partner organizations participated in the comparative evaluation of eight modeling tasks covering a wide range of spatial and temporal scales, geological formations, and coupled processes. This Virtual Special Issue on DECOVALEX-2023 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes. While primarily focused on radioactive waste, much of the work included here has wider application to many geoengineering topics.

Coupled processes↗
DOE OSTI2021

The Low-Lying Electronic States of NO 2 : Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum

Nitrogen dioxide, NO 2 , is a free radical composed of the two most abundant elements in Earth’s atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner–Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N( 4 S u ) and O( 3 P g ) atoms. In this work, three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four ($\tilde X$, $\tilde A$, $\tilde B$, and $\tilde C$) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground $\tilde X$ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm –1 above the zero-point energy) have a root-mean-squared error of 16.5 cm –1 . In addition, dipole moments for each of the lowest four electronic states—and the transition dipoles between them—were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ~1.5% with no empirical adjustments.

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

Do-calculus enables estimation of causal effects in partially observed biomolecular pathways

Abstract Motivation Estimating causal queries, such as changes in protein abundance in response to a perturbation, is a fundamental task in the analysis of biomolecular pathways. The estimation requires experimental measurements on the pathway components. However, in practice many pathway components are left unobserved (latent) because they are either unknown, or difficult to measure. Latent variable models (LVMs) are well-suited for such estimation. Unfortunately, LVM-based estimation of causal queries can be inaccurate when parameters of the latent variables are not uniquely identified, or when the number of latent variables is misspecified. This has limited the use of LVMs for causal inference in biomolecular pathways. Results In this article, we propose a general and practical approach for LVM-based estimation of causal queries. We prove that, despite the challenges above, LVM-based estimators of causal queries are accurate if the queries are identifiable according to Pearl’s do-calculus and describe an algorithm for its estimation. We illustrate the breadth and the practical utility of this approach for estimating causal queries in four synthetic and two experimental case studies, where structures of biomolecular pathways challenge the existing methods for causal query estimation. Availability and implementation The code and the data documenting all the case studies are available at https://github.com/srtaheri/LVMwithDoCalculus. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Representation Learning via Quantum Neural Tangent Kernels

Variational quantum circuits are used in quantum machine learning and variational quantum simulation tasks. Designing good variational circuits or predicting how well they perform for given learning or optimization tasks is still unclear. Here we discuss these problems, analyzing variational quantum circuits using the theory of neural tangent kernels. We define quantum neural tangent kernels, and derive dynamical equations for their associated loss function in optimization and learning tasks. We analytically solve the dynamics in the frozen limit, or lazy training regime, where variational angles change slowly and a linear perturbation is good enough. We extend the analysis to a dynamical setting, including quadratic corrections in the variational angles. We then consider a hybrid quantum classical architecture and define a large-width limit for hybrid kernels, showing that a hybrid quantum classical neural network can be approximately Gaussian. The results presented here show limits for which analytical understandings of the training dynamics for variational quantum circuits, used for quantum machine learning and optimization problems, are possible. These analytical results are supported by numerical simulations of quantum machine-learning experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

CSPlib

This computational singular perturbation library (CSPlib) software is useful for analysis of chemical kinetic systems. The code is used to analyze and describe the dynamical landscape underlying the time evolution of a chemical system. Given a database of state vectors, it allows the analysis of the source term and Jacobian matrix of the chemical system, the identification and decoupling of fast and slow processes, and the analysis of cause-and-effect relationships in the coupled system dynamics. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-2391 O

Najm, HabibN.↗
DOE OSTI2022

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
NASA NTRS1964

Guidance analysis.

Unmanned spacecraft guidance policy and orbit determination, using linear perturbation theory and minimum variance estimation

SPACECRAFT GUIDANCE↗
NASA NTRS1974

Ephemeris data and error analysis in support of a Comet Encke intercept mission

Utilizing an orbit determination based upon 65 observations over the 1961 - 1973 interval, ephemeris data were generated for the 1976-77, 1980-81 and 1983-84 apparitions of short period comet Encke. For the 1980-81 apparition, results from a statistical error analysis are outlined. All ephemeris and error analysis computations include the effects of planetary perturbations as well as the nongravitational accelerations introduced by the outgassing cometary nucleus. In 1980, excellent observing conditions and a close approach of comet Encke to the earth permit relatively small uncertainties in the cometary position errors and provide an excellent opportunity for a close flyby of a physically interesting comet.

Yeomans, D. K.↗
NASA NTRS1974

Tidal parameters from the variation of inclination of GEOS-1 and GEOS-2

Analysis of the luni-solar tidal perturbations of the inclination of GEOS-1 and GEOS-2 has yielded the values 0.22 and 0.31 respectively for the apparent second degree Love number. For GEOS-1 a new purely numerical method involving osculating elements was employed. For GEOS-2 it was necessary to analyze the variations of the mean elements because of the very long period (450 days) of the dominant solar tidal perturbation. The disparate values indicate that the simple second degree zonal harmonic model of the tidal potential is accommodating other effects in addition to those caused by the solid earth tides. A recent paper by Lambeck et al. (1973) indicates that ocean tide effects have significant perturbations on satellite orbits and cannot be neglected.

Douglas, B. C.↗