Engineering Papers⌕ Search

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.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Subprograms for integrating the equations of motion of satellites. FORTRAN 4

The subprograms for the formation of the right members of the equations of motion of artificial Earth satellites (AES), integration of systems of differential equations by Adams' method, and the calculation of the values of various functions from the AES parameters of motion are described. These subprograms are written in the FORTRAN 4 language and constitute an essential part of the package of applied programs for the calculation of navigational parameters AES.

Prokhorenko, V. I.↗
NASA NTRS1983

Spinning mode acoustic radiation from the flight inlet

A mathematical model was developed for spinning mode acoustic radiation from a thick wall duct without flow. This model is based on a series of experiments (with and without flow). A nearly pure azimuthal spinning mode was isolated and then reflection coefficients and far field pressure (amplitude and phase) were measured. In our model the governing boundary value problem for the Helmholtz equation is first converted into an integral equation for the unknown acoustic pressure over a disk, S1, near the mouth of the duct and over the exterior surface, S2, of the duct. Assuming a pure azimuthal mode excitation, the azimuthal dependence is integrated out which yields an integral equation over the generator C1 of S1 and the generator C2 of S2. The sound pressure on C1 was approximated by a truncated modal expansion of the interior acoustic pressure. Piecewise linear spline approximation on C2 was used.

Moss, W. F.↗
DOE OSTI2024

Level 2 Milestone: Develop and Incorporate Novel Code Verification of Fundamental Equations in Gemma and Set Up Appropriate Tests

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementations of the electric-, magnetic-, and combinedfield integral equations, are frequently used to solve Maxwell’s equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this report, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches in Gemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Antenna pattern control using impedance surfaces

A hybrid numerical technique is developed for electrically large pyramidal horn antennas radiating in free space. A stepped-waveguide method is used to analyze the interior surfaces of the horn transition. The Electric Field Integral Equation (EFIE) is employed on the outer surfaces of the pyramidal horn including the radiating aperture. Meanwhile, the Magnetic Field Integral Equation (MFIE) is used on the aperture to relate the aperture fields and those in the horn transition. The resultant hybrid field integral equation (HFIE) is solved numerically by the method of moments. This formulation is both accurate and numerically stable so that high-gain microwave pyramidal horns can be analyzed rigorously. Far-field radiation patterns, both computed and measured, are presented for three electrically-large x-band horn antennas. The comparisons demonstrate that this method is accurate enough to predict the fine pattern structure at wide angles and in the back region. Computed far-field patterns and aperture field distribution of two smaller x-band horns are also presented along with a discussion on the validity of the approximate aperture field distributions routinely used in the analysis and design of pyramidal horns.

Balanis, Constantine A.↗
NASA NTRS1993

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗
NASA NTRS1974

Application of unsteady lifting surface theory to propellers in forward flight

Development of a theory for determining the aerodynamic forces for unsteady, compressible subsonic flow on a propeller in forward flight. The acceleration potential method is used in developing the basic downwash integral equation which governs the flow. This integral equation is solved by the doublet-lattice method, which consists of placing 'load' lines at certain locations on the chord and satisfying the downwash condition at other selected positions. The examples presented include the spanwise and chordwise loading on a rotating propeller for incompressible flow, an example of compressible flow calculations, and, finally, a calculation illustrating the loss of aerodynamic damping of a propeller blade due to the passage of the blade over its own wake.

Hammond, C. E.↗
NASA NTRS1973

Two bonded half planes with a crack going through the interface

The plane problem of two bonded elastic half planes containing a finite crack perpendicular to and going through the interface is considered. The problem is formulated as a system of singular integral equations with generalized Cauchy kernels. Even though the system has three irregular points, it is shown that the unknown functions are algebraically related at the irregular point on the interface and the integral equations can be solved by a method developed previously. The system of integral equations is shown to yield the same characteristic equation as that for two bonded quarter planes in the general case of the through crack, and the characteristic equation for a crack tip terminating at the interface in the special case. The numerical results given in the paper include the stress intensity factors at the crack tips, the normal and shear components of the stress intensity factors at the singular point on the interface, and the crack surface displacements.

Erdogan, F.↗
NASA NTRS1986

The crack problem in a reinforced cylindrical shell

In this paper a partially reinforced cylinder containing an axial through crack is considered. The reinforcement is assumed to be fully bonded to the main cylinder. The composite cylinder is thus modelled by a nonhomogeneous shell having a step change in the elastic properties at the z=0 plane, z being the axial coordinate. Using a Reissner type transverse shear theory the problem is reduced to a pair of singular integral equations. In the special case of a crack tip touching the bimaterial interface it is shown that the dominant parts of the kernels of the integral equations associated with both membrane loading and bending of the shell reduce to the generalized Cauchy kernel obtained for the corresponding plane stress case. The integral equations are solved and the stress intensity factors are given for various crack and shell dimensions. A bonded fiberglass reinforcement which may serve as a crack arrestor is used as an example.

Yahsi, O. S.↗
NASA NTRS1986

The crack problem in a reinforced cylindrical shell

A partially reinforced cylinder containing an axial through crack is considered. The reinforcement is assumed to be fully bonded to the main cylinder. The composite cylinder is thus modelled by a nonhomogeneous shell having a step change in the elastic properties at the z = 0 plane, z being the axial coordinate. Using a Reissner type transverse shear theory the problem is reduced to a pair of singular integral equations. In the special case of a crack tip touching the bimaterial interface it is shown that the dominant parts of the kernels of the integral equations associated with both membrane loading and bending of the shell reduce to the generalized Cauchy kernel obtained for the corresponding plane stress case. The integral equations are solved and the stress intensity factors are given for various crack and shell dimensions. A bonded fiberglass reinforcement which may serve as a crack arrestor is used as an example.

