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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

End wall flow characteristics and overall performance of an axial flow compressor stage

This review indicates the possible future directions for research on endwall flows in axial flow compressors. Theoretical investigations on the rotor blade endwall flows in axial flow compressors reported here include the secondary flow calculation and the development of the momentum integral equations for the prediction of the annulus wall boundary layer. The equations for secondary vorticity at the rotor exit are solved analytically. The solution includes the effects of rotation and the viscosity. The momentum integral equations derived include the effect of the blade boundary layers. The axial flow compressor facility of the Department of Aerospace Engineering at The Pennsylvania State University, which is used for the experimental investigations of the endwall flows, is described in some detail. The overall performance and other preliminary experimental results are presented. Extensive radial flow surveys are carried out at the design and various off design conditions. These are presented and interpreted in this report. The following experimental investigations of the blade endwall flows are carried out. (1) Rotor blade endwall flows: The following measurements are carried out at four flow coefficients. (a) The rotor blade static pressures at various axial and radial stations (with special emphasis near the blade tips). (b) The hub wall static pressures inside the rotor blade passage at various axial and tangential stations. (2) IGV endwall flows: The following measurements are carried out at the design flow coefficient. (a) The boundary layer profiles at various axial and tangential stations inside the blade passage and at the blade exit. (b) Casing static pressures and limiting streamline angles inside the blade passage.

Sitaram, N.↗
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

Reaction formulation for radiation and scattering from plates, corner reflectors and dielectric-coated cylinders

The reaction concept is employed to formulate an integral equation for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders. The surface-current density on the conducting surface is expanded with subsectional bases. The dielectric layer is modeled with polarization currents radiating in free space. Maxwell's equation and the boundary conditions are employed to express the polarization-current distribution in terms of the surface-current density on the conducting surface. By enforcing reaction tests with an array of electric test sources, the moment method is employed to reduce the integral equation to a matrix equation. Inversion of the matrix equation yields the current distribution, and the scattered field is then obtained by integrating the current distribution. The theory, computer program and numerical results are presented for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders.

Wang, N. N.↗
NASA NTRS2020

An Overview of the Null-Field Method. I. Formulation and Basic Results

In this paper we revisit the fundamentals of the null-field method with discrete sources. We prove the unique solvability of the null-field equations for the total field inside the particle and the internal field outside the particle, at all wavenumbers. For this purpose, we use the equivalence between the null-field and surface integral equations methods. Furthermore, we discuss the completeness property of different systems of discrete sources for a surface field approximation, and derive an infinite set of integral equations for the surface fields in a variety of discrete sources. Finally, we formulate the null-field scheme as an approach aiming to construct an approximate solution to the scattering problem. The way in which we introduce the matrix of a particle is different from the standard approach relying on the assumption that the system of regular vector spherical wave functions for the interior problem is a basis.

Null-field method↗
NASA NTRS1978

The thermoelastic patch contact problem for large Peclet number

An equation in the form of a Fourier series was derived in an earlier study to give displacements of the originally straight edge of a plate as the result of a small heat input patch moving at speed c along that edge. That unwieldly formulation is reduced to a simple integral equation, which is shown to be manageable in an example calculation. The integral equation is shown to be closely related to that of Ling and Mow, although it is much simpler to use. Comparison of the results of the simplification to results of the original series formulation suggest that the equation proposed for 'high Peclet number' may be a valid approximation for Peclet number as low as 2.5.

Kilaparti, S. R.↗
NASA NTRS1986

Aerodynamics via acoustics - Application of acoustic formulas for aerodynamic calculations

Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.

Farassat, F.↗
NASA NTRS1986

Aerodynamics Via Acoustics: Application of Acoustic Formulas for Aerodynamic Calculations

Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.

Farassat, F.↗
NASA NTRS2012

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method↗
NASA NTRS1986

Analysis of mixed-mode crack propagation using the boundary integral method

Crack propagation in a rotating inner raceway of a high speed roller bearing is analyzed using the boundary integral equation method. The method consists of an edge crack in a plate under tension, upon which varying Hertzian stress fields are superimposed. A computer program for the boundary integral equation method was written using quadratic elements to determine the stress and displacement fields for discrete roller positions. Mode I and Mode II stress intensity factors and crack extension forces G sub 00 (energy release rate due to tensile opening mode) and G sub r0 (energy release rate due to shear displacement mode) were computed. These calculations permit determination of that crack growth angle for which the change in the crack extension forces is maximum. The crack driving force was found to be the alternating mixed-mode loading that occurs with each passage of the most heavily loaded roller. The crack is predicted to propagate in a step-like fashion alternating between radial and inclined segments, and this pattern was observed experimentally. The maximum changes DeltaG sub 00 and DeltaG sub r0 of the crack extension forces are found to be good measures of the crack propagation rate and direction.

