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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗

NASA NTRS1985

A first-order time-domain Green's function approach to supersonic unsteady flow

A time-domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. The Green's function method is employed in order to convert the potential-flow differential equation into an integral one. This integral equation is then discretized, in space through standard finite-element technique, and in time through finite-difference, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. Consistent with the spatial linear (first-order) finite-element approximations, the potential and co-normalwash are assumed to vary linearly in time. The long range goal of our research is to develop a comprehensive theory for unsteady supersonic potential aerodynamics which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗

NASA NTRS1986

Compliance matrices for cracked bodies

An algorithm is developed to construct the compliance matrix for a cracked solid in the integral-equation formulation of two-dimensional linear-elastic fracture mechanics. The integral equation is reduced to a system of algebraic equations for unknown values of the dislocation-density function at discrete points on the interval from -1 to 1, using the numerical procedure described by Gerasoulis (1982). Sample numerical results are presented, and it is suggested that the algorithm is especially useful in cases where iterative solutions are required; e.g., models of fiber-reinforced concrete, rocks, or ceramics where microcracking, fiber bridging, and other nonlinear effects are treated as nonlinear springs along the crack surfaces (Ballarini et al., 1984).

Ballarini, R.↗

NASA NTRS1993

Shear and compression buckling analysis for anisotropic panels with centrally located elliptical cutouts

An approximate analysis for buckling of biaxial- and shear-loaded anisotropic panels with centrally located elliptical cutouts is presented in the present paper. The analysis is composed of two parts, a prebuckling analysis and a buckling analysis. The prebuckling solution is determined using Lekhnitskii's complex variable equations of plane elastostatics combined with a Laurent series approximation and a boundary collocation method. The buckling solution is obtained using the principle of minimum potential energy. A by-product of the minimum potential energy equation is an integral equation which is solved using Gaussian quadrature. Comparisons with documented experimental results and finite element analyses indicate that the approximate analysis accurately predicts the buckling loads of square biaxial- and shear-loaded panels having elliptical cutouts with major axes up to sixty percent of the panel width. Results of a parametric study are presented for shear- and compression-loaded rectangular anisotropic panels with elliptical cutouts. The effects of panel aspect ratio, cutout shape, cutout size, cutout orientation, laminate anisotropy, and combined loading on the buckling load are examined.

Britt, V. O.↗

NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

NASA NTRS1982

Inverse scattering for an exterior Dirichlet problem

Scattering caused by a metallic cylinder in the field of a wire carrying a periodic current is studied, with a view to determining the location and shape of the cylinder in light of far field measurements between the cylinder and the wire. The associated direct problem is the exterior Dirichlet problem for the Helmholtz equation in two dimensions, and an improved low frequency estimate for its solution by integral equation methods is shown by inverse scattering calculations to be accurate to this estimate. The far field measurements are related to low frequency boundary integral equations whose solutions may be expressed in terms of a mapping function for the exterior of the unknown curve onto the exterior of a unit disk. The conformal transformation's Laurent expansion coefficients can be related to those of the far field, the first of which leads to the calculation of the distance between the source and the cylinder, while the other coefficients are determined by placing the source in a different location.

Hariharan, S. I.↗

NASA NTRS1979

Theoretical determination of nozzle admittances using a finite element method

A finite element method (FEM) is used to predict the admittances of axisymmetric nozzles. The flow in the nozzle is assumed to be isentropic and the disturbances are assumed to be small so that linear analyses apply. An approximate two dimensional compressible flow model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear acoustic wave equation. This partial differential wave equation is transformed to an integral equation using Galerkin's method and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. A two dimensional finite element method using linear triangular elements is used to solve the integral acoustic equation. A one dimensional FEM is used to solve the reduced nozzle acoustic equation developed by Crocco and the solution is used to verify the sufficiency of the boundary residual formation. It is shown that agreement between predicted values of the admittance and experimental data is quite good.

Sigman, R. K.↗

NASA NTRS2007

Using EIGER for Antenna Design and Analysis

