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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily-shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrary channel. The scheme employs a mesh which is a nearly-conformal 'O' mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, J. C., Jr.↗

NASA NTRS1985

Energy management of three-dimensional minimum-time intercept

A real-time computer algorithm to control and optimize aircraft flight profiles is described and applied to a three-dimensional minimum-time intercept mission. The proposed scheme has roots in two well known techniques: singular perturbations and neighboring-optimal guidance. Use of singular-perturbation ideas is made in terms of the assumed trajectory-family structure. A heading/energy family of prestored point-mass-model state-Euler solutions is used as the baseline in this scheme. The next step is to generate a near-optimal guidance law that will transfer the aircraft to the vicinity of this reference family. The control commands fed to the autopilot (bank angle and load factor) consist of the reference controls plus correction terms which are linear combinations of the altitude and path-angle deviations from reference values, weighted by a set of precalculated gains. In this respect the proposed scheme resembles neighboring-optimal guidance. However, in contrast to the neighboring-optimal guidance scheme, the reference control and state variables as well as the feedback gains are stored as functions of energy and heading in the present approach. Some numerical results comparing open-loop optimal and approximate feedback solutions are presented.

Kelley, H. J.↗

NASA NTRS1985

Theoretical analysis of linearized acoustics and aerodynamics of advanced supersonic propellers

The derivation of a formula for prediction of the noise of supersonic propellers using time domain analysis is presented. This formula is a solution of the Ffowcs Williams-Hawkings equation and does not have the Doppler singularity of some other formulations. The result presented involves some surface integrals over the blade and line integrals over the leading and trailing edges. The blade geometry, motion and surface pressure are needed for noise calculation. To obtain the blade surface pressure, the observer is moved onto the blade surface and a linear singular integral equation is derived which can be solved numerically. Two examples of acoustic calculations using a computer program are currently under development.

Farassat, F.↗

NASA NTRS1985

Improving the Accuracy of the Boundary Integral Method Based on the Helmholtz Integral

Several recent papers in the literature have been based on various forms of the Helmholtz integral to compute the radiation fields of vibrating bodies. The surface integral form is given. The symbols of P,R micron, rho,G,R,V, and S micron are acoustic pressure, source coordinate, angular frequency, fluid density, Green function, field coordinate, surface velocity and body surface respectively. A discretized form of the surface integral is also given. Solutions to the surface integral are complicated with the singularity of the Green function at R=R micron and with the uniqueness problem at interior eigen frequencies of the enclosed space. The use of the interior integral circumvents the singularity problem since the field points are chosen in the interior space of the vibrating body where a zero pressure condition exists. The interior integral form is given. The method to improve the accuracy is detailed. Examples of the method is presented for a variety of radiators.

Koopmann, G. H.↗

NASA NTRS1985

A research program in active control/aeroelasticity

In this research period a synthesis methodology for multifunctional robust aeroservoelastic systems was developed. The development consisted of the following stages: (1) development of an universal diagram to determine phase and gain margins of a multi-input multi-output (MIMO) system using singular value based stability margin criteria; (2) determination of singular value gradients with respect to design parameters and their application to improve stability margins of multiloop system; and (3) application of constrained optimization techniques to synthesize a low order robust controller for a high order MIMO aeroservoelastic system while satisfying several design constraints on the dynamic loads and responses and stability margins at the plant input and output. Overall research activities are summarized.

Source record↗

NASA NTRS1986

Theoretical three-and four-axis gimbal robot wrists

In high-performance flight simulations, a four-axis gimbal system allows all possible rotations with acceptable gimbal angle rates while it avoids the so-callled 'gimbal lock' that occurs when gimbal rotational axes are colinear. In this paper, pertinent equations (including quaternions) are assembled for a hypothetical robot wrist, functionally equivalent to this four-axis gimbal system, and also for a true three-axis gimbal robot wrist. These equations are used to simulate the rotation of a robot hand by the robot wrist in response to operator rotational velocity commands to the robot hand. Near gimbal lock (wrist singularity), excessive rotational rates occur. Scaling the rates, which is necessary for the three-gimbal robot wrist to prevent rate limiting, introduces an undesirable time delay in the robot hand rotation with respect to the commanded rotation. However, the merit of the four-gimbal robot wrist is that the fourth gimbal angle keeps the robot wrist away from the singularity so that the robot hand moves exactly as commanded. It appears that in a 'worst-type' maneuver of the robot hand, the fourth gimbal angle can be defined so that none of the gimbal angle rates exceed about twice the commanded rates.