Mendelson, A.↗
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 NTRS1987

On time-dependent radiative transfer

An integral equation is developed for application to time-dependent laboratory experiments in which partial redistribution effects are important. The equation of transport with the Heasly-Kneer (1976) emission coefficient and the equation of statistical equilibrium lead to a time-dependent redistribution function containing an absorption - reemission term which decays exponentially in time and a scattering term which is instantaneous. This integral equation does not agree with an equation written by Payne et al. (1974) that has been used to compare theory with experiments. The difference between the Payne equation and the equation developed here needs to be examined in detail, since it might under some circumstances be on the same order as the difference between partial and complete redistribution.

Streater, A.↗
NASA NTRS2000

Electromagnetic Field Penetration Studies

A numerical method is presented to determine electromagnetic shielding effectiveness of rectangular enclosure with apertures on its wall used for input and output connections, control panels, visual-access windows, ventilation panels, etc. Expressing EM fields in terms of cavity Green's function inside the enclosure and the free space Green's function outside the enclosure, integral equations with aperture tangential electric fields as unknown variables are obtained by enforcing the continuity of tangential electric and magnetic fields across the apertures. Using the Method of Moments, the integral equations are solved for unknown aperture fields. From these aperture fields, the EM field inside a rectangular enclosure due to external electromagnetic sources are determined. Numerical results on electric field shielding of a rectangular cavity with a thin rectangular slot obtained using the present method are compared with the results obtained using simple transmission line technique for code validation. The present technique is applied to determine field penetration inside a Boeing-757 by approximating its passenger cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows by rectangular apertures. Preliminary results for, two windows, one on each side of fuselage were considered. Numerical results for Boeing-757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are presented.

Deshpande, M.D.↗
NASA NTRS2003

Stochastic Mixing Model with Power Law Decay of Variance

Here we present a simple stochastic mixing model based on the law of large numbers (LLN). The reason why the LLN is involved in our formulation of the mixing problem is that the random conserved scalar c = c(t,x(t)) appears to behave as a sample mean. It converges to the mean value mu, while the variance sigma(sup 2)(sub c) (t) decays approximately as t(exp -1). Since the variance of the scalar decays faster than a sample mean (typically is greater than unity), we will introduce some non-linear modifications into the corresponding pdf-equation. The main idea is to develop a robust model which is independent from restrictive assumptions about the shape of the pdf. The remainder of this paper is organized as follows. In Section 2 we derive the integral equation from a stochastic difference equation describing the evolution of the pdf of a passive scalar in time. The stochastic difference equation introduces an exchange rate gamma(sub n) which we model in a first step as a deterministic function. In a second step, we generalize gamma(sub n) as a stochastic variable taking fluctuations in the inhomogeneous environment into account. In Section 3 we solve the non-linear integral equation numerically and analyze the influence of the different parameters on the decay rate. The paper finishes with a conclusion.

Fedotov, S.↗
NASA NTRS1991

A statistical approach to the brittle fracture of a multi-phase solid

A stochastic damage model is proposed to quantify the inherent statistical distribution of the fracture toughness of a brittle, multi-phase solid. The model, based on the macrocrack-microcrack interaction, incorporates uncertainties in locations and orientations of microcracks. Due to the high concentration of microcracks near the macro-tip, a higher order analysis based on traction boundary integral equations is formulated first for an arbitrary array of cracks. The effects of uncertainties in locations and orientations of microcracks at a macro-tip are analyzed quantitatively by using the boundary integral equations method in conjunction with the computer simulation of the random microcrack array. The short range interactions resulting from surrounding microcracks closet to the main crack tip are investigated. The effects of microcrack density parameter are also explored in the present study. The validity of the present model is demonstrated by comparing its statistical output with the Neville distribution function, which gives correct fits to sets of experimental data from multi-phase solids.

Liu, W. K.↗
NASA NTRS1992

A comparison between backscattering models for rough surfaces

The ease of applicability of three scattering models is examined. This is accomplished by considering the time taken to numerically evaluate these models and by comparing their predictions as a function of surface roughness, frequency, incident angle and polarization with the moment-method solution in two dimensions. The complexity of the analytic models in three dimensions and their analytic reduction to high- and low-frequency regions are also compared. The selected models are an integral equation model (IEM), a full wave model (FWM), and the phase perturbation model (PPM). It is noted that in three dimensions, the full-wave model requires an evaluation of a tenfold integral and the phase perturbation model requires a fourfold and a twofold integral, whereas the integral equation model is an algebraic equation in like polarization under single scattering conditions. It is found that both the IEM and PPM agree with the moment-method solution from low to high frequencies numerically.

Chen, K. S.↗
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.↗