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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Non-Hamiltonian perturbations of integrable systems and resonance trapping

This paper studies general, non-Hamiltonian perturbations of integrable systems with two degrees of freedom and derives conditions for temporary and permanent resonance trapping. The analysis involves a noncanonical transformation of variables near the resonant manifold and averaging with respect to the fast phase to investigate oscillatory behavior on the intermediate timescale. The resulting reduced system is Hamiltonian to leading order and permits, after averaging on the intermediate, or libration, timescale, a canonical transformation to action-angle variables in the oscillation zone. The final system so obtained reveals the possible existence of two- and three-dimensional invariant tori in the vicinity of the resonant manifold. An explicit divergence condition for general perturbations to be dissipative on the slow timescale follows from the analysis. An application of this approach to the problem of resonant trapping and escape is outlined for the restricted problem of three bodies subject to dissipative perturbations with a radial symmetry.

Ghil, M.↗
DOE OSTI2024

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The physics of spontaneous parity-time symmetry breaking in the Kelvin–Helmholtz instability

Abstract We show that the dynamical system of an inviscid fluid with velocity shear admits parity-time (PT) symmetry, which provides a physical explanation of the well-known observation that the spectrum of the perturbation eigenmodes of the system is symmetric with respect to the real axis. It is found that the Kelvin–Helmholtz instability is triggered when and only when the PT symmetry is spontaneously broken. The analysis of PT symmetry also reveals that the relative phase between parallel velocity and pressure perturbations needs to be locked at π /2 when the instability is suppressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Nonlinear instability of accelerating shock waves with application to supernovae

We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.

Luo, Ding↗
DOE OSTI2020

NMDQi Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST)

The Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST) is developed to support the thermal-hydraulic aspects for the design of optimal irradiation vehicles for the Nuclear Materials Discovery and Qualification initiative (NMDQi). TOAST is currently a MATLAB-based tool that utilizes a simplified steady-state analytical model where a radial thermal resistance circuit is utilized to account for conduction through up to 2 capsules, a gas gap, a thermal bond, a sample cladding, and a sample, as well as the convection from a water coolant outside the capsules. The geometry is discretized axially along the user-defined height of a basket with a user-defined number of points/nodes, essentially turning TOAST into a semi-2D heat transfer analysis utility for cylindrical irradiation vehicles with practically any configuration, solving for temperatures radially at multiple axial locations. TOAST utilizes a Graphical User Interface (GUI) in-which a user can define the geometrical layout, the materials utilized for each component, the coolant characteristics, and the required solution. The user can define the geometry by inputting the diameters, thicknesses, and heights for each component, whereas the materials are defined via the user-provided constant or variable thermal conductivities. The user can select the coolant's inlet temperature, pressure, and the pressure drop across the height of the problem (which is used to calculate the velocity of the flow). Finally, the user can choose to solve for a sample heat generation rate limit by inputting temperature limits for the sample and capsules or can choose to purely solve for the axial temperature distributions of each component by inputting a pre-selected heat generation rate. Either way, TOAST provides the user with axial temperature distributions of all the components including the coolant, and the heat generation rate limit as well as the maximum outlet coolant temperature. The user can also choose to do one of 6 sensitivity analyses in TOAST. The sensitivity analyses yield plots of the sensitivity of the heat generation rate limit, the maximum coolant temperature, and the pressure limit for the annular components due to perturbations in 1-2 unknown variables based on the selected sensitivity analysis. Benchmarks between computations in TOAST and equivalent 2D axis symmetric ABAQUS models are presented, showing that TOAST results are within less than 3% of ABAQUS results in most cases, and a maximum of 8% difference in some cases. The benchmarks also revealed that this uncertainty is tied to the selection of an appropriate Nusselt number correlation and appropriate thermal conductivities, which the user can do from the GUI. Regardless, TOAST is demonstrated as a computationally efficient, highly accessible, and accurate utility for optimization and scoping s of studies different irradiation vehicles.

42 ENGINEERING↗
DOE OSTI2020

Crossing the Streams – Sampler and the TemplateEngine [Slides]

This presentation discusses Sampler, which is a versatile UQ and parametric study tool that can be applied to any SCALE Sequence. Sampler can perturb any quantity in any SCALE input. Recent work at ORNL has developed new types of covariance data that allow Sampler UQ to be applied to nearly all SCALE applications, including reactor depletion, UNF fuel characterization, source term analysis, and decay heat calculation. In SCALE 6.2 releases, CE data in transport cannot be perturbed. Sampler was originally designed for stochastic sampling with any sequence within SCALE and Parametric capability added in SCALE 6.2.2. Sampler can be used for uncertainty quantification, including sample data in static or depletion calculations and sample inputs for uncertainties in compositions and dimensions. The SCALE TemplateEngine allows for expanding templates to full inputs and the combination provides a powerful UQ tool.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Synthesis of a modern/classical controller for large space structures