Yahsi, O. S.↗
NASA NTRS1976

A nonlinear finite-element analysis of wings in steady incompressible flows with wake roll-up

The problem of lifting surfaces and complex aircraft configurations in steady incompressible flow is considered. For lifting surfaces the problem is formulated in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential on the lifting surface. For complex configurations the problem is formulated in terms of an integral equation relating the potential to its normal derivative on the surface of the aircraft. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential (or the potential discontinuity) and its normal derivative to be constant within each element. The wake geometry is obtained by iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines.

Suciu, E. O.↗
NASA NTRS1992

Wavelets and electromagnetics

Wavelets are an exciting new topic in applied mathematics and signal processing. This paper will provide a brief review of wavelets which are also known as families of functions with an emphasis on interpretation rather than rigor. We will derive an indirect use of wavelets for the solution of integral equations based techniques adapted from image processing. Examples for resistive strips will be given illustrating the effect of these techniques as well as their promise in reducing dramatically the requirement in order to solve an integral equation for large bodies. We also will present a direct implementation of wavelets to solve an integral equation. Both methods suggest future research topics and may hold promise for a variety of uses in computational electromagnetics.

Kempel, Leo C.↗
NASA NTRS1973

The problem of edge cracks in an infinite strip

The elastostatic plane problem of an infinite strip containing two symmetrically located internal cracks perpendicular to the boundary is formulated in terms of a singular integraL equation with the derivative of the crack surface displacement as the density function. The solution of the problem is obtained for various crack geometries and for uniaxial tension applied to the strip away from the crack region. The limiting case of the edge cracks is then considered in some detail. The fundamental function of the integral equation is obtained and a numerical technique for solving the singular integral equations with this particular type of fundamental function which is characteristic of the edge cracks is described. The stress intensity factor for the complete range of net ligament-to-width ratio is calculated. The results also include the solution of the edge crack problem in an elastic half plane.

Gupta, G. D.↗
NASA NTRS1974

The problem of edge cracks in an infinite strip

The elastostatic plane problem of an infinite strip containing two symmetrically located internal cracks perpendicular to the boundary is formulated in terms of a singular integral equation with the derivative of the crack surface displacement as the density function. The solution of the problem is obtained for various crack geometries and for uniaxial tension applied to the strip away from the crack region. The limiting case of the edge cracks is then considered in some detail. The fundamental function of the integral equation is obtained and a numerical technique for solving the singular integral equations with this particular type of fundamental function which is characteristic of the edge cracks is described. The stress-intensity factor for the complete range of net ligament-to-width ratio is calculated. The results also include the solution of the edge crack problem in an elastic half plane.

Gupta, G. D.↗
NASA NTRS1987

Full potential integral solution for transonic flows with and without embedded Euler domains

Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.

Kandil, Osama A.↗
NASA NTRS1995

Electromagnetic Scattering From a Rectangular Cavity Recessed in a 3-D Conducting Surface

The problem of electromagnetic (EM) scattering from an aperture backed by a rectangular cavity recessed in a three-dimensional conducting body is analyzed using the coupled field integral equation approach. Using the free-space Green's function, EM fields scattered outside the cavity are determined in terms of (1) an equivalent electric surface current density flowing on the three-dimensional conducting surface of the object including the cavity aperture and (2) an equivalent magnetic surface current density flowing over the aperture only. The EM fields inside the cavity are determined using the waveguide modal expansion functions. Making the total tangential electric and magnetic fields across the aperture continuous and subjecting the total tangential electric field on the outer conducting three-dimensional surface of the object to zero, a set of coupled integral equations is obtained. The equivalent electric and magnetic surface currents are then obtained by solving the coupled integral equation using the Method of Moments. The numerical results on scattering from rectangular cavities embedded in various three-dimensional objects are compared with the results obtained by other numerical techniques.

Deshpande, M. D.↗
NASA NTRS1997

Electromagnetic Scattering from Arbitrarily Shaped Aperture Backed by Rectangular Cavity Recessed in Infinite Ground Plane

The electromagnetic scattering from an arbitrarily shaped aperture backed by a rectangular cavity recessed in an infinite ground plane is analyzed by the integral equation approach. In this approach, the problem is split into two parts: exterior and interior. The electromagnetic fields in the exterior part are obtained from an equivalent magnetic surface current density assumed to be flowing over the aperture and backed by an infinite ground plane. The electromagnetic fields in the interior part are obtained in terms of rectangular cavity modal expansion functions. The modal amplitudes of cavity modes are determined by enforcing the continuity of the electric field across the aperture. The integral equation with the aperture magnetic current density as an unknown is obtained by enforcing the continuity of magnetic fields across the aperture. The integral equation is then solved for the magnetic current density by the method of moments. The electromagnetic scattering properties of an aperture backed by a rectangular cavity are determined from the magnetic current density. Numerical results on the backscatter radar cross-section (RCS) patterns of rectangular apertures backed by rectangular cavities are compared with earlier published results. Also numerical results on the backscatter RCS patterns of a circular aperture backed by a rectangular cavity are presented.

Cockrell, C. R.↗
DOE OSTI2024

Difference equations and integral families for Witten diagrams

We show that tree-level and one-loop Mellin space correlators in anti-de Sitter space obey certain difference equations, which are the direct analog to the differential equations for Feynman loop integrals in the flat space. Finite-difference relations, which we refer to as “summation-by-parts relations”, in parallel with the integration-by-parts relations for Feynman loop integrals, are derived to reduce the integrals to a basis. We illustrate the general methodology by explicitly deriving the difference equations and summation-by-parts relations for various tree-level and one-loop Witten diagrams up to the four-point bubble level.

AdS-CFT Correspondence↗