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 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Modeling material failure with a vectorized routine

The computational aspects of modelling material failure in structural wood members are presented with particular reference to vector processing aspects. Wood members are considered to be highly orthotropic, inhomogeneous, and discontinuous due to the complex microstructure of wood material and the presence of natural growth characteristics such as knots, cracks and cross grain in wood members. The simulation of strength behavior of wood members is accomplished through the use of a special purpose finite element/fracture mechanics routine, program STARW (Strength Analysis Routine for Wood). Program STARW employs quadratic finite elements combined with singular crack tip elements in a finite element mesh. Vector processing techniques are employed in mesh generation, stiffness matrix formation, simultaneous equation solution, and material failure calculations. The paper addresses these techniques along with the time and effort requirements needed to convert existing finite element code to a vectorized version. Comparisons in execution time between vectorized and nonvectorized routines are provided.

Cramer, S. M.↗

NASA NTRS1984

A method for computing the leading-edge suction in a higher-order panel method

Experimental data show that the phenomenon of a separation induced leading edge vortex is influenced by the wing thickness and the shape of the leading edge. Both thickness and leading edge shape (rounded rather than point) delay the formation of a vortex. Existing computer programs used to predict the effect of a leading edge vortex do not include a procedure for determining whether or not a vortex actually exists. Studies under NASA Contract NAS1-15678 have shown that the vortex development can be predicted by using the relationship between the leading edge suction coefficient and the parabolic nose drag. The linear theory FLEXSTAB was used to calculate the leading edge suction coefficient. This report describes the development of a method for calculating leading edge suction using the capabilities of the higher order panel methods (exact boundary conditions). For a two dimensional case, numerical methods were developed using the double strength and downwash distribution along the chord. A Gaussian quadrature formula that directly incorporates the logarithmic singularity in the downwash distribution, at all panel edges, was found to be the best method.

Ehlers, F. E.↗

NASA NTRS1984

Feedback laws for fuel minimization for transport aircraft

The Theoretical Mechanics Branch has as one of its long-range goals to work toward solving real-time trajectory optimization problems on board an aircraft. This is a generic problem that has application to all aspects of aviation from general aviation through commercial to military. Overall interest is in the generic problem, but specific problems to achieve concrete results are examined. The problem is to develop control laws that generate approximately optimal trajectories with respect to some criteria such as minimum time, minimum fuel, or some combination of the two. These laws must be simple enough to be implemented on a computer that is flown on board an aircraft, which implies a major simplification from the two point boundary value problem generated by a standard trajectory optimization problem. In addition, the control laws allow for changes in end conditions during the flight, and changes in weather along a planned flight path. Therefore, a feedback control law that generates commands based on the current state rather than a precomputed open-loop control law is desired. This requirement, along with the need for order reduction, argues for the application of singular perturbation techniques.

Price, D. B.↗

NASA NTRS1984

Stress intensity factors in a reinforced thick-walled cylinder

An elastic thick-walled cylinder containing a radial crack is considered. It is assumed that the cylinder is reinforced by an elastic membrane on its inner surface. The model is intended to simulate pressure vessels with cladding. The formulation of the problem is reduced to a singular integral equation. Various special cases including that of a crack terminating at the cylinder-reinforcement interface are investigated and numerical examples are given. Results indicate that in the case of the crack touching the interface the crack surface displacement derivative is finite and consequently the stress state around the corresponding crack tip is bounded; and generally, for realistic values of the stiffness parameter, the effect of the reinforcement is not very significant.

Tang, R.↗

NASA NTRS1984

Development of an impact- and solvent-resistant thermoplastic composite matrix

Synthesis, moldability and chloroform, acetone and tricresyl phosphate resistance of 16 polymer compositions are described. These aliphatic heterocyclic polymers include polyimides, polybenzimidazoles, and N-arylenepolybenzimidazoles. A solution condensation (cresol) method to prepare imidized aliphaic polyimides is described. Two polyimides and one polybenzimidazole demonstrate no crazing or cracking during 500 hr exposure to the cited solvents under stress. Modification of one aliphatic polyimide with several aromatic amines suggests that m-phenylenediamine is singular in its behavior to improve the chloroform resistance of that class of polyimides.

Delano, C. B.↗

NASA NTRS1984

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.↗

NASA NTRS1983

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗

NASA NTRS1984

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗

NASA NTRS1984

Minimum-fuel, 3-dimensional flightpath guidance of transfer jets

Minimum fuel, three dimensional flightpaths for commercial jet aircraft are discussed. The theoretical development is divided into two sections. In both sections, the necessary conditions of optimal control, including singular arcs and state constraints, are used. One section treats the initial and final portions (below 10,000 ft) of long optimal flightpaths. Here all possible paths can be derived by generating fields of extremals. Another section treats the complete intermediate length, three dimensional terminal area flightpaths. Here only representative sample flightpaths can be computed. Sufficient detail is provided to give the student of optimal control a complex example of a useful application of optimal control theory.

