Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary 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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

I-BIEM, an iterative boundary integral equation method for computer solutions of current distribution problems with complex boundaries: A new algorithm. I - Theoretical

A new algorithm for an iterative computation of solutions of Laplace's or Poisson's equations in two dimensions, using Green's second identity, is presented. This algorithm converges strongly and geometrically and can be applied to curved, irregular, or moving boundaries with nonlinear and/or discontinuous boundary conditions. It has been implemented in Pascal on a number of micro- and minicomputers and applied to several geometries. Cases with known analytic solutions have been tested. Convergence to within 0.1 percent to 0.01 percent of the theoretical values are obtained in a few minutes on a microcomputer.

Cahan, B. D.↗
NASA NTRS1971

Formulation of boundary integral equations for three-dimensional elasto-plastic flow.

A general theory of elastoplastic flow is formulated for work-hardening materials that may be both anisotropic and compressible. Because the theory is quasi-linear, it may be cast in terms of integral equations and the result is an extended form of Somigliana's identity. When these relations are evaluated on the boundary of a solid, their dimensionality is reduced. Previous experience with simpler materials shows that arbitrary problems may be solved in a direct manner.

Swedlow, J. L.↗
NASA NTRS1982

Experimental boundary integral equation applications in speckle interferometry

A complete data analysis system utilizing laser speckle interferometry and Experimental Integral Equation Method (EBIE) is described. A host computer provides the optical data analysis and serves as an input device to a PDP11/VAX computer for numerical analysis. The basic theory of pointwise filtering and digital correlation is described as experimental data input to the BIE method.

Fraley, J. E.↗
NASA NTRS1996

Singularity Preserving Numerical Methods for Boundary Integral Equations

In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the following five projects: a note on the finite element method with singular basis functions; numerical quadrature for weakly singular integrals; superconvergence of degenerate kernel method; superconvergence of the iterated collocation method for Hammersteion equations; and singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. This final report consists of five papers describing these projects. Each project is preceeded by a brief abstract.

Kaneko, Hideaki↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
NASA NTRS1982

Numerical analysis of Weyl's method for integrating boundary layer equations

A fast method for accurate numerical integration of Blasius equation is proposed. It is based on the limit interchange in Weyl's fixed point method formulated as an iterated limit process. Each inner limit represents convergence to a discrete solution. It is shown that the error in a discrete solution admits asymptotic expansion in even powers of step size. An extrapolation process is set up to operate on a sequence of discrete solutions to reach the outer limit. Finally, this method is extended to related boundary layer equations.

Najfeld, I.↗
NASA NTRS1996

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗
NASA NTRS1986

Analysis of mixed-mode crack propagation using the boundary integral method

Crack propagation in a rotating inner raceway of a high speed roller bearing is analyzed using the boundary integral equation method. The method consists of an edge crack in a plate under tension, upon which varying Hertzian stress fields are superimposed. A computer program for the boundary integral equation method was written using quadratic elements to determine the stress and displacement fields for discrete roller positions. Mode I and Mode II stress intensity factors and crack extension forces G sub 00 (energy release rate due to tensile opening mode) and G sub r0 (energy release rate due to shear displacement mode) were computed. These calculations permit determination of that crack growth angle for which the change in the crack extension forces is maximum. The crack driving force was found to be the alternating mixed-mode loading that occurs with each passage of the most heavily loaded roller. The crack is predicted to propagate in a step-like fashion alternating between radial and inclined segments, and this pattern was observed experimentally. The maximum changes DeltaG sub 00 and DeltaG sub r0 of the crack extension forces are found to be good measures of the crack propagation rate and direction.

Mendelson, A.↗
NASA NTRS1991

A statistical approach to the brittle fracture of a multi-phase solid

A stochastic damage model is proposed to quantify the inherent statistical distribution of the fracture toughness of a brittle, multi-phase solid. The model, based on the macrocrack-microcrack interaction, incorporates uncertainties in locations and orientations of microcracks. Due to the high concentration of microcracks near the macro-tip, a higher order analysis based on traction boundary integral equations is formulated first for an arbitrary array of cracks. The effects of uncertainties in locations and orientations of microcracks at a macro-tip are analyzed quantitatively by using the boundary integral equations method in conjunction with the computer simulation of the random microcrack array. The short range interactions resulting from surrounding microcracks closet to the main crack tip are investigated. The effects of microcrack density parameter are also explored in the present study. The validity of the present model is demonstrated by comparing its statistical output with the Neville distribution function, which gives correct fits to sets of experimental data from multi-phase solids.

Liu, W. K.↗
NASA NTRS1993

A displacement gradient BEM for accurate stress computation near boundaries in 2-D anisotropic problems

A displacement gradient method for 2D anisotropic elasticity problems is presented, which effectively minimizes the boundary layer effect through a two-step procedure. First, the boundary integral equations are solved for the unknown boundary displacements and tractions. Second, a direct integral equation for displacement gradients is developed in terms of boundary tractions. Three methods based on different evaluation procedures and locations for determining the displacement gradients are proposed. In the first method the displacement gradients are averaged at nodes common to adjacent elements. The second method stores the gradients element-wise. In the third method, the gradients are evaluated at the nodes of discontinuous elements. The three methods are applied to near-isotropic plates with circular and elliptic cutouts. It is concluded that all three methods can yield accurate stress distributions.

Sistla, R.↗
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Mager, Arthur↗
NASA NTRS1974

Application of boundary integral method to elastoplastic analysis of v-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problems. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single edge notch and subjected to pure bending. The boundary integral equation method reduced the non-homogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with the method of successive elastic solutions.

Rzasnicki, W.↗
NASA NTRS1975

Application of boundary integral method to elastoplastic analysis of V-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problems. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single edge notch and subjected to pure bending. The boundary integral equation method reduced the nonhomogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with the method of successive elastic solutions.

Rzasnicki, W.↗
NASA NTRS1974

Application of boundary integral method to elastoplastic analysis of V-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problem. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single-edge notch and subjected to pure bending. The boundary integral equation method reduced the inhomogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with a method of successive elastic solutions.

Rzasnicki, W.↗
NASA NTRS1990

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.↗