Search NASA⌕ Search

SEARCH · Search NASA

Results for “comparative modeling”

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 253 records · Page 14

Energy coupling in the magnetospheres of earth and Mercury

The mechanisms involved in the dissipation of solar-wind energy during magnetospheric substorms are considered theoretically, comparing models for earth and Mercury. In the model for terrestrial substorms, IMF lines interconnect with terrestrial field lines near the front of the magnetosphere and are dragged back, carrying plasma and energy, to form tail lobes; a magnetic neutral region is then formed by reconnection of the open lines as the plasma sheet thins, and reconnective heating and acceleration of tail plasma lead to plasma inflow at the poles and formation of a plasmoid flowing down the tail at high velocity. Analogous phenomena on Mercury could produce precipitation of particles carrying 10-1000 GW of power into 'auroral zones' on the dark side of the planet. The feasibility of remote or in situ observations to detect such processes is discussed.

Baker, D. N.↗

Matrix fatigue crack development in a notched continuous fiber SCS-6/Ti-15-3 composite

In this study the extensive matrix fatigue cracking that has been observed in notched SCS-6/Ti-15-3 composites is investigated. Away from the notch a uniform spacing of the fatigue cracks develops. Closer to the notch, fiber-matrix debonding which occurs increases the crack spacing. Crack spacing and debond length determined from shear-lag cylinder models compare favorably with experimental observations. Scanning electron microscope (SEM) fractography showed that the principal fatigue crack initiation occurred around the zero degree fibers. Interface failure in the 90 degree plies does not lead to the development of the primary fatigue cracking.

Hillberry, B. M.↗

A model for the space shuttle main engine high pressure oxidizer turbopump shaft seal system

A simple static model is presented which solves for the flow properties of pressure, temperature, and mass flow in the Space Shuttle Main Engine pressure Oxidizer Turbopump Shaft Seal Systems. This system includes the primary and secondary turbine seals, the primary and secondary turbine drains, the helium purge seals and feed line, the primary oxygen drain, and the slinger/labyrinth oxygen seal pair. The model predicts the changes in flow variables that occur during and after failures of the various seals. Such information would be particularly useful in a post flight situation where processing of sensor information using this model could identify a particular seal that had experienced excessive wear. Most of the seals in the system are modeled using simple one dimensional equations which can be applied to almost any seal provided that the fluid is gaseous. A failure is modeled as an increase in the clearance between the shaft and the seal. Thus, the model does not attempt to predict how the failure process actually occurs (e.g., wear, seal crack initiation). The results presented were obtained using a FORTRAN implementation of the model running on a VAX computer. Solution for the seal system properties is obtained iteratively; however, a further simplified implementation (which does not include the slinger/labyrinth combination) was also developed which provides fast and reasonable results for most engine operating conditions. Results from the model compare favorably with the limited redline data available.

Paxson, Daniel E.↗

A continuous proportional low-thrust propulsion system

Helium gas flows continuously from a low-temperature dewar supplying propellant for satellite thrusters. Thrust is obtained by movement of a spool located between two opposing nozzles, thus differentially restricting the helium flow. The design procedure and test results of a high bandwidth, low power, linear electromagnetic actuator for the pool movement are presented. Two mathematical models are used to describe the helium flow through the thruster under different flow conditions. The experimental data measured in the vacuum chamber and the computational results from the mathematical models are compared. The models agree with the experiment to 3 percent except near a Knudsen number of 1.

Chen, J.-H.↗

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗

Numerical, micro-mechanical prediction of crack growth resistance in a fibre-reinforced/brittle matrix composite

Micromechanics fracture models are incorporated into three distinct fracture process zones which contribute to the crack growth resistance of fibrous composites. The frontal process zone includes microcracking, fiber debonding, and some fiber failure. The elastic process zone is related only to the linear elastic creation of new matrix and fiber fracture surfaces. The wake process zone includes fiber bridging, fiber pullout, and fiber breakage. The R-curve predictions of the model compare well with empirical results for a unidirectional, continuous fiber C/C composite. Separating the contributions of each process zone reveals the wake region to contain the dominant crack growth resistance mechanisms. Fractography showed the effects of the micromechanisms on the macroscopic fracture behavior.

Jenkins, Michael G.↗

Occultation of a compact radio source by Venus

