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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

The dependence of the shape and stability of captive rotating drops on multiple parameters

Asymptotic and numerical techniques in bifurcation theory are applied to the Young-Laplace equation governing meniscus shape in order to analyze the dependence of the shape and stability of rigidly rotating drops held captive between corotating solid faces on multiple parameters. Asymptotic analysis of the evolution of drop shape from the cylindrical as a function of distance between the solid faces, drop volume, rotational Bond number and gravitational Bond number shows that some shape bifurcations from cylinders to wavy, axisymmetric menisci are ruptured by small changes in drop volume or gravity. Computer calculations of axisymmetric drop shapes based on a finite element representation of the interface and numerical algorithms for tracking shape families and singular points are then used to map drop stability for the four-dimensional parameter space. The results of the asymptotic and numerical analyses are shown to agree well within the limited range of parameters where the asymptotic analysis is valid.

Ungar, L. H.↗
NASA NTRS1983

The aerodynamics of propellers and rotors using an acoustic formulation in the time domain

The aerodynamics of propellers and rotors is especially complicated because of the highly three-dimensional and compressible nature of the flow field. However, in linearized theory the problem is governed by the wave equation, and a numerically-efficient integral formulation can be derived. This reduces the problem from one in space to one over a surface. Many such formulations exist in the aeroacoustics literature, but these become singular integral equations if one naively tries to use them to predict surface pressures, i.e., for aerodynamics. The present paper illustrates how one must interpret these equations in order to obtain nonambiguous results. After the regularized form of the integral equation is derived, a method for solving it numerically is described. This preliminary computer code uses Legendre-Gaussian quadrature to solve the equation. Numerical results are compared to experimental results for ellipsoids, wings, and rotors, including effects due to lift. Compressibility and the farfield boundary conditions are satisfied automatically using this method.

Long, L. N.↗
NASA NTRS1983

Three-dimensional algebraic grid generation

One of the major approaches to numerical grid generation is the explicit algebraic expression of a physical grid as a function of a uniform grid in a rectangular computational coordinate system. The algebraic methods are based on mathematical interpolation, and the primary advantages are speed and directness. The relation between interpolation and grid generation is described. For three-dimensional grid generation, transfinite interpolation using the coordinate control processes developed in the multisurface method and two-boundary technique are advocated. Grid singularities encountered in three dimensions are discussed, and the exploration of multiple overlapping grids is proposed. Some aspects of interactive algebraic grid computation in three dimensions are discussed.

Smith, R. E.↗
NASA NTRS1983

Stress intensity factors in two bonded elastic layers containing cracks perpendicular to and on the interface. I Analysis. II - Solution and results

The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled systems of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered. Previously announced in STAR as N80-18428 and N80-18429

Lu, M.-C.↗
NASA NTRS1983

Numerical evaluation of the radiation from unbaffled, finite plates using the FFT

An iteration technique is described which numerically evaluates the acoustic pressure and velocity on and near unbaffled, finite, thin plates vibrating in air. The technique is based on Rayleigh's integral formula and its inverse. These formulas are written in their angular spectrum form so that the fast Fourier transform (FFT) algorithm may be used to evaluate them. As an example of the technique the pressure on the surface of a vibrating, unbaffled disk is computed and shown to be in excellent agreement with the exact solution using oblate spheroidal functions. Furthermore, the computed velocity field outside the disk shows the well-known singularity at the rim of the disk. The radiated fields from unbaffled flat sources of any geometry with prescribed surface velocity may be evaluated using this technique. The use of the FFT to perform the integrations in Rayleigh's formulas provides a great savings in computation time compared with standard integration algorithms, especially when an array processor can be used to implement the FFT.

Williams, E. G.↗
NASA NTRS1983

Separation of time scales in aircraft trajectory optimization

Two methods for analyzing the time-scale properties of aircraft trajectory optimization problems are presented. Time-scale properties must be identified before solutions can be obtained by using singular perturbation methods. Both methods only require a knowledge of the state equations, the aircraft characteristics, and the bounds on the state and control variables. Although these methods give only rough estimates of time-scale separation, they do not require that an 'exact' optimal trajectory be known, as do the more rigorous methods, and they are an improvement on the ad hoc methods currently in use. The two methods are applied to an example problem for a high performance aircraft.

Ardema, M. D.↗
NASA NTRS1983

Fracture mechanics for delamination problems in composite materials