Barker, L. K.↗

NASA NTRS1986

Calculations of axisymmetric vortex sheet roll-up using a panel and a filament model

A method for calculating the self-induced motion of a vortex sheet using discrete vortex elements is presented. Vortex panels and vortex filaments are used to simulate two-dimensional and axisymmetric vortex sheet roll-up. A straight forward application using vortex elements to simulate the motion of a disk of vorticity with an elliptic circulation distribution yields unsatisfactroy results where the vortex elements move in a chaotic manner. The difficulty is assumed to be due to the inability of a finite number of discrete vortex elements to model the singularity at the sheet edge and due to large velocity calculation errors which result from uneven sheet stretching. A model of the inner portion of the spiral is introduced to eliminate the difficulty with the sheet edge singularity. The model replaces the outermost portion of the sheet with a single vortex of equivalent circulation and a number of higher order terms which account for the asymmetry of the spiral. The resulting discrete vortex model is applied to both two-dimensional and axisymmetric sheets. The two-dimensional roll-up is compared to the solution for a semi-infinite sheet with good results.

Kantelis, J. P.↗

NASA NTRS1984

High accuracy heat capacity measurements through the lambda transition of helium with very high temperature resolution

A measurement of the heat capacity singularity of helium at the lambda transition was performed with the aim of improving tests of the Renormalization Group (RG) predictions for the static thermodynamic behavior near the singularity. The goal was to approach as closely as possible to the lambda-point while making heat capacity measurements of high accuracy. To do this, a new temperature sensor capable of unprecedented resolution near the lambda-point, and two thermal control systems were used. A short description of the theoretical background and motivation is given. The initial apparatus and results are also described.

Fairbanks, W. M.↗

NASA NTRS1985

Comparison of uniform perturbation and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T.-C.↗

NASA NTRS1985

Finite difference solutions of the Euler equations in the vicinity of sharp edges

Attempts have been made to explain why finite difference solutions of the Euler equations can describe flows with large vortical structures around sharp-edged bodies. The present paper is concerned with the influence of a singular sharp edge on the truncation error for a set of discretized Euler equations. An analysis is conducted of the distribution of the truncation error of one finite difference approximation of the Euler equations near a sharp edge of a thin plate. The analysis leads to a determination of the size of the region of the neighborhood of such a singularity. Attention is given to the consistency of a discretization of the Euler equations, and numerical experiments.

Hartwich, P.-M.↗

NASA NTRS1985

The crack problem for bonded nonhomogeneous materials under antiplane shear loading

The singular nature of the crack tip stress field in a nonhomogeneous medium having a shear modulus with a discontinuous derivative was investigated. The problem is considered for the simplest possible loading and geometry, namely the antiplane shear loading of two bonded half spaces in which the crack is perpendicular to the interface. It is shown that the square-root singularity of the crack tip stress field is unaffected by the discontinuity in the derivative of the shear modulus. The problem is solved for a finite crack and extensive results are given for the stress intensity factors.

Erdogan, F.↗

NASA NTRS1986

Recent developments in rotary-wing aerodynamic theory

Current progress in the computational analysis of rotary-wing flowfields is surveyed, and some typical results are presented in graphs. Topics examined include potential theory, rotating coordinate systems, lifting-surface theory (moving singularity, fixed wing, and rotary wing), panel methods (surface singularity representations, integral equations, and compressible flows), transonic theory (the small-disturbance equation), wake analysis (hovering rotor-wake models and transonic blade-vortex interaction), limitations on computational aerodynamics, and viscous-flow methods (dynamic-stall theories and lifting-line theory). It is suggested that the present algorithms and advanced computers make it possible to begin working toward the ultimate goal of turbulent Navier-Stokes calculations for an entire rotorcraft.

Johnson, W.↗

NASA NTRS1987

Convergence of strain energy release rate components for edge-delaminated composite laminates

