Search NASA⌕ Search

SEARCH · Search NASA

Results for “drag force”

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 145 records · Page 8

Assessment of Particle Drag Models Using L2K and GBK Experimental Data

Properly assessing the drag force is an important element to computing the trajectories of particles traveling through a shock layer. This study will use experimental data in the German Aerospace Center (DLR) L2K arcjet to evaluate the ability of various particle drag models to reproduce the experimental particle velocity data.

particle flows↗

Absolute plate motions by boundary velocity minimizations

The main interaction of the earth's interior with the lithosphere is as a material source and sink. An absolute reference frame defined by minimizing the translational motion of tectonic plate boundaries differs by 0.6 cm/year from a frame defined by hot spot traces and by 0.4 cm/year from the frame defined by the most plausible model of drag forces on the plates. The rms absolute translational velocities are about 2 cm/year for ocean-ocean plate boundaries and 1.5 cm/year for ocean-continent plate boundaries. The close agreement between the source and sink and the drag-dependent definitions suggests that the lithosphere, as a stress guide, to some extent controls the locations of its sources and sinks.

Kaula, W. M.↗

Attitude Control System Design for CubeSats Configured with Exo-Brake Parachute

This paper develops a novel attitude control strategy for an Earth orbiting CubeSat spacecraft by utilizing the exo-brake parachute to modulate the atmospheric drag forces as a source of attitude control authority, enabling orbital exo-sail maneuvers. In particular, the spacecraft attitude controls can be realized through the two dimensional exo-sail maneuvers in pitch and yaw directions. The uncertain atmospheric drag induced disturbance torque is estimated through an adaptive parameter estimation process which makes use of the adaptive least-squares minimization techniques. The covariance updating law with a variable forgetting factor is adopted and it can be shown that the convergent rate for the estimation errors can be chosen at the same level as the forgetting factor, in order to meet the design needs. The proposed approach is best suited for Earth orbiting micronano-satellite applications, which are configured with exo-brake parachute. With integration of exo-sail actuation mechanism and disturbance estimation, we demonstrate through simulations that exo-sail induced control torque for CubeSat attitude maneuver is feasible.

Disturbance estimation↗

Colorimetric qualification of shear sensitive liquid crystal coatings

The work that has been done to date on the Shear Sensitive Liquid Crystal Project demonstrated that cholesteric liquid crystal coatings respond to both the direction and magnitude of a shearing force. The response of the coating is to selectively scatter incident white light into a spectrum of colors. Discernible color changes at a fixed angle of observation and illumination are the result of an applied shear stress. The intention was to be able to convert these observable color patterns from a flow visualization technique into a quantitative tool. One of the earlier intentions was to be able to use liquid crystals in dynamic flow fields. This was assumed possible because liquid crystals had made it possible to visualize transients in surface shear forces. Although the transients were visualized by color changes to an order one micro second, the time response of a coating to align to a shearing force is dependent on the magnitude of the change between its initial and final states. Unfortunately, the response is not instantaneous. It is for this reason any future attempt at quantifying the magnitude and directions of a shearing force are limited to surface shear stress vector fields in three dimensional steady state flows. This limitation does not significantly detract from the utility of liquid crystal coatings. The measurement of skin friction in the study of transition on wings, prediction of drag forces, performance assessment, and the investigation of boundary layer behavior is of great importance in aerodynamics. There exist numerous examples of techniques for the measurement of surface shear stress. Most techniques require arduous calibrations and necessitate extensive preparation of the receiving surfaces. However, the main draw back of instruments such as Preston tubes, hot films, buried wire gages, and floating element balances is that they only provide a point measurement. The advantages of capturing global shear data would be appreciable when compared with conventional point measurement sensors. It has yet to be determined if a repeatable correlation exists between the measured color of a liquid crystal coating and the magnitude/directional components of a shear vector imposed onto it.

Muratore, Joseph J., Jr.↗

