Search NASA⌕ Search

SEARCH · Search NASA

Results for “momentum conservation”

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.

At least 19 records

A discrete momentum-conserving explicit algorithm for multibody dynamics analysis

A discrete momentum-conserving, explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

A discrete momentum-conserving explicit algorithm for rigid body dynamics analysis

A discrete momentum-conserving explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

Direct Observation of Momentum Conservation at the Au/Si Interface Using BEEM (abstract)

The verification of parallel momentum conservation at a metal/semiconductor interface is a fundamental issue in interface transport. Ballistic-electron-microscopy (BEEM) is a recently developed method for probing interfaces with nanometer resolution. In this talk, BEEM spectroscopy on Au/Si(111) structures as a function of Au thickness and temperature will be described.

metal/semiconductor interface ballistic-electron-e↗

Evidence of Momentum Conservation at a Nonepitaxial Metal/Semiconductor Interface Using Ballistic Electron Emission Microscopy

Ballistic-Electron-Emission Microscopy (BEEM) spectroscopy has been performed on Au/Si(111) structures as a function of Au thickness and temperature. At 77 K a direct signature of parallel momentum conservation at the Au/Si interface is observed in the BEEM spectra. The variation in spectral shape with both Au thickness and temperature places restrictions on allowable values of inelastic and elastic mean-free paths in the metal, and also requires the presence of multiple electron passes within the Au layer. An independent indication of multiple reflections is directly observed in the attenuation of BEEM current with Au thickness.

Bell, L. D.↗

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗

IUE high dispersion spectra of luminous stars in symmetric nebulae

The objects are + 60 deg 2522 in NGC 7635, Sigma Ori AB and Sigma Ori E centered in a faint elliptical nebula, and AG Car in its nebula. Definite spectra are identified and classified according to photospheric, P Cygni, or interstellar types. Estimates of equivalent width or FWHM are included in the tables. Velocities are given and discussed. The extreme ultraviolet edges of P Cygni profiles imply terminal wind speeds that improve to extend earlier velocity data for use in the theories of energy conserving or momentum conserving modes of interaction of stellar winds with ambient interstellar media. The respective modes lead to mass loss rates of 1-10 x 10 to the -7th solar masses per yr for + 60 deg 2522 and 6-14 x 10 to the -6th solar masses per yr for AG Car. The expansion velocity for the Sigma Ori nebula is not yet known for application of the interaction theory.

Johnson, H. M.↗

Momentum Management Tool for Low-Thrust Missions

A momentum management tool was designed for the Dawn low-thrust interplanetary spacecraft en route to the asteroids Vesta and Ceres, in an effort to better understand the early creation of the solar system. Momentum must be managed to ensure the spacecraft has enough control authority to perform necessary turns and hold a fixed inertial attitude against external torques. Along with torques from solar pressure and gravity-gradients, ion-propulsion engines produce a torque about the thrust axis that must be countered by the four reaction wheel assemblies (RWA). MomProf is a ground operations tool built to address these concerns. The momentum management tool was developed during initial checkout and early cruise, and has been refined to accommodate a wide range of momentum-management issues. With every activity or sequence, wheel speeds and momentum state must be checked to avoid undesirable conditions and use of consumables. MomProf was developed to operate in the MATLAB environment. All data are loaded into MATLAB as a structure to provide consistent access to all inputs by individual functions within the tool. Used in its most basic application, the Dawn momentum tool uses the basic principle of angular momentum conservation, computing momentum in the body frame, and RWA wheel speeds, for all given orientations in the input file. MomProf was designed specifically to be able to handle the changing external torques and frequent de - saturations. Incorporating significant external torques adds complexity since there are various external torques that act under different operational modes.

Swenka, Edward R.↗

Redundancy of space manipulator on free-flying vehicle and its nonholonomic path planning

The nonholonomic mechanical structure of space robots and path planning is discussed. The angular momentum conservation works as a nonholonomic constraint while the linear momentum conservation is a holonomic one. Thus, a vehicle with a 6 d.o.f. manipulator is described as a 9 variable system with 6 inputs. This implies the possibility of controlling the vehicle orientation and the joint variables of the manipulator by actuating the joint variables, but only if the trajectory is carefully planned; however, both of them cannot be controlled independently. It means that by assuming feasible-path planning, a system that consists of a vehicle and a 6 d.o.f. manipulator can be utilized as 9 d.o.f. system. Initially, the nonholonomic mechanical structure of space vehicle/manipulator system is shown. Then a path planning scheme for nonholonomic systems is proposed using Lyapunov functions.

