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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Digital imaging techniques in experimental stress analysis

Digital imaging techniques are utilized as a measure of surface displacement components in laser speckle metrology. An image scanner which is interfaced to a computer records and stores in memory the laser speckle patterns of an object in a reference and deformed configuration. Subsets of the deformed images are numerically correlated with the references as a measure of surface displacements. Discrete values are determined around a closed contour for plane problems which then become input into a boundary integral equation method in order to calculate surface traction in the contour. Stresses are then calculated within this boundary. The solution procedure is illustrated by a numerical example of a case of uniform tension.

Peters, W. H.↗

NASA NTRS1998

The Radiated Field Generated by a Monopole Source in a Short, Rigid, Rectangular Duct

This thesis develops a method to model the acoustic field generated by a monopole source placed in a moving rectangular duct. The walls of the duct are assumed to be infinitesimally thin and the source is placed at the center of the duct. The total acoustic pressure is written in terms of the free-space pressure, or incident pressure, and the scattered pressure. The scattered pressure is the augmentation to the incident pressure due to the presence of the duct. It satisfies a homogeneous wave equation and is discontinuous across the duct walls. Utilizing an integral representation of the scattered pressure, a set of singular boundary integral equations governing the unknown jump in scattered pressure is derived. This equation is solved by the method of collocation after representing the jump in pressure as a double series of shape functions. The solution obtained is then substituted back into the integral representation to determine the scattered pressure, and the total acoustic pressure at any point in the field. A few examples are included to illustrate the influence of various geometric and kinematic parameters on the radiated sound field.

Lakota, Barbara Anne↗

DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

Long-distance sound propagation over discontinuous impedances

A calculation method is presented for sound propagation over an impedance discontinuity in flat ground with a homogeneous, still atmosphere. The method is based on an approximate solution to a two dimensional boundary integral equation formulation of the problem, which expresses the wave field as the solution for homogeneous ground plus an integral over half of the boundary. Through recognizing this integral as a generalized Fourier integral, asymptotic methods are applied to evaluate the part of the integral most expensive to compute by numerical quadrature. Single frequency excess attenuation results for propagation from a point source above rigid ground to a receiver above absorbing ground are discussed. The results are applied, with air attenuation and A-weighting, to a notional jet engine noise source. Simple trends are noted.

Chandler-Wilde, Simon N.↗

NASA NTRS2000

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗

NASA NTRS1991

Hydrodynamic behavior and electrochemical impedance of the Hanging Meniscus Rotating Disk (HMRD) electrode. I - Meniscus shape under rotation. II - I-BIEM calculations of frequency dispersion and minimization

The shape equations for an HMRD in static and rotating configurations are developed and solved numerically. A rationale for the applicability of the standard Levich equations to the rotating case is given. The region of stability of the HMRD is examined, and the observed small negative intercept for a Levich plot is explained. The iterative boundary integral equation method is applied to the problem of frequency dispersion at an HMRD electrode. It is shown that a range of disk sizes and heights can be chosen to give almost uniform primary and secondary current distribution and minimal frequency disperison.

Cahan, Boris D.↗

NASA NTRS1988

Parameter estimation technique for boundary value problems by spline collocation method

A parameter-estimation technique for boundary-integral equations of the second kind is developed. The output least-squares identification technique using the spline collocation method is considered. The convergence analysis for the numerical method is discussed. The results are applied to boundary parameter estimations for two-dimensional Laplace and Helmholtz equations.

Kojima, Fumio↗

NASA NTRS1992

Boundary formulations for shape sensitivity of temperature dependent conductivity problems

Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.

Kane, James H.↗

NASA NTRS1991

I-BIEM calculations of the frequency dispersion and AC current distribution at disk and ring-disk electrodes

The Iterative Boundary Integral Equation Method (I-BIEM) has been applied to the problem of frequency dispersion at a disk electrode in a finite geometry. The I-BIEM permits the direct evaluation of the AC potential (a complex variable) using complex boundary conditions. The point spacing was made highly nonuniform, to give extremely high resolution in those regions where the variables change most rapidly, i.e., in the vicinity of the edge of the disk. Results are analyzed with respect to IR correction, equipotential surfaces, and reference electrode placement. The current distribution is also examined for a ring-disk configuration, with the ring and the disk at the same AC potential. It is shown that the apparent impedance of the disk is inductive at higher frequencies. The results are compared to analytic calculations from the literature, and usually agree to better than 0.001%.