Use of potential flow theory to evaluate subsonic inlet data from a simulator-powered nacelle at cruise conditions

Incompressible potential flow theory corrected for compressibility effects, using the Lieblein-Stockman compressibility correction, was used to predict surface and flow field static pressures for a subsonic inlet at cruise conditions. The calculated internal and external surface static pressures were in good agreement with data at most conditions. The analysis was used to determine the capture stream-tube location and static-pressure distribution. Additive drag coefficients obtained from these results were consistently higher than those obtained using one-dimensional compressible flow theory. Increasing the distance between the inlet and boattail increased the cowl drag force. The effect of the boundary layer on internal and external surface static-pressure distributions was small at the design cruise condition. The analytical results may be used as an aid to data reduction and for predicting inlet mass flow, stagnation point location, and inlet additive drag.

Bober, L. J.↗

Modelling the Huge, Herschel-Resolved Debris Ring around HD 207129

Debris disks, which are inferred from the observed infrared excess to be ensembles of dust, rocks. and probably planetesimals, arc common features of stellar systems. As the mechanisms of their fonnation and evolution are linked to those of planetary bodies, they provide valuable infonnation. The few well-resolved debris disks are even more valuable because they can serve as modelling benchmarks and help resolve degeneracies in modelling aspects such as typical grain sizes and distances. Here, we present an analysis of the HD 207129 debris disk, based on its well-covered spectral energy distribution and Herschel/PACS images obtained in the framework of the DUNES (DUst around NEarby Stars) programme. We use an empirical power-law approach to the distribution of dust and we then model the production and removal of dust by means of collisions, direct radiation pressure, and drag forces. The resulting best-fit model contains a total of nearly 10(sup -2) Earth masses in dust, with typical grain sizes in the planetesimal beh ranging from 4 to 7 micrometers. We constrain the dynamical excitation to be low, which results in very long collisional lifetimes and a drag that notably fills the inner gap, especially at 70 micrometers. The radial distribution stretches from well within 100 AU in an unusual, outward-riSing slope towards a rather sharp outer edge at about 170-190 AU. The inner edge is therefore smoother than that reported for Fomalhaut, but the contribution from the extended halo of barely bound grains is similarly small. Both slowly self-stirring and planetary perturbations could potentially have formed and shaped this disk.

Lohne, T.↗

X-57 Ground Dynamics Modeling and Analysis

This paper presents simulation work to assess the ground handling behavior of the X-57 airplane, prior to its first flight under electric propulsion. A tire deformation model was fit to manufacturer ground testing data to develop a detailed simulation model of previously unmodeled rolling drag forces on the tires. A special modeling effort corrected for a difference between the runway surface in the manufacturer data (grass) and in intended X-57 surface operations (asphalt/concrete). Analysis based on variation in the test conditions of the manufacturer test data showed that uncertainty in the new estimated rolling drag model should lead to no more than 1.5-percent uncertainty in the takeoff acceleration. Another batch simulation series swept across ground operating conditions and airplane configurations, characterizing turning radii and roll angles experienced for different open-loop steering commands in these various situations. This simulation allowed the identification of ground handling characteristics even beyond the regime of normal operations. As part of this task, a new, modified variant of the Fibonacci search algorithm identified the correct steering inputs for maintaining a straight path as a function of the operating condition. Finally, this paper describes a fixed-base piloted simulator assessment. This assessment supplements the previous simulated ground dynamics data with pilot-in-the-loop data and pilot feedback and confirms the overall acceptability of the ground handling qualities of the X-57 airplane prior to first flight.

Loren J. Newton↗

Navigation of the Space Shuttle

