LOW-DENSITY SPHERE DRAG WITH EQUILIBRIUM AND NONEQUILIBRIUM WALL TEMPERATURE
Rarefied gas dynamics - flow regime between continuum flow and free molecular flow
CONTINUUM FLOW↗
SEARCH · Engineering Papers
Results for “Flow dynamics”
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.
(abstract) Three Dimensional Ice-Flow Velocity Estimation Using Satellite Radar Interferometry
An understanding of the flow dynamics of an ice sheet's outlet glaciers and ice streams requires knowledge of their flow velocity and strain rates (i.e., velocity gradients). Prior to the recent advent of satellite radar interferometry, it was not possible to measure detailed ice-flow velocity over the vast featureless areas that comprise most of the ice sheets. Since the launch of ERS-1, the use of satellite radar interferometry data for making densly sampled ice-flow velocity measurements has been firmly established by several studies. We have combined data from nonparallel orbits with surface slope information to make vector ice-flow measurements for the Ryder Glacier, Greenland. Our results for the Ryder are promising and indicate that repeat-pass interferometric data can be used to make vector measurements of ice velocity.
ice sheets radar interferometry glaciers flow dyna↗
Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification
The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.
Gupta, A. K.↗
Embedded cavity drag in steady and unsteady flows
The numerical solution of the laminar boundary-layer flow over an embedded cavity is studied. The purpose is to examine the relevant drag characteristics of laminar cavity flow. The solution field is obtained in terms of velocity and vorticity variables, with the stream function and pressure derivable from the directly computed variables. An analysis and comparison is made among four square cavities, ranging in size from 0.25 to 1.00 boundary-layer thicknesses deep. The dominant flow features are examined in the vicinity of the cavity by means of the stream function and iso-vorticity contours. The dominant physics in the overall drag characteristics of the flow is examined by an analysis of the pressure and wall shear stress distributions in the cavity, and upstream and downstream of the cavity. Pressure forces and frictional forces in, and in the vicinity of, the cavity are determined. Stress relaxation distances, both upstream and downstream of the cavity, are calculated and analyzed. The flow dynamics of the boundary-layer flow over an embedded cavity is summarized. Finally, the relevance of the results to the control of flow separation in such flows is discussed.
Gatski, T. B.↗
Embedded cavity drag in steady and unsteady flows
The numerical solution of the laminar boundary-layer flow over an embedded cavity is studied. The purpose of the study is to examine the relevant drag characteristics of laminar cavity flow. The solution field is obtained in terms of velocity and vorticity variables, with the stream function and pressure derivable from the directly computed variables. An analysis and comparison is made among four square cavities, ranging in size from 0.25 to 1.00 boundary-layer thicknesses deep. The dominant flow features are examined in the vicinity of the cavity by means of the stream function and iso-vorticity contours. The dominant physics in the overall drag characteristics of the flow is examined by an analysis of the pressure and wall shear stress distributions in the cavity, and upstream and downstream of the cavity. Pressure forces and frictional forces in, and in the vicinity of, the cavity are determined. Stress relaxation distances, both upstream and downstream of the cavity, are calculated and analyzed. The flow dynamics of the boundary-layer flow over an embedded cavity is summarized. Finally, the relevance of the present results to the control of flow separation in such flows is discussed.
Gatski, T. B.↗
Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production
In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).
In situ Resource Utilization↗
A Study of Bubble and Slug Gas-Liquid Flow in a Microgravity Environment
The influence of gravity on the two-phase flow dynamics is obvious.As the gravity level is reduced,there is a new balance between inertial and interfacial forces, altering the behavior of the flow. In bubbly flow,the absence of drift velocity leads to spherical-shaped bubbles with a rectilinear trajectory.Slug flow is a succession of long bubbles and liquid slug carrying a few bubbles. There is no flow reversal in the thin liquid film as the long bubble and liquid slug pass over the film. Although the flow structure seems to be simpler than in normal gravity conditions,the models developed for the prediction of flow behavior in normal gravity and extended to reduced gravity flow are unable to predict the flow behavior correctly.An additional benefit of conducting studies in microgravity flows is that these studies aide the development of understanding for normal gravity flow behavior by removing the effects of buoyancy on the shape of the interface and density driven shear flows between the gas and the liquid phases. The proposal calls to study specifically the following: 1) The dynamics of isolated bubbles in microgravity liquid flows will be analyzed: Both the dynamics of spherical isolated bubbles and their dispersion by turbulence, their interaction with the pipe wall,the behavior of the bubbles in accelerated or decelerated flows,and the dynamics of isolated cylindrical bubbles, their deformation in accelerated/decelerated flows (in converging or diverging channels), and bubble/bubble interaction. Experiments will consist of the use of Particle Image Velocimetry (PIV) and Laser Doppler Velocimeters (LDV) to study single spherical bubble and single and two cylindrical bubble behavior with respect to their influence on the turbulence of the surrounding liquid and on the wall 2) The dynamics of bubbly and slug flow in microgravity will be analyzed especially for the role of the coalescence in the transition from bubbly to slug flow (effect of fluid properties and surfactant), to identify clusters that promote coalescence and transition the void fraction distribution in bubbly and slug flow,to measure the wall friction in bubbly flow. These experiments will consist of multiple bubbles type flows and will utilize hot wire and film anemometers to measure liquid velocity and wall shear stress respectively and double fiber optic probes to measure bubble size and velocity as a function of tube radius and axial location.
McQuillen, J.↗
Prediction of prompt NO(x) in hydrocarbon air flames
The gas turbine industry is directing particular attention to very low NOx combustors, whether for aircraft or land based CCGT systems. These low NOx combustors frequently use liquid fuels or natural gas burning under very lean premixed conditions with air or under rich-lean conditions, although only the first case is studied here. In land based systems, diluted steam or nitrogen are sometimes injected into the combustion chambers to reduce flame temperature. The NOx emissions from such systems are the product of three chemical mechanisms which are interrelated: the hydrocarbon prompt NO, the thermal NO (extended Zeldovich mechanism), and the nitrous oxide route to NO. Formation of NO2 from NO also occurs, as well as emission of carbon monoxide and unburnt hydrocarbons. When the fuel-oxidant proportion decreases towards leaner conditions, flame temperatures are lowered, resulting in the total NO being reduced and the thermal-NO contribution greatly diminished to the benefit of the remaining two mechanisms of NO formation. While knowledge of the elementary reactions and their chemical kinetics concerning methane and simple hydrocarbons combustion has existed for a number of years, its use for computer modeling is limited to simple flow dynamics configurations. Nevertheless, understanding of such combusting flows under a wide range of experimental conditions allows for analogies or speculations with more complex actual systems. Such understanding can be achieved by means of one dimensional laminar premixed flame modeling, with a full chemical mechanism which incorporates the three routes of NO formation. Complementary to this understanding is the modeling of the actual combustion system using a full description of the fluid dynamics coupled with a reduced chemical scheme, which is then compared against the first model. The objective of this investigation is to evaluate the relative importance of the three mechanisms of NO formation in lean premixed methane-air combustion with increasing pressure using the one dimensional plug flow package PREMIX, and to test the validity of a three dimensional model with a global chemical mechanism against the one dimensional model in the atmospheric pressure case. Methane is chosen because it is the only mechanism which is reasonably well known and is a good guide to the behavior of other hydrocarbons. The mixture ratio chosen is richer than that in lean gas turbines, but the combustion of this mixture with the low preheat gives realistic gas turbine final flame temperatures. Conditions of NO2 formation are also analyzed in the one dimensional model and results extrapolated to the case of gas turbines.
Valerie Dupont↗
Synchronization trigger control system for flow visualization
The use of cinematography or holographic interferometry for dynamic flow visualization in an internal combustion engine requires a control device that globally synchronizes camera and light source timing at a predefined shaft encoder angle. The device is capable of 0.35 deg resolution for rotational speeds of up to 73 240 rpm. This was achieved by implementing the shaft encoder signal addressed look-up table (LUT) and appropriate latches. The developed digital signal processing technique achieves 25 nsec of high speed triggering angle detection by using direct parallel bit comparison of the shaft encoder digital code with a simulated angle reference code, instead of using angle value comparison which involves more complicated computation steps. In order to establish synchronization to an AC reference signal whose magnitude is variant with the rotating speed, a dynamic peak followup synchronization technique has been devised. This method scrutinizes the reference signal and provides the right timing within 40 nsec. Two application examples are described.
Chun, K. S.↗
Aspirator increases relief valve poppet stroke
Addition of an aspirator to a relief valve increases the valve poppet stroke under dynamic flow conditions. The aspirator allows poppet inlet dynamic forces to overcome relief valve spring force. It reduces the fluid pressure in the skirt cavity by providing a low pressure sense probe.
Biddle, M. E.↗
Generation of Martian chaos and channels by debris flows
A debris flow mechanism is proposed to account for the formation of chaos and the large channels debouching into Crysae Planitia from the adjacent southern uplands of Mars. Based on considerations of the juxtaposition of individual channel environments, the morphological assemblages within each environment and flow dynamics, it is suggested that the debris flows were triggered by the large-scale failure of subsurface sediments, possibly initiated by a seismic event. During the initial, slow-moving phase of the flow, the debris would have formed gently sinuous channels with multiple side-wall slumps, grooves and ridges, and elongate erosional remnants. The flow would have gained mobility as the debris moved downslope, producing travel distances greatly in excess of those characteristic of terrestrial examples, and eroded, streamlined remnants at the distal reaches of the channel. Finally, due to internal and boundary friction, the flow would have been slowed down once it entered the Chryse plains, resulting in a thin debris blanket with no depositional relief.
Nummedal, D.↗
Visualization of a Sweeping Jet by Laser Speckle Retro-reflective Background Oriented Schlieren
Sweeping jet actuators are currently being studied as a method of active flow control to increase the performance of wings and tails. Computational Fluid Dynamics simulations show that the jet does not sweep at a constant rate, instead spending more time at the minimum and maximum of its extent. As part of a test done to provide dynamic flow measurements for Computation Fluid Dynamics validation in NASA Ames Fluid Mechanics Lab, Particle Image Velocimetry (PIV) was performed on a single sweeping jet actuator designed by NASA and Boeing that was inset in a small airfoil. However, even at low jet pressures, the PIV could not be acquired at high enough frequencies to both spatially and temporally resolve the stickiness of the jet motion. Therefore, a small portion of the PIV laser was picked off and used to acquire flow visualization at higher frequencies through a new technique, laser speckle Retroreflective Background Oriented Schlieren (RBOS).
sweeping jet actuator↗
A subsonic transonic and supersonic nozzle flow by the inverse technique
Inverse solution of two dimensional gas dynamic flow fields of rotational or irrotational character
Norton, David J.↗
Local dynamic subgrid-scale models in channel flow
The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.
Cabot, William H.↗
Transport processes in boiling and two-phase systems, including near-critical fluids
Aspects of pool boiling are considered, taking into account nucleate boiling, the nucleate boiling mechanism, film boiling, and the transition between nucleate and film boiling. The characteristics of two-phase flow are also investigated, giving attention to two-phase flow parameters and equations, the flow pattern in two-phase flow, the pressure drop in two-phase flow, heat transfer in two-phase flow, two-phase flow dynamics, the boiling crisis in two-phase flow, the critical flow rate, the propagation of the pressure pulse and the sonic velocity in two-phase media, instrumentation for two-phase flow, and geometry and field effects on boiling and two-phase flow. Near-critical fluids are also considered.
Hsu, Y.-Y.↗
Wind and Rotation Enhanced Escape from the Early Terrestrial Atmospheres
The earliest atmospheres of the terrestrial planets are thought to have been hotter, have stronger winds and rotate faster than atmospheres of today. Since these primitive atmospheres were weakly bound, they evolved rapidly because atmospheric escape was very strong, often referred to as "blowoff." Such escape has been treated as hydrodynamic, transonic flow; similar to solar wind flow dynamics. However, in many cases, although the outward flow is hydrodynamic at low altitudes, it becomes collisionless at higher altitudes, before sonic speeds are ever attained. Recent models dealing with the transition from fluid to kinetic flow have applied the Jeans escape flux at the exobase. This approach leads to escape rates that are too low, because thermospheric winds and planetary rotation are known to increase the escape flux above the corresponding Jeans flux. In particular, for a given density and temperature at the exobase, the escape flux increases as the wind speed and/or the rotation rate increase. Also, for a given wind speed and rotation rate, the escape flux enhancement over the Jeans flux increases as the mass of an escaping constituent increases, an important factor in isotope fractionation, especially the enrichment of deuterium on Mars. Accounting for a range of possible temperatures, thermospheric wind speeds and planetary rotation rates in the primitive atmospheres of the terrestrial planets, estimates are made of light constituent escape flux increases over the corresponding Jeans fluxes.
Hartle, Richard E.↗
CFD Simulations of the Supersonic Inflatable Aerodynamic Decelerator (SIAD) Ballistic Range Tests
A series of ballistic range tests were performed on a scaled model of the Supersonic Flight Demonstration Test (SFDT) intended to test the Supersonic Inflatable Aerodynamic Decelerator (SIAD) geometry. The purpose of these experiments were to provide aerodynamic coefficients of the vehicle to aid in mission and vehicle design. The experimental data spans the moderate Mach number range, $3.8-2.0$, with a total angle of attack ($alpha_T$) range, $10o-20o$. These conditions are intended to span the Mach-$alpha$ space for the majority of the SFDT experiment. In an effort to validate the predictive capabilities of Computational Fluid Dynamics (CFD) for free-flight aerodynamic behavior, numerical simulations of the ballistic range experiment are performed using the unstructured finite volume Navier-Stokes solver, US3D. Comparisons to raw vehicle attitude, and post-processed aerodynamic coefficients are made between simulated results and experimental data. The resulting comparisons for both raw model attitude and derived aerodynamic coefficients show good agreement with experimental results. Additionally, near body pressure field values for each trajectory simulated are investigated. Extracted surface and wake pressure data gives further insights into dynamic flow coupling leading to a potential mechanism for dynamic instability.
CFD modeling↗
Wind-US Unstructured Flow Solutions for a Transonic Diffuser
The Wind-US Computational Fluid Dynamics flow solver computed flow solutions for a transonic diffusing duct. The calculations used an unstructured (hexahedral) grid. The Spalart-Allmaras turbulence model was used. Static pressures along the upper and lower wall agreed well with experiment, as did velocity profiles. The effect of the smoothing input parameters on convergence and solution accuracy was investigated. The meaning and proper use of these parameters are discussed for the benefit of Wind-US users. Finally, the unstructured solver is compared to the structured solver in terms of run times and solution accuracy.
Mohler, Stanley R., Jr.↗