The position of Venus in the radio reference frame was measured on the basis of an observation of an occultation of the compact radio source P0507 + 17 by Venus on July 19, 1988. By comparing model light curves to the observed light curves for both ingress and egress, the frame-tie offsets Delta(lambda) and Delta(beta) (ecliptic longitude and latitude) were determined. The solution for both Delta(lambda) and Delta(beta) gave results consistent with zero offsets, but with an error of about 0.2 arcsec. The egress light curve alone yielded higher accuracy for a linear combination of Delta(lambda) and Delta(beta): Delta(lambda) + 0.51 Delta(beta) = 0.026 +/-0.04 arcsec.

Linfield, R. P.↗

Matrix fatigue crack development in a notched continuous fiber SCS-6/Ti-15-3 composite

In this study the extensive matrix fatigue cracking that has been observed in notched SCS-6/Ti-15-3 composites is investigated. Away from the notch uniform spacing of the fatigue cracks develops. Closer to the notch, fiber-matrix debonding which occurs increases the crack spacing. Crack spacing and debond length determined from shear-lag cylinder models compare favorably with experimental observations. Scanning electron microscope (SEM) fractography showed that the principal fatigue crack initiation occurred around the zero degree fibers. Interface failure in the 90 degree plies does not lead to the development of the primary fatigue cracking.

Hillberry, B. M.↗

A comparison of the predictive capabilities of several turbulence models using upwind and central-difference computer codes

Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.

Rumsey, Christopher L.↗

Gravity and topography

The paper summarizes the fundamental gravity field constants for Mars and a brief historical review of early determinations and current-day accurate estimates. These include the planetary gravitational constant, global figure, dynamical oblateness, mean density, and rotational period. Topographic results from data acquired from the 1967 opposition to the most recent, 1988, opposition are presented. Both global and selected local topographic variations and features are discussed. The inertia tensor and the nonhydrostatic component of Mars are examined in detail. The dimensionless moment of inertia about the rotational axis is 0.4 for a body of uniform density and 0.37621 if Mars were in hydrostatic equilibrium. By comparing models of both gravity and topography, inferences are made about the degree and depth of compensation in the interior and stresses in the lithosphere.

Esposito, P. B.↗

Large eddy simulation of shock turbulence interaction

A nonconservative formulation of the energy equation (solving for internal energy) was used to perform large eddy simulations of compressible turbulence in Moin et al. due to its simplicity in implementing SGS models compared to the conservative formulation (solving for total energy). In problems with shocks in the domain, however, the total energy formulation is preferred due to its conservative nature. A conservative set of equations for the LES were derived from the nonconservative equations derived by Moin et al. Performance of the conservative formulation was compared with the experiment on decaying grid-generated turbulence as well as with the filtered DNS field. Various shock-capturing schemes were tested, and an ENO shock-capturing scheme of Shu and Osher was chosen for the simulation of shock/turbulence interaction. The scheme was tested and validated against the data base generated by DNS of weak shock waves. The results obtained with the essentially nonoscillatory (ENO) scheme were within 5 percent from the DNS results, and used less than 25 percent of the CPU time used in the DNS.

Lee, Sangsan↗

A simple analytical aerodynamic model of Langley Winged-Cone Aerospace Plane concept

A simple three DOF analytical aerodynamic model of the Langley Winged-Coned Aerospace Plane concept is presented in a form suitable for simulation, trajectory optimization, and guidance and control studies. The analytical model is especially suitable for methods based on variational calculus. Analytical expressions are presented for lift, drag, and pitching moment coefficients from subsonic to hypersonic Mach numbers and angles of attack up to +/- 20 deg. This analytical model has break points at Mach numbers of 1.0, 1.4, 4.0, and 6.0. Across these Mach number break points, the lift, drag, and pitching moment coefficients are made continuous but their derivatives are not. There are no break points in angle of attack. The effect of control surface deflection is not considered. The present analytical model compares well with the APAS calculations and wind tunnel test data for most angles of attack and Mach numbers.

Pamadi, Bandu N.↗

The NASA/MSFC Global Reference Atmospheric Model-1995 version (GRAM-95)