The paper is concerned with wide-sense stationary random vibrations of a representative large space structure subject to nonwhite excitations, and with synthesis of its controller amid modeling uncertainties due to truncated modes. Dynamics is analyzed by residue calculus in frequency domain. Steady state variance of a modal coordinate is shown to depend on, among other things, power spectral density of excitation at the modal frequency. Modal cost analysis is performed for an order reduction with a performance objective of attitude, attitude rate and energy control. Optimal and suboptimal output feedback controllers are compared. By treating deleted modes as an additive and a multiplicative perturbation in an idealized transfer function, the low frequency modes omitted according to modal cost analysis are found to cause instability. Finally, the size of the idealized model is shown to be governed by light damping and stringency in performance and stability specifications.

Hablani, H. B.↗
NASA NTRS1979

Origin of longitudinal grooving in Tiu Vallis, Mars - Isolation of responsible fluid-types

Fluid instability analysis developed for the growth or decay of large, longitudinal roller vortex perturbations imposed at the base of fluids of exponentially stratified viscosities is presented. The origin of longitudinal grooving and fluting in a region of converging flow in Tiu Vallis, Mars, is analyzed, and the propensity of fluids, such as water, mudflow/debris flow, and lava, to develop longitudinal roll patterns as part of their flow field is tested. It is found that the grooving pattern in Tiu Vallis was probably formed by vortex motion in a mudflow or debris-laden slurry, but exact representation of sediment/debris concentration or sediment transport machanism is not yet possible. It is concluded that longitudinal vortex motion modelled in this analysis can develop in water, lava, or ice, but that the specific flow characteristics of these models which include the vortex motion are unreasonable for flow in Tiu Vallis.

Thompson, D. E.↗
DOE OSTI2023

Resurgence, conformal blocks, and the sum over geometries in quantum gravity

In two dimensional conformal field theories the limit of large central charge plays the role of a semi-classical limit. Certain universal observables, such as conformal blocks involving the exchange of the identity operator, can be expanded around this classical limit in powers of the central charge c. This expansion is an asymptotic series, so — via the same resurgence analysis familiar from quantum mechanics — necessitates the existence of non-perturbative effects. In the case of identity conformal blocks, these new effects have a simple interpretation: the CFT must possess new primary operators with dimension of order the central charge. This constrains the data of CFTs with large central charge in a way that is similar to (but distinct from) the conformal bootstrap. We study this phenomenon in three ways: numerically, analytically using Zamolodchikov’s recursion relations, and by considering non-unitary minimal models with large (negative) central charge. In the holographic dual to a CFT2, the expansion in powers of c is the perturbative loop expansion in powers of ћ. So our results imply that the graviton loop expansion is an asymptotic series, whose cure requires the inclusion of new saddle points in the gravitational path integral. In certain cases these saddle points have a simple interpretation: they are conical excesses, particle-like states with negative mass which are not in the physical spectrum but nevertheless appear as non-manifold saddle points that control the asymptotic behaviour of the loop expansion. This phenomenon also has an interpretation in SL(2, R) Chern-Simons theory, where the non-perturbative effects are associated with the non-Teichmüller component of the moduli space of flat connections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗
DOE OSTI2022

Vacancy-engineered nodal-line semimetals

Symmetry-enforced nodal-line semimetals are immune to perturbations that preserve the underlying symmetries. This intrinsic robustness enables investigations of fundamental phenomena and applications utilizing diverse materials design techniques. The drawback of symmetry-enforced nodal-line semimetals is that the crossings of energy bands are constrained to symmetry-invariant momenta in the Brillouin zone. On the other end are accidental nodal-line semimetals whose band crossings, not being enforced by symmetry, are easily destroyed by perturbations. Some accidental nodal-line semimetals have, however, the advantage that their band crossings can occur in generic locations in the Brillouin zone, and thus can be repositioned to tailor material properties. We show that lattice engineering with periodic distributions of vacancies yields a hybrid type of nodal-line semimetals which possess symmetry-enforced nodal lines and accidental nodal lines, with the latter endowed with an enhanced robustness to perturbations. Both types of nodal lines are explained by a symmetry analysis of an effective model which captures the relevant characteristics of the proposed materials, and are verified by first-principles calculations of vacancy-engineered borophene polymorphs. Our findings offer an alternative path to relying on complicated compounds to design robust nodal-line semimetals; one can instead remove atoms from a common monoatomic material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-resolution infrared spectroscopy of supersonically cooled singlet carbenes: Bromomethylene (HCBr) in the CH stretch region

