Integral-differential equation of plates whose material behaves according to linear rheological relationships
Integral differential equations of plates whose materials exhibit linear rheological relationships
INTEGRAL EQUATION↗
SEARCH · Engineering Papers
Results for “integral equations”
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.
Laplace's equation and the Dirichlet-Neumann map in multiply connected domains
A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.
Greenbaum, A.↗
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.↗
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↗
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.↗
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.↗
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.↗
Deconvolution of wide-field-of-view satellite radiometer measurements of reflected solar radiation
An approach is presented by which the deconvolution of wide-field-of-view (WFOV) measurements of solar radiation reflected from earth is reduced from the solution of a two-dimensional integral equation to the solution of a set of one-dimensional integral equations. This reduces the storage needed for the required matrices by an order of magnitude and the computations by two orders of magnitude. Also, the theoretical and operational difficulties associated with the solution of ill-posed problems are greatly simplified by working with one-dimensional integral equations rather than with a two-dimensional integral equation.
Smith, G. L.↗
A Hybrid Symmetric FEM/MOM Formulation Applied to Scattering by Inhomogeneous Bodies of Revolution
A new symmetric formulation of the Hybrid Finite Element Method (HFEM) is described which combines elements of the Electric Field Integral Equation (EFIE) and the Magnetic Field Integral Equation (MFIE) for the exterior region along with the finite element solution for the interior problem. The formulation is applied to scattering by inhomogeneous bodies of revolution. To avoid spurious modes in the interior region a combination of vector and nodal based finite elements are used. Integral equations in the exterior region are used to enforce the Sommerfield radiation condition by matching both the tangential electric and magnetic fields between interior and exterior regions.
Lee, J-F.↗
Opening of an interface flaw in a layered elastic half-plane under compressive loading
A static analysis is given of the problem of an elastic layer perfectly bonded, except for a frictionless interface crack, to a dissimilar elastic half-plane. The free surface of the layer is loaded by a finite pressure distribution directly over the crack. The problem is formulated using the two dimensional linear elasticity equations. Using Fourier transforms, the governing equations are converted to a pair of coupled singular integral equations. The integral equations are reduced to a set of simultaneous algebraic equations by expanding the unknown functions in a series of Jacobi polynomials and then evaluating the singular Cauchy-type integrals. The resulting equations are found to be ill-conditioned and, consequently, are solved in the least-squares sense. Results from the analysis show that, under a normal pressure distribution on the free surface of the layer and depending on the combination of geometric and material parameters, the ends of the crack can open. The resulting stresses at the crack-tips are singular, implying that crack growth is possible. The extent of the opening and the crack-top stress intensity factors depend on the width of the pressure distribution zone, the layer thickness, and the relative material properties of the layer and half-plane.
Kennedy, J. M.↗
Methods of Attenuation Correction for Dual-Wavelength and Dual-Polarization Weather Radar Data
In writing the integral equations for the median mass diameter and number concentration, or comparable parameters of the raindrop size distribution, it is apparent that the forms of the equations for dual-polarization and dual-wavelength radar data are identical when attenuation effects are included. The differential backscattering and extinction coefficients appear in both sets of equations: for the dual-polarization equations, the differences are taken with respect to polarization at a fixed frequency while for the dual-wavelength equations, the differences are taken with respect to frequency at a fixed polarization. An alternative to the integral equation formulation is that based on the k-Z (attenuation coefficient-radar reflectivity factor) parameterization. This-technique was originally developed for attenuating single-wavelength radars, a variation of which has been applied to the TRMM Precipitation Radar data (PR). Extensions of this method have also been applied to dual-polarization data. In fact, it is not difficult to show that nearly identical equations are applicable as well to dualwavelength radar data. In this case, the equations for median mass diameter and number concentration take the form of coupled, but non-integral equations. Differences between this and the integral equation formulation are a consequence of the different ways in which attenuation correction is performed under the two formulations. For both techniques, the equations can be solved either forward from the radar outward or backward from the final range gate toward the radar. Although the forward-going solutions tend to be unstable as the attenuation out to the range of interest becomes large in some sense, an independent estimate of path attenuation is not required. This is analogous to the case of an attenuating single-wavelength radar where the forward solution to the Hitschfeld-Bordan equation becomes unstable as the attenuation increases. To circumvent this problem, the equations can be expressed in the form of a final-value problem so that the recursion begins at the far range gate and proceeds inward towards the radar. Solving the problem in this way traditionally requires estimates of path attenuation to the final gate: in the case of orthogonal linear polarizations, the attenuations at horizontal and vertical polarizations (same frequency) are required while in the dual-wavelength case, attenuations at the two frequencies (same polarization) are required.
Meneghini, R.↗
Green's function solution to heat transfer of a transparent gas through a tube
A heat transfer analysis of a transparent gas flowing through a circular tube of finite thickness is presented. This study includes the effects of wall conduction, internal radiative exchange, and convective heat transfer. The natural mathematical formulation produces a nonlinear, integrodifferential equation governing the wall temperature and an ordinary differential equation describing the gas temperature. This investigation proposes to convert the original system of equations into an equivalent system of integral equations. The Green's function method permits the conversion of an integrodifferential equation into a pure integral equation. The proposed integral formulation and subsequent computational procedure are shown to be stable and accurate.
Frankel, J. I.↗
General airfoil theory
On the assumption of infinitely small disturbances the author develops a generalized integral equation of airfoil theory which is applicable to any motion and compressible fluid. Successive specializations yield various simpler integral equations, such as Possio's, Birnbaum's, and Prandtl's integral equations, as well as new ones for the wing of infinite span with periodic downwash distribution and for the oscillating wing with high aspect ratio. Lastly, several solutions and methods for solving these integral equations are give.
Kussner, H G↗
The k-space formulation of the n-dimensional scattering problem
The n-dimensional scattering problem is solved by means of a k-space formulation of the field equations, thereby replacing the conventional integral equation formulation by a set of two algebraic equations in two unknowns in two spaces (the constitutive equation being an algebraic equation in x-space). These equations are solved by an iterative method with the aid of the fast Fourier transform (FFT) algorithm connecting the two spaces, requiring very simple initial approximations. Since algebraic and FFT equations are used, the number of arithmetic multiple-add operations and storage allocations required for a numerical solution are reduced from the order of N sq (for solving the matrix equations resulting from the conventional integral equations) to the order of N(log base 2 of N) and N, respectively (where N is the number of data points required for the specification of the problem). The advantage gained in speed and storage is thus of the order of N/log base 2 of N and N, respectively. This method is thus considerably more efficient than the conventional matrix method, and permits exact numerical solutions for much larger problems. Arguments are presented toward the view that the field equations are more fundamental in k-space. The details and some numerical results of the application of this method to the three-dimensional electromagnetic scattering problems are presented as an example.
Bojarski, N. N.↗
A cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet
The stress intensity factor is determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer. Since the stringer is modeled as a semiinfinite sheet, the solution is most appropriate for a crack tip located near a stringer edge. Both adherends are treated as homogeneous, orthotropic media. It is assumed they are in plane stress and the adhesive is in pure shear. From Green's functions and the complex variable theory of orthotropic elasticity developed by Lekhnitskiy, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations which is solved to obtain the shear stress distribution in the adhesive layer. With these stresses, equations for the stress intensity factors at both crack tips are found. A parametric study is conducted to determine the sensitivity of the system to material properties and specimen configuration. Unless the crack tip is very close to or under the stringer the stress intensity factor is approximately that of an unstiffened sheet. However, as the crack propagates beneath the stringer, the stress intensity factor decreases significantly. Increasing the stiffness of the stringer or the adhesive also results in a decrease in the stress intensity factor.
Bigelow, C. A.↗
Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet
The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.
Bigelow, C. A.↗
Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet
The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.
Bigelow, C. A.↗
Theoretical Prediction of Pressure Distributions on Nonlifting Airfoils at High Subsonic Speeds
Theoretical pressure distributions on nonlifting circular-arc airfoils in two-dimensional flows with high subsonic free-stream velocity are found by determining approximate solutions, through an iteration process, of an integral equation for transonic flow proposed by Oswatitsch. The integral equation stems directly from the small-disturbance theory for transonic flow. This method of analysis possesses the advantage of remaining in the physical, rather than the hodograph, variable and can be applied in airfoils having curved surfaces. After discussion of the derivation of the integral equation and qualitative aspects of the solution, results of calculations carried out for circular-arc airfoils in flows with free-stream Mach numbers up to unity are described. These results indicate most of the principal phenomena observed in experimental studies.
Spreiter, John R↗