Navigational systems and operations for the Space Shuttle are described. All navigational instrumentation is controlled from within the pressurized main cabin. Measurements of the state vector and the attitude are made with an inertial measurement unit (IMU), which uses data initialized at the moment of take-off. Orbital location is calculated in approximations using the initial propulsion conditions, models of the gravity field, and aerodynamic drag forces. Updates are periodically received from ground tracking stations. IMU continues attitude information, and additional references are made with an automated startracker device. Information can also be gathered by optical alignment, and future systems will include radar tracking in an approach mode. Deorbit is accompanied by IMU altitude measurements as well as calculations of altitude based on drag measurements. Barometric measurements begin at about 80,000 ft altitude. Signals are received from TACAN beginning at 145,000 ft, and the microwave scanning beam landing system is started at 20,000 ft. Various radionavigation systems are also employed in all flight phases.

Edwards, A., Jr.↗

Experimental Investigation of a Point Design Optimized Arrow Wing HSCT Configuration

The M2.4-7A Arrow Wing HSCT configuration was optimized for straight and level cruise at a Mach number of 2.4 and a lift coefficient of 0.10. A quasi-Newton optimization scheme maximized the lift-to-drag ratio (by minimizing drag-to-lift) using Euler solutions from FL067 to estimate the lift and drag forces. A 1.675% wind-tunnel model of the Opt5 HSCT configuration was built to validate the design methodology. Experimental data gathered at the NASA Langley Unitary Plan Wind Tunnel (UPWT) section #2 facility verified CFL3D Euler and Navier-Stokes predictions of the Opt5 performance at the design point. In turn, CFL3D confirmed the improvement in the lift-to-drag ratio obtained during the optimization, thus validating the design procedure. A data base at off-design conditions was obtained during three wind-tunnel tests. The entry into NASA Langley UPWT section #2 obtained data at a free stream Mach number, M(sub infinity), of 2.55 as well as the design Mach number, M(sub infinity)=2.4. Data from a Mach number range of 1.8 to 2.4 was taken at UPWT section #1. Transonic and low supersonic Mach numbers, M(sub infinity)=0.6 to 1.2, was gathered at the NASA Langley 16 ft. Transonic Wind Tunnel (TWT). In addition to good agreement between CFD and experimental data, highlights from the wind-tunnel tests include a trip dot study suggesting a linear relationship between trip dot drag and Mach number, an aeroelastic study that measured the outboard wing deflection and twist, and a flap scheduling study that identifies the possibility of only one leading-edge and trailing-edge flap setting for transonic cruise and another for low supersonic acceleration.

Narducci, Robert P.↗

A Plasma Aerocapture and Entry System for Manned Missions and Planetary Deep Space Orbiters

A rapid mission to Mars requires a large change in vehicle velocity upon arrival to establish a stable orbit. This demand is even greater for a Neptune science, requiring many kilometers per second of V. It is clear from past mission studies that a manned Mars mission and deep space planetary orbiters require aerobraking and aerocapture which use aerodynamic drag forces to slow the spacecraft. Aerocapture would enable long term studies of the outer planets and moons that would not be possible with existing braking methodologies. While the ability to utilize these atmospheres to slow down and capture spacecraft would dramatically reduce the cost, launch mass, and travel time, currently planned approaches require significant additional spacecraft mass and risk as the spacecraft must descend deep into the planetary atmosphere in order to produce significant drag on a relatively small aeroshell. The plasma based Magnetoshell being developed in this program holds the potential to perform the desired braking with significantly increased drag and control while dramatically reducing mass. Most importantly, this technology significantly lowers the risk involved with aerocapture thereby making manned planetary missions possible. The fundamental physics of the Magnetoshell is based on demonstrated experimental results. Successful implementation will dramatically decrease radiation exposures, mission risk, launch cost, and launch mass. Implementation of aerobraking by employing a solid deflector or aeroshell as a method for orbit insertion and circularization has been successfully demonstrated in the past, resulting in launch mass savings greater than 50%. In order to reduce the effect of frictional heating and dynamic pressure on the typically fragile aeroshell, or worse solar panels, the braking must be distributed over many orbital passes at a high altitude in the less dense regions of the atmosphere. It can thus take several months for a meter-scale, 1000 kg, aeroshell to execute the many elliptic orbital passes through the atmosphere to achieve the required V. This rather slow method of braking not only reduces frictional heating and dynamic forces, but also avoids unpredictable dynamic behavior due to turbulence, as well as unknown and seasonally variable atmospheric composition and temperature which has led to dangerous, mission-critical events. For exploration class missions such as DRA 5.0 aerocapture and Thermal Protection Systems (TPS) are proposed for breaking at Mars for cargo missions. Traditional aerocapture is considered too risky for manned missions. Even with the enormous mass savings that aerocapture allows, it still requires 80 tons of aeroshell and significantly increased launch mass and propellant. As will be shown, by using Magnetoshells for aerobraking, the DRA 5.0 mission will save 224 metric tons (MT) and greater than $2 B in launch costs. Beyond the dramatic savings for existing mission architectures, a low-mass, risk-free aerocapture system would allow much more rapid missions to Mars and deep space orbiters by allowing direct, faster trajectories. As will be shown, the plasma Magnetoshell Aerobraking, Aerocapture, and Entry System (AAES) not only reduces mass and cost while enabling significant new mission architectures, but also significantly reduces radiation exposures by decreasing trip times.

