A fundamental investigation of the phenomena that characterize liquid-film cooling
Mass transfer net rate /evaporation plus entrainment/ from thin liquid film, using methanol, water, butanol and RP-1 as coolants
SEARCH · Engineering Papers
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.
Mass transfer net rate /evaporation plus entrainment/ from thin liquid film, using methanol, water, butanol and RP-1 as coolants
Real clouds and clouds systems are inherently three dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud system with large horizontal domains at the National Center for Atmospheric Research. The results indicate that surface precipitation and latent heating profiles are very similar between the 2D and 3D simulations of these same cases. The reason for the strong similarity between the 2D and 3D CRM simulations is that the observed large-scale advective tendencies of potential temperature, water vapor mixing ratio, and horizontal momentum were used as the main forcing in both the 2D and 3D models. Interestingly, the 2D and 3D versions of the CRM used in CSU and U.K. Met Office showed significant differences in the rainfall and cloud statistics for three ARM cases. The major objectives of this project are to calculate and axamine: (1)the surface energy and water budgets, (2) the precipitation processes in the convective and stratiform regions, (3) the cloud upward and downward mass fluxes in the convective and stratiform regions; (4) cloud characteristics such as size, updraft intensity and lifetime, and (5) the entrainment and detrainment rates associated with clouds and cloud systems that developed in TOGA COARE, GATE, SCSMEX, ARM and KWAJEX. Of special note is that the analyzed (model generated) data sets are all produced by the same current version of the GCE model, i.e. consistent model physics and configurations. Trajectory analyse and inert tracer calculation will be conducted to identify the differences and similarities in the organization of convection between simulated 2D and 3D cloud systems.
This study presents a detailed analysis of a sweeping jet (fluidic oscillator) actuator. The sweeping jet actuator promises to be a viable flow control actuator candidate due to its simple, no moving part structure and its high momentum, spatially oscillating flow output. Hot-wire anemometer and particle image velocimetry measurements were carried out with an emphasis on understanding the actuator flow field in a quiescent environment. The time averaged, fluctuating, and instantaneous velocity measurements are provided. A modified actuator concept that incorporates high-speed solenoid valves to control the frequency of oscillation enabled phase averaged measurements of the oscillating jet. These measurements reveal that in a given oscillation cycle, the oscillating jet spends more time on each of the Coanda surfaces. In addition, the modified actuator generates four different types of flow fields, namely: a non oscillating downward jet, a non oscillating upward jet, a non oscillating straight jet, and an oscillating jet. The switching from an upward jet to a downward jet is accomplished by providing a single pulse from the solenoid valve. Once the flow is switched, the flow stays there until another pulse is received. The oscillating jet is compared with a non oscillating straight jet, which is a typical planar turbulent jet. The results indicate that the oscillating jet has a higher (5 times) spreading rate, more flow entrainment, and higher velocity fluctuations (equal to the mean velocity).
The correct representation of lowlevel midlatitude clouds found in the wake of cold fronts strongly relies on the representation of planetary boundary layer (PBL) and convection processes, which are typically parameterized separately in numerical models. Using the Weather Research and Forecasting Model (WRF), this study investigates how distinct pairs of PBL and convection parameterization schemes represent cloud fraction in the postcoldfrontal region. The simulations focus on the region of the DOEARM Eastern North Atlantic observation site in the Azores Islands in the wake of a cold front that passed on December 25, 2015. Different PBL and convection schemes are combined to create 12 distinct configurations. The main differences between the selected physical parameterizations are the strength of vertical mixing and the entrainment. The simulations produce a wide range of cloud fractions, where some configurations significantly underestimate while others clearly overestimate satellite and surfacebased estimates of cloud fraction. A skill score is used to quantitatively assess the performance of each configuration with respect to groundbased radar data. The key processes that are found to significantly impact the cloud fraction distribution are the strength of the PBL decoupling, the vertical wind shear, entrainment and detrainment rates in shallow convection and the occurrence of drizzle. This indicates that to successfully simulate postcoldfrontal clouds, modeled physics must balance strong internal vertical mixing and weak exchange with the free troposphere. For this case study, cloud fraction was more sensitive to the choice of convection scheme than PBL scheme.
It is well known that buoyancy has a major influence on the flow structure of turbulent nonpremixed jet flames. Buoyancy acts by inducing baroclinic torques, which generate large-scale vortical structures that can significantly modify the flow field. Furthermore, some suggest that buoyancy can substantially influence the large-scale structure of even nominally momentum-dominated flames, since the low velocity flow outside of the flame will be more susceptible to buoyancy effects. Even subtle buoyancy effects may be important because changes in the large-scale structure affects the local entrainment and fluctuating strain rate, and hence the structure of the flame. Previous studies that have compared the structure of normal- and micro-gravity nonpremixed jet flames note that flames in microgravity are longer and wider than in normal-gravity. This trend was observed for jet flames ranging from laminar to turbulent regimes. Furthermore, imaging of the flames has shown possible evidence of helical instabilities and disturbances starting from the base of the flame in microgravity. In contrast, these characteristics were not observed in normal-gravity. The objective of the present study is to further advance our knowledge of the effects of weak levels of buoyancy on the structure of transitional and turbulent nonpremixed jet flames. In later studies we will utilize the drop tower facilities at NASA Glenn Research Center (GRC), but the preliminary work described in this paper was conducted using the 1.25-second drop tower located at the University of Texas at Austin. A more detailed description of these experiments can be found in Idicheria et al.
It has been proposed that a vestibular reflex originating in the otolith organs and other body graviceptors modulates sympathetic activity during changes in posture with regard to gravity. To test this hypothesis, we selectively stimulated otolith and body graviceptors sinusoidally along different head axes in the coronal plane with off-vertical axis rotation (OVAR) and recorded sympathetic efferent activity in the peroneal nerve (muscle sympathetic nerve activity, MSNA), blood pressure, heart rate, and respiratory rate. All parameters were entrained during OVAR at the frequency of rotation, with MSNA increasing in nose-up positions during forward linear acceleration and decreasing when nose-down. MSNA was correlated closely with blood pressure when subjects were within +/-90 degrees of nose-down positions with a delay of 1.4 s, the normal latency of baroreflex-driven changes in MSNA. Thus, in the nose-down position, MSNA was probably driven by baroreflex afferents. In contrast, when subjects were within +/-45 degrees of the nose-up position, i.e., when positive linear acceleration was maximal along the naso-ocipital axis, MSNA was closely related to gravitational acceleration at a latency of 0.4 s. This delay is too short for MSNA changes to be mediated by the baroreflex, but it is compatible with the delay of a response originating in the vestibular system. We postulate that a vestibulosympathetic reflex, probably originating mainly in the otolith organs, contributes to blood pressure maintenance during forward linear acceleration. Because of its short latency, this reflex may be one of the earliest mechanisms to sustain blood pressure upon standing.
The spreading rates of convective thermals are linked to their net entrainment, and previous literature has suggested differences in spreading rates between moist and dry thermals. Here, growth rates of idealized numerically simulated axisymmetric dry and moist convective thermals are directly compared. In an environment with dry-neutral stratification, the increase of thermal radius with height dR / dz is a larger by a factor of 1.7 for dry thermals relative to moist thermals. The fractional change in thermal volume ε is also greater for dry thermals within a distance of ~4 radii from the initial thermal height. Values of dR / dz are nearly constant with height for both moist and dry thermals consistent with classical theory based on dimensional analysis. The simulations are also consistent with theory relating impulse, circulation, and spreading rate for dry thermals proposed in previous papers and extended here to moist thermals, predicting they will spread less than dry thermals. Tests adding heating to dry thermals, either spread uniformly across the thermal volume or concentrated in the inner core, indicate that dR / dz and ε are smaller for moist thermals because latent heating is confined mostly to their cores. Additional axisymmetric moist simulations using modified lapse rates and large-eddy simulations support this analysis. Overall, these results indicate that slow spreading rates are a fundamental feature of moist thermals caused by latent heating that alters the spatial distribution of buoyancy within them relative to dry thermals.
Fluid flow inside a beam blank mold fed through a three-port SEN (two lateral ports and one bottom port), positioned at the center of the mold, has been investigated. Literature survey shows that this kind of configuration is not frequently used. It had been shown that it is possible to get a symmetrical flow with this configuration. The downward inclination of the ports should influence the transient fluid flow. The greater the downward inclination, the greater is the instability of the jet and of the interface between the immiscible fluids that simulate the slag/metal behavior. The flow characterization was made by dye dispersion, PIV technique, measurement of the meniscus oscillation using ultrasonic sensors, and CFD simulations. The slag/metal behavior was simulated using water and oils with different physical properties. The results from PIV as well as from observations of the water–oil interface have been used to validate the CFD simulations. Oils with density close to water resulted in more instability at the interface with entrainment starting from fluid flow rate of 125 L/min (equivalent to casting speed of 1 m/min). Decreasing the viscosity of slag (oils) and increasing the casting velocity (water flow rate) result in reduction of the interfacial stability.
A method is developed which uses numerical tracers to make accurate diagnoses of entraimnent and detrainment rates and of the properties of the entrained and detrained air in numerically simulated clouds. The numerical advection scheme is modified to make it nondispersive, as required by the use of the tracers. Tests of the new method are made, and an appropriate definition of clouds is selected. Distributions of mixing fractions in the model consistently show maximums at the end points, for nearly undilute environmental air or nearly undilute cloud air, with a uniform distribution between. The cumulonimbus clouds simulated here entrain air that had been substantially changed by the clouds, and detrained air that is not necessarily representative of the cloud air at the same level.
Gas flow rates in Directed Energy Deposition (DED) Additive Manufacturing (AM) can significantly affect the quality of built parts by altering melt pool geometry. Using a DED process replicator and in situ synchrotron radiography, together with analogous experiments in an industrial DED machine, we investigate the impact of carrier gas and shield gas flow rates on build quality. The results reveal that there is a critical shield gas flow rate above which melt pools are flattened, tracks widen, and thus layer thickness decreases. The reduction in layer thickness is most prominent in conditions with low carrier gas flow rate, as the highly turbulent shield gas flow may divert slow moving powder particles away from the melt pool, decreasing capture efficiency. Very high flow rates increase internal porosity, as fast-moving particles impacting the melt pool surface can entrain chamber gas behind them. High gas flow rates also cool the melt pool, creating shallower melt pools with increased thermal gradients near the solidification front, increasing pore entrapment in the solidified track.
In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.
The results of direct numerical simulations of chemically reacting turbulent mixing layers are presented. This is an extension of earlier work to a more detailed study of previous three dimensional simulations of cold reacting flows plus the development, validation, and use of codes to simulate chemically reacting shear layers with heat release. Additional analysis of earlier simulations showed good agreement with self similarity theory and laboratory data. Simulations with a two dimensional code including the effects of heat release showed that the rate of chemical product formation, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release. Subsequent three dimensional simulations showed similar behavior, in agreement with laboratory observations. Baroclinic torques and thermal expansion in the mixing layer were found to produce changes in the flame vortex structure that act to diffuse the pairing vortices, resulting in a net reduction in vorticity. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers were shown to result from vorticity generation by baroclinic torques.
The effect of low-speed mixed convection (forced plus buoyant) on concurrent flow flame spread over a thin solid is examined. Computations are carried out using an existing model. Results indicate that seemingly small levels of gravity can significantly alter flame spread rates. Starting with a purely forced flow condition, as gravity is added, the entrainment due to buoyancy lengthens the flame and increases the spread rate significantly. Buoyancy has an influence on the extinction limits. At low speed, the presence of a small gravity level widens the flammability limit.
Here, this work reports on the design and performance of a jet pump system to pump slurries. The jet pump uses a high-velocity water jet to entrain and pump the slurry. This pumping action requires no valves and does not use moving solid parts. Hence, the system can be readily sealed, and the pumping action creates no observable mechanical damage to slurry particles. We tested carbon black beads suspended in water at concentrations of 0.0, 0.2, and 2.0 wt% at primary jet flow rates of 50 and 100 mL/min. Our measurements show an approximately linear decrease in entrainment ratio (slurry-to-jet flow rate) and absolute slurry flow rate with increasing adverse (load) pressure difference. The entrainment ratio is a strong function of the absolute jet flow rate. Pump performance was insensitive to slurry composition. Micron-resolution particle image velocimetry (micro-PIV) revealed qualitative changes in the pump’s internal flow due to changes in absolute jet flow rate, including evidence of the suppression of a vortex at higher flow rates. These findings establish key operational parameters for the jet pump and highlight its potential in a variety of slurry flow pumping applications. The pump is particularly suited for flow electrode slurries (e.g. as used for energy storage and water treatment) as it maintains particle integrity and can be used with standard tubing, tube fittings, and sealing methods.
Increases in aerosol concentration are well known to influence the microphysical processes and radiative properties of clouds. By reducing droplet size, an increase in aerosol can lessen collision efficiency and increase liquid water path (LWP) in precipitating clouds or enhance evaporation rate and decrease LWP in non-precipitating clouds. We utilize large eddy simulations to further investigate these aerosol indirect effects in Arctic mixed-phase clouds and find, in agreement with previous studies, precipitating clouds to experience an increase in LWP and non-precipitating clouds a decrease in LWP. Most importantly however, our results reveal a different explanation for why such an LWP decrease occurs in decoupled, non-precipitating clouds. We find enhanced evaporation near cloud top to be driven primarily by a strengthening of maximum radiative cooling rate with aerosol concentration which drives stronger entrainment, an effect that holds true even in clouds that are optically thick.
Bubble cavities formed by air entrainment and attached to a rotating shaft in an oil reservoir were studied. The cavities appear to the unaided eye as toroidal. High speed photography, however, reveals the individuality of the bubble cavities and their near solid body rotational characteristics. The cavities are distorted by the rotation effects but remain attached and tend to merge because of edge effects in the axial direction. The flow field within the reservoir is influenced by the unusual character of the two phase fluid found there; the vorticity is readily visualized. Other examples of vapor entrapment at the inlet of an eccentric rotor are also discussed. A simplified analytical method is provided, and a numerical analysis is being investigated. Vapor (void) entrainment and generation can significantly alter leakage rates and stability of seals, bearings, and dampers. Recognition of these effects in the component design systems will result only after detailed studies of the above phenomena.
When a closed body or a duct envelope moves through the atmosphere, air pressure and temperature rises occur ahead of the body or, under ram conditions, within the duct. If cloud water droplets are encountered, droplet evaporation will result because of the air-temperature rise and the relative velocity between the droplet and stagnating air. It is shown that the solution of the steady-state psychrometric equation provides evaporation rates which are the maximum possible when droplets are entrained in air moving along stagnation lines under such conditions. Calculations are made for a wide variety of water droplet diameters, ambient conditions, and flight Mach numbers. Droplet diameter, body size, and Mach number effects are found to predominate, whereas wide variation in ambient conditions are of relatively small significance in the determination of evaporation rates. The results are essentially exact for the case of movement of droplets having diameters smaller than about 30 microns along relatively long ducts (length at least several feet) or toward large obstacles (wings), since disequilibrium effects are then of little significance. Mass losses in the case of movement within ducts will often be significant fractions (one-fifth to one-half) of original droplet masses, while very small droplets within ducts will often disappear even though the entraining air is not fully stagnated. Wing-approach evaporation losses will usually be of the order of several percent of original droplet masses. Two numerical examples are given of the determination of local evaporation rates and total mass losses in cases involving cloud droplets approaching circular cylinders along stagnation lines. The cylinders chosen were of 3.95-inch (10.0+ cm) diameter and 39.5-inch 100+ cm) diameter. The smaller is representative of icing-rate measurement cylinders, while with the larger will be associated an air-flow field similar to that ahead of an airfoil having a leading-edge radius comparable with that of the cylinder. It is found that the losses are less than 5 percent. It is concluded that such losses are, in general, very small (less than 1 percent) in the case of smaller obstacles (of icing-rate measurement- cylinder size); the motional dynamics are such, however, that exceptions will occur by reason of failure of very small droplets (moving along stagnation lines) to impinge upon obstacle surfaces. In such cases, the droplets will evaporate completely.
A unified theory of entrainment into the planetary boundary layer is presented. It is assumed that the rates of buoyant and shear production of turbulence kinetic energy can be determined in terms of the entrainment mass flux. An expression is derived from the conservation law for turbulence kinetic energy, which, with the introduction of an empirical parameter, can be used together with a second relation between turbulence kinetic energy and the turbulence velocity scale to obtain the mass entrainment flux. The theory provides descriptions of storage-limited entrainment, buoyancy-limited entrainment into a clear mixed layer, and shallowing. It has been incorporated into a simulation of Day 33 of the Wangara experiment using a simple mixed layer model.