Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Topics in structural dynamics: Nonlinear unsteady transonic flows and Monte Carlo methods in acoustics

The results are reported of two unrelated studies. The first was an investigation of the formulation of the equations for non-uniform unsteady flows, by perturbation of an irrotational flow to obtain the linear Green's equation. The resulting integral equation was found to contain a kernel which could be expressed as the solution of the adjoint flow equation, a linear equation for small perturbations, but with non-constant coefficients determined by the steady flow conditions. It is believed that the non-uniform flow effects may prove important in transonic flutter, and that in such cases, the use of doublet type solutions of the wave equation would then prove to be erroneous. The second task covered an initial investigation into the use of the Monte Carlo method for solution of acoustical field problems. Computed results are given for a rectangular room problem, and for a problem involving a circular duct with a source located at the closed end.

Haviland, J. K.↗
NASA NTRS1975

General relativity and satellite orbits

The general relativistic correction to the position of a satellite is found by retaining Newtonian physics for an observer on the satellite and introducing a potential. The potential is expanded in terms of the Keplerian elements of the orbit and substituted in Lagrange's equations. Integration of the equations shows that a typical earth satellite with small orbital eccentricity is displaced by about 17 cm. from its unperturbed position after a single orbit, while the periodic displacement over the orbit reaches a maximum of about 3 cm. The moon is displaced by about the same amounts. Application of the equations to Mercury gives a total displacement of about 58 km. after one orbit and a maximum periodic displacement of about 12 km.

Rubincam, D. P.↗
NASA NTRS1974

A note on a broken layer in an orthotropic laminate composite

An orthotropic laminate composite containing a completely broken layer is considered. The problem is formulated in terms of integral transforms and then reduced to a singular integral equation which is solved numerically. The strength of stress singularity at the crack tip is determined from a characteristic equation which is obtained by studying the dominant part of the singular integral equation near the end points. The stress intensity factors are given for various material properties.

Arin, K.↗
NASA NTRS1977

A note on a broken layer in an orthotropic laminate composite

An orthotropic laminate composite containing a completely broken layer is considered. The problem is formulated in terms of integral transforms and then reduced to a singular integral equation which is solved numerically. The strength of stress singularity at the crack tip is determined from a characteristic equation which is obtained by studying the dominant part of the singular integral equation near the end points. The stress intensity factors are given for various material properties.

Arin, K.↗
NASA NTRS1981

New approach to the solution of large, full matrix equations

A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.

Clark, R. W.↗
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

Ignition of a combustible half space

A half space of combustible material is subjected to an arbitrary energy flux at the boundary where convection heat loss is also allowed. An asymptotic analysis of the temperature growth reveals two conditions necessary for ignition to occur. Cases of both large and order unity Lewis number are shown to lead to a nonlinear integral equation governing the thermal runaway. Some global and asymptotic properties of the integral equation are obtained.

Olmstead, W. E.↗
NASA NTRS1983

Estimates of the seasonal mean vertical velocity fields of the extratropical Northern Hemisphere

Indirect methods are employed to estimate the wintertime and summertime mean vertical velocity fields of the extratropical Northern Hemisphere and intercomparisons are made, together with comparisons with mean seasonal patterns of cloudiness and precipitation. Twice-daily NMC operational analyses produced general circulation statistics for 11 winters and 12 summers, permitting calculation of the seasonal NMC averages for 6 hr forecasts, solution of the omega equation, integration of continuity equation downward from 100 mb, and solution of the thermodynamic energy equation in the absence of diabatic heating. The methods all yielded similar vertical velocity patterns; however, the magnitude of the vertical velocities could not be calculated with great accuracy. Orography was concluded to have less of an effect in summer than in winter, when winds are stronger.

White, G. H.↗
NASA NTRS1983

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

Numerical solution of three-dimensional unsteady transonic flow over wings including inviscid/viscous interactions

A numerical procedure is presented for computing the unsteady transonic flow field about three dimensional swept wings undergoing general time dependent motion. The outer inviscid portion of the flow is assumed to be governed by the modified unsteady transonic small disturbance potential equation which is integrated in the time domain by means of an efficient alternating direction implicit approximate factorization algorithm. Gross dominant effects of the shock boundary layer interaction are accounted for by a simple empirically defined model. Viscous flow regions adjacent to the wing surface and in the trailing wake are described by a set of integral equations appropriate for compressible turbulent shear layers. The two dimensional boundary layer equations are applied quasi-statically stripwise across the span. Coupling with the outer inviscid flow is implemented through use of the displacement thickness concept within the limitations of small disturbance theory. Validity of the assumptions underlying the method is established by comparison with experimental data for the flow about a high aspect ratio transport wing having an advanced airfoil section.

Rizzetta, D. P.↗
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.↗