Strain energy release rates for edge delaminated composite laminates were obtained using quasi 3 dimensional finite element analysis. The problem of edge delamination at the -35/90 interfaces of an 8-ply composite laminate subjected to uniform axial strain was studied. The individual components of the strain energy release rates did not show convergence as the delamination tip elements were made smaller. In contrast, the total strain energy release rate converged and remained unchanged as the delamination tip elements were made smaller and agreed with that calculated using a classical laminated plate theory. The studies of the near field solutions for a delamination at an interface between two dissimilar isotropic or orthotropic plates showed that the imaginary part of the singularity is the cause of the nonconvergent behavior of the individual components. To evaluate the accuracy of the results, an 8-ply laminate with the delamination modeled in a thin resin layer, that exists between the -35 and 90 plies, was analyzed. Because the delamination exists in a homogeneous isotropic material, the oscillatory component of the singularity vanishes.

Raju, I. S.↗

NASA NTRS1986

Plates and shells containing a surface crack under general loading conditions

The severity of the underlying assumptions of the line-spring model (LSM) are such that verification with three-dimensional solutions is necessary. Such comparisons show that the model is quite accurate, and therefore, its use in extensive parameter studies is justified. Investigations into the endpoint behavior of the line-spring model have led to important conclusions about the ability of the model to predict stresses in front of the crack tip. An important application of the LSM was to solve the contact plate bending problem. Here the flexibility of the model to allow for any crack shape is exploited. The use of displacement quantities as unknowns in the formulation of the problem leads to strongly singular integral equations, rather than singular integral equations which result from using displacement derivatives. The collocation method of solving the integral equations was found to be better and more convenient than the quadrature technique. Orthogonal polynomials should be used as fitting functions when using the LSM as opposed to simpler functions such as power series.

Joseph, Paul F.↗

NASA NTRS1987

A hybrid perturbation Galerkin technique with applications to slender body theory

A two step hybrid perturbation-Galerkin method to solve a variety of applied mathematics problems which involve a small parameter is presented. The method consists of: (1) the use of a regular or singular perturbation method to determine the asymptotic expansion of the solution in terms of the small parameter; (2) construction of an approximate solution in the form of a sum of the perturbation coefficient functions multiplied by (unknown) amplitudes (gauge functions); and (3) the use of the classical Bubnov-Galerkin method to determine these amplitudes. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is applied to some singular perturbation problems in slender body theory. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the degree of applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

NASA NTRS1987

An eigensystem realization algorithm using data correlations (ERA/DC) for modal parameter identification

A modification to the Eigensystem Realization Algorithm (ERA) for modal parameter identification is presented in this paper. The ERA minimum order realization approach using singular value decomposition is combined with the philosophy of the Correlation Fit method in state space form such that response data correlations rather than actual response values are used for modal parameter identification. This new method, the ERA using data correlations (ERA/DC), reduces bias errors due to noise corruption significantly without the need for model overspecification. This method is tested using simulated five-degree-of-freedom system responses corrupted by measurement noise. It is found for this case that, when model overspecification is permitted and a minimum order solution obtained via singular value truncation, the results from the two methods are of similar quality.

Juang, Jer-Nan↗

NASA NTRS1987

Viscous damping of Alfven normal modes in non-uniform plasmas

The Alfven normal mode in a nonuniform, viscous plasma is investigated. Because the ideal hydromagnetic equation is singular in a nonuniform magnetic field, viscosity is included in order to regularize the equation, analogous to the removal of the singular point by resistivity (Mok and Einaudi, 1985). The eigenvalue equation is then solved numerically for a particular magnetic configuration. The real part of the frequency of this viscous normal mode is found to be similar to the one in the resistive case, while the damping is shown to be comparable, and sometimes to exceed, the resistive effect under certain conditions in which the Lundquist number is sufficiently large. The damping rate is evaluated for various plasma conditions corresponding to different parts of the solar atmosphere. The effects of viscosity are found to dominate resistivity in the quiet sun corona and solar wind.

Mok, Y.↗

NASA NTRS1987

Perturbation analysis of internal balancing for lightly damped mechanical systems with gyroscopic and circulatory forces

Approximate expressions are developed for internally balanced singular values corresponding to the modes of mechanical systems with gyroscopic forces, light damping, and small circulatory forces. A brief overview is first given of the balanced realization model reduction method, including a discussion of recent work. The models considered are defined, and a perturbation analysis is used to show that the modal representation becomes asymptotically balanced as damping reduces to zero. The approximate balanced singular values are calculated, and a simple example of a flexible, dual-spin spacecraft is given as an illustration of the results.

Blelloch, P. A.↗