Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Application of singular perturbation techniques /SPT/ and continuation methods for on-line aircraft trajectory optimization

The singular perturbation theory (SPT) approximation of optimal feedback control laws is presented and methods for on-line application of these approximations are discussed. It is demonstrated that SPT control laws break down when the current state is near the terminal target state. The use of continuation methods to improve the accuracy of the SPT approximation and to obtain global solutions of two-point boundary value problems is also discussed. As an illustration, consideration is given to the minimum-time control of a supersonic aircraft for a three-dimensional intercept problem.

Washburn, R. B.↗
NASA NTRS1982

Subsonic panel method for the efficient analysis of multiple geometry perturbations

An accurate and efficient method was developed for the aerodynamic analysis of a series of arbitrary small geometry perturbations to a given baseline configuration. The method is appropriate for wing-fuselage configurations in incompressible potential flow. Mathematical formulations are presented for three computer programs that are employed. The first program is a conventional surface panel method for calculating the baseline singularity distribution. The second program calculates a partial derivative matrix. Each element of the matrix is the rate of change of singularity strength at one point with respect to a surface coordinate of a different point. For each baseline configuration, the calculated quantities from the first two programs establish an input file for the third. The third program calculates the surface pressure distribution and forces and moments for a series of geometry perturbations.

Bristow, D. R.↗
NASA NTRS1984

Sensitivity analysis in computational aerodynamics

Information on sensitivity analysis in computational aerodynamics is given in outline, graphical, and chart form. The prediction accuracy if the MCAERO program, a perturbation analysis method, is discussed. A procedure for calculating perturbation matrix, baseline wing paneling for perturbation analysis test cases and applications of an inviscid sensitivity matrix are among the topics covered.

Bristow, D. R.↗
NASA NTRS1973

A semi-analytical method of computation of oceanic tidal perturbations in the motion of artificial satellites

The method of expansion of the satellite's perturbations, as caused by the oceanic tides, into Fourier series is discussed. The coefficients of the expansion are purely numerical and peculiar to each particular satellite. Such a method is termed as semi-analytical in celestial mechanics. Gaussian form of the differential equations for variation of elements, with the right hand sides averaged over the orbit of the satellite, is convenient to use with the semi-analytical expansion.

Musen, P.↗
DOE OSTI2021

Interference in nonlinear Compton scattering using a Schrödinger-equation approach

The interference between Compton scattering and nonlinear Compton scattering from a two-color field in the X-ray regime is theoretically examined for bound electrons. Here, the underlying phase shifts are analyzed using a perturbative approach in the incoming classical field. The perturbative approach is bench marked with a non-perturbative approach in the classical field. The interference for different combinations of linear polarization of the two fields is examined when the Compton and the nonlinear Compton scattered waves have the same wave vector and polarization. Only two cases exhibit interference. When there is interference, the calculations reveal an intrinsic phase difference between the Compton scattered wave function and the nonlinear Compton scattered wave function of either 0 or pi depending on the scattering angle.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulations for x-ray imaging of wave-packet dynamics

Previous work on imaging wave-packet dynamics with x-ray scattering revealed that the scattering patterns deviate substantially from the notion of instantaneous momentum density of the wave packet. Here we show that scattering patterns can provide clear insights into the electron wave-packet dynamics if the final state of the scattered electron and the scattered photon momentum are determined simultaneously. The scattering probability is shown to be proportional to the modulus square of the Fourier transform of the instantaneous electronic spatial wave function weighted by the final state of the electron. Several cases for the choice of final state of the electron are explored. First, the case where the final state can be measured up to a given principal quantum number n and orbital angular momentum l are presented. Next, the case where the final states can only be determined up to a given energy is discussed. Lastly, the case of an initial wave packet consisting of a large amount of a known stationary state and a small amount of an unknown stationary state is examined. The scattering profile is used to determine the properties of the unknown state in the wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Easy_PERT: a Python tool for writing PERT cards and parsing PERT card results [Slides]

