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 91 records · Page 5

Analysis and testing for rotordynamic coefficients of turbulent annular seals with different, directionally-homogeneous surface-roughness treatment for rotor and stator elements

A combined analytical-computational method has been developed to calculate the transient pressure field and rotordynamic coefficients for high-pressure annular seals. The solution procedure applies to constant-clearance or the convergent tapered geometries of multistage centrifugal pumps which may have different surface roughness conditions on the stator and rotor seal elements. In experimental calculations with the method, the turbulent equations of Hirs (1973) are modified slightly to account for different surface roughness conditions, and a perturbation analysis is employed to develop zeroth and first order perturbation equations. Zeroth equations are also used to define both the leakage and circumferential flow due to shear stresses around the stator and rotor surfaces. The solution to analytical equations for four different surface roughnesses confirm the predicted net damping for the seals. A round-holed stator pattern yielded the highest net damping and lowest leakage of all the seals tested.

Childs, D. W.↗

Analysis for leakage and rotordynamic coefficients of surface roughened tapered annular gas seals

In order to soften the effects of rub, the smooth stators of turbine gas seals are sometimes replaced by a honeycomb surface. This deliberately roughened stator and smooth rotor combination retards the seal leakage and may lead to enhanced rotor stability. However, many factors determine the rotordynamic coefficients and little is known as to the effectiveness of these honeycomb seals under various changes in the independent seal parameters. An analytical-computational method to solve for the rotordynamic coefficients of this type of compressible-flow seal is developed. The governing equations for surface roughned tapered annular gas seals are based on a modified Hirs' turbulent bulk flow model. A perturbation analysis is employed to develop zeroth and first-order perturbation equations. These equations are numerically integrated to solve for the leakage, pressure, density, and velocity for small motion of the shaft about the centered position. The resulting pressure distribution is then integrated to find the corresponding rotor-dynamic coefficients.

Nelson, C. C.↗

Aerodynamics of high-performance Ludwieg tubes.

A small perturbation analysis of the flow in a Ludwieg tube supply is described which includes both the growing turbulent boundary-layer and axial pipe nonuniformities. Measured pressure time histories are used to check the theory and its extension to high pipe Mach numbers and large flow perturbations. The nonsteady coupling of this flow to a nozzle is investigated. Calculations are made of the effects on the nozzle boundary layer, while a separate analysis is carried out for the case of complete mixing in a nozzle plenum. The double expansion nozzle is shown to extend Ludwieg-tube capability. Analytical and experimental studies are made of the aerodynamic problems that can arise in these nozzles.

Russell, D. A.↗

Evaluation of Rotordynamic Coefficients of Look-through Labyrinths by Means of a Three Volume Bulk Flow Model

To describe the compressible, turbulent flow in a labyrinth seal, a three volume bulk flow model is presented. The conservation equations for mass, momentum, and energy are established in every control volume. A perturbation analysis is performed, yielding zeroth order equations for centric rotor position and first order equations describing the flow field for small rotor motions around the seal center. The equations are integrated numerically. From perturbation pressure, the forces on the shaft and the dynamic coefficients are calculated.

Nordmann, R.↗

The passage of a distorted velocity field through a cascade of airfoils

An analysis has been developed to predict the unsteady force and moment generated by the passage of a timewise periodic total pressure distortion through an arbitrary cascade of airfoils. The mathematical formulation of this analysis is based on the assumption that the magnitudes of the timewise fluctuations of the variables which describe the flow field are small compared to their time average values. This assumption permits the development of a linear unsteady perturbation analysis about a steady flow field. In addition to this linearization assumption the fluid medium is assumed to be incompressible and inviscid. The mathematical development begins by decomposing the velocity field surrounding an infinite cascade of airfoils into its irrotational and rotational components. The rotational component is associated with an upstream unsteady total pressure distortion and is defined in terms of the vorticity field associated with the distortion pattern. The irrotational component is further decomposed into a steady and unsteady part. A combined analytical and numerical procedure has been developed to solve the field equations which govern the rotational and irrotational velocity fields. Results of this analysis show a strong influence of mean loading on the unsteady force generated by the passage of a one dimensional gust through a cascade of compressor blades.

Adamczyk, J. J.↗

Photochemical-radiative damping and instability in the stratosphere. II Numerical results

In a paper by Strobel (1977), opacity effects on the joint photochemical-radiative relaxation in the stratosphere were evaluated by a linear perturbation analysis. Significant reductions in the photochemical acceleration of the thermal relaxation rate and the ozone photochemical relaxation rate were obtained in the upper stratosphere as a result of opacity changes. For very large vertical scale height ozone density perturbations, amplifying solutions were generated by large opacity changes in the O3 dissociation and heating rates. Such solutions could indicate that the stratosphere does not possess an equilibrium ozone concentration. In the present paper, some illustrative numerical calculations are presented to demonstrate the properties of the amplifying solutions. The results indicate that ungeophysically large disturbances are required and that they cease growth before attainment of geophysically interesting amplitudes, and decay to the unperturbed state.

Strobel, D. F.↗

Nonlinear zero-sum differential game analysis by singular perturbation methods