First high-resolution spectra of cold (~35 K) singlet bromomethylene HCBr in the CH stretching (v 1 ) region from 2770 to 2850 cm-1 are reported using near quantum shot-noise limited laser absorption methods in a slit jet supersonic discharge expansion source. Three rovibrational bands are identified at high S/N (20:1 - 40:1) and rotationally assigned to i) the CH stretch fundamental (v 1 ) band $\tilde{X}$(1,0,0) ← $\tilde{X}$ (0,0,0) and ii) vibrational hot bands ($\tilde{X}$(1,1,0) ← $\tilde{X}$(0,1,0) and $\tilde{X}$(1,0,1) ← $\tilde{X}$(0,0,1)) arising from vibrationally excited HCBr populated in the discharge with single quanta in either the H–C–Br bend (v 2 ) or C–Br stretch (v 3 ) modes. Precision rotational constants are reported for a total of six states, with an experimentally determined CH stretch vibrational frequency (2799.38 cm -1 ) in good agreement with previous low-resolution fluorescence studies. Detailed analysis of the fundamental v 1 band highlights the presence of perturbations in the $\tilde{X}$(1,0,0) level, which we tentatively attribute to arise from the nearby triplet state $\tilde{a}$(0,0,1) through spin-orbit interaction or the multiple quanta $\tilde{X}$(0,2,1) singlet state via c-type Coriolis coupling. Reduced-Doppler resolution (60 MHz) in the slit-jet IR spectrometer permits clear observation of nuclear spin hyperfine structure, with experimental line shapes well reproduced by nuclear quadrupole/spin-rotation coupling constants from microwave studies. Lastly, the a-type to b-type transition intensity ratio for the fundamental CH stretch band is notably larger than predicted by a bond-dipole model, which from high level ab initio quantum calculations (CCSD(T)/PVQZ) can be attributed to vibrationally induced “charge-sloshing” of electron density along the polar C–Br bond.

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

Theory of mean E × B shear in a stochastic magnetic field: ambipolarity breaking and radial current

The mean E $ \times $ B shear in a stochastic magnetic field is calculated, using the radial force balance relation and transport equations. This analysis is relevant to the L → H transition with resonant magnetic perturbations, and special focus is placed upon the physics of non-ambipolar transport and radial current. The key physical process is the flow of fluctuating currents along wandering magnetic fields. The increments in poloidal and toroidal rotation, density and ion pressure are calculated. The radial envelope of the magnetic perturbations inside the plasma defines a new scale ${\ell _{{\text{env}}}}$, which is the characteristic scale of the magnetic fluctuation intensity profile. Here, the net particle outflow due to stochastic magnetic fields is calculated and is determined by the net radial current through the separatrix. Implications for the L → H transition are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Improving Robustness of Spectrogram Classifiers with Neural Stochastic Differential Equations

Signal analysis and classification is fraught with high levels of noise and perturbation. Computer-vision-based deep learning models applied to spectrograms have proven useful in the field of signal classification and detection; however, these methods aren't designed to handle the low signal-to-noise ratios inherent within non-vision signal processing tasks. While they are powerful, they are currently not the method of choice in the inherently noisy and dynamic critical infrastructure domain, such as smart-grid sensing, anomaly detection, and non-intrusive load monitoring. Currently, these models can be brittle, which makes them susceptible to noisy input. This also means they have sub-optimal stability of explanation outputs. Experts and technicians using these models to make decisions in real world scenarios need assurance that a model is performing as it is supposed to. The classification or prediction outputs it generates should be sound and grounded, not likely to change in the presence of shifting noise landscapes. In this work, we explore the idea of Neural Stochastic Differential Equations (NSDE's) to improve the robustness of models trained to classify time series data and the effect of NSDE's on the explainability of outputs. We then test the effectiveness of these approaches by applying them to a non-intrusive load monitoring (NILM) dataset that consists of simulated harmonic signals injected into a real building.

Brogan, Joel↗
DOE OSTI2020

Machine Learning Characterization of Alfvénic and Sub-Alfvénic Chirping and Correlation With Fast-Ion Loss at NSTX

Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behaviour of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. Analysis allows for comparison between different mode character (such as quiescent, fixedfrequency, chirping, avalanching) and plasma parameters obtained from the TRANSP code such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v inj. /v A ), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (β beam,i /β). In agreement with previous work by Fredrickson et al., in this work, we find correlation between β beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to non-quiescent behaviour for magnetic fluctuations in the 50 - 200 kHz frequency band is observed along the boundary v Φ ≲ (1/4) (v inj. - 3v A ) where v Φ is the rotation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory and applications of parton pseudodistributions

We review the basic theory of the parton pseudodistributions approach and its applications to lattice extractions of parton distribution functions. The crucial idea of the approach is the realization that the correlator $M(z,p)$ of the parton fields is a function $M(ν,-z^2)$ of Lorentz invariants $ν=-(zp)$, the Ioffe time, and the invariant interval $z^2$. This observation allows to extract the Ioffe-time distribution $M(ν,-z^2)$ from Euclidean separations z accessible on the lattice. Another basic feature is the use of the ratio $M(ν,-z^2)≡M(ν,-z^2)/M(0,-z^2)$, that allows to eliminate artificial ultraviolet divergence generated by the gauge link for spacelike intervals. The remaining $z^2$-dependence of the reduced Ioffe-time distribution $M(ν,-z^2)$ corresponds to perturbative evolution, and can be converted into the scale-dependence of parton distributions $f(x,μ^2)$ using matching relations. The ν-dependence of $M(ν,-z^2)$ governs the x-dependence of parton densities $f(x,μ^2)$. The perturbative evolution was successfully observed in exploratory quenched lattice calculation. The analysis of its precise data provides a framework for extraction of parton densities using the pseudodistributions approach. It was used in the recently performed calculations of the nucleon and pion valence quark distributions. We also discuss matching conditions for the pion distribution amplitude and generalized parton distributions, the lattice studies of which are now in progress.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