Manned Missions↗

Momentum Flux Estimates for South Georgia Island Mountain Waves in the Stratosphere Observed via Satellite

We show high-resolution satellite observations of mountain wave events in the stratosphere above South Georgia Island in the remote southern Atlantic Ocean and compute the wave momentum fluxes for these events. The fluxes are large, and they imply important drag forces on the circulation. Small island orography is generally neglected in mountain wave parameterizations used in global climate models because limited model resolution treats the grid cell containing the island as ocean rather than land. Our results show that satellite observations can be used to quantitatively constrain mountain wave momentum fluxes, and they suggest that mountain waves from island topography may be an important missing source of drag on the atmospheric circulation.

SOUTH GEORGIA ISLANDS↗

Neutral motions in the polar thermosphere for northward Interplanetary Magnetic Field

During the rare occurrence of November 24, 1982, when the north-south B(z) component of the IMF became positive for a period of about 11 hours and reached a steady value of about 25 nT, the DE-2 spacecraft simultaneously measured the ionic and neutral species of the high latitude F-region along the track of the polar-orbiting satellite. Data from two northern winter polar passes of DE-2 illustrate the response of the neutral F-region to ion drag forcing that arises from a configuration of ion convection characteristics of strongly northward IMF. The multicellular ion drift pattern associated with positive B(z) is observed to drive a similar, but less structured and weaker, neutral wind configuration in the winter polar cap. While major features of the ion drift pattern are mimicked by the neutral circulation, smaller scale ion flow structures are not. These results demonstrate that the sunward flow of neutral gas can be established and maintained by ion drag in the central polar cap, for positive B(z).

Killeen, T. L.↗

Flight Dynamics of an Aeroshell Using an Attached Inflatable Aerodynamic Decelerator

An aeroelastic analysis of the behavior of an entry vehicle utilizing an attached inflatable aerodynamic decelerator during supersonic flight is presented. The analysis consists of a planar, four degree of freedom simulation. The aeroshell and the IAD are assumed to be separate, rigid bodies connected with a spring-damper at an interface point constraining the relative motion of the two bodies. Aerodynamic forces and moments are modeled using modified Newtonian aerodynamics. The analysis includes the contribution of static aerodynamic forces and moments as well as pitch damping. Two cases are considered in the analysis: constant velocity flight and planar free flight. For the constant velocity and free flight cases with neutral pitch damping, configurations with highly-stiff interfaces exhibit statically stable but dynamically unstable aeroshell angle of attack. Moderately stiff interfaces exhibit static and dynamic stability of aeroshell angle of attack due to damping induced by the pitch angle rate lag between the aeroshell and IAD. For the free-flight case, low values of both the interface stiffness and damping cause divergence of the aeroshell angle of attack due to the offset of the IAD drag force with respect to the aeroshell center of mass. The presence of dynamic aerodynamic moments was found to influence the stability characteristics of the vehicle. The effect of gravity on the aeroshell angle of attack stability characteristics was determined to be negligible for the cases investigated.