The latest version of the Global Reference Atmospheric Model (GRAM-95) is presented and discussed. GRAM-95 uses the new Global Upper Air Climatic Atlas (GUACA) CD-ROM data set, for 0- to 27-km altitudes. As with earlier versions, GRAM-95 provides complete geographical and altitude coverage for each month of the year. Individual years 1985 to 1991 and a period-of-record (1980 to 1991) can be simulated for the GUACA height range. GRAM-95 uses a specially developed data set, based on Middle Atmosphere Program (MAP) data, for the 20- to 120-km height range, and the NASA Marshall Engineering Thermosphere (MET) model for heights above 90 km. Fairing techniques assure a smooth transition in the overlap height ranges (20 to 27 km and 90 to 120 km). In addition to the traditional GRAM variables of pressure, density, temperature and wind components, GRAM-95 now includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He, and H). A new, variable-scale perturbation model provides both large-scale and small-scale deviations from mean values for the thermodynamic variables and horizontal and vertical wind components. The perturbation model includes new features that simulate intermittency (patchiness) in turbulence and small-scale perturbation fields. The density perturbations and density gradients (density shears) computed by the new model compare favorably in their statistical characteristics with observed density perturbations and density shears from 32 space shuttle reentry profiles. GRAM-95 provides considerable improvement in wind estimates from the new GUACA data set, compared to winds calculated from the geostrophic wind relations previously used in the 0- to 25-km height range. The GRAM-95 code has been put into a more modular form, easier to incorporate as subroutines in other programs (e.g., trajectory codes). A complete user's guide for running the program, plus sample input and output, is provided.

Justus, C. G.↗

Results of the Quasi-Steady Acceleration Environment from the STS-62 Missions

One of the clear benefits of conducting scientific research in space is to take advantage of the reduced acceleration environment. Many accelerometer packages have proven to accurately measure the acceleration environment at frequency levels above one Hz. However, for particular classes of experiments the quality of science returns is a direct function of the extremely low frequency (less than 0.01 Hz), quasi-steady acceleration environment. One class particularly interested in this acceleration regime is the group of crystal growth experimenters. These scientists are primarily interested in knowing the resultant quasi-steady acceleration vector at their respective crystal growth locations. The objective of many of these scientists is to minimize the amount of convective flow acting in a direction perpendicular to the growth axis of the crystal. Convective flow within the crystal can be induced by the direction and magnitude of the quasi-steady acceleration vector. Convective flows acting perpendicular to the growth axis of the crystal can cause nonuniformity within the crystal, thus reducing the quality of the results. The Orbital Acceleration Research Experiment (OARE), an accelerometer package hardmounted to the bottom of the payload bay of the orbiter Columbia (OV-102), has the capability of monitoring and recording the quasi-steady acceleration environment. This paper will describe the components that make up the on-orbit quasi-steady acceleration environment, detail how results from the OARE device were achieved, and compare modelled acceleration results with actual on-orbit OARE results from the STS-62 and STS-65 flights. A summary of the results will be provided along with possible recommendations of how to combine modelled and realtime quasi-steady accelerometer data for future Shuttle flights.

Matisak, Brian↗

Infalling Envelopes and Pre-Main Sequence Disks

The goal of this project is to understand the observed infrared emission of young stellar objects, and explore the implications of this emission for the evolution of dusty envelopes and circumstellar disks. We are using sophisticated radiative transfer methods to compare models with observations, thereby making critical tests of the standard picture of low-mass star formation.

Hartmann, Lee W.↗

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↗

Bursting Bubbles from Combustion of Thermoplastic Materials in Microgravity

Many thermoplastic materials in common use for a wide range of applications, including spacecraft, develop bubbles internally as they burn due to chemical reactions taking place within the bulk. These bubbles grow and migrate until they burst at the surface, forceably ejecting volatile gases and, occasionally, molten fuel. In experiments in normal gravity, Kashiwagi and Ohlemiller observed vapor jets extending a few centimeters from the surface of a radiatively heated polymethylmethacrylate (PMMA) sample, with some molten material ejected into the gas phase. These physical phenomena complicated the combustion process considerably. In addition to the non-steady release of volatiles, the depth of the surface layer affected by oxygen was increased, attributed to the roughening of the surface by bursting events. The ejection of burning droplets in random directions presents a potential fire hazard unique to microgravity. In microgravity combustion experiments on nylon Velcro fasteners and on polyethylene wire insulation, the presence of bursting fuel vapor bubbles was associated with the ejection of small particles of molten fuel as well as pulsations of the flame. For the nylon fasteners, particle velocities were higher than 30 cm/sec. The droplets burned robustly until all fuel was consumed, demonstrating the potential for the spread of fire in random directions over an extended distance. The sequence of events for a bursting bubble has been photographed by Newitt et al.. As the bubble reaches the fluid surface, the outer surface forms a dome while the internal bubble pressure maintains a depression at the inner interface. Liquid drains from the dome until it breaks into a cloud of droplets on the order of a few microns in size. The bubble gases are released rapidly, generating vortices in the quiescent surroundings and transporting the tiny droplets. The depression left by the escaping gases collapses into a central jet, which rises with a high velocity and may break up, releasing one or more relatively large drops (on the order of a millimeter in these experiments). A better understanding of bubble development and bursting processes, the effects of bursting behavior on burning rate of the bulk material, and the circumstances under which large droplets are expelled, as well as their trajectories, sizes, and burning rates, is sought through computer modeling compared with experiment.

