SEARCH · Engineering Papers
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.
A Fredholm integral equation is solved in closed form by reducing it to a standard boundary value problem for a sectionally analytic function. The procedure maybe extended to more general classes of functions whose Fourier transforms only exist along certain lines in the complex k-plane. Examples are used to illustrate the results.
The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The solution involves both ideal-MHD and resistive-MHD processes, and the solution shows an evolution which is remarkably suggestive of the preflare, impulsive, and main phases of the flare-cycle.
For a flow over an airfoil with laminar separation, a feedback cycle may exist whereby a Kelvin-Helmholtz instability wave emanating from the separation point on the airfoil surface grows along the shear layer and is diffracted as it interacts with the sharp trailing edge of the airfoil, causing acoustic radiation which , in turn, propagates upstream and regenerates the initial instability wave. The analysis is restricted to the high frequency limit. Solutions to the boundary-value problem are obtained using the slowly varying approximation and the method of matched asymptotic expansions. Resonant solutions exist for certain discrete values of the Reynolds and Strouhal numbers. The results are discussed and compared with available data.
It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.
This paper concerns a numerical model of a regenerator running at very low temperatures. The model consists of the usual three equations for a compressible fluid with an additional equation for a matrix temperature. The main difficulty with the model is the very low Mach number (approximately 1.E-3). The divergence of the velocity is not small, the pressure divergence is small, and the pressure fluctuation in time is not small. An asymptotic expansion based on the bounded derivative method of Kreiss is used to give a reduced model which eliminates acoustic waves. The velocity is then determined by a two-point boundary value problem which does not contain a time derivative. The solution obtained from the reduced system is compared with the numerical solution of the original system.
In this paper a certain class of singular integral equations that may arise from the mixed boundary value problems in nonhomogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernel has strong singularities of the form (t-x) sup-2, x sup n-2 (t+x) sup n, (n or = 2, 0x,tb). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.
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.
A certain class of singular integral equations that may arise from the mixed boundary value problems in nonhonogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernal has strong singularities of the form (t-x)(-2), x(n-2) (t+x)(n), (n is = or 2, 0 x, t b). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.
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.
A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.
The analytic properties of the Jost functions for the inverse scattering transform associated with the forced Toda lattice are shown to determine the time evolution of this particular boundary value problem. It is suggested that inverse scattering methods may be used generally to analyze forced integrable systems. Thus an extension of the applicability of the inverse scattering transform is indicated.
The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.
Steady state harmonic stress waves in an isotropic elastic plate excited on one face by a circular transducer are analyzed theoretically. The transmitting transducer transforms an electrical voltage into a uniform normal stress at the top of the plate. To solve the boundary value problem,the radiation into a half-space is considered. The receiving transducer produces an electrical voltage proportional to the average spatially integrated normal stress over its face due to an incident wave. A numerical procedure is given to evaluate the frequency response at a receiving point due to a multiply reflected wave in the near field. Its stability and convergence are discussed. Parameterization plots which determine the particular wave whose frequency response has maximum magnitude compared with other multiple reflected waves are given for a range of values of dimensionless parameters. The effects of changes in the values of the parameters are discussed.
The purpose of the paper is to demonstrate that constitutive equations for moisture sorption and swelling can be obtained with an internal variable method which is rooted in the theory of irreversible thermodynamics. Realism is injected into the model through the concept of free volume. The power of the method is illustrated by the solution of the formidable nonlinear boundary value problem for diffusion with swelling-stress dependent diffusivity.
A mathematical model is developed to describe the effect of aircraft-engine inlet geometry on the reflected and radiated acoustic field without flow, as studied experimentally using a spinning-mode synthesizer by Silcox (1983). The acoustic pressure in the inlet interior and exterior is modeled by a pure cylindrical azimuthal mode for the Helmholtz equation with hardwall boundary and by the Helmholtz equation and the radiation condition at infinity, respectively. The analytical approach to the solution of the resulting boundary-value problem and the program implementation are explained; numerical results are presented in tables and graphs; and the uniqueness of the problem is demonstrated.
The standard method of multiple shooting for a system of n first-order differential equations with k unknown initial conditions requires the integration of k sets of variational equations on the first shot and n sets of variational equations on every shot thereafter. This paper describes a variant of multiple shooting that requires the solution of k sets of variational equations on every shot. The technique applies to both linear and nonlinear boundary-value problems. Techniques to deal with difficulties unique to the solution of nonlinear problems are suggested.
Report discusses computer simulations of 300- to 900-GHz triplers using nonideal GaAs Schottky diodes operating in varistor mode. Nonlinear boundary-value problem posed by diode state equations not solved in closed form. Consequently, such computer simulations needed to optimize tripler configuration and operating conditions.
Recent work has shown that the use of separable symmetric functions of the principal stretches can adequately describe the response of certain propellant materials and, further, that a data reduction scheme gives a convenient way of obtaining the values of the functions from experimental data. Based on representation of the energy, a computational scheme was developed that allows finite element analysis of boundary value problems of arbitrary shape and loading. The computational procedure was implemental in a three-dimensional finite element code, TEXLESP-S, which is documented herein.