SEARCH · Search NASA

Results for “SINGULAR INTEGRAL EQUATION”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Interface crack in a nonhomogeneous elastic medium

The linear elasticity problem for an interface crack between two bonded half planes is reconsidered. It is assumed that one of the half planes is homogeneous and the second is nonhomogeneous in such a way that the elastic properties are continuous throughout the plane and have discontinuous derivatives along the interface. The problem is formulated in terms of a system of integral equations and the asymptotic behavior of the stress state near the crack tip is determined. The results lead to the conclusion that the singular behavior of stresses in the nonhomogeneous medium is identical to that in a homogeneous material provided the spacial distribution of material properties is continuous near and at the crack tip. The problem is solved for various values of the nonhomogeneity parameter and for four different sets of crack surface tractions, and the corresponding stress intensity factors are tabulated.

Delale, F.↗

NASA NTRS1976

Arrangement of vortex lattices

A new method is developed for solving the lifting-surface equation for thin wings. The solution requires the downwash equation to be in the form of Cauchy integrals which can be interpreted as a vortex lattice with the positions of the vortices and control points dictated by the finite sum used to approximate the integrals involved. Lan's continuous loading method is employed for the chordwise integral since it properly accounts for the leading-edge singularity, Cauchy singularity, and Kutta condition. Unlike Lan, the spanwise loading is also continuous and the Cauchy singularity in the spanwise integral is also properly accounted for by using the midpoint trapezoidal rule and the theory of Chebyshev polynomials. This technique yields the exact classical solution to Prandtl's lifting-line equation.The solution to the lifting-surface equation for rectangular wings was found to compare well with other continuous loading methods, but with much smaller computational times, and to converge faster than other vortex lattice methods.

Dejarnette, F. R.↗

NASA NTRS1987

Three dimensional numerical scattering from ocean-like surfaces

A new method of calculating the electric field scattered from two-dimensional conducting rough surfaces is presented. This method uses the operator conjugate gradients technique and the extended boundary condition equation. This avoids the build-up of round-off errors encountered in the inversion of large matrices. It also avoids the singularities of the integral kernels usually encountered in the method of moments. Techniques for improving the rate of convergence by applying physical constraints are discussed. This method is presently being used to calculate the polarization signatures from simulated ocean-like surfaces.

Rodriguez, Ernesto↗

NASA NTRS1987

Application of a sensitivity analysis technique to high-order digital flight control systems

A sensitivity analysis technique for multiloop flight control systems is studied. This technique uses the scaled singular values of the return difference matrix as a measure of the relative stability of a control system. It then uses the gradients of these singular values with respect to system and controller parameters to judge sensitivity. The sensitivity analysis technique is first reviewed; then it is extended to include digital systems, through the derivation of singular-value gradient equations. Gradients with respect to parameters which do not appear explicitly as control-system matrix elements are also derived, so that high-order systems can be studied. A complete review of the integrated technique is given by way of a simple example: the inverted pendulum problem. The technique is then demonstrated on the X-29 control laws. Results show linear models of real systems can be analyzed by this sensitivity technique, if it is applied with care. A computer program called SVA was written to accomplish the singular-value sensitivity analysis techniques. Thus computational methods and considerations form an integral part of many of the discussions. A user's guide to the program is included. The SVA is a fully public domain program, running on the NASA/Dryden Elxsi computer.

Paduano, James D.↗

NASA NTRS1983

Derivation of the fundamental equation of sound generated by moving aerodynamic surfaces

Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.

Aggarwal, H. R.↗

NASA NTRS1973

An improved method for the aerodynamic analysis of wing-body-tail configurations in subsonic and supersonic flow. Part 1: Theory and application

A new method has been developed for calculating the pressure distribution and aerodynamic characteristics of wing-body-tail combinations in subsonic and supersonic potential flow. A computer program has been developed to perform the numerical calculations. The configuration surface is subdivided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. The normal components of velocity induced at specified control points by each singularity distribution are calculated and make up the coefficients of a system of linear equations relating the strengths of the singularities to the magnitude of the normal velocities. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations using an interactive procedure. Once the singularity strengths are known, the pressure coefficients are calculated, and the forces and moments acting on the configuration determined by numerical integration.

Woodward, F. A.↗

NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗

NASA NTRS1972

A method for integrating the boundary-layer equations through a region of reverse flow.

Development of a procedure for numerically integrating the boundary-layer equations through a region of reverse flow which takes downstream influence into account. This method is applied to the problem of uniform flow past a parallel flat plate of finite length whose surface has a constant velocity directed opposite to that of the main stream. Although singularities occur at both the point of detachment (x sub s) and reattachment (x sub r) of the psi = 0 streamline, this integration technique provides a solution which ceases to apply only in the close proximity of these singular points. From this solution it is evident that, throughout a large portion of the separated region, the flow is strongly affected by conditions near x sub r, thereby demonstrating the importance of allowing information to be transmitted upstream in a region of backflow. Near x sub s, however, it is found that, in spite of the presence of reverse flow, the solution has a self-similar form in this particular example.

Klemp, J. B.↗

NASA NTRS1977

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗

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 NTRS1972

Stresses in bonded materials with a crack perpendicular to the interface.

The problem of two elastic bonded half planes containing a crack perpendicular to the interface is considered. First the solution for the semi-infinite crack under concentrated wedge loading is given. Then the problem of a finite crack fully imbedded in one of the half planes of terminating at the interface is considered. The Mellin transform in conjunction with the dislocations is used to formulate the problem and to derive the integral equation. The integral equation is solved; the stress intensity factors, the crack surface displacement, and the stresses around the crack tip terminating at the interface are obtained. Noting that the power of the singularity for the interface tip of the crack is not -1/2, a tentative fracture criterion dealing with the initiation of fracture propagation is proposed.

Cook, T. S.↗

NASA NTRS1996

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗

NASA NTRS1987

A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach

A new formulation of electromagnetic wave scattering by convex, two-dimensional conducting bodies is reported. This formulation, called the on-surface radiation condition (OSRC) approach, is based upon an expansion of the radiation condition applied directly on the surface of a scatterer. It is now shown that application of a suitable radiation condition directly on the surface of a convex conducting scatterer can lead to substantial simplification of the frequency-domain integral equation for the scattered field, which is reduced to just a line integral. For the transverse magnetic case, the integrand is known explicitly. For the transverse electric case, the integrand can be easily constructed by solving an ordinary differential equation around the scatterer surface contour. Examples are provided which show that OSRC yields computed near and far fields which approach the exact results for canonical shapes such as the circular cylinder, square cylinder, and strip. Electrical sizes for the examples are ka = 5 and ka = 10. The new OSRC formulation of scattering may present a useful alternative to present integral equation and uniform high-frequency approaches for convex cylinders larger than ka = 1. Structures with edges or corners can also be analyzed, although more work is needed to incorporate the physics of singular currents at these discontinuities. Convex dielectric structures can also be treated using OSRC.

Kriegsmann, Gregory A.↗

NASA NTRS1982

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

The crack problem for a nonhomogeneous plane

Delale, F.↗

NASA NTRS1992

Confining potential in momentum space

A method is presented for the solution in momentum space of the bound state problem with a linear potential in r space. The potential is unbounded at large r leading to a singularity at small q. The singularity is integrable, when regulated by exponentially screening the r-space potential, and is removed by a subtraction technique. The limit of zero screening is taken analytically, and the numerical solution of the subtracted integral equation gives eigenvalues and wave functions in good agreement with position space calculations.

Norbury, John W.↗

NASA NTRS1992

Confining potential in momentum space

Norbury, John W.↗

NASA NTRS1975

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