Butler, K. B.↗

Bubble Generation in a Flowing Liquid Medium and Resulting Two-Phase Flow in Microgravity

The present investigation reports a study of bubble generation under reduced gravity conditions, using both a co-flow and a cross-flow configuration. This study may be used in the conceptual design of a space-based thermal management system. Ensuing two-phase flow void fraction can be accurately monitored using a single nozzle gas injection system within a continuous liquid flow conduit, as utilized in the present investigation. Accurate monitoring of void fraction leads to precise control of heat and mass transfer coefficients related to a thermal management system; hence providing an efficient and highly effective means of removing heat aboard spacecraft or space stations. Our experiments are performed in parabolic flight aboard the modified DC-9 Reduced Gravity Research Aircraft at NASA Lewis Research Center, using an air-water system. For the purpose of bubble dispersion in a flowing liquid, we use both a co-flow and a cross-flow configuration. In the co-flow geometry, air is introduced through a nozzle in the same direction with the liquid flow. On the other hand, in the cross-flow configuration, air is injected perpendicular to the direction of water flow, via a nozzle protruding inside the two-phase flow conduit. Three different flow conduit (pipe) diameters are used, namely, 1.27 cm, 1.9 cm and 2.54 cm. Two different ratios of nozzle to pipe diameter (D(sub N))sup * are considered, namely (D(sub N))sup * = 0.1 and 0.2, while superficial liquid velocities are varied from 8 to 70 cm/s depending on flow conduit diameter. It is experimentally observed that by holding all other flow conditions and geometry constant, generated bubbles decrease in size with increase in superficial liquid velocity. Detached bubble diameter is shown to increase with air injection nozzle diameter. Likewise, generated bubbles grow in size with increasing pipe diameter. Along the same lines, it is shown that bubble frequency of formation increases and hence the time to detachment of a forming bubble decreases, as the superficial liquid velocity is in-creased. Furthermore, it is shown that the void fraction of the resulting two-phase flow increases with volumetric gas flow rate Q(sub d), pipe diameter and gas injection nozzle diameter, while they decrease with surrounding liquid flow. The important role played by flowing liquid in detaching bubbles in a reduced gravity environment is thus emphasized. We observe that the void fraction can be accurately controlled by using single nozzle gas injection, rather than by employing multiple port injection, since the later system gives rise to unpredictable coalescence of adjacent bubbles. It is of interest to note that empirical bubble size and corresponding void fraction are somewhat smaller for the co-flow geometry than the cross-flow configuration at similar flow conditions with similar pipe and nozzle diameters. In order to supplement the empirical data, a theoretical model is employed to study single bubble generation in the dynamic (Q(sub d) = 1 - 1000 cu cm/s) and bubbly flow regime within the framework of the co-flow configuration. This theoretical model is based on an overall force balance acting on the bubble during the two stages of generation, namely the expansion and the detachment stage. Two sets of forces, one aiding and the other inhibiting bubble detachment are identified. Under conditions of reduced gravity, gas momentum flux enhances, while the surface tension force at the air injection nozzle tip inhibits bubble detachment. In parallel, liquid drag and inertia can act as both attaching and detaching forces, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with our experimental results. However, at higher superficial liquid velocities, as the bubble loses its spherical form, empirical bubble size no longer matches the theoretical predictions. In summary, we have developed a combined experimental and theoretical work, which describes the complex process of bubble generation and resulting two-phase flow in a microgravity environment. Results of the present study can be used in a wide range of space-based applications, such as thermal energy and power generation, propulsion, cryogenic storage and long duration life support systems, necessary for programs such as NASA's Human Exploration for the Development of Space (HEDS).

Pais, S. C.↗