SEARCH · Search NASA

Results for “SINGULAR 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗

NASA NTRS1979

Reduction of computer usage costs in predicting unsteady aerodynamic loadings caused by control surface motions: Computer program description

A digital computer program was developed to calculate unsteady loadings caused by motions of lifting surfaces with leading edge and trailing edge controls based on the subsonic kernel function approach. The pressure singularities at hinge line and side edges were extracted analytically as a preliminary step to solving the integral equation of collocation. The program calculates generalized aerodynamic forces for user supplied deflection modes. Optional intermediate output includes pressure at an array of points, and sectional generalized forces. From one to six controls on the half span can be accomodated.

Petrarca, J. R.↗

NASA NTRS1975

Prediction of unsteady aerodynamic loadings caused by leading edge and trailing edge control surface motions in subsonic compressible flow: Computer program description

A digital computer program has been developed to calculate unsteady loadings caused by motions of lifting surfaces with leading edge or trailing edge controls based on the subsonic kernel function approach. The pressure singularities at hinge line and side edges have been extracted analytically as a preliminary step to solving the integral equation by collocation. The program calculates generalized aerodynamic forces for user supplied deflection modes. Optional intermediate output includes pressure at an array of points, and sectional generalized forces. From one to six controls on the half span can be accommodated.

Redman, M. C.↗

NASA NTRS1991

RWF rotor-wake-fuselage code software reference guide

The RWF (Rotor-Wake-Fuselage) code was developed from first principles to compute the aerodynamics associated with the complex flow field of helicopter configurations. The code is sized for a single, multi-bladed main rotor and any configuration of non-lifting fuselage. The mathematical model for the RWF code is based on the integration of the momentum equations and Green's theorem. The unknowns in the problem are the strengths of prescribed singularity distributions on the boundaries of the flow. For the body (fuselage) a surface of constant strength source panels is used. For the rotor blades and rotor wake a surface of constant strength doublet panels is used. The mean camber line of the rotor airfoil is partitioned into surface panels. The no-flow boundary condition at the panel centroids is modified at each azimuthal step to account for rotor blade cyclic pitch variation. The geometry of the rotor wake is computers at each time step of the solution. The code produces rotor and fuselage surface pressures, as well as the complex geometry of the evolving rotor wake.

Berry, John D.↗

NASA NTRS1974

Computation of solar perturbations with Poisson series

Description of a project for computing first-order perturbations of natural or artificial satellites by integrating the equations of motion on a computer with automatic Poisson series expansions. A basic feature of the method of solution is that the classical variation-of-parameters formulation is used rather than rectangular coordinates. However, the variation-of-parameters formulation uses the three rectangular components of the disturbing force rather than the classical disturbing function, so that there is no problem in expanding the disturbing function in series. Another characteristic of the variation-of-parameters formulation employed is that six rather unusual variables are used in order to avoid singularities at the zero eccentricity and zero (or 90 deg) inclination. The integration process starts by assuming that all the orbit elements present on the right-hand sides of the equations of motion are constants. These right-hand sides are then simple Poisson series which can be obtained with the use of the Bessel expansions of the two-body problem in conjunction with certain interation methods. These Poisson series can then be integrated term by term, and a first-order solution is obtained.

Broucke, R.↗

NASA NTRS1975

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗

NASA NTRS1983

A computational method for subsonic compressible flow in diffusers

Integral methods are used to model subsonic flow in planar diffusers operating with or without stall. Turbulent boundary layers are computed with a lag-entrainment method employing new correlations for shape factor and skin friction, which are based on a Coles-Van Driest velocity profile. The inviscid flow is modeled as one dimensional and is easily coupled with the boundary layer method to yield equations that can be solved by integrating in the downstream direction. Analysis of the equations describing flow in a diffuser shows that the boundary layer separation singularity may not be eliminated by viscous-inviscid interaction; a modified interaction scheme is used to guarantee a non-singular equation set. The nature of this singularity offers an explanation for the apparent universality of asymmetric stall in planar diffusers. Results for several flows demonstrate the accuracy and range of applicability of the method.

Childs, R. E.↗

NASA NTRS1976

Multi-star processing and gyro filtering for the video inertial pointing system

The video inertial pointing (VIP) system is being developed to satisfy the acquisition and pointing requirements of astronomical telescopes. The VIP system uses a single video sensor to provide star position information that can be used to generate three-axis pointing error signals (multi-star processing) and for input to a cathode ray tube (CRT) display of the star field. The pointing error signals are used to update the telescope's gyro stabilization system (gyro filtering). The CRT display facilitates target acquisition and positioning of the telescope by a remote operator. Linearized small angle equations are used for the multistar processing and a consideration of error performance and singularities lead to star pair location restrictions and equation selection criteria. A discrete steady-state Kalman filter which uses the integration of the gyros is developed and analyzed. The filter includes unit time delays representing asynchronous operations of the VIP microprocessor and video sensor. A digital simulation of a typical gyro stabilized gimbal is developed and used to validate the approach to the gyro filtering.

