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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An entropy and viscosity corrected potential method for rotor performance prediction

An unsteady Full-Potential Rotor code (FPR) has been enhanced with modifications directed at improving its drag prediction capability. The shock generated entropy has been included to provide solutions comparable to the Euler equations. A weakly interacted integral boundary layer has also been coupled to FPR in order to estimate skin-friction drag. Pressure distributions, shock positions, and drag comparisons are made with various data sets derived from two-dimensional airfoil, hovering, and advancing high speed rotor tests. In all these comparisons, the effect of the nonisentropic modification improves (i.e., weakens) the shock strength and wave drag. In addition, the boundary layer method yields reasonable estimates of skin-friction drag. Airfoil drag and hover torque data comparisons are excellent, as are predicted shock strength and positions for a high speed advancing rotor.

Bridgeman, John O.↗
NASA NTRS1989

A general dynamic model of flexible robot arms for control

Hamilton's principle is used to derive the dynamic model for a large class of flexible robot arms. The resultant dynamic model consists of a system of partial differential-integral equations and the dynamic boundary conditions associated with it. Some properties of the model are observed, and its application to control is discussed. This model represents an infinite-dimensional nonlinear dynamic system and yet can be turned into a finite-dimensional system that could be obtained by modal expansion, if it is desired. This provides more flexibility for control purposes as well as for the analysis of the system.

Ding, X.↗
NASA NTRS1993

Analysis of microstrip patch antennas with nonzero surface resistance

The scattering properties of a microstrip patch antenna with nonzero surface impedance are examined. The electric field integral equation for a current element on a grounded dielectric slab is developed for a rectangular geometry by using Galerkin's technique with subdomain piecewise linear basis functions. The integral equation includes a resistive boundary condition on the surface of the patch. The incident field on the patch is expressed as a function of incidence angle. The resulting system of equations is then solved for the unknown current modes on the patch, and the radar cross section is calculated for a given scattering angle. Theoretical results in the form of radar cross section as a function of frequency are compared with results measured at the NASA Langley Research Center.

Shively, David G.↗
NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu↗
NASA NTRS1995

Viscous aerodynamic design using the adjoint variable approach

The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.

DeRise, George↗
NASA NTRS1997

Variational Methods in Design Optimization and Sensitivity Analysis for Two-Dimensional Euler Equations

Variational methods (VM) sensitivity analysis employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗
NASA NTRS1973

Two-dimensional laminar incompressible separated flow past airfoils

A method is proposed to treat the problem of steady, two-dimensional, laminar, incompressible high Reynolds number separated flow past thin airfoils. An integral form of the boundary layer equations with interaction is used and the interaction between the inviscid and viscous flow fields is provided for by use of a thin-airfoil integral. Documentation of the attempts at obtaining a solution is presented. A survey of the current state-of-the-art of problems involving viscous-inviscid interactions in flow fields with separation is given.

Plotkin, A.↗
NASA NTRS1987

A finite-difference outer layer and integral inner layer method for the solution of the turbulent boundary layer equations

A new turbulent boundary-layer method is developed which models the inner region with the law of the wall while the outer region uses Clauser's eddy viscosity in Matsuno's finite-difference method. The match point between the inner and outer regions as well as the wall shear stress are determined at each marching step during the computation. Results obtained for incompressible, two-dimensional flow over flat plates and ellipses are compared with solutions from a baseline method which uses a finite-difference method for the entire boundary layer. Since the present method used the finite-difference method in the outer region only, the number of grid points required was about half that needed for the baseline method. Accurate displacement and momentum thicknesses were predicted for all cases. Skin friction was predicted well for the flat plate, but the accuracy decreased significantly for the ellipses. Adding a wake functions to the law of the wall allows some of the pressure gradient effect to be taken into account thereby increasing the accuracy of the method.

Barnwell, R. W.↗
NASA NTRS1990

Accurate calculation and instability of supersonic wake flows

This study is concerned with the computation and linear stability of a class of laminar compressible wake flows. The emphasis is on correct basic flow profiles that satisfy the steady equations of motion, and to this end the unperturbed state is obtained through numerical integration of the compressible boundary-layer equations. The linear stability of the flow is examined via the Rayleigh equation that describes evolution of inviscid disturbances. Analytical results are given for short- and long-wavelength disturbances and some numerical results of the general eigenvalue problem are also reported.

Papageorgiou, Demetrius T.↗
NASA NTRS2002

Higher Order Bases in a 2D Hybrid BEM/FEM Formulation