Cruz, Juan R.↗

Ion-neutral coupling in the high latitude thermosphere, part 2

On the 24th November, 1982, The North-South (Bz) component of the Interplanetary Magnetic Field (IMF) became positive for a period of about 11 hours reaching a relatively large and steady value of approximately 25 nT. During this rare occurrence, the Dynamics Explorer-2 (DE-2) spacecraft was in a configuration that enabled the dynamics of both ionic and neutral species of the high latitude F region to be measured simultaneously along the track of the polar orbiting satellite. Results from two Northern (winter) polar passes of DE-2, extracted from a larger data set, are shown to illustrate the response of the neutral F region to ion drag forcing arising from a configuration of ion convection characteristics of strongly northward IMF. The measured neutral winds differ appreciably from those more commonly observed for periods of southward IMF. The multi-cellular ion drift pattern associated with positive Bz is observed to drive a similar but less structured and weaker neutral wind configuration in the winter polar cap. Major features of the ion drift pattern are mimicked by the neutral circulation but smaller scale and more irregular sturctures of ion flow are not. This is ascribed to the relatively long time constant (few hours) for momentum exchange between the ion and neutral gases. The results demonstrate that sunward flow of neutral gas can be established and maintained by ion drag in the central polar cap for positive Bz.

Killeen, T. L.↗

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds

Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.

Farrugia, C. J.↗

Demonstration of an Aerocapture GN and C System Through Hardware-in-the-Loop Simulations

Aerocapture is an orbit insertion maneuver in which a spacecraft flies through a planetary atmosphere one time using drag force to decelerate and effect a hyperbolic to elliptical orbit change. Aerocapture employs a feedback Guidance, Navigation, and Control (GN&C) system to deliver the spacecraft into a precise postatmospheric orbit despite the uncertainties inherent in planetary atmosphere knowledge, entry targeting and aerodynamic predictions. Only small amounts of propellant are required for attitude control and orbit adjustments, thereby providing mass savings of hundreds to thousands of kilograms over conventional all-propulsive techniques. The Analytic Predictor Corrector (APC) guidance algorithm has been developed to steer the vehicle through the aerocapture maneuver using bank angle control. Through funding provided by NASA's In-Space Propulsion Technology Program, the operation of an aerocapture GN&C system has been demonstrated in high-fidelity simulations that include real-time hardware in the loop, thus increasing the Technology Readiness Level (TRL) of aerocapture GN&C. First, a non-real-time (NRT), 6-DOF trajectory simulation was developed for the aerocapture trajectory. The simulation included vehicle dynamics, gravity model, atmosphere model, aerodynamics model, inertial measurement unit (IMU) model, attitude control thruster torque models, and GN&C algorithms (including the APC aerocapture guidance). The simulation used the vehicle and mission parameters from the ST-9 mission. A 2000 case Monte Carlo simulation was performed and results show an aerocapture success rate of greater than 99.7%, greater than 95% of total delta-V required for orbit insertion is provided by aerodynamic drag, and post-aerocapture orbit plane wedge angle error is less than 0.5 deg (3-sigma). Then a real-time (RT), 6-DOF simulation for the aerocapture trajectory was developed which demonstrated the guidance software executing on a flight-like computer, interfacing with a simulated IMU and simulated thrusters, with vehicle dynamics provided by an external simulator. Five cases from the NRT simulations were run in the RT simulation environment. The results compare well to those of the NRT simulation thus verifying the RT simulation configuration. The results of the above described simulations show the aerocapture maneuver using the APC algorithm can be accomplished reliably and the algorithm is now at TRL-6. Flight validation is the next step for aerocapture technology development.

Masciarelli, James↗

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