Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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 1,477 records · Page 82

Cryogenic hydrogen-induced air liquefaction technologies

Extensively utilizing a special advanced airbreathing propulsion archives database, as well as direct contacts with individuals who were active in the field in previous years, a technical assessment of cryogenic hydrogen-induced air liquefaction, as a prospective onboard aerospace vehicle process, was performed and documented. The resulting assessment report is summarized. Technical findings are presented relating the status of air liquefaction technology, both as a singular technical area, and also that of a cluster of collateral technical areas including: compact lightweight cryogenic heat exchangers; heat exchanger atmospheric constituents fouling alleviation; para/ortho hydrogen shift conversion catalysts; hydrogen turbine expanders, cryogenic air compressors and liquid air pumps; hydrogen recycling using slush hydrogen as heat sink; liquid hydrogen/liquid air rocket-type combustion devices; air collection and enrichment systems (ACES); and technically related engine concepts.

Escher, William J. D.↗

Probabilistic finite elements for fracture and fatigue analysis

The fusion of the probabilistic finite element method (PFEM) and reliability analysis for probabilistic fracture mechanics (PFM) is presented. A comprehensive method for determining the probability of fatigue failure for curved crack growth was developed. The criterion for failure or performance function is stated as: the fatigue life of a component must exceed the service life of the component; otherwise failure will occur. An enriched element that has the near-crack-tip singular strain field embedded in the element is used to formulate the equilibrium equation and solve for the stress intensity factors at the crack-tip. Performance and accuracy of the method is demonstrated on a classical mode 1 fatigue problem.

Liu, W. K.↗

RWF rotor-wake-fuselage code software reference guide

The RWF (Rotor-Wake-Fuselage) code was developed from first principles to compute the aerodynamics associated with the complex flow field of helicopter configurations. The code is sized for a single, multi-bladed main rotor and any configuration of non-lifting fuselage. The mathematical model for the RWF code is based on the integration of the momentum equations and Green's theorem. The unknowns in the problem are the strengths of prescribed singularity distributions on the boundaries of the flow. For the body (fuselage) a surface of constant strength source panels is used. For the rotor blades and rotor wake a surface of constant strength doublet panels is used. The mean camber line of the rotor airfoil is partitioned into surface panels. The no-flow boundary condition at the panel centroids is modified at each azimuthal step to account for rotor blade cyclic pitch variation. The geometry of the rotor wake is computers at each time step of the solution. The code produces rotor and fuselage surface pressures, as well as the complex geometry of the evolving rotor wake.

Berry, John D.↗

Progress in multirate digital control system design

A new methodology for multirate sampled-data control design based on a new generalized control law structure, two new parameter-optimization-based control law synthesis methods, and a new singular-value-based robustness analysis method are described. The control law structure can represent multirate sampled-data control laws of arbitrary structure and dynamic order, with arbitrarily prescribed sampling rates for all sensors and update rates for all processor states and actuators. The two control law synthesis methods employ numerical optimization to determine values for the control law parameters. The robustness analysis method is based on the multivariable Nyquist criterion applied to the loop transfer function for the sampling period equal to the period of repetition of the system's complete sampling/update schedule. The complete methodology is demonstrated by application to the design of a combination yaw damper and modal suppression system for a commercial aircraft.

Berg, Martin C.↗

The prediction of pseudo-resonance positions in the Schwinger variational principle

The Schwinger Variational Principle is appled to s-wave electron-hydrogen atom scattering. It is shown computationally, that, consistent with a paper by B. Apagyi, P. Levay, and K. Ladanyi, there are pseudo-resonances at the static exchange level of approximation, but not at the static level. The T-matrix was employed, as well as, the K-matrix version of the Schwinger Principle, with a real Slater basis, and obtained the same results in both. The origin of the pseudo-resonances as resulting from singularities in the separable potential that is effectively employed in the Lippmann-Schwinger equation from which the Schwinger Variational Principle can be derived. The determination of the pseudo-resonance parameters from the separable potential is computationally inexpensive and may be used to predict the pseudo-resonance parameters for the scattering calculations so that they may be avoided.

Weatherford, Charles A.↗

A comparison of numerical methods for the Rayleigh equation in unbounded domains

