SEARCH · Search NASA
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.
The theory of consistent step-by-step methods for solving Volterra integral equations is extended to nonsingular Volterra integro-differential equations. It is shown that standard step-by-step algorithms for these more general equations are convergent. Several numerical examples are included.
Numerical solution of boundary value problem for nonlinear integro-differential equation system of hypersonic boundary layer
The response of a clamped panel to supersonically convected turbulence is considered. A theoretical model in the form of an integro-differential equation is employed that takes into account the coupling between the panel motion and the surrounding acoustic medium. The kernels of the integrals, which represent induced pressures due to the panel motion, are Green's functions for sound radiations under various moving and stationary sources. An approximate analysis is made by following a finite-element Ritz-Galerkin procedure. Preliminary numerical results, in agreement with experimental findings, indicate that the acoustic damping is the controlling mechanism of the response.
Explicit difference equations are presented for the solution of a signal of arbitrary waveform propagating in an ohmic dielectric, a cold plasma, a Debye model dielectric, and a Lorentz model dielectric. These difference equations are derived from the governing time-dependent integro-differential equations for the electric fields by a finite difference method. A special difference equation is derived for the grid point at the boundary of two different media. Employing this difference equation, transient signal propagation in an inhomogeneous media can be solved provided that the medium is approximated in a step-wise fashion. The solutions are generated simply by marching on in time. It is concluded that while the classical transform methods will remain useful in certain cases, with the development of the finite difference methods described, an extensive class of problems of transient signal propagating in stratified dispersive media can be effectively solved by numerical methods.
The failure of most viscous-inviscid interaction methods at strong interactions is attributed to the presence of a normal pressure gradient. A new theory is proposed for supersonic laminar boundary layers that can generate normal pressure gradients. The Navier-Stokes equations are reexamined by an order of magnitude analysis and all first and second order terms are retained. The approximation is found to be dependent not only on the boundary layer thickness but also on the ratio of the dimensionless viscosity and density. The equations are transformed into two quasi-similar, nonlinear, third order, ordinary integro-differential equations for the velocity and pressure as functions of a single transverse variable. The properties of the equations at the boundaries are discussed.
Radiative transfer of partially polarized radiation in an anisotropically scattering, inhomogeneous atmosphere containing an arbitrary polydispersion of particles is described using Jones's (1941) amplitude vectors and matrices. This approach exploits the close analogy between the quantum-mechanical states of spin 1/2 systems and the polarization states of electromagnetic radiation described by Jones's vector. The complete equivalence between the transport equation for Jones's vectors and the classical radiative-transfer equation for Stokes's intensity vectors is demonstrated in two independent ways after deriving the transport equations for the polarization coherency matrices and for the quaternions corresponding to the Jones's vectors. A compact-operator formulation of the theory is provided and used to derive the necessary equations for both a local and a global description of the transport of Jones's vectors. The integro-differential equations for the amplitude reflection and transmission matrices are derived and related to the usual corresponding equations.
The physical properties of the upper atmosphere are determined mainly by heat conduction, the heat sources and the barometric law. An analysis of the integro-differential equation which describes these physical processes has been carried out. It is found that heating of the thermosphere due to absorption of the solar extreme ultraviolet (EUV) radiation alone cannot explain the observed diurnal variation of density and temperature, since it would yield a maximum of these properties at about 17h local time, instead of 14h where it is observed. Secondly, if the EUV flux is adjusted to give the observed average temperature, then the diurnal variation in density would be much too large compared with the observed amplitude. Thirdly, it would require an extremely high efficiency for the conversion of EUV radiation into heat, if we compare the required flux with Hinteregger's measurements of the EUV flux. Thus, it is necessary to have another heat source in addition to the heating due to absorption of EUV radiation. If an additional beat source is used, which has a maximum at about 9h local time and a flux of 1 erg cm−1 sec−1, a time-dependent model of the upper atmosphere is obtained that is in good agreement with the observed densities. There is evidence that this additional heat source derives its energy ultimately from the solar corpuscular radiation. In this paper we present the results of calculations for a model in the equatorial and temperature zones of the earth, for those times when the average solar activity corresponds to a solar radiation flux of 200×110−22 Wm−2 (cps)−1 at 10.7-cm wavelength. The physical properties (temperature, density, pressure, scale height, mean molecular weight and the number densities of N2, O2, O, He and H) are given as a function of local time and for the altitudes between 120 km and 2050 km.
Integro-differential equation describing cosmic ray interaction with solar wind
An explicit definition of the radiation-stream concept is stated and applied to approximate the integro-differential equation of radiative transfer with a set of twelve coupled differential equations. Computational efficiency is enhanced by distributing the corresponding streams in three-dimensional space in a totally symmetric way. Polarization is then incorporated in this model. A computer program based on the model is briefly compared with a Monte Carlo program for simulation of horizon scans of the earth's atmosphere. It is found to be considerably faster.
Discussion of various zeroth-order approximations to s-wave scattering of electrons by hydrogen atoms below the first excitation threshold. The formalism previously developed by the author (1967, 1968) is applied to Feshbach operators to derive integro-differential equations, with the optical-potential set equal to zero, for the singlet and triplet cases. Phase shifts of s-wave scattering are computed in the zeroth-order approximation of the Feshbach operator method and in the static-exchange approximation. It is found that the convergence of numerical computations is faster in the former approximation than in the latter.
Final stretch criterion of failure is applied to the problem of quasi-static extension of a crack embedded in an elastic-plastic or viscoelastic-plastic matrix. The slow growth under subcritical conditions in a rate-sensitive Tresca solid is shown to be a superposition of creep rupture and McClintock's ductile growth. This type of growth occurs at subcritical magnitude of the imposed K-factor and can be accounted for only through a recognition of inelastic properties of solids. In the subcritical range there is no unique value for K sub c independent of geometrical configuration and flaw size. Not only the produced states of stress and strain are dependent on the loading path, but also the material resistance to fracture turns out to be a function of the history of loading that precedes catastrophic failure. A nonlinear integro-differential equation of motion is derived for a crack progressing through a viscoelastic medium with some limited ability to plastic flow. Examples of numerical integration are given incorporating both monotonic and cyclic loading programs.
A general theory of subsonic potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the velocity potential is obtained for both supersonic and subsonic flow. Under the small perturbation assumption, the potential at any point in the field depends only upon the values of the potential and its normal derivative on the surface of the body. On the surface of the body, this representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The theory is applied to finite-thickness wings in subsonic steady and oscillatory flows.
This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.
The general theory of potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the potential is obtained for both supersonic and subsonic flow. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface, sigma, of the body. Hence, if the point P approaches the surface of the body, the representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface sigma. For the important practical case of small harmonic oscillation around a rest position, the equation reduces to a two-dimensional Fredholm integral equation of second-type. It is shown that this equation reduces properly to the lifting surface theories as well as other classical mathematical formulas. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings.
In order to accommodate a quadrature amplitude-shift-keyed (QASK) signal, Simon and Smith (1974) have modified the decision-feedback loop which tracks a quadrature phase-shift-keyed (QPSK). In the investigation reported approaches are considered to modify the loops in such a way that offset QASK signals can be tracked, giving attention to the special case of an offset QPSK. The development of the stochastic integro-differential equation of operation for a decision-feedback offset QASK loop is discussed along with the probability density function of the phase error process.
The propagation of transient scalar waves in a three-dimensional random medium is considered. The analysis is based on the smoothing method. An integro-differential equation for the coherent (or average) wave is derived and solved for the case of a statistically homogeneous and isotropic medium and a delta-function source. This yields the coherent Green's function of the medium. It is found that the waveform of the coherent wave depends generally on the distance from the source measured in terms of a certain dimensionless parameter. Based on the magnitude of this parameter, three propagation zones, called the near zone, the far zone, and the intermediate zone, are defined.
A simple probabilistic model was constructed for the average value of a cosmogenic nuclide as a function of depth in a regolith. An arbitrary function was chosen for the size distribution of craters. The resulting integro-differential equation was found to reduce in limiting cases to the marching equation with a characteristic residence time and to the diffusion equation. The regolith diffusion constant is shown to be a simple integral of the cratering rate weighted by geometrical terms. This formal treatment provides a direct and general connection between cosmogenic nuclides and cratering rates and crater population in a simple analytical form. The validity of this model remains to be tested.
Subcritical crack growth in a linear viscoelastic material subjected to cyclic loading is investigated starting with the thermodynamic power balance. Physically, it is supposed that the subcritical slow crack growth is due to local weakening of the material in the neighborhood of the crack tip. Under fatigue loading it is assumed that all energy dissipation goes into heat and that this heat build-up is the dominant irreversible process governing crack growth. The cycle averaged temperature distribution around the crack tip is obtained from local application of the first law of thermodynamics (i.e., conservation of energy). The analysis leads to the solution of a nonlinear integro-differential equation for crack length as a function of time which is coupled with the local energy equation. A regular perturbation technique is used to obtain an analytical solution which compares very well with experimental results.