The advantages of using higher order, interpolatory basis functions are examined in the analysis of transverse electric (TE) plane wave scattering by homogeneous, dielectric cylinders. A boundary-element/finite-element (BEM/FEM) hybrid formulation is employed in which the interior dielectric region is modeled with the vector Helmholtz equation, and a radiation boundary condition is supplied by an Electric Field Integral Equation (EFIE). An efficient method of handling the singular self-term arising in the EFIE is presented. The iterative solution of the partially dense system of equations is obtained using the Quasi-Minimal Residual (QMR) algorithm with an Incomplete LU Threshold (ILUT) preconditioner. Numerical results are shown for the case of an incident wave impinging upon a square dielectric cylinder. The convergence of the solution is shown versus the number of unknowns as a function of the completeness order of the basis functions.

Fink, Patrick W.↗
NASA NTRS1973

Unsteady subsonic compressible flow around finite thickness wings.

A general formulation for the unsteady subsonic compressible potential flow around aircraft having arbitrary configurations is presented. An integral representation of the velocity potential is obtained. From this a linear integral equation relating the perturbation potential and its normal derivative (which is known from the boundary conditions) is derived. For the numerical solution of the integral equation, the surface of the aircraft is divided into small elements and the potential is assumed to be constant within each element. Numerical results are obtained for an oscillating finite-thickness wing and indicate good convergence and excellent agreement with existing lifting surface solutions.

Morino, L.↗
NASA NTRS1974

Parametric study of relaminarization of turbulent boundary layers on nozzle walls

By means of comparisons between theoretical predictions and experimental data, the accuracy of a boundary procedure to predict the effect of large streamwise accelerations upon initially turbulent boundary layers is assessed. The boundary layer procedure is based upon simultaneous solution of the boundary layer partial differential equations and the integral turbulence kinetic energy equation. The results of the present investigation show the ability of the procedure to accurately predict properties of boundary layers subjected to large streamwise accelerations. The procedure is used to conduct a parametric study of the effect of free stream turbulence, heat transfer, Reynolds number, acceleration, and Mach number on boundary layers in supersonic nozzles to assist in the design of a quiet tunnel. Results of the investigation show that, even in the presence of moderate free-stream turbulence levels, the boundary layer in the approach section of the quiet tunnel nozzle relaminarizes and becomes thin enough to be removed by a small slot in the nozzle wall. Furthermore, the calculations indicate that it should be possible to maintain a laminar boundary layer for the entire length of the supersonic portion of the quiet tunnel nozzle.

Kreskovsky, J. P.↗
NASA NTRS1978

Mixed boundary-value problems in mechanics

Definitions in the case of multiple series equations and multiple integral equations are examined. In considering the solution of a given mixed boundary value problem perhaps the simplest technique is the direct application of the method of complex potentials provided the problem admits such potentials and the domain and the boundary conditions are suitable for such an application. The direct application of complex potentials is described with the aid of examples, taking into account a problem in potential theory, the case of periodic cuts, and an elasticity problem for a nonhomogeneous plane. The reduction to singular integral equations is discussed along with the numerical solution of singular integral equations of the first kind, integral equations with generalized Cauchy kernels, and singular integral equations of the second kind.

Erdogan, F.↗
NASA NTRS1992

A finite element-boundary integral method for conformal antenna arrays on a circular cylinder

Conformal antenna arrays offer many cost and weight advantages over conventional antenna systems. In the past, antenna designers have had to resort to expensive measurements in order to develop a conformal array design. This is due to the lack of rigorous mathematical models for conformal antenna arrays, and as a result the design of conformal arrays is primarily based on planar antenna design concepts. Recently, we have found the finite element-boundary integral method to be very successful in modeling large planar arrays of arbitrary composition in a metallic plane. Herewith we shall extend this formulation for conformal arrays on large metallic cylinders. In this we develop the mathematical formulation. In particular we discuss the finite element equations, the shape elements, and the boundary integral evaluation, and it is shown how this formulation can be applied with minimal computation and memory requirements. The implementation shall be discussed in a later report.

Kempel, Leo C.↗
NASA NTRS1992

A finite element-boundary integral method for conformal antenna arrays on a circular cylinder

Conformal antenna arrays offer many cost and weight advantages over conventional antenna systems. In the past, antenna designers have had to resort to expensive measurements in order to develop a conformal array design. This was due to the lack of rigorous mathematical models for conformal antenna arrays. As a result, the design of conformal arrays was primarily based on planar antenna design concepts. Recently, we have found the finite element-boundary integral method to be very successful in modeling large planar arrays of arbitrary composition in a metallic plane. We are extending this formulation to conformal arrays on large metallic cylinders. In doing so, we will develop a mathematical formulation. In particular, we discuss the finite element equations, the shape elements, and the boundary integral evaluation. It is shown how this formulation can be applied with minimal computation and memory requirements.

Kempel, Leo C.↗
NASA NTRS1993

An efficient formulation of the coupled finite element-integral equation technique for solving large 3D scattering problems

It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.

Cwik, T.↗