This presentation begins by providing an overview of the PERT card. The PERT card uses differential operator method to compute first- and second-order tally variations due to density, composition, and reaction cross-sections. It is possible to have multiple PERT cards in one MCNP input deck to study tally variations for several sets of nuclides, reactions, and energy ranges. Furthermore, the METHOD option tells MCNP to calculate either the perturbed tally (METHOD=-1, -2, -3) or the change in the unperturbed tally (METHOD=1, 2, 3). In summation, a powerful use-case for the MCNP code PERT card is that it facilitates calculating tally sensitivities to nuclear data. Writing PERT card entries and parsing output MCTAL files is tedious and error prone. however, Easy_PERT makes use of existing tools (Faust and MCNPTools) to handle writing PERT card entries and parsing the output MCTAL files. The PERT card is early in the development process and planned upcoming capabilities include calculating sensitivities and combining MCTAL files from separate runs into one JSON file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1979

Perturbation solutions of combustion instability problems

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Googerdy, A.↗
DOE OSTI2021

Strobed imaging as a method for the determination and diagnosis of local seeing

ABSTRACT The image quality budget of many telescopes can have substantial contributions from local seeing, both ‘mirror’ and ‘dome’, which arise from turbulence and temperature variations that are difficult to quantify, measure directly, and ameliorate. We describe a method to determine the ‘local’ seeing degradation due to wavefront perturbations within the final few metres of the optical path from celestial sources to the focal plane of a ground-based telescope, using the primary instrument and along the same path taken by light from celestial sources. The concept involves placing strobed emitters along the light path to produce images on the main focal plane that ‘freeze’ different realizations of index perturbations. This method has the advantage of measuring directly the image motion and scintillation imparted by the dynamic spatial and temporal structure of local perturbations in the index of refraction along the light path, with a clean separation from seeing induced in the atmosphere above the dome. The strobed-source approach allows for rapid image motion and scintillation to be measured directly on the focal plane, even for large-aperture telescopes with wide-field instruments and slow shutters, such as that being constructed for the Rubin Observatory. A conceptual design is presented that uses the ‘guider’ CCDs in the Rubin telescope focal plane to make local-seeing measurements on demand, perhaps even during science exposures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Automated Direct Perturbation Calculations with SCALE TSUNAMI [Abstract]

In nuclear criticality safety analysis, the sensitivity of the eigenvalue keff to uncertainties in nuclear data and its evaluation are crucial. The TSUNAMI sequences within the SCALE code system offer users various options with both multigroup (MG) and continuous-energy (CE) 3D Monte Carlo (MC) transport capabilities for calculating keff sensitivity coefficients and storing them in a sensitivity data file (SDF). Each methodology available in TSUNAMI offers distinct advantages and limitations, and its effectiveness can vary based on the specific problem being solved. As a best practice, practitioners typically use the direct perturbation (DP) method as a confirmatory step alongside their sensitivity calculations to verify the accuracy of the sensitivity data generated. In this process, DP calculations are usually performed on select nuclides, those considered most important for validating their total sensitivities. However, because of code limitations, analysts use a workaround method when conducting DP calculations for a single nuclide: rather than perturbing the nuclide's microscopic cross section, an equivalent number density for this nuclide is calculated to reflect the effect of a change in the macroscopic cross section due to a perturbation in the microscopic cross section. The current approach requires rerunning the CSAS criticality calculation several times with model changes. Although this method can yield results with acceptable accuracy, it is labor-intensive and prone to errors.

AZURE: SAMMY↗
DOE OSTI2022

Differential methods for assessing sensitivity in biological models

Differential sensitivity analysis is indispensable in fitting parameters, understanding uncertainty, and forecasting the results of both thought and lab experiments. Although there are many methods currently available for performing differential sensitivity analysis of biological models, it can be difficult to determine which method is best suited for a particular model. In this paper, we explain a variety of differential sensitivity methods and assess their value in some typical biological models. First, we explain the mathematical basis for three numerical methods: adjoint sensitivity analysis, complex perturbation sensitivity analysis, and forward mode sensitivity analysis. We then carry out four instructive case studies. (a) The CARRGO model for tumor-immune interaction highlights the additional information that differential sensitivity analysis provides beyond traditional naive sensitivity methods, (b) the deterministic SIR model demonstrates the value of using second-order sensitivity in refining model predictions, (c) the stochastic SIR model shows how differential sensitivity can be attacked in stochastic modeling, and (d) a discrete birth-death-migration model illustrates how the complex perturbation method of differential sensitivity can be generalized to a broader range of biological models. Finally, we compare the speed, accuracy, and ease of use of these methods. We find that forward mode automatic differentiation has the quickest computational time, while the complex perturbation method is the simplest to implement and the most generalizable.

59 BASIC BIOLOGICAL SCIENCES↗