A second-order finite difference and two spectral methods, including a Chebyshev tau and a Chebyshev collocation method were implemented to determine the linear hydrodynamic stability of an unbounded shear flow. The velocity profile of the basic flow in the stability analysis mimicks that of a one-stream free mixing layer. Local and global eigenvalue solution methods are used to determine individual eigenvalues and the eigenvalue spectrum, respectively. The calculated eigenvalue spectrum includes a discrete mode, a continuous spectrum associated with the equation singularity and a continuous spectrum associated with the domain unboundedness. The efficiency and the accuracy of these discretization methods in the prediction of the eigensolutions of the discrete mode were evaluated by comparison with a conventional shooting procedure. Their capabilities in mapping out the continuous eigenvalue spectra are also discussed.

Liou, W. W.↗

A hybrid perturbation-Galerkin technique that combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be used to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

Quasilinear saturation of forced current sheet tearing modes

Numerical studies of tearing modes in a nearly singular forced current sheet equilibrium (Liewer and Payne, 1990) show that the modes saturate quasilinearly when the width of the magnetic island formed by the reconnection is on the order of several times the linear mode width which scales as approximately (kS) exp -2/5, where S is the Lundquist number and k is the wavenumber. The modes saturate quasilinearly by flattening the current profile, converting magnetic energy into plasma energy. The longer wavelength modes, which saturate at higher levels, release the most energy. These modes may, nonlinearly, play a role in coronal heating when sharp current sheets form as a result of global magnetic stresses.

Liewer, Paulett C.↗

A boundary integral approach to the scattering of nonplanar acoustic waves by rigid bodies

The acoustic scattering of an incident wave by a rigid body can be described by a singular Fredholm integral equation of the second kind. This equation is derived by solving the wave equation using generalized function theory, Green's function for the wave equation in unbounded space, and the acoustic boundary condition for a perfectly rigid body. This paper will discuss the derivation of the wave equation, its reformulation as a boundary integral equation, and the solution of the integral equation by the Galerkin method. The accuracy of the Galerkin method can be assessed by applying the technique outlined in the paper to reproduce the known pressure fields that are due to various point sources. From the analysis of these simpler cases, the accuracy of the Galerkin solution can be inferred for the scattered pressure field caused by the incidence of a dipole field on a rigid sphere. The solution by the Galerkin technique can then be applied to such problems as a dipole model of a propeller whose pressure field is incident on a rigid cylinder. This is the groundwork for modeling the scattering of rotating blade noise by airplane fuselages.

Gallman, Judith M.↗

A parallel algorithm for nonlinear convection-diffusion equations

A parallel algorithm for the efficient solution of nonlinear time-dependent convection-diffusion equations with small parameter on the diffusion term is presented. The method is based on a physically motivated domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. The method is suitable for the solution of problems arising in the simulation of fluid dynamics. Experimental results for a nonlinear equation in two-dimensions are presented.

Scroggs, Jeffrey S.↗

Three-dimensional stress intensity factor analysis of a surface crack in a high-speed bearing

The boundary element method is applied to calculate the stress intensity factors of a surface crack in the rotating inner raceway of a high-speed roller bearing. The three-dimensional model consists of an axially stressed surface cracked plate subjected to a moving Hertzian contact loading. A multidomain formulation and singular crack-tip elements were employed to calculate the stress intensity factors accurately and efficiently for a wide range of configuration parameters. The results can provide the basis for crack growth calculations and fatigue life predictions of high-performance rolling element bearings that are used in aircraft engines.

Ballarini, Roberto↗

Three-dimensional kinematic reconnection of plasmoids

The kinematic reconnection model of Lau and Finn (1990) is applied in a theoretical investigation of three-dimensional plasmoid morphology, as seen in the solar corona and earth magnetotail. The derivation of the governing equations is outlined; long, short, and periodic plasmoid models are developed; and the evolution of the stable and unstable manifolds in these models is presented graphically. It is inferred that sheet currents and tangential discontinuities can form on surfaces topologically identical to those where Delta(phi) about equal to Delta(z) singularities occur in the kinematic reconnection model, and can be broadened in a similar way by nonideal effects. Two such surfaces exist in long plasmoids, and also in short plasmoids in the presence of finite resistivity; the intertwined surfaces characteristic of the periodic plasmoid form a fractal set but merge in the presence of finite resistivity, producing structures similar to those proposed by the sheet-current theory of Parker (1983).