Murphy, J. P.↗

NASA NTRS1987

Some numerical calculations by using linear classical sonic theories approached from sub- or supersonic speeds

The unsteady aerodynamics of a two-dimensional wing at sonic speed are studied by using so-called classical sonic theories (linear), approached from supersonic flow (M=1+0) or subsonic flow (M=1-0). In the former approach, the exact expressions of lift and lift distribution are obtained in terms of Fresnel integrals, while in the latter approach an integral equation must be solved, the kernel function of which is obtained from the subsonic Possio's equation and has a root singularity. The discrete analysis is adopted on the basis of the semicircle method (SCM) and the weighting function for subsonic-flow-Gauss-quadrature, as well as modified characteristics obtained from both approaches agree quite well with each other. The results obtained by the present computations are compared with those of DLM-C (subsonic 2D code) developed by ANDO et al, and are found to give a reasonable outer boundary for subsonic unsteady aerodynamics.

Yamamoto, Y.↗

NASA NTRS1981

General relaxation schemes in multigrid algorithms for higher order singularity methods

Relaxation schemes based on approximate and incomplete factorization technique (AF) are described. The AF schemes allow construction of a fast multigrid method for solving integral equations of the second and first kind. The smoothing factors for integral equations of the first kind, and comparison with similar results from the second kind of equations are a novel item. Application of the MD algorithm shows convergence to the level of truncation error of a second order accurate panel method.

Oskam, B.↗

NASA NTRS1955

On the Kernel function of the integral equation relating the lift and downwash distributions of oscillating finite wings in subsonic flow

This report treats the Kernel function of an integral equation that relates a known prescribed downwash distribution to an unknown lift distribution for a harmonically oscillating finite wing in compressible subsonic flow. The Kernel function is reduced to a form that can be accurately evaluated by separating the Kernel function into two parts: a part in which the singularities are isolated and analytically expressed and a nonsingular part which may be tabulated. The form of the Kernel function for the sonic case (Mach number 1) is treated separately. In addition, results for the special cases of Mach number of 0 (incompressible case) and frequency of 0 (steady case) are given. The derivation of the integral equation which involves this Kernel function is reproduced as an appendix. Another appendix gives the reduction of the form of the Kernel function obtained herein for the three-dimensional case to a known result of Possio for two-dimensional flow. A third appendix contains some remarks on the evaluation of the Kernel function, and a fourth appendix presents an alternate form of expression for the Kernel function.

Watkins, Charles E↗

NASA NTRS1972

Variational solutions by the use of stepwise constant functions. I - Linear case.

Demonstration that accurate solutions to the integral equation or to the extremal function can be easily obtained by the variational method in many cases when a stepwise constant function is used for the trial function. The evaluation procedure is simple and straightforward; integrals of the kernel function can be evaluated analytically; the method provides, when the solution involves singularities, the best mean value across the singularity; the results are accurate in both the detailed physical quantities and their averages; the resultant solution can be further integrated analytically over the parameters associated with the problem; and the method can be readily applied to nonlinear integral equations.

Yoshikawa, K. K.↗

NASA NTRS1977

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗

NASA NTRS2000

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗

NASA NTRS1984

The crack-inclusion interaction problem

The general plane elastostatic problem of interaction between a crack and an inclusion is considered. The Green's functions for a pair of dislocations and a pair of concentrated body forces are used to generate the crack and the inclusion. Integral equations are obtained for a line crack and an elastic line inclusion having an arbitrary relative orientation and size. The nature of stress singularity around the end points of rigid and elastic inclusions is described and three special cases of this intersection problem are studied. The problem is solved for an arbitrary uniform stress state away from the crack-inclusion region. The nonintersecting crack-inclusion problem is considered for various relative size, orientation, and stiffness parameters, and the stress intensity factors at the ends of the inclusion and the crack are calculated. For the crack-inclusion intersection case, special stress intensity factors are defined and are calculated for various values of the parameters defining the relative size and orientation of the crack and the inclusion and the stiffness of the inclusion.

Xue-Hui, L.↗

NASA NTRS1986

The crack-inclusion interaction problem

Liu, X.-H.↗

NASA NTRS1987

A new approach to the effect of sound on vortex dynamics

Analytical results are presented on the effect of acoustic radiation on three-dimensional vortex motions in a homogeneous, slightly compressible, inviscid fluid. The flow is considered as linear and irrotational everywhere except inside a very thin cylindrical core region around the vortex filament. In the outside region, a velocity potential is introduced that must be multivalued, and it is shown how to compute this scalar potential if the motion of the vortex filament is prescribed. To find the motion of this singularity in an external potential flow, a variational principle involving a volume integral that must exclude the singular region is considered. A functional of the external potential and vortex filament position is obtained whose extrema give equations to determine the sought-after evolution. Thus, a generalization of the Biot-Savart law to flows with constant sound speed at low Mach number is obtained.

Lund, Fernando↗