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.
Optimal temperature control design using IBM 7094 digital computer program for thermal environment simulation of stabilized multisurfaced earth- orbiting spacecraft structure
As Pioneer 10 approaches the planet Jupiter, its trajectory will carry it behind the planet, as seen from the earth. The microwave telemetry signal from the spacecraft will then pass through the Jovian atmosphere on its way to the earth, and analysis of the effects of the atmosphere on the signal should yield information on the atmospheric structure. To help define the capabilities of this experiment and to aid the analysis of the real data, we have simulated the occultation experiment on a computer by using model Jovian atmospheres and a predicted trajectory. The results indicate that the useful depth of penetration of the microwave signal will be governed by ammonia absorption rather than superrefractivity. The experiment is apparently far more sensitive to variations in cloud-top temperature than to variations in the hydrogen-helium ratio. Characteristic 'flashes' in the attenuation record of the signal may yield information on the location of the top of the Jovian convective zone.
The formulation and employment of a computer code designed to simulate the directional solidification of lead-rich Pb-Sn alloys in the form of an ingot with a uniform and circular cross-section are described. The formulation is for steady-state solidification in which convection in the all-liquid zone is ignored. Particular attention was given to designing a code to simulate the effect of a subtle variation of temperature in the radial direction. This is important because a very small temperature difference between the center and the surface of the ingot (e.g., less than 0.5 C ) is enough to cause substantial convection within the mushy-zone when the solidification rate is approximately 0.001 to 0.0001 cm/s.
The goal of this project is to develop a magneto-hydrodynamic (MHD) computer code for simulation of the VASIMR engine. This code is designed be easy to modify and use. We achieve this using the Cactus framework, a system originally developed for research in numerical relativity. Since its release, Cactus has become an extremely powerful and flexible open source framework. The development of the code will be done in stages, starting with a basic fluid dynamic simulation and working towards a more complex MHD code. Once developed, this code can be used by students and researchers in order to further test and improve the VASIMR engine.
A combination of computational fluid dynamic analysis and analytical solutions is being used to characterize the dominant modes in liquid rocket engines in conjunction with laboratory experiments. The analytical solutions are based on simplified geometries and flow conditions and are used for careful validation of the numerical formulation. The validated computational model is then extended to realistic geometries and flow conditions to test the effects of various parameters on chamber modes, to guide and interpret companion laboratory experiments in simplified combustors, and to scale the measurements to engine operating conditions. In turn, the experiments are used to validate and improve the model. The present paper gives an overview of the numerical and analytical techniques along with comparisons illustrating the accuracy of the computations as a function of grid resolution. A representative parametric study of the effect of combustor mean flow Mach number and combustor aspect ratio on the chamber modes is then presented for both transverse and longitudinal modes. The results show that higher mean flow Mach numbers drive the modes to lower frequencies. Estimates of transverse wave mechanics in a high aspect ratio combustor are then contrasted with longitudinal modes in a long and narrow combustor to provide understanding of potential experimental simulations.
No abstract available
A finite element model of cardiac conduction was used to simulate two experimental protocols: 1) fibrillation threshold measurements and 2) clinical electrophysiologic (EP) testing procedures. The model consisted of a cylindrical lattice whose properties were determined by four parameters: element length, conduction velocity, mean refractory period, and standard deviation of refractory periods. Different stimulation patterns were applied to the lattice under a given set of lattice parameter values and the response of the model was observed through a simulated electrocardiogram. The studies confirm that the model can account for observations made in experimental fibrillation threshold measurements and in clinical EP testing protocols.
Equipment available in the primary training facility for the Space Transportation System (STS) flight crews includes the Fixed Base Simulator, the Motion Base Simulator, the Spacelab Simulator, and the Guidance and Navigation Simulator. The Shuttle Mission Simulator (SMS) consists of the Fixed Base Simulator and the Motion Base Simulator. The SMS utilizes four visual Computer Generated Image (CGI) systems. The Motion Base Simulator has a forward crew station with six-degrees of freedom motion simulation. Operation of the Spacelab Simulator is planned for the spring of 1983. The Guidance and Navigation Simulator went into operation in 1982. Aspects of orbital visual simulation are discussed, taking into account the earth scene, payload simulation, the generation and display of 1079 stars, the simulation of sun glare, and Reaction Control System jet firing plumes. Attention is also given to landing site visual simulation, and night launch and landing simulation.
The simulation of converter-controller combinations by means of a flexible digital computer program which produces output to a graphic display is discussed. The procedure is an alternative to mathematical analysis of converter systems. The types of computer programming involved in the simulation are described. Schematic diagrams, state equations, and output equations are displayed for four basic forms of inductor-energy-storage dc to dc converters. Mathematical models are developed to show the relationship of the parameters.
The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).
A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.
The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.
A computer program has been developed to solve the compressible flow equation for the velocity potential. The exterior of the airfoil is mapped onto the unit circle and the flow is computed on a grid in the circle plane. A relaxation method using backward differencing in the flow direction at supersonic points permits solutions for large supersonic areas. The pressure distribution resulting from the flow becomes the input to the von Karman momentum equation which when integrated gives the displacement thickness. This displacement thickness is smoothed and added to the airfoil to account for the turbulent boundary layer. The boundary layer correction is computed iteratively with the flow. Results from this program and test data agree well.
A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.
A computational study has been conducted to gain a better understanding of the role of turbulence in the merger of corotating vortices and in the decay of a single vortex. The flow field is simulated by using an inviscid flow model and zero-, one-, and two-equation turbulent-flow models. An implicit finite-difference procedure is employed to integrate the unsteady, two-dimensional, incompressible flow equations in a cross-plane. In the zero- and one-equation models, a spatially varying mixing-length model that simulates the streamline curvature effect is used to define an eddy viscosity. The computational merging distance criterion is in good agreement with the experimental one. In the case of the two-equation model, the turbulent kinetic energy and a modified rate of dissipation equation incorporating a curvature correction term are solved. A comparative study of the computational results corresponding to the various flow models and the available experimental data is presented.
Far infrared collision induced absorption spectra of gaseous CO2 were computed using molecular dynamics simulations. The quadrupole and hexadecapole multipolar induction, through the trace, and the anisotropy of the molecular polarizability were found to be insufficient to represent properly the dipole induction mechanism. For a detailed analysis of the induction process the spectra obtained were decomposed into components resulting from different terms of the induced dipole. Based on this decomposition, all additional overlap contribution for each term was proposed. When spectra were recomputed including such overlap, good agreement between experiment and simulation was achieved over the temperature range at which measurements exist (233-400 K). The use of an anisotropic intermolecular potential was found to be of critical importance for obtaining the right shape of the far wings of the spectra.
The automation of a small, solid propellant mixer was studied. Temperature control is under investigation. A numerical simulation of the system is under development and will be tested using different control options. Control system hardware is currently being put into place. The construction of mathematical models and simulation techniques for understanding various engineering processes is also studied. Computer graphics packages were utilized for better visualization of the simulation results. The mechanical mixing of propellants is examined. Simulation of the mixing process is being done to study how one can control for chaotic behavior to meet specified mixing requirements. An experimental mixing chamber is also being built. It will allow visual tracking of particles under mixing. The experimental unit will be used to test ideas from chaos theory, as well as to verify simulation results. This project has applications to extraterrestrial propellant quality and reliability.