Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 775 records · Page 43

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence both Duffing's equation and the van der Pol equation are solved using each method.

Ahmed, A.↗

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence Duffing's equation, the van der Pol equation and the oscillator with quadratic damping problem are solved using each method.

Ahmed, A. H.↗

Depletion GPT-free sensitivity analysis for reactor eigenvalue problems

This manuscript introduces a novel approach to solving depletion perturbation theory problems without the need to set up or solve the generalized perturbation theory (GPT) equations. The approach, hereinafter denoted generalized perturbation theory free (GPT-Free), constructs a reduced order model (ROM) using methods based in perturbation theory and computes response sensitivity profiles in a manner that is independent of the number or type of responses, allowing for an efficient computation of sensitivities when many responses are required. Moreover, the reduction error from using the ROM is quantified in the GPT-Free approach by means of a Wilks' order statistics error metric denoted the K-metric. Traditional GPT has been recognized as the most computationally efficient approach for performing sensitivity analyses of models with many input parameters, e.g. when forward sensitivity analyses are computationally intractable. However, most neutronics codes that can solve the fundamental (homogenous) adjoint eigenvalue problem do not have GPT capabilities unless envisioned during code development. The GPT-Free approach addresses this limitation by requiring only the ability to compute the fundamental adjoint. This manuscript demonstrates the GPT-Free approach for depletion reactor calculations performed in SCALE6 using the 7x7 UAM assembly model. A ROM is developed for the assembly over a time horizon of 990 days. The approach both calculates the reduction error over the lifetime of the simulation using the K-metric and benchmarks the obtained sensitivities using sample calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational analysis of a stability robustness margin for structured real-parameter perturbations

An efficient computational method is presented for stability robustness analysis with structured real-parameter perturbations. A generic model of a class of uncertain dynamical systems is used as an example. The parameter uncertainty is characterized by a real scalar, epsilon. Multilinearity of the closed-loop characteristic polynomial is exploited to permit application of the mapping theorem to calculate the stability robustness margin. It is found that sensitive geometry of the stability boundary in the epsilon, omega-plane renders problematic the calculation of the minimum epsilon as a function of omega. This difficulty is avoided by calculating the minimum distance to the image of the uncertainty domain over omega as a function of epsilon. It is also shown that a certain class of uncertain dynamical systems has the required multilinearity property and are thus amenable to the proposed technique.

Wedell, Evan↗

Effects of gravity-gradient torque on the rotational motion of a triaxial satellite in a precessing elliptic orbit.

A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the effects of gravity-gradient torque on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is free to precess and spin. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.

Cochran, J. E.↗

A hybrid perturbation-Galerkin technique for partial differential equations

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Geer, James F.↗

Long-term motion in a restricted problem of rotational motion

A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the long-term effects of gravity-gradient torque and orbital evolution on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is constrained to precess and spin with constant rates. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite to second order in n/omega sub 0, where n is the orbital mean motion of the center of mass and omega sub 0 is a reference value of the magnitude of the satellite's rotational angular velocity. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.

Cochran, J. E.↗

Second-order estimates in singular perturbation of output regulator problems

A method of order reduction in singular perturbation problems due to Chow and Kokotovic is studied. The method is used for an output regulator problem in control theory. The quadratic term in the asymptotic expansion of Chow and Kokotovic is calculated and it is shown by some examples how this term can give evidence of the accuracy of the procedure.

Rublein, G.↗

Rapid approximate determination of nonlinear solutions - Application to aerodynamic flows and design/optimization problems

Stahara et al. (1978) have considered the use of an approximation technique which employs two or more nonlinear base solutions determined by the full computational method to predict entire families of related nonlinear solutions. The present investigation provides results for several applications of that method which demonstrate both its accuracy and its utility for engineering applications. Attention is given to the perturbation concept and methods, aspects of coordinate straining, aspects of analytical formulation, and an application to surface properties. In a discussion of the results, single and multiple parameter perturbations are considered along with a combination of the approximation method with optimization procedures. The results show that it is possible to combine in certain cases large savings in computational cost with improved optimization.

Stahara, S. S.↗

High-resolution electron spectroscopy of lanthanide (Ce, Pr, and Nd) complexes of cyclooctatetraene: The role of 4f electrons

