Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary integral equation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integral boundary conditions in phase field models

Modeling the chemical, electric and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential equations often with integral boundary conditions. Employing the familiar Poisson equation describing the electric potential evolution in a material exhibiting insulator to metal transitions, we exploit a special property of such an integral boundary condition, and we properly formulate the variational problem and establish its well-posedness. Next, we compare our method with the commonly-used Lagrange multiplier method that can also handle such boundary conditions. Numerical experiments demonstrate that our new method achieves optimal convergence rate in contrast to the conventional Lagrange multiplier method. Furthermore, the linear system derived from our method is symmetric positive definite, and can be efficiently solved by Conjugate Gradient method with algebraic multigrid preconditioning.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

On the solution of integral equations with strong ly singular kernels

In this paper some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m, m or = 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t,x) sup-m, terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗
NASA NTRS1986

On the solution of integral equations with strongly singular kernels

Some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m ,m greater than or equal 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t-x) sup -m , terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗
NASA NTRS1987

On the solution of integral equations with strongly singular kernels

Some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m, m greater than or equal 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t-x) sup-m, terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗
NASA NTRS1982

Construction of solutions for some nonlinear two-point boundary value problems

Constructive existence and uniqueness results for boundary value problems associated with some simple special cases of the second order equation y'' = f(x,y,y') 0 or = x or = 1, are sought. The approach considered is to convert the differential equation and boundary conditions to an integral equation via Green's functions, and then to apply fixed point and contraction map principles to a sequence of successive approximations. The approach is tested on several applied problems. Difficulties in trying to prove general theorems are discussed.

Pennline, J. A.↗
NASA NTRS1995

A Global Interpolation Function (GIF) boundary element code for viscous flows

Using global interpolation functions (GIF's), boundary element solutions are obtained for two- and three-dimensional viscous flows. The solution is obtained in the form of a boundary integral plus a series of global basis functions. The unknown coefficients of the GIF's are determined to ensure the satisfaction of the governing equations at selected collocation points. The values of the coefficients involved in the boundary integral equations are determined by enforcing the boundary conditions. Both primitive variable and vorticity-velocity formulations are examined.

Reddy, D. R.↗
NASA NTRS1993

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

Analyzing Transonic Flow Over A Wing And Fuselage

TAWFIVE computer program calculates transonic flow over transport-type wing and fuselage. Combines multigrid acceleration technique in iterative solution of potential equation with use of integral-form boundary-layer equations to provide software tool computationally efficient, yet accurate enough for use in design. Simplifies solution process by breaking problem into loosely coupled set of modified equations. Inviscid-flow method, using standard inviscid-flow equations, valid in "outer" region away from wing, whereas boundary-layer equations valid in thin region near solid surface of wing. Written in FORTRAN 77.

Melson, N. Duane↗
NASA NTRS1949

Two-Dimensional Potential Flows

Contents include the following: Characteristic differential equations - initial and boundary conditions. Integration of the second characteristic differential equations. Direct application of Meyer's characteristic hodograph table for construction of two-dimensional potential flows. Prandtl-Busemann method. Development of the pressure variation for small deflection angles. Numerical table: relation between deflection, pressure, velocity, mach number and mach angle for isentropic changes of state according to Prandtl-Meyer for air (k = 1.405). References.

Schaefer, Manfred↗
NASA NTRS2003

Electromagnetic Interactions GenERalized (EIGER): Applications at Johnson Space Center

This slide presentation reviews the software framework for the analysis and design of complex electromagnetic systems. The system called Electromagnetic Interactions GenERalized (EIGER) is a hybrid finite element solution for wave equations and deru=ives boundary element solutions for integral equation formulations. Written in Fortran 90 Eiger uses object-oriented design methods to abstract key analysis components. EIGER has been applied to higher order modeling for analysis of antennas to assist in the design of the antennas of a mini-AERCam. EIGER might also be used to model the Ion Cyclotron Resonant heating stage for the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR).

Dobbins, Justin A.↗
NASA NTRS1988

Three-dimensional unsteady transonic viscous-inviscid interaction using the Euler and boundary-layer equations

The objective of this study is the development of a numerical technique which can provide three-dimensional, time-accurate, compressible, turbulent flow solutions in a practical and relatively economical way. The approach taken is that of the method of viscous-inviscid interaction. The Euler equations are assumed to govern the outer inviscid portion of the flow, and the viscous layer close to the solid wall is described by a set of integral boundary-layer equations. The viscous solutions are obtained in a direct fashion with a weighted-average phase error scheme. The method of equivalent sources is used for viscous-inviscid coupling. Steady-state and unsteady computations for an AGARD airfoil and a wing show that satisfactory engineering solutions can be obtained for attached, high Reynolds number flows using this method. Quasi-unsteady interactions are shown to produce similar results to those provided by true-unsteady interactions. Considerable computer resources can be saved for some cases by using quasi-unsteady interactions.

Whitfield, David L.↗
NASA NTRS1986

Prediction of sound fields in cavities using boundary element methods

Two boundary element formulations of acoustical behavior exist; the Direct Boundary Element Methods (DBEM) based on the Helmholtz Integral Equations and the Indirect Boundary Element Methods (IBEM) based on Huygen's Principle. In this investigation, both methods are implemented utilizing a simple linear superparametric element. In addition the IBEM is studied using a quadratic isoparametric element. The accuracy and relative efficiency of the various techniques are examined. In order to properly model aircraft interior cavities the additional capability to model wall treatments and internal point sources is added to the methods. The procedures are verified for several well-understood cavity problems. The relative merits of each boundary element method and the finite element method are examined.

Bernhard, R. J.↗
NASA NTRS1988

Dendritic growth in the presence of convection

The motion of the freezing front between a dendritic crystal and a supercooled liquid is studied using an interface evolution equation derived from a boundary integral transformation of the transient convective-diffusion equation. A new steady-state theory is introduced that incorporates the effects of convection in dendritic growth. It is shown that in the absence of capillary effects the shape of the crystal-melt interface is a paraboloid of revolution, similar to that found in situations where diffusion is the sole heat transfer mechanism. A relation between the supercooling, the product of the tip velocity and tip radius, and the strength of the flow is derived which reduces to the well-known Ivantsov theory in the absence of convection. A non-linear interface-tracking algorithm is developed and used to study the temporal and spatial evolution of the dendritic interface. The important role of capillarity and convection on the interface dynamics is established and the response of the interface to finite amplitude disturbances is examined for the first time. Tip splitting is identified as the dominant destabilization mechanism in the limit of zero surface tension. Finite surface tension leads to interface stabilization, irrespective of the magnitude and structure of the external perturbations. Finally, convection significantly decreases the magnitude of the freezing velocity.

Beaghton, Pantelis John↗
NASA NTRS1977

Transient wave propagation in composite plates due to impact

Stress-wave propagation in a composite plate due to an impact has been examined using anisotropic elasticity theory. The plate is modelled as a number of identical anisotropic layers whose elastic moduli are given. Mindlin's (1959) approximate theory of plates is then applied to each layer to obtain a set of difference-differential equations of motion by using the interlaminar stresses and displacements as explicit variables. Dispersion relationships for harmonic waves are found when traction-free boundary conditions are applied to both surfaces of the plate. The difference-differential equations are reduced to difference equations via integral transforms. With given impact boundary conditions these equations are solved for an arbitrary number of layers in the plate, and the transient propagation of a stress wave is calculated by means of a fast Fourier transform algorithm.

Sung Kim, B.↗
NASA NTRS1951

Analysis of turbulent free-convection boundary layer on flat plate

With the use of Karman's integrated momentum equation for the boundary layer and data on the wall-shearing stress and heat transfer in forced-convection flow, a calculation was carried out for the flow and heat transfer in the turbulent free-convection boundary layer on a vertical flat plate. The calculation is for a fluid with a Prandtl number that is close to 1. A formula was derived for the heat-transfer coefficient that was in good agreement with experimental data in the range of Grashof numbers from 10sup10 to 10sup12. Because of the good agreement between the theoretical formula and the experimental data, the formula may be used to obtain data for high Grashof numbers. The calculation also yielded formulas for the maximum velocity in the boundary layer and for boundary-layer thickness.

Eckert, E R G↗