A fracture mechanics approach to the well-known delamination problem in composite materials is presented. Based on the theory of anisotropic laminate elasticity and interlaminar fracture mechanics concepts, the composite delamination problem is formulated and solved. The exact order of the delamination crack-tip stress singularity is determined. Asymptotic stress and displacement fields for an interlaminar crack are obtained. Fracture mechanics parameters such as mixed-mode stress intensity factors, KI, KII, KIII, and the energy release rate, G, for composite delamination problems are defined. To illustrate the fundamental nature of the delamination crack behavior, solutions for edge-delaminated graphite-epoxy composites under uniform axial extension are presented. Effects of fiber orientation, ply thickness, and delamination length on the interlaminar fracture are examined.

Wang, S. S.↗
NASA NTRS1983

The stress intensity factor for the double cantilever beam

Fourier transforms and the Wiener-Hopf technique are used in conjunction with plane elastostatics to examine the singular crack tip stress field in the double cantilever beam (DCB) specimen. In place of the Dirac delta function, a family of functions which duplicates the important features of the concentrated forces without introducing unmanageable mathematical complexities is used as a loading function. With terms of order h-squared/a-squared retained in the series expansion, the dimensionless stress intensity factor is found to be K (h to the 1/2)/P = 12 to the 1/2 (a/h + 0.6728 + 0.0377 h-squared/a-squared), in which P is the magnitude of the concentrated forces per unit thickness, a is the distance from the crack tip to the points of load application, and h is the height of each cantilever beam. The result is similar to that obtained by Gross and Srawley by fitting a line to discrete results from their boundary collocation analysis.

Fichter, W. B.↗
NASA NTRS1982

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three-dimensional flows having thin shear layers are analyzed, using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamic properties, such as pressure, shear stress, and vorticity on the body surface, are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is discussed and demonstrated.

Deiwert, G. S.↗
NASA NTRS1983

The crack problem for a nonhomogeneous plane

The plane elasticity problem for a nonhomogeneous medium containing a crack is considered. It is assumed that the Poisson's ratio of the medium is constant and the Young's modulus E varies exponentially with the coordinate parallel to the crack. First the half plane problem is formulated and the solution is given for arbitrary tractions along the boundary. Then the integral equation for the crack problem is derived. It is shown that the integral equation having the derivative of the crack surface displacement as the density function has a simple Cauchy type kernel. Hence, its solution and the stresses around the crack tips have the conventional square root singularity. The solution is given for various loading conditions. The results show that the effect of the Poisson's ratio and consequently that of the thickness constraint on the stress intensity factors are rather negligible.

Delale, F.↗
NASA NTRS1983

Rain estimation from satellites: An examination of the Griffith-Woodley technique

The Griffith-Woodley Technique (GWT) is an approach to estimating precipitation using infrared observations of clouds from geosynchronous satellites. It is examined in three ways: an analysis of the terms in the GWT equations; a case study of infrared imagery portraying convective development over Florida; and the comparison of a simplified equation set and resultant rain map to results using the GWT. The objective is to determine the dominant factors in the calculation of GWT rain estimates. Analysis of a single day's convection over Florida produced a number of significant insights into various terms in the GWT rainfall equations. Due to the definition of clouds by a threshold isotherm the majority of clouds on this day did not go through an idealized life cycle before losing their identity through merger, splitting, etc. As a result, 85% of the clouds had a defined life of 0.5 or 1 h. For these clouds the terms in the GWT which are dependent on cloud life history become essentially constant. The empirically derived ratio of radar echo area to cloud area is given a singular value (0.02) for 43% of the sample, while the rainrate term is 20.7 mmh-1 for 61% of the sample. For 55% of the sampled clouds the temperature weighting term is identically 1.0. Cloud area itself is highly correlated (r=0.88) with GWT computed rain volume. An important, discriminating parameter in the GWT is the temperature defining the coldest 10% cloud area. The analysis further shows that the two dominant parameters in rainfall estimation are the existence of cold cloud and the duration of cloud over a point.

Negri, A. J.↗
NASA NTRS1983

Kinematic equations for resolved-rate control of an industrial robot arm

An operator can use kinematic, resolved-rate equations to dynamically control a robot arm by watching its response to commanded inputs. Known resolved-rate equations for the control of a particular six-degree-of-freedom industrial robot arm and proceeds to simplify the equations for faster computations are derived. Methods for controlling the robot arm in regions which normally cause mathematical singularities in the resolved-rate equations are discussed.

