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.
The theory of homogeneous condensation is reviewed and equations describing this process are presented. Numerical computer solutions to transient problems in nucleation (relaxation to steady state) are presented and compared to a prior computation.
The problem of designing optimal inputs in the identification of linear systems with unknown random parameters is considered using a Bayesian approach. The information matrix, which is positive definite for the class of systems analyzed, gives a measure of performance for the system inputs. The computation of the optimal closed-loop input mappings is shown to be a nontrivial exercise in adaptive control. Deterministic optimal inputs are shown to be easily computable. Numerical examples are given. A Kalman filter is used to estimate the parameters. A necessary condition for the Kalman filter not to diverge when applying linear feedback is also given.
The thermal relaxation process for a spatially uniform two dimensional plasma in a uniform dc magnetic field is simulated numerically. Thermal relaxation times are defined in terms of the time necessary for the numerically computer Boltzman H-function to decrease through a given part of the distance to its minimum value. Dependence of relaxation time on two parameters is studied: number of particles per Debye square and ratio of gyrofrequency to plasma frequency.
Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published companion paper on finite-element appendage vibration equations, this paper provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.
Four W Ursae Majoris-type binaries were investigated using the simple model of blackbody solid spheres. Numerical computations were carried out for the rotationally broadened profiles of spectral lines Sr II (4077.7) and Ca I (4226.7) and the mean radial velocities. These mean velocities at different phases outside eclipses were compared with sinusoidal values of the radial velocity curve for point masses. Within the present accuracy of observations, the reflection effect is small and may be neglected in the analysis of radial velocity curves.
Review of recent work in recasting the energy type of approximation to aircraft flight in terms of singular perturbation theory and in extending this theory to three-dimensional maneuvers. Singular perturbations for differential equations arising in optimal control are first examined for a system of fairly general form but low order. The attitude dynamics for optimal flight of a rocket in vacuum is then studied as an introductory, and fairly transparent, example. The question of the choice of variables is then discussed. Finally, optimal aircraft flight in various reduced-order approximations is investigated. In particular, the problem of three-dimensional aircraft flight is formulated for singular perturbation treatment, and possibilities for decoupling into several lower-order problems are illustrated. The use of reduced-order approximation facilitates numerical computations by reducing the number of multiplier initial values that must be determined simultaneously and by improving the conditioning of the differential equations.
In order to investigate the dipolar and nondipolar contributions of spacecraft magnetic fields, a simple magnetic field model is proposed. This model consists of a number N of randomly oriented dipoles of strength M sub k in a given volume. Two sets of formulas are presented that give the rms multipole field components, first, for isotropic orientations of the dipoles at given positions and second, for isotropic orientations of the dipoles distributed uniformly throughout a cube or sphere. The statistical results for a cube (8 cu m in size) together with individual examples computed numerically show the following features. Beyond about 2- to 3-m distance from the center of the cube the field is dominated by an equivalent dipole. The magnitude of the magnetic moment of the dipolar part is approximated by N to the 1/2 power times M for equal magnetic moments or generally by the Pythagorean sum of the dipole moments. The radial component tend to be greater than either of the transverse components for the dipole portion as well as for the nondipolar field contributions.
The unsteady, inviscid flowfield that results when a supersonic vehicle strikes a planar oblique shock wave, though difficult to simulate experimentally, is quite easy to model and compute numerically. The complicated flowfield, which contains multiple shock wave interactions, is determined using a second-order, shock-capturing, finite-difference approach which solves the time-dependent Euler equations under a self-similar transformation. A series of numerical results for a simple two-dimensional wedge is presented which describes the entire disturbed region, including the wave structure, and shows good agreement with the available two- and three-dimensional experimental data.-
Results from numerical computations of the diffuse sky radiances and degree of polarization are given. These results were obtained by solving the equation of radiative transfer in the formulation of Eschelbach. In these computations, multiple scattering and absorption by aerosol particles were considered. In addition, results from experimental measurements of the complex refractive index of aerosol samples at different relative humidities are presented. From the radiation flux divergences which were computed based on the determined properties of atmospheric aerosol particles, atmospheric heating rates were derived which were found to be comparable to the heating rates by water vapor.
A numerical computing procedure was developed for calculating vaporization histories of evaporating drops in a combustor in which travelling transverse oscillations occurred. The liquid drop was assumed to have a finite thermal conductivity. The system of equations was solved by using a finite difference method programmed for solution on a high speed digital computer. Oscillations in the ratio of vaporization of an array of repetitivity injected drops in the combustor were obtained from summation of individual drop histories. A nonlinear in-phase frequency response factor for the entire vaporization process to oscillations in pressure was evaluated. A nonlinear out-of-phase response factor, in-phase and out-of-phase harmonic response factors, and a Princeton type 'n' and 'tau' were determined. The resulting data was correlated and is presented in graphical format. Qualitative agreement with the open literature is obtained in the behavior of the in-phase response factor. Quantitatively the results of the present finite conductivity spray analysis do not correlate with the results of a single drop model.
A terminal guidance and navigation scheme developed in earlier work was modified and evaluated for a solar electric propulsion rendezvous mission to comet Encke. The scheme is intended for autonomous, on-board use. The guidance algorithm is based on optimal control theory and minimizes the time integrated square of thrust acceleration. The navigation algorithm employs a modified Kalman filter set in measurement variables. Random sequences were generated to simulate measurement errors, and the evaluation was conducted with detailed numerical computations which include actual motions of spacecraft and comet. The evaluations showed that the scheme attains rendezvous and maintains station after rendezvous within less than 10 km for estimated best measurements and within less than 100 km for estimated worst measurements. The measurements required are angles, range, and range rate. Angles and range appear to be absolutely necessary; range rate is not as strong a measurement type, and further modifications of the filter will allow a scheme that does not require the rate measurements.
Several expressions for the determination of the acoustic field of moving bodies are presented. The analysis is based on the Ffowcs Williams-Hawkings equation. Applying some proposed criteria, one of these expressions is singled out for numerical computation of acoustic pressure signature. The compactness of sources is not assumed and the main results are not restricted by the observer position. The distinction between compact and noncompact sources on moving surfaces is discussed. Some thickness noise calculations of helicopter rotors and comparison with experiments are included which suggest this mechanism as the source of high-speed blade slap of rotors.
The CAD/CAM interactive computer graphics system was described; uses to which it has been put were shown, and current developments of the system were outlined. The system supports batch, time sharing, and fully interactive graphic processing. Engineers using the system may switch between these methods of data processing and problem solving to make the best use of the available resources. It is concluded that the introduction of on-line computing in the form of teletypes, storage tubes, and fully interactive graphics has resulted in large increases in productivity and reduced timescales in the geometric computing, numerical lofting and part programming areas, together with a greater utilization of the system in the technical departments.
A tracking filter that includes an output that is logarithmically proportional to the amplitude of the filtered signal is used to measure the damping exponent of a transient decay. This method is analogous to the digital technique referred to as the moving-block or peak-plot method. The method is simple to use and quite accurate, avoids the numerical computations associated with similar digital techniques, but is constrained by the poor time-domain response characteristics of commercial tracking filters presently available.
The use of the R sub 0 criterion for modulation system design is investigated for channels with non-white Gaussian noise. A signal space representation of the waveform channel is developed, and the cut-off rate R sub 0 for vector channels with additive nonwhite Gaussian noise and unquantized demodulation is derived. When the signal unput to the channel is a continuous random vector, maximization of R sub 0 with constrained average signal energy leads to a water-filling interpretation of optimal energy distribution in signal space. The necessary condition for a finite signal set to maximize R sub 0 with constrained energy and an equally likely probability assignment of signal vectors is presented, and an algorithm is outlined for numerically computing the optimum signal set. A necessary condition on a constrained energy, finite signal set is found which maximizes a Taylor series approximation of R sub 0. This signal set is compared with the finite signal set which has the water-filling average energy distribution.
A lifting surface theory based on a parallel shear flow model is presented for steady, incompressible flows. The theory is intended to account approximately for the presence of a boundary layer. The method of Fourier transforms is used to calculate the pressure on a surface of infinite extent and arbitrary contour. Immediately above the surface is a region of sheared flow (the boundary layer), outside of which the flow velocity is constant. The Fourier transform of the pressure on this surface is used to derive the shear flow equivalent to the kernel function of classical potential flow lifting surface theory. The kernel function provides an integral relation between the upwash at a given point on the surface and the pressure everywhere on the surface. This relation is treated as an integral equation for the pressure, and is solved numerically. Computations are presented for the lift and pitching moment on a flat plate in two-dimensional flow, and for flat, rectangular wings of aspect ratio 1, 2, and 5. As expected, the shear layer decreases the lift curve slope; however, the shear layer (whose thickness is constant along the wing chord) has little effect on the center of pressure.
Using the dissipation-fluctuation approach, the brightness temperature of a stratified medium with inhomogeneous permittivities and nonuniform temperature profiles is solved. The solution is expressed in closed form ready for machine computation. Numerical results are illustrated and compared with closed-form analytical solutions and results obtained from the WKB method for simple profiles.
Comparison of the Jones-Launder, Ng-Spalding, Saffman-Wilcox, and Wilcox-Traci two-equation turbulence models has been conducted. It was shown that the Saffman-Wilcox and Wilcox-Traci dissipation-rate formulations admit straightforward integration through the viscous sublayer, whereas integration through the viscous sublayer is a more difficult issue with the Jones-Launder dissipation-function and the Ng-Spalding length-scale formulations. Numerical computations were conducted in which the models were applied to four equilibrium boundary layer flows including adverse, zero, and favorable pressure gradients. Computations of zero pressure gradient flow over a convex wall composed the final part of the comparison.