EIGER (Electromagnetic Interactions GenERalized) is a frequency-domain electromagnetics software package that is built upon a flexible framework, designed using object-oriented techniques. The analysis methods used include moment method solutions of integral equations, finite element solutions of partial differential equations, and combinations thereof. The framework design permits new analysis techniques (boundary conditions, Green#s functions, etc.) to be added to the software suite with a sensible effort. The code has been designed to execute (in serial or parallel) on a wide variety of platforms from Intel-based PCs and Unix-based workstations. Recently, new potential integration scheme s that avoid singularity extraction techniques have been added for integral equation analysis. These new integration schemes are required for facilitating the use of higher-order elements and basis functions. Higher-order elements are better able to model geometrical curvature using fewer elements than when using linear elements. Higher-order basis functions are beneficial for simulating structures with rapidly varying fields or currents. Results presented here will demonstrate curren t and future capabilities of EIGER with respect to analysis of installed antenna system performance in support of NASA#s mission of exploration. Examples include antenna coupling within an enclosed environment and antenna analysis on electrically large manned space vehicles.

Champagne, Nathan J.↗

NASA NTRS1981

Inverse scattering for an exterior Dirichlet program

Scattering due to a metallic cylinder which is in the field of a wire carrying a periodic current is considered. The location and shape of the cylinder is obtained with a far field measurement in between the wire and the cylinder. The same analysis is applicable in acoustics in the situation that the cylinder is a soft wall body and the wire is a line source. The associated direct problem in this situation is an exterior Dirichlet problem for the Helmholtz equation in two dimensions. An improved low frequency estimate for the solution of this problem using integral equation methods is presented. The far field measurements are related to the solutions of boundary integral equations in the low frequency situation. These solutions are expressed in terms of mapping function which maps the exterior of the unknown curve onto the exterior of a unit disk. The coefficients of the Laurent expansion of the conformal transformations are related to the far field coefficients. The first far field coefficient leads to the calculation of the distance between the source and the cylinder.

Hariharan, S. I.↗

NASA NTRS1984

Random equations in aerodynamics

Literature was reviewed to identify aerodynamic models which might be treated by probablistic methods. The numerical solution of some integral equations that arise in aerodynamical problems were investigated. On the basis of the numerical studies a qualitative theory of random integral equations was developed to provide information on the behavior of the solutions of these equations (in particular, boundary and asymptotic behavior, and stability) and their statistical properties without actually obtaining explicit solutions of the equations.

Bharucha-Reid, A. T.↗

NASA NTRS1969

Remote probing of the optical strength of atmospheric turbulence and of wind velocity

A procedure for determining the optical strength of turbulence of the atmosphere and the wind velocity at various altitudes by measuring the spatial and temporal covariance of scintillation is developed. Emphasis is placed on the development of the formal relationships that have to be inverted to obtain the desired results. For determination of optical strength of turbulence, it is a linear integral equation that is developed. However, for determination of remote wind velocity, a nonlinear integral equation is obtained. A computer approach for solving each of the equations is suggested. The configuration and performance requirements of the measurement apparatus are discussed.

Fried, D. L.↗

NASA NTRS1979

Shape integral method for magnetospheric shapes

A method is developed for calculating the shape of any magnetopause to arbitrarily high precision. The method uses an integral equation which is evaluated for a trial shape. The resulting values of the integral equation as a function of auxiliary variables indicate how close one is to the desired solution. A variational method can then be used to improve the trial shape. Some potential applications are briefly mentioned.

Michel, F. C.↗

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 NTRS2019

Electromagnetic Scattering by Discrete Random Media. III: The Vector Radiative Transfer Equation

A vector radiative transfer equation with an additional source term typical of dense media is obtained. The analysis includes (i) the derivation of an integral equation for the correlation matrix of the exciting field coefficients accounting for the correlation between the particles, (ii) the derivation of an integral representation for the specific coherency dyadicin terms of this matrix, and (iii) the simplification of the integral equation for the correlation matrix and of the integral representation for the specific coherency dyadic by employing a series of approximations which are characteristic of sparse media.

Adrian Doicu↗

NASA NTRS1973

The problem of a finite strip compressed between two rough rigid stamps

A finite strip compressed between two rough rigid stamps is considered. The elastostatic problem is formulated in terms of a singular integral equation from which the proper stress singularities at the corners are determined. The singular integral equation is solved numerically to determine the stresses along the fixed ends of the strip. The effect of material properties and strip geometry on the stress intensity factor is presented graphically.

Gupta, G. D.↗

NASA NTRS1975

The problem of a finite strip compressed between two rough rigid stamps

Gupta, G. D.↗

NASA NTRS1984

Determination of equivalent model geometry for tunnel wall interference assessment/correction

A formula for the determination of equivalent model geometry with two variables measured at the interface is derived, based on two dimensional subsonic flow. This predicted model profile is a reasonable initial estimate for transonic flow as long as the sonic region does not reach the interface. A general formula is given in two forms. One is in terms of complex variable functions and the other is an integral equation. The complex-function formula has the advantage of using analytic expressions. The integral equation form requires a numerical solution after assuming the model geometry as a polynomial function. Examples are given to illustrate the application of the formulas.

Lo, C. F.↗

NASA NTRS1993

Scattering by a groove in an impedance plane

An analysis of two-dimensional scattering from a narrow groove in an impedance plane is presented. The groove is represented by a impedance surface and the problem reduces to that of scattering from an impedance strip in an otherwise uniform impedance plane. On the basis of this model, appropriate integral equations are constructed using a form of the impedance plane Green's functions involving rapidly convergent integrals. The integral equations are solved by introducing a single basis representation of the equivalent current on the narrow impedance insert. Both transverse electric (TE) and transverse magnetic (TM) polarizations are treated. The resulting solution is validated by comparison with results from the standard boundary integral method (BIM) and a high frequency solution. It is found that the presented solution for narrow impedance inserts can be used in conjunction with the high frequency solution for the characterization of impedance inserts of any given width.

Bindiganavale, Sunil↗