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 325 records · Page 18

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗

A domain decomposition preconditioning for an inverse volume scattering problem

In this work, we propose domain decomposition preconditioners for the solution of an integral equation formulation of the acoustic forward and inverse scattering problems. We study both forward and inverse volume problems and propose preconditioning techniques to accelerate the iterative solvers. For the forward scattering problem, we extend the domain decomposition based preconditioning techniques presented for partial differential equations in Cai and Sarkis (1999 SIAM J. Sci. Comput. 21 792–797), to integral equations. We combine this domain decomposition preconditioner with a low-rank correction, which is easy to construct, forming a new preconditioner. For the inverse scattering problem, we use the forward problem preconditioner as a building block for constructing a preconditioner for the Gauss–Newton Hessian. We present numerical results that demonstrate the performance of both preconditioning strategies.

97 MATHEMATICS AND COMPUTING↗

Partial-wave projection of relativistic three-body amplitudes

We derive the integral equations for partial-wave projected three-body scattering amplitudes, starting from the integral equations for three-body amplitudes developed for lattice QCD analyses. The results, which hold for generic three-body systems of spinless particles, build upon the recently derived partial-wave projected one-particle exchange, a primary component of the relativistic framework proven to satisfy 𝑆 matrix unitarity. We derive simplified expressions for factorizable short-distance interactions, 𝒦 3 , in two equivalent formalisms—one symmetric under particle interchange and one asymmetric. For the asymmetric case, we offer parametrizations useful for amplitude analysis. Finally, we examine toy models for 3⁢𝜋 systems at unphysically heavy pion masses with total isospins 0, 1, and 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗

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.↗

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.↗

Scalar field in Reissner–Nordström spacetime: Bound state and scattering state (with appendix on eliminating oscillation in partial sum approximation of periodic function)

Highlights: • Bound-state and scattering-state solutions of massive scalar fields in R-N spacetime. • Bound-state wave functions and eigenvalues of massive scalar fields in R-N spacetime. • Solving explicit expressions of scattering phase shifts by integral equation methods. • Introducing tortoise coordinates for R-N spacetime. • Eliminating oscillations in the partial sum approximation of periodic functions. In this paper, we solve the massive scalar field in the Reissner–Nordström spacetime. The scalar field in the Reissner–Nordström spacetime has both bound states and scattering states. For bound states, we solve the bound-state wave function and the eigenvalue spectrum. For scattering states, we solve the scattering wave function and give an explicit expression for scattering phase shift by the integral equation method. Especially, we introduce the tortoise coordinate for the Reissner–Nordström spacetime. Moreover, in the calculation of scattering cross sections, we encounter a difficulty in partial sum approximation. If the sum of partial waves cannot be performed exactly, one has to turn to the partial sum approximation which approximates a function by the first several terms of the series. However, in the partial sum approximation there exists an incorrect oscillation which cannot be eliminated by keeping more terms. In the appendix we suggest an approach for eliminating such oscillations in the partial sum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

The effect of adhesive layer on crack propagation in laminates

The effect of the adhesive layer on crack propagation in composite materials is investigated. The composite medium consists of parallel load carrying laminates and buffer strips arranged periodically and bonded with thin adhesive layers. The strips, assumed to be isotropic and linearly elastic, contain symmetric cracks of arbitrary lengths located normal to the interfaces. Two problems are considered: (1) thin adhesive layers are approximated by uncoupled tension and shear springs distributed along the interfaces of the strips for which only the case of internal cracks can be treated rigorously; (2) broken laminates and the true singular behavior in the presence of the adhesive layer are studied. The adhesive is then treated as an isotropic, linearly elastic continuum. General expressions for field quantities are obtained in terms of infinite Fourier integrals. These expressions give a system of singular integral equations in terms of the crack surface displacement derivatives. By using appropriate quadrature formulas, the integral equations reduce to a system of linear algebraic equations which are solved numerically.

Gecit, M. R.↗

PANEL Code for Planar Cascades

Plane cascade flow solved using improved surface-singularity methods. Current external aerodynamic integral-equation techniques adapted for use in internal flow calculation. Inherent computational speed and flexibility of integral equation solutions make them useful for design calculations.

Mcfarland, E. R.↗

Microstress Analysis of Periodic Composites

Local elastic fields in the unit cell of a periodic composite are examined numerically with an integral equation approach. Techniques of Fourier series and Green's functions are used to construct the integral equations. Numerical solutions are obtained using the Fourier series approach with rectangular subvolume elements. Specific results are given for a tungsten/copper metal matrix composite.

Walker, Kevin P.↗

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 1: Theory document (version 3.0)

An outline of the derivation of the differential equation governing linear subsonic and supersonic potential flow is given. The use of Green's Theorem to obtain an integral equation over the boundary surface is discussed. The engineering techniques incorporated in the Panel Aerodynamics (PAN AIR) program (a discretization method which solves the integral equation for arbitrary first order boundary conditions) are then discussed in detail. Items discussed include the construction of the compressibility transformation, splining techniques, imposition of the boundary conditions, influence coefficient computation (including the concept of the finite part of an integral), computation of pressure coefficients, and computation of forces and moments. Principal revisions to version 3.0 are the following: (1) appendices H and K more fully describe the Aerodynamic Influence Coefficient (AIC) construction; (2) appendix L now provides a complete description of the AIC solution process; (3) appendix P is new and discusses the theory for the new FDP module (which calculates streamlines and offbody points); and (4) numerous small corrections and revisions reflecting the MAG module rewrite.

Epton, Michael A.↗