A class of nonlinear, zero-sum differential games, exhibiting time-scale separation properties, can be analyzed by singular-perturbation techniques. The merits of such an analysis, leading to an approximate game solution, as well as the 'well-posedness' of the formulation, are discussed. This approach is shown to be attractive for investigating pursuit-evasion problems; the original multidimensional differential game is decomposed to a 'simple pursuit' (free-stream) game and two independent (boundary-layer) optimal-control problems. Using multiple time-scale boundary-layer models results in a pair of uniformly valid zero-order composite feedback strategies. The dependence of suboptimal strategies on relative geometry and own-state measurements is demonstrated by a three dimensional, constant-speed example. For game analysis with realistic vehicle dynamics, the technique of forced singular perturbations and a variable modeling approach is proposed. Accuracy of the analysis is evaluated by comparison with the numerical solution of a time-optimal, variable-speed 'game of two cars' in the horizontal plane.

Sinar, J.↗

A Quantum Mechanical MP2 Study of the Electronic Effect of Nonplanarity on the Carbon Pyramidalization of Fullerene C 60

Among C 60 ’s diverse functionalities, its potential application in CO 2 sequestration has gained increasing interest. However, the processes involved are sensitive to the molecule’s electronic structure, aspects of which remain debated and require greater precision. To address this, we performed structural optimization of fullerene C 60 using the QM MP2/6–31G* method. The nonplanarity of the optimized icosahedron is characterized by two types of dihedral angles: 138° and 143°. The 120 dihedrals of 138° occur between two hexagons intersecting at C–C bonds of 1.42 Å, while the 60 dihedrals of 143° are observed between hexagons and pentagons at C–C bonds of 1.47 Å. NBO analysis reveals less pyramidal sp 1.78 hybridization for carbons at the 1.42 Å bonds and more pyramidal sp 2.13 hybridization for the 1.47 Å bonds. Electrostatic potential charges range from -0.04 a.u. to 0.04 a.u. on the carbon atoms. Second-order perturbation analysis indicates that delocalization interactions in the C–C bonds of 1.42 Å (143.70 kcal/mol) and 1.47 Å (34.98 kcal/mol) are 22% and 38% higher, respectively, than those in benzene. MP2/Def2SVP calculations yield a correlation energy of 13.49 kcal/mol per electron for C 60 , slightly higher than the 11.68 kcal/mol for benzene. However, the results from HOMO-LUMO calculations should be interpreted with caution. This study may assist in the rational design of fullerene C 60 derivatives for CO 2 reduction systems.

36 MATERIALS SCIENCE↗

Basic Mechanisms of Solar Variability

In order to simulate the behavior of a changing Sun in a realistic way, researchers have used a perturbation analysis (the results are summarized in Endal et al. 1985). In this approach, they use a standard solar model and then vary several of the model parameters to mimic the sudden or gradual change of some physical property within the Sun. The evolution following the perturbation is followed in a physically self-consistent way, that is, hydrostatic and thermal processes occur in their normal timescales. This treatment has allowed us to determine the sensitivity of the various global parameters to physical processes affecting the solar interior, and it has guided us in defining the least complex global solar model which can address the question of climatically significant variability. Not surprisingly, magnetic fields play the crucial role. The current status of the model, what researchers have already learned from it, and the future prospects are discussed.

Sofia, Sabatino↗

A near optimal guidance algorithm for aero-assisted orbit transfer

The paper presents a near optimal guidance algorithm for aero-assited orbit plane change, based on minimizing the energy loss during the atmospheric portion of the maneuver. The guidance algorithm makes use of recent results obtained from energy state approximations and singular perturbation analysis of optimal heading change for a hypersonic gliding vehicle. This earlier work ignored the terminal constraint on altitude needed to insure that the vehicle exits that atmosphere. Thus, the resulting guidance algorithm was only appropriate for maneuvering reentry vehicle guidance. In the context of singular perturbation theory, a constraint on final altitude gives rise to a difficult terminal boundary layer problem, which cannot be solved in closed form. This paper will demonstrate the near optimality of a predictive/corrective guidance algorithm for the terminal maneuver. Comparisons are made to numerically optimized trajectories for a range or orbit plane angles.

Calise, Anthony J.↗

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of computerized symbolic manipulation in rotorcraft dynamics analysis

The potential role of symbolic manipulation programs in development and solution of the governing equations for rotorcraft dynamics problems is discussed and illustrated. Nonlinear equations of motion for a helicopter rotor blade represented by a rotating beam are developed making use of the computerized symbolic manipulation program MACSYMA. The use of computerized symbolic manipulation allows the analyst to concentrate on more meaningful tasks, such as establishment of physical assumptions, without being sidetracked by the tedious and trivial details of the algebraic manipulations. Furthermore, the resulting equations can be produced, if necessary, in a format suitable for numerical solution. A perturbation-type solution for the resulting dynamical equations is shown to be possible with a combination of symbolic manipulation and standard numerical techniques. This should ultimately lead to a greater physical understanding of the behavior of the solution than is possible with purely numerical techniques. The perturbation analysis of the flapping motion of a rigid rotor blade in forward flight is presented, for illustrative purposes, via computerized symbolic manipulation with a method that bypasses Floquet theory.

Crespo Da Silva, Marcelo R. M.↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