Barker, L. K.↗
NASA NTRS1983

Solution of Poisson equations for 3-dimensional grid generations

A method for generating three dimensional, finite difference grids about complicated geometries by using Poisson equations is developed. The inhomogenous terms are automatically chosen such that orthogonality and spacing restrictions at the body surface are satisfied. Spherical variables are used to avoid the axis singularity, and an alternating-direction-implicit (ADI) solution scheme is used to accelerate the computations. Computed results are presented that show the capability of the method. Since most of the results presented have been used as grids for flow-field computations, this is indicative that the method is a useful tool for generating three-dimensional grids about complicated geometries.

Fujii, K.↗
NASA NTRS1983

The compressible aerodynamics of rotating blades based on an acoustic formulation

An acoustic formula derived for the calculation of the noise of moving bodies is applied to aerodynamic problems. The acoustic formulation is a time domain result suitable for slender wings and bodies moving at subsonic speeds. A singular integral equation is derived in terms of the surface pressure which must then be solved numerically for aerodynamic purposes. However, as the 'observer' is moved onto the body surface, the divergent integrals in the acoustic formulation are semiconvergent. The procedure for regularization (or taking principal values of divergent integrals) is explained, and some numerical examples for ellipsoids, wings, and lifting rotors are presented. The numerical results show good agreement with available measured surface pressure data.

Long, L. N.↗
NASA NTRS1983

On-board near-optimal climb-dash energy management

On-board real time flight control is studied in order to develop algorithms which are simple enough to be used in practice, for a variety of missions involving three dimensional flight. The intercept mission in symmetric flight is emphasized. Extensive computation is required on the ground prior to the mission but the ensuing on-board exploitation is extremely simple. The scheme takes advantage of the boundary layer structure common in singular perturbations, arising with the multiple time scales appropriate to aircraft dynamics. Energy modelling of aircraft is used as the starting point for the analysis. In the symmetric case, a nominal path is generated which fairs into the dash or cruise state.

Weston, A. R.↗
NASA NTRS1983

Optimal symmetric flight with an intermediate vehicle model

Optimal flight in the vertical plane with a vehicle model intermediate in complexity between the point-mass and energy models is studied. Flight-path angle takes on the role of a control variable. Range-open problems feature subarcs of vertical flight and singular subarcs. The class of altitude-speed-range-time optimization problems with fuel expenditure unspecified is investigated and some interesting phenomena uncovered. The maximum-lift-to-drag glide appears as part of the family, final-time-open, with appropriate initial and terminal transient exceeding level-flight drag, some members exhibiting oscillations. Oscillatory paths generally fail the Jacobi test for durations exceeding a period and furnish a minimum only for short-duration problems.

Menon, P. K. A.↗
NASA NTRS1983

A cylindrical shell with a stress-free end which contains an axial part-through or through crack

The interaction problem of a through or a part through crack with a stress free boundary in a semi-infinite cylindrical shell is considered. It is assumed that the crack lies in a meridional plane which is a plane of symmetry with respect to the external loads as well as the geometry. The circular boundary of the semi-infinite cylinder is assumed to be stress free. By using a transverse shear theory the problem is formulated in terms of a system of singular integral equations. The line spring model is used to treat the part through crack problem. In the case of a through crack the interaction between the perturbed stress fields due to the crack and the free boundary is quite strong and there is a considerable increase in the stress intensity factors caused by the interaction. On the other hand in the problem of a surface crack the interaction appears to be much weaker and consequently the magnification in the stress intensity factors is much less significant.

Erdogan, F.↗
NASA NTRS1984

The crack-inclusion interaction problem

The general plane elastostatic problem of interaction between a crack and an inclusion is considered. The Green's functions for a pair of dislocations and a pair of concentrated body forces are used to generate the crack and the inclusion. Integral equations are obtained for a line crack and an elastic line inclusion having an arbitrary relative orientation and size. The nature of stress singularity around the end points of rigid and elastic inclusions is described and three special cases of this intersection problem are studied. The problem is solved for an arbitrary uniform stress state away from the crack-inclusion region. The nonintersecting crack-inclusion problem is considered for various relative size, orientation, and stiffness parameters, and the stress intensity factors at the ends of the inclusion and the crack are calculated. For the crack-inclusion intersection case, special stress intensity factors are defined and are calculated for various values of the parameters defining the relative size and orientation of the crack and the inclusion and the stiffness of the inclusion.

Xue-Hui, L.↗