A Method for Estimating Global Subgrid‐Scale Orographic Gravity‐Wave Temperature Perturbations in Chemistry‐Climate Models
Explore the source record for details and available documents.
M. Weimer↗
SEARCH · Engineering Papers
Results for “Perturbative methods”
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.
Optimal Transport for Super Resolution Applied to Astronomy Imaging
Super resolution is an essential tool in optics, especially on interstellar scales, due to physical laws restricting possible imaging resolution. We propose using optimal transport and entropy for super resolution applications. We prove that the reconstruction is accurate when sparsity is known and noise or distortion is small enough. We prove that the optimizer is stable and robust to noise and perturbations. We compare this method to a state of the art convolutional neural network and get similar results for much less computational cost and greater methodological flexibility.
Rawson, Michael G.↗
Resonant frequency calculations using a hybrid perturbation-Galerkin technique
A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.
Geer, James F.↗
Resonant frequency calculations using a hybrid perturbation-Galerkin technique
A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degrees of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.
Geer, James F.↗
Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling Wave Tube Structures
Conventional methods used to measure the cold- test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. With the advent of accurate three-dimensional (3-D) helical circuit models, these standard approximations can be fully Investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the 3-D simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using 3-D computational methods rather than performing costly and time consuming experimental cold-test measurements.
Kory, Carol L.↗
The calculation of upper and lower bounds of energy eigenvalues in perturbation theory by means of partitioning techniques.
Quantum mechanical perturbation theory calculation of upper and lower bounds of energy eigenvalues using partitioning methods
PERTURBATION THEORY↗
Real-time coupled-cluster approach for the cumulant Green's function
Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.
Vila, Fernando D.↗
Effects of resonant tesseral gravity coefficients on Viking-type orbits.
A special perturbation technique has been developed for the long-period and secular motion of a Viking-type orbiter. The technique uses a method of singly averaging the perturbations over the mean anomaly to determine changes in orbital elements due to atmospheric drag, solar radiation pressure, solar gravity, and the asymmetrical Mars gravitational field. The technique has been specialized to treat resonant tesseral gravity effects resulting from the commensurability of the satellite orbit period with the rotational period of the primary. Results are presented which indicate that orbiter lifetimes can be significantly reduced by resonant tesserals for orbits near the critical inclination. Conversely, it is shown that lifetimes can be extended via capture in resonant orbits which are stable against small drag forces.
Compton, H. R.↗
VARI3D & PERSENT: Perturbation and Sensitivity Analysis
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Evolution of solar magnetic fields - A new approach to MHD initial-boundary value problems by the method of nearcharacteristics
A method of analysis for the MHD initial-boundary problem is presented in which the model's formulation is based on the method of nearcharacteristics developed by Werner (1968) and modified by Shin and Kot (1978). With this method, the physical causality relationship can be traced from the perturbation to the response as in the method of characteristics, while achieving the advantage of a considerable reduction in mathematical procedures. The method offers the advantage of examining not only the evolution of nonforce free fields, but also the changes of physical conditions in the atmosphere accompanying the evolution of magnetic fields. The physical validity of the method is demonstrated with examples, and their significance in interpreting observations is discussed.
Nakagawa, Y.↗
Asymptotic Energies and QED Shifts for Rydberg States of Helium
This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.
Drake, G.W.F.↗
Advances for QCD and the standard model: color-confining light-front holography and the principle of maximum conformality
Here, I review how the application of superconformal quantum mechanics and light-front holography leads to new insights into the physics of color confinement, the spectroscopy and dynamics of hadrons, as well as surprising supersymmetric relations between the masses of mesons, baryons, and tetraquarks. Spontaneous chiral symmetry breaking is automatically fulfilled by supersymmetric Light-Front QCD. The light-front holographic approach (HLFQCD) also predicts the behavior of the QCD running coupling and other observables from the nonperturbative color-confining domain to the perturbative domain. One can determine the QCD running coupling to high precision from the data of just a single experiment over the entire perturbative regime by using the Principle of Maximum Conformality (PMC). The PMC, which generalizes the conventional Gell-Mann-Low method for scale-setting in perturbative QED to non-Abelian QCD, provides a rigorous method for achieving unambiguous scheme-independent, fixed-order Standard Model predictions, consistent with the principles of the renormalization group. I also briefly review a novel feature of hadronic physics predicted by QCD: intrinsic heavy quarks.
Brodsky, Stanley J. [SLAC National Accelerator Lab↗
Application of Krylov-Bogoliubov method to the solution of the stellar three body problem.
Higher order perturbations of elements in three- body problem determined using methods of Krylov- Bogoliubov and Poincare
THREE-BODY PROBLEM↗
Loss Landscape Analysis for Reliable Quantized ML Models for Scientific Sensing
In this paper, we propose a method to perform empirical analysis of the loss landscape of machine learning (ML) models. The method is applied to two ML models for scientific sensing, which necessitates quantization to be deployed and are subject to noise and perturbations due to experimental conditions. Our method allows assessing the robustness of ML models to such effects as a function of quantization precision and under different regularization techniques -- two crucial concerns that remained underexplored so far. By investigating the interplay between performance, efficiency, and robustness by means of loss landscape analysis, we both established a strong correlation between gently-shaped landscapes and robustness to input and weight perturbations and observed other intriguing and non-obvious phenomena. Our method allows a systematic exploration of such trade-offs a priori, i.e., without training and testing multiple models, leading to more efficient development workflows. This work also highlights the importance of incorporating robustness into the Pareto optimization of ML models, enabling more reliable and adaptive scientific sensing systems.
Baldi, Tommaso [Pisa, Scuola Normale Superiore]↗
Estimation of geopotential from satellite-to-satellite range rate data: Numerical results
A technique for high-resolution geopotential field estimation by recovering the harmonic coefficients from satellite-to-satellite range rate data is presented and tested against both a controlled analytical simulation of a one-day satellite mission (maximum degree and order 8) and then against a Cowell method simulation of a 32-day mission (maximum degree and order 180). Innovations include: (1) a new frequency-domain observation equation based on kinetic energy perturbations which avoids much of the complication of the usual Keplerian element perturbation approaches; (2) a new method for computing the normalized inclination functions which unlike previous methods is both efficient and numerically stable even for large harmonic degrees and orders; (3) the application of a mass storage FFT to the entire mission range rate history; (4) the exploitation of newly discovered symmetries in the block diagonal observation matrix which reduce each block to the product of (a) a real diagonal matrix factor, (b) a real trapezoidal factor with half the number of rows as before, and (c) a complex diagonal factor; (5) a block-by-block least-squares solution of the observation equation by means of a custom-designed Givens orthogonal rotation method which is both numerically stable and tailored to the trapezoidal matrix structure for fast execution.
Thobe, Glenn E.↗
Numerical computation of steady-state acoustic disturbances in flow
Two time domain methods for computing two dimensional steady-state acoustic disturbances propagating through internal subsonic viscous flow fields in the presence of variable area are investigated. The first method solves the Navier-Stokes equations for the combined steady and acoustic field together and subtracts the steady flow to obtain the acoustic field. The second method solves a system of perturbation equations to obtain the acoustic disturbances, making use of a separate steady flow computation as input to the system. In each case the periodic steady-state acoustic fluctuations are obtained numerically on a supercomputer using a second order unsplit explicit MacCormack predictor-corrector method. Results show that the first method is not very effective for computing acoustic disturbances of even moderate amplitude. It appears that more accurate steady flow algorithms are required for this method to succeed. On the other hand, linear and nonlinear acoustic disturbances extracted from the perturbation approach are shown to exhibit expected behavior for the problems considered. It is also found that inflow boundary conditions for an equivalent uniform duct can be successfully applied to a nonuniform duct to obtain steady-state acoustic disturbances.
Watson, W. R.↗
Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond
While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Cosmological perturbation theory for large scale structure in phase space
We develop a framework for Large Scale Structure (LSS) perturbation theory, that solves the Vlasov-Poisson system of equations for the distribution function in full phase space. This approach relaxes the usual apriori assumption of negligible velocity dispersion underlying the Standard Perturbation Theory (SPT). We apply the new method to rederive the usual SPT kernels up to third order in the perturbative expansion. We also show that a counterterm, identical to the one introduced by standard Effective Field Theory (EFT) methods, naturally arises within our framework. We finish by making a precise connection to EFT techniques, which reveals the necessity of the EFTofLSS to self-consistently model the long-wavelength fluid, and illustrates the importance of having theoretical control over short distance fluctuations.
Cosmological perturbation theory in GR and beyond↗