Nakamura, Yoshihiko↗

Numerical Modeling of Flow Distribution in Micro-Fluidics Systems

This paper describes an application of a general purpose computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels. GFSSP employs a finite volume formulation of mass and momentum conservation equations in a network consisting of nodes and branches. Mass conservation equation is solved for pressures at the nodes while the momentum conservation equation is solved at the branches to calculate flowrate. The system of equations describing the fluid network is solved by a numerical method that is a combination of the Newton-Raphson and successive substitution methods. The numerical results have been compared with test data and detailed CFD (computational Fluid Dynamics) calculations. The agreement between test data and predictions is satisfactory. The discrepancies between the predictions and test data can be attributed to the frictional correlation which does not include the effect of surface tension or electro-kinetic effect.

Majumdar, Alok↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Human Outer Solar System Exploration via Q-Thruster Technology

Propulsion technology development efforts at the NASA Johnson Space Center continue to advance the understanding of the quantum vacuum plasma thruster (QThruster), a form of electric propulsion. Through the use of electric and magnetic fields, a Q-thruster pushes quantum particles (electrons/positrons) in one direction, while the Qthruster recoils to conserve momentum. This principle is similar to how a submarine uses its propeller to push water in one direction, while the submarine recoils to conserve momentum. Based on laboratory results, it appears that continuous specific thrust levels of 0.4 - 4.0 N/kWe are achievable with essentially no onboard propellant consumption. To evaluate the potential of this technology, a mission analysis tool was developed utilizing the Generalized Reduced Gradient non-linear parameter optimization engine contained in the Microsoft Excel® platform. This tool allowed very rapid assessments of "Q-Ship" minimum time transfers from earth to the outer planets and back utilizing parametric variations in thrust acceleration while enforcing constraints on planetary phase angles and minimum heliocentric distances. A conservative Q-Thruster specific thrust assumption (0.4 N/kWe) combined with "moderate" levels of space nuclear power (1 - 2 MWe) and vehicle specific mass (45 - 55 kg/kWe) results in continuous milli-g thrust acceleration, opening up realms of human spaceflight performance completely unattainable by any current systems or near-term proposed technologies. Minimum flight times to Mars are predicted to be as low as 75 days, but perhaps more importantly new "retro-phase" and "gravity-augmented" trajectory shaping techniques were revealed which overcome adverse planetary phasing and allow virtually unrestricted departure and return opportunities. Even more impressively, the Jovian and Saturnian systems would be opened up to human exploration with round-trip times of 21 and 32 months respectively including 6 to 12 months of exploration at the destinations. Finally, interstellar trip times are assessed at milli-g acceleration levels.

Joosten, B. Kent↗

Causal compensated perturbations in cosmology

A theoretical framework is developed to calculate linear perturbations in the gravitational and matter fields which arise causally in response to the presence of stiff matter sources in a FRW cosmology. It is shown that, in order to satisfy energy and momentum conservation, the gravitational fields of the source must be compensated by perturbations in the matter and gravitational fields, and the role of such compensation in containing the initial inhomogeneities in their subsequent evolution is discussed. A complete formal solution is derived in terms of Green functions for the perturbations produced by an arbitrary source in a flat universe containing cold dark matter. Approximate Green function solutions are derived for the late-time density perturbations and late-time gravitational waves in a universe containing a radiation fluid. A cosmological energy-momentum pseudotensor is defined to clarify the nature of energy and momentum conservation in the expanding universe.

Veeraraghavan, Shoba↗

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

Phase space factors in multiparticle processes

General phase space theorems are discussed for the cases (A) with only energy conservation applied and (B) with energy and momentum conservation applied. It is shown that in the non-relativistic limit for N particles there is a very close relationship between the multiparticle phase space integral in case B and that for case A and N-1 particles.

Gould, R. J.↗

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