Gravitational experiments with a collisionless two-dimensional computer model
Two dimensional computer model for collisionless self gravitating galactic system with gravitational field obtained from Poisson equation
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.
Two dimensional computer model for collisionless self gravitating galactic system with gravitational field obtained from Poisson equation
Discussion of the interactive man-computer system (WIMDCO) for measuring cloud motion from ATS pictures. The multipicture model of ATS navigation is used for consecutive ATS digital picture alignment to provide the required precision in cloud motion measurements by automated tracking techniques. Cloud motion is measured by tracking clouds between two digital pictures with the aid of two-dimensional cross correlation analysis. The fast Fourier transform method and other programming techniques are used for time and cost minimization. Cloud motion vectors were computed over three time intervals in July, 1969. The EW and NS components were reproducible within roughly 2 knots.
A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.
Attention is drawn to recent advancements in the fluctuating line-reversal temperature measurement, the development of the two-dimensional drag-sensing probe into a three-dimensional drag-sensing probe, and the fluctuating density gradient cross beam Schlieren technique. An experimental apparatus is explained whereby the temperature fluctuation in a pulsating air-fuel ratio Bunsen burner is measured by means of sodium D-line reversal methods with a new photoelectric circuit to obtain both fluctuating and mean temperature data. The three-D drag probe is made possible by a newly invented differential and total magnetic sensing system that separates signals due to three orthogonal movements so that the vectorial momentum fluctuation can be measured. In the case of the Schlieren technique, the effect of anisotropic density gradients with respect to the angle between the polarization plane and the knife edge of a laser Schlieren system is studied.
The structure of two-dimensional turbulent jets and wakes is studied using two eddy viscosity models. The turbulent energy equation is used together with the mean momentum equations, and the system is closed by introducing eddy coefficients. The fine structure turbulence is obtained by applying proper boundary conditions at the mean turbulent interface position. The results, after a second averaging process, are compared with available experimental data.
A very general method for calculating compressible three-dimensional laminar and turbulent boundary layers on arbitrary wings is described. The method utilizes a nonorthogonal coordinate system for the boundary-layer calculations and includes a geometry package that represents the wing analytically. In the calculations all the geometric parameters of the coordinate system are accounted for. The Reynolds shear-stress terms are modeled by an eddy-viscosity formulation developed by Cebeci. The governing equations are solved by a very efficient two-point finite-difference method used earlier by Keller and Cebeci for two-dimensional flows and later by Cebeci for three-dimensional flows.
The first phases of a fundamental analytical study of STOL ground effects were presented. Ground effects were studied in two dimensions to establish the importance of nonlinear effects, to examine transient aspects of ascent and descent near the ground, and to study the modelling of the jet impingement on the ground. Powered lift system effects were treated using the jet-flap analogy. The status of a three-dimensional jet-wing ground effect method was presented. It was shown, for two-dimensional unblown airfoils, that the transient effects are small and are primarily due to airfoil/freestream/ground orientation rather than to unsteady effects. The three-dimensional study showed phenomena similar to the two-dimensional results. For unblown wings, the wing/freestream/ground orientation effects were shown to be of the same order of magnitude as for unblown airfoils. This may be used to study the nonplanar, nonlinear, jet-wing ground effect.
Flutter of two parallel elastically coupled flat plates was investigated analytically. A closed-form solution including both aerodynamic and structural damping is presented for flutter of flat orthotropic plates coupled by an elastic medium. Both plates are simply supported along the side edges but are supported by deflectional, rotational, and torsional springs of arbitrary stiffness at the leading and trailing edges. Two-dimensional quasi-steady aerodynamics was utilized in the solution. Since the large number of variables present in the problem precludes extensive parametric studies, results are presented to indicate the basic flutter characteristics of coupled two-plate systems and to assess the validity of previously published modal solutions for similar problems.
The results reviewed are divided into two categories: those that relate two-dimensional shell theory to three-dimensional elasticity theory and those concerned with shell theory per se. In the second category results for general elastic systems that carry over, by specialization or analogy, to shells and results that are unique to shell theory itself are considered.
The sensitivity of silica heat shield requirements to gap width, tile edge radius, and heat transfer distribution within tile gaps was investigated. A two-dimensional thermal model was modified and used to determine the effect of two dimensional heat transfer distributions at high temperature reusable surface insulation edges on shuttle thermal protection system (TPS) requirements. The sensitivity of TPS requirements to coating thickness, emissivity, substructure thickness, and changes in gap heating for several locations on shuttle was also studied. An inverse solution technique was applied to temperature data obtained in the Ames 20 MW turbulent duct in order to examine the effect of tile edge radius on TPS requirements. The derived heating values were then used to predict TPS requirements. Results show that increasing tile radius reduces TPS requirements.
An interactive procedure was developed for supersonic viscous flows that can be used for either two-dimensional or axisymmetric configurations. The procedure is directed to supersonic internal flows as well as those supersonic external flows that require consideration of mutual interaction between the outer flow and the boundary layer flow. The flow field is divided into two regions: an inner region which is highly viscous and mostly subsonic and an outer region where the flow is supersonic and in which viscous effects are small but not negligible. For the outer region a numerical solution is obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines. The inner region is treated by a system of equations of the boundary layer type that includes higher order effects such as longitudinal and transverse curvature and normal pressure gradients. These equations are coupled and solved simultaneously in the physical coordinates by using an implicit finite difference scheme. This system can also be used to calculate laminar and turbulent boundary layers using a scalar eddy viscosity concept.
Tests were conducted to demonstrate a rapid scan two dimensional laser velocimeter (LV) measurement technique for aerial applications research. The LV system is capable of simultaneously measuring both vertical and axial flow velocity components in a near or far field vortex system. Velocity profiles were successfully measured in the wake vortex of a representative agricultural aircraft model, with the vortex system rapidly transporting in ground effect. Results indicate that the laser velocimetry technique can provide quantitative information of wake vortex characteristics in ground effect.
An interactive numerical procedure has been developed for supersonic viscous flows (either two-dimensional or axisymmetric configurations). The flow field is divided into two regions: (1) an inner region which is highly viscous and mostly subsonic, and (2) an outer region where the flow is supersonic and in which viscous effects are small, but not negligible. This paper presents a detailed description of: I. Outer Region - numerical solution obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines; II. Inner Region - treated by a system of equations of the boundary-layer type that includes higher order effects, such as longitudinal and transverse curvature and normal pressure gradients (equations are coupled and solved simultaneously in physical coordinates, using an implicit finite-difference scheme); III. Interactive Procedure - in the interaction mode, the two regions are coupled iteratively along a matching line, where the Mach number is of the order of 1.2.
A ceramic/metallic aircraft gas turbine outer gas path seal designed for improved engine performance was studied. Transient temperature and stress profiles in a test seal geometry were determined by numerical analysis. During a simulated engine deceleration cycle from sea-level takeoff to idle conditions, the maximum seal temperature occurred below the seal surface, therefore the top layer of the seal was probably subjected to tensile stresses exceeding the modulus of rupture. In the stress analysis both two- and three-dimensional finite element computer programs were used. Predicted trends of the simpler and more easily usable two-dimensional element programs were borne out by the three-dimensional finite element program results.
A two dimensional oscillating airfoil test apparatus is presented as a method of measuring unsteady aerodynamic forces on an airfoil or rotor blade section. The oscillating airfoil test rig, which is being built for use in an 11 X 11-foot transonic wind tunnel (speed range M = 0.4 - 1.4), will allow determination of unsteady loadings and detailed pressure distributions on representative airfoil sections undergoing simulated pitching and flapping motions. The design details of the motion generating system and supporting structure are presented. This apparatus is now in the construction phase.
The inlet was investigated with both a two-oblique-shock and an isentropic compression ramp. Provisions for a variety of bypass systems were incorporated and none decreased pressure recovery by more than about 1 percent. By utilizing throat boundary-layer control, critical pressure recoveries of 64 and 70 percent were obtained with the two-shock and isentropic ramps, respectively. Diffuser-exit air distortion was about 5 percent at critical operation Theoretical cowl pressure drag was equal to about 10 percent of the net thrust of an assumed engine.
The extension to a compressible fluid of Lagally's theorem on the moment on a body in an incompressible fluid and Poggi's method of treating the flow of compressible fluids are employed for the determination of the effect of compressibility on the moment of an arbitrary body. Only the case of the two-dimensional subsonic flow of an ideal compressible fluid is considered. As examples of the application of the general theory, two well-known systems of profiles are treated; namely the elliptic profile and the symmetrical Joukowski profiles with sharp trailing edges. The effect of compressibility on the position of the center of pressure is also discussed. Several numerical examples of both the elliptic and the Joukowski profiles are given.
For the purpose of testing the theory of an oscillating airfoil of two degrees of freedom, a wing was mounted in the wind tunnel between two walls in such a manner that it could execute vertical (flexural) oscillations as well as torsional oscillations about an arbitrary axis parallel to the span. It was possible to vary the inertia and elasticity parameters and also to increase artificially the negligibly small natural damping of the system. The oscillations were recorded to a strongly magnified scale. The experimentally determined critical (or flutter) velocities fully agree with the theoretical ones of Wagner and Kussner within the limits of computational and measuring accuracy. An extremely narrow wing without end walls (three dimensional problem) showed the same oscillations as one with end walls (two-dimensional problem).