Cerium, praseodymium, and neodymium complexes of 1,3,5,7-cyclooctatetraene (COT) complexes were produced in a laser-vaporization metal cluster source and studied by pulsed-field ionization zero electron kinetic energy spectroscopy and quantum chemical calculations. The computations included the second-order Moller-Plesset perturbation theory, the coupled cluster method with single, double, and perturbative triple excitations, and the state-average complete active space self-consistent field method. The spectrum of each complex exhibits multiple band systems and is assigned to ionization of several low-energy electronic states of the neutral complex. This observation is different from previous studies of M(COT) (M = Sc, Y, La, and Gd), for which a single band system was observed. The presence of the multiple low-energy electronic states is caused by the splitting of the partially filled lanthanide 4f orbitals in the ligand field, and the number of the low-energy states increases rapidly with increasing number of the metal 4f electrons. On the other hand, the 4f electrons have a small effect on the geometries and vibrational frequencies of these lanthanide complexes.

36 MATERIALS SCIENCE↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

Subsonic 3-D surface panel method for rapid analysis of multiple geometry perturbations

For any baseline aircraft configuration, a matrix of partial derivatives of surface velocity potential with respect to geometry coordinates is calculated. Linear extrapolation is then used to analyze the subcritical potential flow corresponding to a series of arbitrary small geometry perturbations. Each perturbation analysis is more than an order of magnitude more efficient than a conventional panel method solution because no influence coefficients are calculated and no large system of linear algebraic equations is solved. Wing and wing-fuselage examples are presented to demonstrate that the predicted pressure distributions are nearly exact for large changes to wing camber, thickness, and leading edge radius.

Bristow, D. R.↗

Measurement of dijet production with a veto on additional central jet activity in pp collisions at $ \sqrt {s} = 7 $ TeV using the ATLAS detector

A measurement of jet activity in the rapidity interval bounded by a dijet system is presented. Events are vetoed if a jet with transverse momentum greater than 20 GeV is found between the two boundary jets. The fraction of dijet events that survive the jet veto is presented for boundary jets that are separated by up to six units of rapidity and with mean transverse momentum 50 < p¯ T < 500 GeV. The mean multiplicity of jets above the veto scale in the rapidity interval bounded by the dijet system is also presented as an alternative method for quantifying perturbative QCD emission. The data are compared to a next-to-leading order plus parton shower prediction from the powheg-box, an all-order resummation using the hej calculation and the pythia, herwig++ and alpgen event generators. In conclusion, the measurement was performed using pp collisions at √s = 7 TeV using data recorded by the ATLAS detector in 2010.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analysis of wind perturbation activity for the NSTS launches

A method for statistical analysis of wind perturbation activity is proposed for application in NSTS prelaunch assessment of the wind environment. Methods for comparing day-of-launch wind perturbation standard deviations with historical data at the launch site are described. A statistical model based on the Gumbel extreme value distribution is developed for the largest vector and component perturbations in a wind profile as a function of high-pass filter cutoff wavelength. Application of the model is proposed for assessing extreme perturbations within prelaunch sequential Jimsphere wind profiles.

Adelfang, Stanley I.↗

Energy levels, weighted oscillator strengths, transition rates, lifetimes, hyperfine structures and isotope shifts of He-like-Sn deduced from the relativistic multiconfiguration Dirac–Hartree–Fock and second-order many-body perturbation theory calculations

We disclose relativistic multiconfiguration Dirac–Hartree–Fock (MCDHF) spectrum calculations for SnXLIX. Energy levels, weighted oscillator strengths, isotope shifts, hyperfine structure and Landé g{sub J} factors are calculated for 127 odd- and even-parity states as well as lifetimes and transition rates between these states. To scrutinize the accuracy of our results, we have implemented parallel calculations using a Flexible Atomic Code (FAC) by introducing the Second-Order Many-Body Perturbation Theory (MBPT) method. Additionally, the Breit interaction and leading quantum electrodynamic effects (QED) are included as perturbations in extensive relativistic configuration interaction (RCI) calculations. We signal that, our calculations for SnXLIX are made for the first time and they provide to date the most accurate and complete atomic data related to this ion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Second-order small disturbance theory for hypersonic flow over power-law bodies

A mathematical method for determining the flow field about power-law bodies in hypersonic flow conditions is developed. The second-order solutions, which reflect the effects of the second-order terms in the equations, are obtained by applying the method of small perturbations in terms of body slenderness parameter to the zeroth-order solutions. The method is applied by writing each flow variable as the sum of a zeroth-order and a perturbation function, each multiplied by the axial variable raised to a power. The similarity solutions are developed for infinite Mach number. All results obtained are for no flow through the body surface (as a boundary condition), but the derivation indicates that small amounts of blowing or suction through the wall can be accommodated.

Townsend, J. C.↗

Accurate ab initio Quartic Force Fields of Cyclic and Bent HC2N Isomers

Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.

Inostroza, Natalia↗