Boris D Cahan↗

NASA NTRS1991

I-BIEM calculations of the frequency dispersion and ac current distribution at disk and ring-disk electrodes

Cahan, Boris D.↗

NASA NTRS1979

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The stress intensity factors were calculated by a three dimensional finite element method. The finite element models employed singularity elements along the crack front and linear strain elements elsewhere. The models had about 6500 degrees of freedom. The stress intensity factors were evaluated from a nodal force method. An equation for the stress intensity factors was obtained form the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary integral equation method were in agreement (+ or - 2 percent) and those from another finite element were in fair agreement (+ or - 8 percent) with the present results.

Newman, J. C., Jr.↗

NASA NTRS1980

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The paper presents stress-intensity factors for a wide range of semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 to to 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method using singularity elements along the crack front and linear-strain elements elsewhere. An equation for the stress-intensity factors was obtained which applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was also made between the results and other analyses of internal surface cracks in cylinders. The results from a boundary-integral equation method were in agreement and those from another finite-element method were in fair agreement (+ or - 8 percent) with the results.

Newman, J. C., Jr.↗

NASA NTRS1992

Shape reanalysis and sensitivities utilizing preconditioned iterative boundary solvers

The computational advantages associated with the utilization of preconditined iterative equation solvers are quantified for the reanalysis of perturbed shapes using continuum structural boundary element analysis (BEA). Both single- and multi-zone three-dimensional problems are examined. Significant reductions in computer time are obtained by making use of previously computed solution vectors and preconditioners in subsequent analyses. The effectiveness of this technique is demonstrated for the computation of shape response sensitivities required in shape optimization. Computer times and accuracies achieved using the preconditioned iterative solvers are compared with those obtained via direct solvers and implicit differentiation of the boundary integral equations. It is concluded that this approach employing preconditioned iterative equation solvers in reanalysis and sensitivity analysis can be competitive with if not superior to those involving direct solvers.

Guru Prasad, K.↗

NASA NTRS1991

Transient thermoelasticity and other body force effects in boundary element shape sensitivity analysis

Structural shape design sensitivity analysis (DSA) is presented for objects subjected to gravity, centrifugal, and general thermal loading. DSA employs implicit differentiation of the governing boundary integral equations. Both a surface integral and a particular integral approach are discussed for computation of these effects in response analysis and sensitivity analysis. It is shown that the particular integral approach is capable of computing the sensitivities of objects subjected to nonharmonic temperature distributions. Attention is also given to several techniques for the modeling of temperature distributions in the DSA process. In particular, an approach that eliminates the need for the incorporation of thermal response analysis and sensitivity analysis is described in the overall shape optimization process. A sharp corner formulation is also presented that can account for the jump in the surface traction components at a sharp corner in an interzone boundary.

Kane, James H.↗

NASA NTRS1977

A study of damage zones or characteristic lengths as related to the fracture behavior of graphite/epoxy laminates

Uniaxial tensile tests conducted on a variety of graphite/epoxy laminates, containing narrow rectangular slits, square or circular holes with various aspect ratios are discussed. The techniques used to study stable crack or damage zone growth--namely, birefringence coatings, COD gages, and microscopic observations are discussed. Initial and final fracture modes are discussed as well as the effect of notch size and shape, and laminate type on the fracture process. Characteristic lengths are calculated and compared to each other using the point, average and inherent flaw theories. Fracture toughnesses are calculated by the same theories and compared to a boundary integral equation technique. Finite width K-calibration factors are also discussed.

Yeow, Y. T.↗

NASA NTRS1980

Measurements of mixed-mode crack surface displacements and comparison with theory

The problem of a finite-width tension specimen containing a crack oriented at various angles to the load axis is attacked from experimental and theoretical viewpoints. Displacements of an electro-machined slot, 12.5 mm long and oriented at angles of 0, 15, 30, 45, 60, and 75 deg, are measured using a laser-based in-plane measuring technique. Various width specimens, ranging from a crack-length/width ratio of 0.167 to 0.794, are tested. A boundary-integral equation method is extended to deal with the presence of a sharp crack. Agreement between the two approaches is generally good except near the tips of the cracks.

Sharpe, W. N., Jr.↗

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

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