Neuman, F.↗

NASA NTRS1983

A cylindrical shell with an arbitrarily oriented crack

The general problem of a shallow shell with constant curvatures is considered. It is assumed that the shell contains an arbitrarily oriented through crack and the material is specially orthotropic. The nonsymmetric problem is solved for arbitrary self equilibrating crack surface tractions, which, added to an appropriate solution for an uncracked shell, would give the result for a cracked shell under most general loading conditions. The problem is reduced to a system to five singular integral equations in a set of unknown functions representing relative displacements and rotations on the crack surfaces. The stress state around the crack tip is asymptotically analyzed and it is shown that the results are identical to those obtained from the two dimensional in plane and antiplane elasticity solutions. The numerical results are given for a cylindrical shell containing an arbitrarily oriented through crack. Some sample results showing the effect of the Poisson's ratio and the material orthotropy are also presented. Previously annunced in STAR as N83-16783

Yahsi, O. S.↗

NASA NTRS1984

Stress analysis of multilayered plates around circular holes

The elasticity problem for a thick plate consisting of two bonded dissimilar homogeneous layers is considered. It is assumed that the plate is infinite, contains a circular hole, and is subjected to axisymmetric external loads. A technique suitable to the formulation of three-dimensional boundary value problems is introduced and applied to the three-dimensional elasticity solution of an infinite plate with a hole. The effect of Poisson's ratio is studied in some detail, and the axisymmetric problem of a laminated plate is solved. The stresses at the interface are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singularity is presented and the stress intensity factors are computed. Some of the results are compared with those found from the finite element method. The agreement between the two sets with regard to trend and magnitude is found to be very satisfactory.

Delale, F.↗

NASA NTRS1982

Yet another method for triangulation and contouring for automated cartography

An algorithm is presented for hierarchical subdivision of a set of three-dimensional surface observations. The data structure used for obtaining the desired triangulation is also singularly appropriate for extracting contours. Some examples are presented, and the results obtained are compared with those given by Delaunay triangulation. The data points selected by the algorithm provide a better approximation to the desired surface than do randomly selected points.

De Floriani, L.↗

NASA NTRS1983

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.

Zebker, H. A.↗

NASA NTRS1983

Very high-resolution heat-capacity measurements near the lambda point of helium

New measurements of the heat capacity of a sample of helium 3-mm high are reported, which extend to within 5 x 10 to the -8th deg of the lambda transition at the vapor pressure. From an analysis of the results allowing for the effect of gravity, the values -0.0127 + or - 0.0026 (2 sigma) for the exponent alpha (= alpha-prime) and 1,058 + or - 0.004 for the leading singularity ratio A/A-prime are obtained. These values are in closer agreement with the theoretical predictions than those reported previously.

Lipa, J. A.↗

NASA NTRS1982

A general algorithm for solving the algebraic Riccati equation

The generalized eigenvalue problem provides a suitable framework for reliable solutions of many system theoretic, control, and estimation problems. A general algorithm for solving the matrix algebraic Riccati equation (ARE) which utilizes a pencil structure is described here. This algorithm avoids unnecessary inversion of cost or transition matrices, making it a numerically sound way to solve for the gains and/or ARE with singular quadratic costs, for cases satisfying detectability and stabilizability conditions. Examples are solution with discrete dead-beat control, noiseless measurements in Kalman filters and time-delays in discrete-time systems, which cause difficulties in the Hamiltonian standard eigenvalue problem formulation. The ARE algorithm implementatiton and numerical examples are shown.

Walker, R. A.↗

NASA NTRS1984

Condensation modes in magnetized cylindrical plasmas

In order to understand the formation and stability of solar prominences, attention is given to the condensation modes of the thermal instability in a cylindriical plasma whose magnetic field has both potential (longitudinal) and nonpotential (poloidal) components forming helically twisted field lines. It is noted that the field line twist has a significant effect on the stability of condensation modes which are unstable if field lines are straight, but become stable as the twist increases. The problem is treated in a fully self-consistent way, to derive a second-order ordinary differential equation for radiative MHD stability. The approximation made by neglecting inertial terms is valid for coronal loop conditions if there is no singular surface in the plasma, and a simplified differential equation is derived. Attention is given to stability for the m = zero, 1, and 2 modes numerically, for various loop parameters.

An, C.-H.↗

NASA NTRS1984

Strong fluctuation theory for electromagnetic wave scattering by a layer of random discrete scatterers

The strong fluctuation theory is applied to the study of electromagnetic wave scattering from a layer of random discrete scatterers. The singularity of the dyadic Green's function is taken into account in the calculation of the effective permittivity functions. The correlation functions for the random medium with different scatterer constituents and size distributions are derived. Applying the dyadic Green's function for a two-layer medium and using the bilocal and distorted Born approximations, the first and the second moments of the fields are then calculated. Both the backscattering and bistatic scattering coefficients are obtained, and the former is shown to match favorably with experimental data obtained from snow fields.

Jin, Y. Q.↗