Lau, Yun-Tung↗

A time-accurate implicit method for chemically reacting flows at all Mach numbers

The objective of this work is to develop a unified solution algorithm capable of treating time-accurate chemically reacting flows at all Mach numbers, ranging from molecular diffusion velocities to supersonic speeds. A rescaled pressure term is used in the momentum equation to circumvent the singular behavior of pressure at low Mach numbers. A dual time-stepping integration procedure is established. The system eigenvalues become well behaved and have the same order of magnitude, even in the very low Mach number regime. The computational efficiency for moderate and high speed flow is competitive with the conventional density-based scheme. The capabilities of the algorithm are demonstrated by applying it to selected model problems including nozzle flows and flame dynamics.

Withington, J. P.↗

Approach to control moment gyroscope steering using feedback linearization

This paper presents an approach for controlling spacecraft equipped with control moment gyroscopes. A technique from feedback linearization theory is used to transform the original nonlinear problem to an equivalent linear form without approximating assumptions. In this form, the spacecraft dynamics appear linearly, and are decoupled from redundancy in the system of gyroscopes. A general approach to distributing control effort among the available actuators is described which includes provisions for redistribution of rotors, explicit bounds in gimbal rates, and guaranteed operation at or near singular configurations. A particular algorithm is developed for systems of double-gimbal devices, and demonstrated in two examples for which existing approaches fail to give adequate performance.

Dzielski, John↗

Prediction of vortex shedding from forebodies with chines

An engineering prediction method to calculate vortex shedding from noncircular forebodies with sharp chine edges in subsonic flow at large incidence angles is presented. The forebody is represented by two- and three-dimensional singularities, and the lee side vortex wake is modeled by discrete vortices in crossflow planes along the body. The computational procedure is described, and comparisons of measured and predicted surface pressure distributions and predicted flow field vectors are presented to illustrate the method.

Mendenhall, Michael R.↗

Trajectory optimization for kinematically redundant arms

A review of local optimization methods for resolving joint configurations in underconstrained manipulation tasks is conducted. A new approach is developed for observing joint limits and avoiding obstacles during the trajectory planning. The methodology is used in a four-link arm example to avoid a workspace singularity and is compared with results using the extended Moore-Penrose technique. An alternative measure of arm 'manipulability' based directly on the rank of the Jacobian is also introduced.

Carignan, Craig R.↗

Constructing second-order structural models with specified mass submatrix condition numbers

This paper presents a simple new algorithm for constructing second-order models for any given flexible structure, with these models satisfying the constraint that each leading principal submatrix of the mass matrix must have a specified condition number. Now, the effects of rounding errors in practical computation generally increase with increasing condition number, i.e. as the matrix considered approaches singularity. Thus, the new models should prove useful for the systematic testing of the numerical properties of many algorithms in the dynamics and control of flexible structures.

Williams, T.↗

Intraseasonal oscillations in the global atmosphere. I - Northern Hemisphere and tropics

Oscillatory modes in the Northern Hemisphere and in the tropics were examined systematically. The 700 mb heights were used to analyze extratropical oscillations, and the outgoing longwave radiation to study tropical oscillations in convection. All datasets were band-pass filtered to focus on the intraseasonal (IS) band of 10-120 days. Leading spatial patterns of variability were obtained by applying empirical orthogonal function analysis to these IS data. The leading principal components were subjected to singular spectrum analysis. In the Northern Hemisphere, there are two important modes of oscillation with periods near 48 and 23 days, respectively. The 48-day mode is the most important of the two. It has both traveling and standing components, and is dominated by a zonal wavenumber two. The 23-day mode has the spatial structure and propagation properties described by Branstator and by Kushnir (1987). In the tropics, the 40-50 day oscillation documented by Madden and Julian (1972), Weickmann (1983), Lau, and their colleagues, dominates the Indian and Pacific oceans from 60 deg E to the date line. From 170 deg W to 90 deg W, however, a 24-28 day oscillation is equally strong. The extratropical modes are often independent of, and sometimes lead, the tropical modes.

Ghil, Michael↗