Entrainment due to a thermal impinging on a stratified interface with and without buoyancy reversal
|
|
- Polly Gaines
- 5 years ago
- Views:
Transcription
1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 105, NO. D12, PAGES 15,457-15,467, JUNE 27, 2000 Entrainment due to a thermal impinging on a stratified interface with and without buoyancy reversal Qing Zhang Department of Mechanical and Industrial Engineering, University of Manitoba, Winnipeg, Manitoba, Canada Aline Corel Department of Civil and Environmental Engineering, University of Michigan, Ann Arbor Abstract. The aim of the present paper is to understand the interaction between a rising thermal and an inversion in the atmosphere and to quantify the turbulent entrainment rate due to a thermal impingement on a stratified interface. This problem was simulated in the laboratory using a water tank. The thermal is created by releasing a small volume of buoyant fluid into a stratified environment composed of two layers of different densities. A thin interface separates the lower layer from the lighter upper layer. The entrainment of upper layer fluid into the thermal is investigated using a passive dye flow visualization technique. The entrainment rate is found to obey a Ri -3/ power law, as predicted by Cotel and Breidenthal [1997]. The effect of simulated evaporative cooling on the entrainment of a thermal impinging on a stratified interface is also investigated experimentally. Evaporative cooling in atmospheri clouds is simulated in the laboratory using alcohol-water mixtures, so that the mixed fluid is denser than either parent parcel. This is realized in the laboratory by releasing a mixture of ethyl alcohol and ethylene glycol in an aqueous solution. It rises first through a relatively dense lower layer fluid and then impinges on a thin stratified interface, above which is a layer of relatively light fluid. The entrainment rate for values of the buoyancy reversal parameter D* between 0 and 0.5 was found to obey a Ri -3/2 power law. The entrainment rate is independent of D* between 0 and 0.5 for a range of Richardson numbers Ri from 3 to 25. This is consistent with the behavior of the buoyancy-reversing thermal in an unstratified environment observed by Johari [1992]. 1. Introduction Most fluid motions in the atmosphere and oceans are strongly influenced by the dynamics of stratified flows. It is therefore essential to understand their dynamics in order to comprehend the environment. This paper focuses on understanding the turbulent entrainment of a thermal impinging on a stratified interface. The goal of the present study is to quantify the entrainment rate due to the impinging thermal with and without buoyancy reversal. A direct application of this work is to study the interaction between clouds and the tropopause. Accurate modeling of mass transport across the tropopause is necessary to understand the concentration of ozone and other chemicals in the stratosphere. Some cumulus clouds contain pollutant air, and their penetration through atmospheric inversions may affect the ozone budget. The stratosphere is just like a "reservoir" for certain types of atmospheric pollution. The quantity of pollutants penetrating into the stratosphere might reduce the equilibrium concentration of ozone, thereby allowing more ultraviolet radiation to reach the Earth's surface [Wallace and Hobbs, 1977]. This is considered to be an important factor in the issue of global warming. It is necessary to gain further insight into stratified entrainment to yield better models of atmospheric dynamics in global climate models as well as Copyright 2000 by the American Geophysical Union. Paper number 2000JD /00/2000JD weather forecasts. The study of a rising thermal will also enhance our understanding of atmospheric convection, an essential part of atmospheric dynamics [Emanuel, 1994]. Among turbulent flows the motion of buoyant fluid known as the "thermal" has been studied theoretically, experimentally, and computationally for several decades. However, only preliminary work has been done to study the interaction between a stratified interface and a thermal, and no entrainment measurements have been performed for this specific case Background on Thermals The term thermal is used to denote a parcel of isolated buoyant fluid suddenly released from rest, which subsequently moves under the action of buoyancy forces alone. An isolated thermal is the ultimate form of convective elements in the atmosphere [Woodward, 1959]. From the motion of cumulus clouds convection, which is produced by buoyancy forces, appears to consist largely of more or less isolated masses of buoyant air rising into and mixing with their surroundings [Scorer, 1957]. A thermal is characterized by its total buoyancy force and does not possess any circulation, kinetic energy, or momentum at the moment of release. The circulation is generated entirely by buoyancy after the thermal is released. Entrainment in thermals is accomplished through engulfment of fluid in the rear of thermals because the induced velocity from the vorticity field causes ambient fluid to enter the backside of the thermal. The motion within a thermal is believed to be composed of a rising motion in the middle and a descending 15,457
2 15,458 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL where Ri = (2) w ' Figure 1. Cloud entrainment. one on the sides, similar to a vortex ring [Woodward, 1959; Turner, 1963]. The advancing front of a thermal has traditionally been used for measurementsince it is a relatively distinctive point [Johari, 1989]. A self-similar thermal has a width of about half the distance the front traveled [Johari, 1992] Background on Stratified Entrainment A. I. Taylor first described the entrainment hypothesis at the Pacific Science Association meeting in 1949; however, Batchelor [1954] gave it wider exposure through a review lecture. The basic form of this hypothesis is that the mean inflow velocity across the edge of a turbulent flow is assumed to be proportional to the local characteristic velocity, for a thermal, the rising speed (Figure 1). The total volume inflow depends on the surface area, the geometry, and the dynamics of the flow [Turner, 1986]. Scorer [1957], Woodward [1959], and Turner [1963] performed tank experiments involving heavy brine parcel sinking through freshwater. They observed that the diameter increases linearly with vertical distance from the source. Their results seemed to confirm the applicability of similarity solutions. Grabowski and Clark [1991, 1993] used two- and threedimensional numerical experiments to simulate thermals in a stratified environment. The subject of their studies, namely, the cloud-environment interface instability, is considered by them a primary entrainment mechanism. In stratified turbulence the interaction between stratified interfaces and turbulent flows is a critical problem. Experiments have shown that in general the entrainment rate declines with increasing stratification. The relationship between the entrainment rate and the Richardson number is observed to obey a power law [Turner, 1973] when the Richardso number Ri is between one and a few hundreds. The relationship is often defined as: g': ---. (3) p Here g' is the buoyancy acceleration due to the density difference A9; the acceleration of gravity is g; A 9 is the relevant density difference across the interface or between the inside and outside of a thermal; 9 is the density of the surrounding fluid; (5 is the length scale of the incident turbulence at the interface; and w i is its characteristic velocity (Figure 2). For Ri < 1 the kinetic energy of the largest eddies is greater than the potential energy investment in engulfing a tongue of fluid, so there is little or no effect of stratification on entrainment. When Ri is very large, the interface becomes flat, and so, potential energy can no longer affect the entrainment. Therefore the entrainment rate becomes independent of Ri [Cotel and Breidenthal, 1997]. Diffusion becomes the only available entrainment mechanism. The entrainment velocity we represents the rate at which the interface is rising because of entrainment of upper layer fluid into the turbulent flow emerging from the lower layer (see Figure 2). The dimensionless entrainment rate is defined as the ratio of the entrainment velocity w e to the characteristic velocity w 1- In the literature the value of the entrainment expo- nent a has been found to be -1/2, -1, -3/2 or, in one case, -2. Linden [1973] concluded from theoretical arguments that a should be -3/2 for a vortex ring. Kumagai [1984] measured a for a plume, finding it close to -3/2. For a jet impinging on a stratified interface [Cotel et al., 1997], a was measured to be -1/2. Turner [1973] also found that the value of a seemed to depend on the Schmidt or the Prandtl number from his stirring grid experiments with heat and salt. The Schmidt number Sc is defined to be SC ---Din, (4) where Dm is the molecular diffusivity and v is the kinematic viscosity. The Prandtl number Pr is Interface Pr = Dr' (5) We w1 where w e is the entrainment velocity across the interface, w i is the characteristic velocity of the impinging turbulent flow, and c and a are dimensionless constants. Ri (the ratio of potential energy to kinetic energy) is defined in terms of the impinging turbulence quantities: Figure 2. Schematic of a thermal impinging on a stratified interface.
3 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL 15,459 where D r is the thermal diffusivity. When Ri < 7, the value of a equals -3/2. If Ri > 7, the value of a remains -3/2 for salt stratification but becomes -1 for heat stratification in the case of stirring grid turbulence [Turner, 1973]. Cotel and Breidenthal [1997] suggested that the persistence of the entraining eddies is a critical parameter in determining the entrainment regime, in addition to the conventional parameters, such as Ri, Reynolds number Re, and Sc. Ri and Sc have been defined previously. The Reynolds number Re is defined by w 8 Re = --, (6) PA PB I Ptl Lid Plate PB>PA>Pt where 15 and w represent the characteristic lengthscale and Cylinder containing the thermal fluid velocity scale of the flow, respectively, and v is the kinematic viscosity. The dimensionless vortex persistence parameter T is defined to be the number of rotations a vortex makes during Figure 3. Sketch of the water tank. the time it moves its own diameter with respect to the inter- [1995]. According to Johari [1992] and Grabowski and Clark face. T is proportional to the ratio of the rotational to the [1993], the entrainment should not be influenced dramatically translational speed of the entraining eddy. The translational by buoyancy reversal as long as it is weak. Moreover, the speed is the speed at which the eddy is moving with respect to entrainment rate should decrease in general with increasing the interface. The rotational speed is unambiguous to any stratification as is the case for ordinary thermals. However, the inertial observer. T is independent of Ri, Re, or Sc since they effect of evaporative cooling has not yet been tested in the case all can be varied separately while T is held constant. According of impinging turbulence on a stratified interface. to Corel and Breidenthal's [1997] investigation a thermal has a A water tank with a thermal release mechanism was used for persistence parameter of unity or less because of its nonstathe experiments. The laboratory experiments and results for tionary character as it impinges on the interface, similar to that the classical thermals impinging on a stratified interface are of a vortex ring. The entrainment rate of the thermal is then presented in section 2. The buoyancy reversal case is discussed predicted to be proportional to Ri -3/2 for a thermal for a in section 3. Finally, the summary and conclusions are prerange ofri > 1 but <Ri = Re /4, at which point the interface becomes flat. Little laboratory work has been done on thermals interacting with a stratified interface. In the early 1960s, Richards [1962] and Saunders [1962] observed the penetration of thermals through an interface between different density fluids. However, they did not quantify the entrainment rate of thermals Background on Evaporative Cooling Evaporative cooling is another factor that affects the dynamics of clouds in the atmosphere. It occurs when dry environmental air is entrained into positively buoyant cumulus clouds, the resulting evaporation reducing or even completely reversing the positive buoyancy. Buoyancy reversal from evaporative cooling has been considered to have a dramatic effect on the global motions of clouds but is believed to have a lesser effect on the entrainment rate and mixing of buoyant fluid [Johari, 1992; Grabowski and Clark, 1993]. Grabowski [1995] suggested that the buoyancy reversal affects the entrainment in two different ways. First, buoyancy reversal directly enhances the kinetic energy of small-scal entraining eddies so the enhanced entrainment is thought to occur when the magnitude of the reversal is comparable to the initial stratification. Second, buoyancy reversal indirectly leads to more entrainment and mixing through large-scal eddies entraining air into cloud. In other words, enhancement of entrainment occurs through the enhancement of global flow structures. Furthermore, Turner and Yang's [1963] early laboratory experiments found that the rate of entrainment at the top of stratocumulus clouds was not enhanced significantly by modest buoyancy reversal. Buoyancy reversal has a stronger effect on the large-scale structure of the flow than on the small-scale dynamics, as observed in the laboratory by Johari [1992] and numerically by Grabowski sented in section Thermals Impinging on a Stratified Interface Without Buoyancy Reversal As discussed in section 1, there is an obvious need to measure the entrainment rate of a thermal impinging on a stratified interface Experimental Apparatus An ideal mechanism for releasing thermals would simulate the atmospheric release process as closely as possible. The dimensions of the tank are 30 x 30 x 65 cm in height. A vertical cylinder, 2.5 cm inside diameter and 2.9 cm high, is located at the bottom of the tank and covered by a thin, horizontal sliding stainles steel lid. For the purpose of generating a thin stratified interface the tank is separated into two chambers by a thin, horizontal sliding stainles steel plate. The material of the tank walls consists of chemically resistant lucite to give an unobstructed view of the flow. Thermals are released from the bottom of the tank as the rising clouds in the atmosphere. The apparatus is sketched in Figure Experimental Techniques Set up. The top and bottom layers are composed of two fluids of different salt concentrations. The density of the lower layer fluid Pa is greater than that of the upper layer fluid 9.4- The cylinder is filled with a mixture of pure water and water soluble dye, and the lid is slid back over the cylinder. Then the lower section of the tank is filled with a fluid of density Pa. The sliding plate is pushed in, and finally, a fluid with density 9.4 is poured into the upper section. The sliding plate and the lid are used to minimize mixing between the
4 15,460 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL Figure 4. Photograph of (a) a classical thermal just after release, (b) a classical thermal before impinging on the stratified interface, (c) a classical thermal impinging on the stratified interface, and (d) a classical thermal after impingement on the stratified interface. different fluids before a run. The density of the thermal p, is always less than both pu and P Experimental procedure. First, the sliding plate, located in the middle of the tank, is slowly retracted. In order to obtain sufficient entrainment and to minimize wall effects, the tank is then slowly drained so that the interface is positioned at a desired level. The initial level of the interface is set at 30 cm; the final position is between 20 and 30 cm from the tank bottom, depending on the conditions. The interface remains thin throughout the draining process. Then the lid above the cylinder is withdrawn to begin a run. Considerable attention is necessary at this stage so that minimal initial circulation is imposed on the thermal Experimental Results The evolution of a thermal rising toward the interface is shown in a sequence of photographs from a single run (Figures 4a-4d). The initial interface is marked with dots in the photographs. The interface is thin compared with the thermal eddy size. The buoyant thermal fluid is released from the cylinder at the bottom of the tank (Figure 4a). In Figure 4b the thermal becomes approximately self-similar as it has advanced 3.5 times the effective diameter of the released fluid [Johari, 1992]. The thermal penetrates through the interface, rebounds from it, and spreads laterally (Figure 4c). At this point in the exper- iment the entraining eddy may not be the largest one because of stratification. In an unstratified incompressiblenvironment with no external forcing applied the largest eddies are the most efficient eddies at entraining and mixing fluid [Roshko, 1976]. When stratification is present, the largest eddies' Ri is now greater than unity, and entrainment by these vortices is impeded. The most efficient entraining eddy is now the largest eddy with a Ri still below unity. This entrainment process causes the interface to rise by a certain amount. Finally, the motion decays, and the interface becomes flat (Figure 4d). Although a single thermal impinging on the interface described above is a single discrete event, by considering the process to be repeated many times we can define a corresponding average entrainment velocity: Mir We : --, (7) tit where H is the change in interface height due to thermal entrainment measured at t - t (where t represents the duration of the entrainment process) from the time the thermal is released until the thermal stops entraining fluid from the upper layer. The entrainment process starts shortly after the release. There is a small error in the time measurement, accounted for in the error bars. A stopwatch is used to record the different times. The thermal velocity is
5 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL 15,461 d tr W1 = tt r, where d tr is a distance of 2 cm just below the interface and t tr is the time the thermal traveled the distance dtr recorded with a stopwatch. Thus w i represents the thermal velocity just before impingement. The height of the interface and the distance traveled by the thermal are measured visually from a ruler on the side of the tank. The width of the thermal 8 is measured just before the thermal impinging on the interface (Figure 2). For all experiments the interface is reached by a thermal in a few seconds after release, but the entrainment process at the interface takes at least 40 s. The dimensionless entrainment rate we/w i versus the Ri is shown in log-log plots (Figure 5). The entrainment rate we/w 1 is the velocity at which the interface is rising normalized by the measured thermal velocity as it approaches the interface. Ri is based on the density difference between the upper and the lower layer fluid, and the density p is the density of the upper layer fluid (equation(2)). We found that the entrainment rate declines with increasing Ri as in (1). The value of the exponent determined by curve fitting is found to be -3/2. There is -_ % error due to lack of accuracy when measuring the time. Re defined in (5), is based on the thermal quantities: its diameter 8 and velocity w 1. The value of Re for all experiments is Discussion The results are similar to those obtained for stirring grid turbulence [Turner, 1973] and for a plume impinging on a stratified interface [Kumagai, 1984]. There is almost no entrainment at high Ri; because of the lack of measurement accuracy for such small values of the entrainment rate, no data are plotted above Ri - 30 in Figure 5. On the other hand, no matter which density difference is used in the definition of Ri, the same results are obtained ø 101 Ri 10 2 Figure 5. Entrainment rate w /w i versus Richardson number Ri for the classical thermal. 3. Evaporative Cooling Effects on Thermal Entrainment Evaporative cooling happens when a positively buoyant cloud parcel becomes negatively buoyant because of a change in density due to the evaporation [Johari, 1992]. For instance, atmospheric moist convection is an example of a buoyancy reversing system because the mixing of positively buoyant cloudy air and dry environmental air usually leads to the formation of negative buoyancy [Grabowski, 1995] Buoyancy Reversal Parameter In order to simulate the density changes accompanying evaporative cooling in a water tank, Johari [1992] used mixtures of ethyl alcohol and ethylene glycol as thermals, which were released in a uniform environment from the bottom of a The results confirm the prediction that the entrainment rate tank. Johari [1992] found that the behavior of thermals deof a thermal impinging on a stratified interface is proportional pends strongly on a dimensionless parameter, called the buoyto Ri -3/2 [Cotel and Breidenthal, 1997]. The present experi- ancy reversal parameter: ments reinforce the importance of the persistence parameter. As we mentioned in section 1, a thermal has a nonpersistent character, moving relatively rapidly away from the interface. p* - p(0) D* p(0) - p(1)' (8) Let us compare a thermal with a vortex ring, the former The density is a function of the volume fraction of released driven by buoyancy and the latter driven by momentum. Asthermal fluid. Thus p(0) is the density of the environment, and suming the Ri, Re, and Sc are the same for both cases, the p(1) is that of the initial thermal fluid. The peak density of a entrainment regime is completely independent of momentum given alcohol-glycol mixture is p*. or buoyancy. For both the vortex ring and the thermal the The same technique was used in our experiments. Several persistence parameter is the same, less than its presumed crit- curves (Figure 6) describe the density of thermals against the ical value of order unity, and the same entrainment law is volume fraction of thermal fluid P. The horizontal line repremeasured in spite of the different origins of the flows. As sents the density of the lower layer fluid initially surrounding further examples, a steady jet is persistent while a plume is not the thermal. The curves above and below this line show the since the plume fluctuates much more than the jet [Dai et al., 1995]. There is a direct correlation between the value of the persistence parameter and the entrainment regime. It does not matter whether the flow originates from buoyancy or momentum. A fully accurate laboratory simulation of atmospheric thermals would include the effect of evaporative cooling and latent heat. We choose to add the effect of buoyancy reversal to a thermal in order to simulate evaporative cooling in the laboratory. Evaporative cooling is more practical to simulate in the laboratory than latent heat release at this stage of our study. thermal with negative and positive buoyancy, respectively. The pure thermal fluid is represented by P = 1, and the ambient fluid is represented by P = 0. When P equals P*, the alcoholglycol mixture reaches its peak density p*. A thermal would then have maximum negative buoyancy. The dimensionless parameter D* (equation(8)) determines the intensity of buoyancy reversal. In the absence of any buoyancy reversal, D* = 0. Thermals rise continuously in a uniform environment where D* < 0.5. Within this regime, thermals do not stop or turn back even though they decelerate. They also leave some fluid
6 15,462 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL ø ø ' I Density of Water (gm/cc)... D*= D*= D*=0.471 _ m _ D*=0 I I i I I I I I I I II I,, I!, I I ø 'l ø Volume Fraction P Figure 6. Density p as a function of volume fraction P. behind in their wake. For the values of D* close to 0.5, thermals detrain a large volume of fluid [Johari, 1992]. In general, the trajectory of thermals with buoyancy reversal depends strongly on D*, and the "evaporating" thermals do not penetrate as far as the thermals without buoyancy reversal. The range of D* for the present experiments is from 0 to 0.5. Within this range, thermals rise toward the top of the tank and do not break into two parts or turn down shortly after having been released, as discussed above. There is a small amount of entrainment across the interface when the range of D* is >0.5, so that no accurate measurements could be performed. We chose to restrict the range of D* so that the thermals produced are representative of the real atmospheri case [Siems et al., 1990]. alcohol varies from 0.04 to 0.17% in this layer for different runs. The lower layer fluid is pure water with density greater than p.. The release mechanism is a cylinder covered with a small sliding lid to prevent the alcohol-glycol mixture from mixing with the surrounding fluid. The mixture is visualized by adding water soluble dye. The volume ratio of alcohol to glycol in the initial thermals is 94 to 71% for different runs. The densities of the two fluids, p. and pb, are always greater than the thermal density lot. The experimental procedure used in the nonevaporative case is repeated for the evaporative case. To ensure that thermals are reaching the interface without wall effects, the interface is lowered by draining the water from the lower chamber. The appropriate position was found to be at a height between 14 and 20 cm from the bottom of the tank for different runs. The draining process must be slow to avoid disturbing the interface. Then the lid is pulled out to begin a run Results The buoyant thermal fluid mixture rising from the cylinder at the bottom of the tank is shown in Figure 7a. When the initial positively buoyant thermal begins to mix with the ambient fluid, the density of the thermal is increased but is still less than that of the lower layer (Figure 7a). A partially negatively buoyanthermal is generated as shown in Figure 7b. It appears that the negative buoyancy decelerates the thermal. The thermal continues to travel upward, because the total density of the thermal is still less than that of the environment, penetrates through the interface, rebounds from it, and spreads laterally along the interface (Figure 7c). The interface rises a certain amount because of turbulent entrainment of the upper layer fluid into the thermal. The thermal collapsed after several seconds. Finally, the thermal entrainment process is finished (Figure 7d). The measurement techniques used in the classical thermal case are repeated here, and the same quantities are measured, i.e., the instantaneous width of the thermal below the interface, the velocity of the thermal, and the entrainment velocity. Johari's [1992] experiment show that the molecular-scale mixing features of these thermals are similar to those of rising thermals without buoyancy reversal as long as these thermals The dimensionless entrainment rate W e/w versus Ri is shown in six log-log plots (Figures 8a-8e), and the value of the exponent is determined from curve fitting to between maintain upward motion, even though the trajectory is strongly and -1.63, following (1). The entrainment rate W e/w is deaffected by D*. This implies that buoyancy reversal does not fined as before, i.e., on the basis of the thermal quantities influence the fundamental mixing and entrainment dynamics below the interface. Ri is based on the same thermal quantiof convective elements; that is, entrain tongues of ambient fluid at each vortex rotation independent of buoyancy reversal. ties. The entrainment rates clearly decline with increasing Ri (Figures 8a-8e) for different values of D*. The relationship The same result was also found in numerical experiments. between the entrainment rate W e/w and the dimensionless Grabowski and Clark [1993] and Grabowski [1995] believe that entrainment is driven by the baroclinic instability and not by parameter D* between 0 and 0.5 is plotted in log-log plots for several values of Ri (Figures 9a-9c). The slope between We/W the buoyancy reversal process. They suggested that buoyancy and D* is from 0.03 to 0.08, obtained by curve fitting. The reversal associated with evaporative cooling affects the overall error is of the order of _+10-20% because of the lack of flow evolution but has small effects on the dynamics of the entraining eddies. In other words, buoyancy reversal has little effect at the initial stages of cloud development, but the entrainment rate is likely to change significantly when the buoyancy reversal becomes very strong. accuracy when measuring the time intervals. Re based on the thermal quantities at the interface is 104. The possibility of buoyancy reversal between the lower and upper layer fluids has been tested. No buoyancy reversal between those two layers was observed Experimental Procedure The apparatus is identical to that used for the experiments of the classical thermals (Figure 3). Two fluids with different densities represent two layers separated by a stainles steel sliding plate in the middle of the tank. The upper layer fluid is a mixture of water and alcohol with density p.. The content of 3.4. Discussion The exponent a of the power law describing the entrainment rate as a function of Ri is similar to the one found for classical thermals. A difference in the magnitude of the dimensionless coefficient c for the entrainment rate with and without buoyancy reversal of a factor of 10 (Figure 10) exists. This can be
7 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL 15,463 Figure 7. (a) a thermal with buoyancy reversal (D* = 0.16) just after release, (b) a thermal with buoyancy reversal (D* = 0.16) before impinging on the stratified interface, (c) a thermal with buoyancy reversal (D* = 0.16) impinging on the stratified interface, and (d) a thermal with buoyancy reversal (D* = 0.16) after impingement on the stratified interface.
8 15,464 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL Figure 7. (continued) explained by the decrease in size and acceleration of the thermal due to buoyancy reversal. As seen in Figure 7c, some of the thermal fluid is left behind as the thermal rises. A net decrease in velocity is also easily observed. The entrainment velocity at large Ri could not be measured because of low measurement accuracy. These results indicate that the relationship between the entrainment rate and Ri follows the same power law for thermals with or without buoyancy reversal, impinging on a stratified interface. The results in Figures 9a-9c demonstrate that the entrain-
9 ... ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL 15,465 a loø2 C 10-2 D*=0.053 D*=0.16 lo a.46 lo l0 ø l01 ] ø 10' 102 gi d 10-2 Ri b 10'2 D*=0.078 D*=0.214 l I i i i,,,,, i,,,,,!,, 10 ø RI ø l0 t 10 2 D*=0.471 Ri 10-4 I ø 101 Ri lo 2 Figure 8. Entrainment rate We/Wl versus Ri for different values ofd*' (a) D* = 0.053, (b) D* = 0.078, (c) D* = 0.16, (d)d* = 0.214, and (e)d* = ment rates can be assumed to be independent of the buoyancy reversal parameter for 0 < D* < 0.5 for the case of thermals impinging on a stratified interface. This is in agreement with Grabowski and Clark [1993], who stipulate that the effects of buoyancy reversal are negligible on the dynamics of entraining thermals when the buoyancy reversal parameter is small. These results are also supported by Johari's [1992] laboratory experiments. The mixing and entrainment dynamics of convective elements do not depend on the presence of the buoyancy reversal phenomenon as long as upward motion is maintained. Therefore that buoyancy reversal does not significantly affect the rate of mixing between convecting fluid and its environ- ment is believed. 4. Summary and Conclusions Laboratory experiments provide a powerful tool for investigating the fundamentals of stratified turbulence, leading to new insight in natural atmospheric phenomena. The present experiments aim to improve the current understanding of transport across the tropopause by making the first measurements of entrainment across an interface due to a rising ther- mal with and without evaporative cooling. For thermals without buoyancy reversal the entrainment rate was found to be proportional to Ri - /2, as predicted by Corel and Breidenthal [1997]. The effect of evaporative cooling on thermal entrainment was investigated in the laboratory. The
10 15,466 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL l0-2 Ri=4 10 ø D ß Ri=6 I i I I I I I II I I I I I,, I ø D*,, I, i, i Ri=10 of 10) between the dimensionless entrainment rate coefficients with and without buoyancy reversal because of the deceleration and the reduced size of the thermal in the buoyancy reversal case. Ri ranged between 2.5 and 25, Re was of the order of 104, and S c was around 600 for the experiments performed. The large differences in Re and Sc between the laboratory and the atmosphere are important issues for geophysical interpretation of laboratory results. Re in the atmosphere ( ) is much larger than Re in the laboratory (of the order of 104). The results from laboratory experiments are expected to be applicable to the atmosphere as long as Re is above a critical value of order 103 called the mixing transition [Breidenthal and Baker, 1985], which is necessary for the development of small-scale turbulence. Furthermore, the entrainment rate of a jet impinging on a stratified interface was found to be independent of Re for Ri between 1 and 100 [Cotel, 1995]. Moreover, Siems et al. [1990] analyzed soundings from atmospheric data on stratocumulus and trade cumulus clouds, taken from the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE) in 1987 and Joint Hawaiian Warm Rain Project experiment in 1985, respectively. They suggesthat D* observed in stratocumulus is often <1.3. The buoyancy reversal parameter in the laboratory is of the same order as that in the atmosphere [Siems et al., 1990]. The laboratory experiments represent a case of weak buoyancy reversal. The experiments reported here represent the first step in quantifying turbulent entrainment of thermals impinging on a stratified interface. The logical next step would be to include latent heat to simulate completely the atmosphericase in the laboratory. Also, the amount of mixing produced by a thermal penetrating in a stratified environment needs to be quantified using laser-induced fluorescence. Particle image velocimetry measurements of the velocity and vorticity fields would be useful in determining the interface dynamics. These future steps would further our understanding of thermals in a strati- fied environment I I i I I i i i I i i i i i i i i D* 10-3 Figure 9. Entrainment rate We/W1 versus buoyancy reversal parameter D* for different values of Ri: (a) Ri = 4, (b) Ri = 6, and (c) Ri = 10. entrainment rate was also found to be proportional to Ri -3/2. In our experiments the value of the entrainment rate exponent does not depend on the buoyancy reversal parameter D* for the range of 0-0.5, which is representative of the atmosphere [Siems et al., 1990]. There is a difference in magnitude (factor 10-5 Thermal with Buoyancy Reversal D*=0.053 i, i i i i I! i, I, i i, 10 ø Ri Figure 10. Comparison of the entrainment rate We/W versus Ri with (D* = 0.053) and without buoyancy reversal.
11 ZHANG AND COTEL: ENTRAINMENT DUE TO AN IMPINGING THERMAL 15,467 Acknowledgments. This work was supported by the Natural Sciences and Engineering Research Council of Canada and a University Linden, P. F., The interaction of a vortex ring with a sharp density interface, J. Fluid Mech., 60, , of Manitoba Research Grant. Richards, J. M., Experiments on the penetration of an interface by References buoyant thermals, J. Fluid Mech., 11, , Roshko, A., Turbulent shear flow: A new look, A/AA J., 14, , Batchelor, G. K., Heat convection and buoyancy effects in fluids, Saunders, P.M., Penetrative convection in stably stratified fluids, Tellus, 14, , Q. J. R. Meteorol. Soc., 80, , Breidenthal, R. E., and M. B. Baker, Convection and entrainment Scorer, R. S., Experiments on convection of isolated masses of buoyant fluid, J. Fluid Mech., 2, , across stratified interfaces, J. Geophys. Res., 90, 13,055-13,062, Cotel, A. J., Entrainment and detrainment of a jet impinging on a Siems, S. T., S. B. Christopher, M. B. Baker, S.S. Shy, and R. E. stratified interface, Ph.D. thesis, Univ. of Wash., Seattle, Breidenthal, Buoyancy reversal and cloud-top entrainment instabil- Cotel, A. J., and R. E. Breidenthal, A model of stratified entrainment ity, Q. J. R. Meteorol. Soc., 116, , Turner, J. S., Model experiments relating to thermals with increasing using vortex persistence, Appl. Sci. Res., 57, , Cotel, A. J., J. A. Gjestvang, N. N. Ramkhelawan, and R. E. Brebuoyancy, Q. J. R. Meteorol. Soc., 89, 62-74, Turner, J. S., Buoyancy Effects in Fluids, 368 pp., Cambridge Univ. identhal, Laboratory experiments of a jet impinging on a stratified Press, New York, interface, Exp. Fluids, 23, , Turner, J. S., Turbulence entrainment: The development of the en- Dai, Z., L. K. Tseng, and G. M. Faeth, Velocity statistics of round, trainment assumption, and its application to geophysical flows, J. fully-developed, buoyant turbulent plumes, J. Heat Transfer, 117, Fluid Mech., 173, , , Turner, J. S., and I. K. Yang, Turbulent mixing at the top of stratocu- Emanuel, K. A., Atmospheric Convection, 580 pp., Oxford Univ. Press, mulus clouds, J. Fluid Mech., 17, , New York, Wallace, J. M., and P. V. Hobbs, Atmospheric science: An introduc- Grabowski, W. W., Entrainment and mixing in buoyancy-reversing tory survey, Academic, San Diego, Calif., convection with applications to cloud-top entrainment instability, Woodward, B., The motion in and around isolated thermals, Q. J. R. Q. J. R. Meteorol. Soc., 121, , Meteorol. Soc., 85, , Grabowski, W. W., and T. L. Clark, Cloud-environment interface instability: Rising thermal calculations in two spatial dimensions, J. Atmos. Sci., 48, , A. Cotel, Department of Civil and Environmental Engineering, Uni- Grabowski, W. W., and T. L. Clark, Cloud-environment interface versity of Michigan, Ann Arbor, MI (acotel@engin. instability, part II, Extension to three spatial dimensions, J. Atmos. umich.edu) Sci., 50, , Q. Zhang, Department of Mechanical and Industrial Engineering, Johari, H., An experimental investigation of mixing in buoyant flows, University of Manitoba, Winnipeg, Manitoba, Canada R3T 2N2. Ph.D. thesis, Univ. of Wash., Seattle, Johari, H., Mixing in thermals with and without buoyancy reversal, J. Atmos. Sci., 49, , Kumagai, M., Turbulent buoyant convection from a source in a confined two-layer region, J. Fluid Mech., 147, , (Received June 23, 1999; revised November 15, 1999; accepted December 18, 1999.)
Study on Fire Plume in Large Spaces Using Ground Heating
Available online at www.sciencedirect.com Procedia Engineering 11 (2011) 226 232 The 5 th Conference on Performance-based Fire and Fire Protection Engineering Study on Fire Plume in Large Spaces Using
More informationDUE TO EXTERNAL FORCES
17B.6 DNS ON GROWTH OF A VERTICAL VORTEX IN CONVECTION DUE TO EXTERNAL FORCES Ryota Iijima* and Tetsuro Tamura Tokyo Institute of Technology, Yokohama, Japan 1. INTRODUCTION Various types of vertical vortices,
More informationChapter 2. Turbulence and the Planetary Boundary Layer
Chapter 2. Turbulence and the Planetary Boundary Layer In the chapter we will first have a qualitative overview of the PBL then learn the concept of Reynolds averaging and derive the Reynolds averaged
More informationINTERACTION BETWEEN WIND-DRIVEN AND BUOYANCY-DRIVEN NATURAL VENTILATION Bo Wang, Foster and Partners, London, UK
INTERACTION BETWEEN WIND-DRIVEN AND BUOYANCY-DRIVEN NATURAL VENTILATION Bo Wang, Foster and Partners, London, UK ABSTRACT Ventilation stacks are becoming increasingly common in the design of naturally
More informationInfluence of rounding corners on unsteady flow and heat transfer around a square cylinder
Influence of rounding corners on unsteady flow and heat transfer around a square cylinder S. K. Singh Deptt. of Mech. Engg., M. B. M. Engg. College / J. N. V. University, Jodhpur, Rajasthan, India Abstract
More informationQuantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges
Quantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges T. Abbas 1 and G. Morgenthal 2 1 PhD candidate, Graduate College 1462, Department of Civil Engineering,
More informationAn experimental study of internal wave generation through evanescent regions
An experimental study of internal wave generation through evanescent regions Allison Lee, Julie Crockett Department of Mechanical Engineering Brigham Young University Abstract Internal waves are a complex
More informationInvestigation of Suction Process of Scroll Compressors
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2006 Investigation of Suction Process of Scroll Compressors Michael M. Cui Trane Jack Sauls
More informationLOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12
LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/ This experiment will introduce you to the kinetic properties of low-pressure gases. You will make observations on the rates with which selected
More informationSuper-parameterization of boundary layer roll vortices in tropical cyclone models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-parameterization of boundary layer roll vortices in tropical cyclone models PI Isaac Ginis Graduate School of Oceanography
More information10.6 The Dynamics of Drainage Flows Developed on a Low Angle Slope in a Large Valley Sharon Zhong 1 and C. David Whiteman 2
10.6 The Dynamics of Drainage Flows Developed on a Low Angle Slope in a Large Valley Sharon Zhong 1 and C. David Whiteman 2 1Department of Geosciences, University of Houston, Houston, TX 2Pacific Northwest
More information3.3 USING A SIMPLE PARCEL MODEL TO INVESTIGATE THE HAINES INDEX
3.3 USING A SIMPLE PARCEL MODEL TO INVESTIGATE THE HAINES INDEX Mary Ann Jenkins 1 Steven K. Krueger 2 and Ruiyu Sun 2 1 York University, Toronto, Canada 2 University of Utah, Salt Lake City, Utah 1. INTRODUCTION
More informationScott Denning CSU CMMAP 1
Thermodynamics, Buoyancy, and Vertical Motion Temperature, Pressure, and Density Buoyancy and Static Stability Adiabatic Lapse Rates Dry and Moist Convective Motions Present Atmospheric Composition What
More informationUnsteady Wave-Driven Circulation Cells Relevant to Rip Currents and Coastal Engineering
Unsteady Wave-Driven Circulation Cells Relevant to Rip Currents and Coastal Engineering Andrew Kennedy Dept of Civil and Coastal Engineering 365 Weil Hall University of Florida Gainesville, FL 32611 phone:
More informationNumerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water
Numerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water Hong Xu, Chokri Guetari ANSYS INC. Abstract Transient numerical simulations of the rise of a train of gas bubbles in a liquid
More informationThe effect of back spin on a table tennis ball moving in a viscous fluid.
How can planes fly? The phenomenon of lift can be produced in an ideal (non-viscous) fluid by the addition of a free vortex (circulation) around a cylinder in a rectilinear flow stream. This is known as
More informationIntegration of Impact Factors of Gas-Liquid Transfer Rate
Zhiyong Duan and James L. Martin Department of Civil Engineering Mississippi State University P.O. Box 9546 Mississippi State, MS 39762-9546 662-325-2902 E-mail: zd9@msstate.edu ABSTRACT The gas transfer
More informationGeophysical Fluid Dynamics of the Earth. Jeffrey B. Weiss University of Colorado, Boulder
Geophysical Fluid Dynamics of the Earth Jeffrey B. Weiss University of Colorado, Boulder The Earth is a spinning sphere Coriolis force depends on latitude solar flux depends on latitude Michael Ritter,
More informationVertical Motion and Atmospheric Stability
Lesson 4 Vertical Motion and Atmospheric Stability This lesson describes the vertical structure of the atmosphere, atmospheric stability and the corresponding vertical motion. Adiabatic diagrams are introduced
More informationLOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary
ADH 1/7/014 LOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary This experiment will introduce you to the kinetic properties of low-pressure gases. You will make observations on the
More informationForest Winds in Complex Terrain
Forest Winds in Complex Terrain Ilda Albuquerque 1 Contents Project Description Motivation Forest Complex Terrain Forested Complex Terrain 2 Project Description WAUDIT (Wind Resource Assessment Audit and
More informationThe Effect of Von Karman Vortex Street on Building Ventilation
The Effect of Von Karman Vortex Street on Building Ventilation P.Praveen Kumar Abstract This paper deals with the utilisation of the von Karman vortex street principle to maximise air flow into buildings.
More informationFluid Flow. Link. Flow» P 1 P 2 Figure 1. Flow Model
Fluid Flow Equipment: Water reservoir, output tubes of various dimensions (length, diameter), beaker, electronic scale for each table. Computer and Logger Pro software. Lots of ice.temperature probe on
More informationNumerical Investigation of Air Bubbles Evolution and Coalescence from Submerged Orifices Based on OpenFOAM
Numerical Investigation of Air Bubbles Evolution and Coalescence from Submerged Orifices Based on OpenFOAM Pan Feng, He Ying, Li-zhong Mu 2018-7-6 Dalian University of Technology, China Outline Background
More informationWinds and Ocean Circulations
Winds and Ocean Circulations AT 351 Lab 5 February 20, 2008 Sea Surface Temperatures 1 Temperature Structure of the Ocean Ocean Currents 2 What causes ocean circulation? The direction of most ocean currents
More informationChapter 3 Atmospheric Thermodynamics
Chapter 3 Atmospheric Thermodynamics Spring 2017 Partial Pressure and Dalton Dalton's law of partial pressure: total pressure exerted by a mixture of gases which do not interact chemically is equal to
More informationENVIRONMENTAL PHYSICS
ENVIRONMENTAL PHYSICS Atmospheric Stability An understanding of why and how air moves in the atmosphere is fundamental to the prediction of weather and climate. What happens to air as it moves up and down
More informationMeteorology & Air Pollution. Dr. Wesam Al Madhoun
Meteorology & Air Pollution Dr. Wesam Al Madhoun Dispersion = Advection (Transport) + Dilution (Diffusion) Source Transport Receptor Re-entrainment Fick s law of diffusion J= - D * D C/Dx Where, J= Mass
More informationATMS 310 Tropical Dynamics
ATMS 310 Tropical Dynamics Introduction Throughout the semester we have focused on mid-latitude dynamics. This is not to say that the dynamics of other parts of the world, such as the tropics, are any
More informationNatural ventilation of an enclosure containing two buoyancy sources
J. Fluid Mech. (1996), vol. 311,pp. 153-176 Copyright 0 1996 Cambridge University Press 153 Natural ventilation of an enclosure containing two buoyancy sources By P. COOPER' AND P. F. LINDEN2 Department
More informationA Study on the Effects of Wind on the Drift Loss of a Cooling Tower
A Study on the Effects of Wind on the Drift Loss of a Cooling Tower Wanchai Asvapoositkul 1* 1 Department of Mechanical Engineering, Faculty of Engineering, King Mongkut s University of Technology Thonburi
More informationLarge-eddy simulation of a turbulent buoyant helium plume
Center for Turbulence Research Annual Research Briefs 8 45 Large-eddy simulation of a turbulent buoyant helium plume By G. Blanquart AND H. Pitsch. Motivation and objectives The numerical simulation of
More informationNumerical and Experimental Investigation of the Possibility of Forming the Wake Flow of Large Ships by Using the Vortex Generators
Second International Symposium on Marine Propulsors smp 11, Hamburg, Germany, June 2011 Numerical and Experimental Investigation of the Possibility of Forming the Wake Flow of Large Ships by Using the
More informationZIN Technologies PHi Engineering Support. PHi-RPT CFD Analysis of Large Bubble Mixing. June 26, 2006
ZIN Technologies PHi Engineering Support PHi-RPT-0002 CFD Analysis of Large Bubble Mixing Proprietary ZIN Technologies, Inc. For nearly five decades, ZIN Technologies has provided integrated products and
More informationPlane Turbulent Wall Jets in Limited Tailwater Depth
International Journal of Engineering & Technology IJET-IJENS Vol: 11 No: 159 Plane Turbulent Wall Jets in Limited Tailwater Depth Shazy A. Shabayek 1 Abstract This paper presents laboratory study of plane
More informationAir Pollution Dispersion
Air Pollution Dispersion Dispersion Processes Convective Dispersion Air Parcel Dynamics Adiabatic Process Lapse Rate Equilibrium and Stability Atmospheric Stability Stability and Dispersion Temperature
More informationThe total precipitation (P) is determined by the average rainfall rate (R) and the duration (D),
Orographic precipitation Common ingredients of heavy orographic precipitation The total precipitation (P) is determined by the average rainfall rate (R) and the duration (D), P = RD. (1) The rainfall rate
More informationEXPERIMENTAL STUDY ON THE DISCHARGE CHARACTERISTICS OF SLUICE FOR TIDAL POWER PLANT
EXPERIMENTAL STUDY ON THE DISCHARGE CHARACTERISTICS OF SLUICE FOR TIDAL POWER PLANT Sang-Ho Oh 1, Kwang Soo Lee 1 and Dal Soo Lee 1 The discharge characteristics of sluice caisson for tidal power plant
More informationMountain Forced Flows
Mountain Forced Flows Jeremy A. Gibbs University of Oklahoma gibbz@ou.edu February 3, 2015 1 / 45 Overview Orographic Precipitation Common Ingredients of Heavy Orographic Precipitation Formation and Enhancement
More informationThe Monsoon and Its Variability Prof. Sulochana Gadgil Centre for Atmospheric & Oceanic Sciences Indian Institute of Science Bangalore
The Monsoon and Its Variability Prof. Sulochana Gadgil Centre for Atmospheric & Oceanic Sciences Indian Institute of Science Bangalore Lecture - 04 Background about the Atmosphere and Rotating Systems
More informationYasuyuki Hirose 1. Abstract
Study on Tsunami force for PC box girder Yasuyuki Hirose 1 Abstract In this study, a waterway experiment was performed in order to understand the influence of tsunami forms on tsunami forces acting on
More informationWIND-INDUCED LOADS OVER DOUBLE CANTILEVER BRIDGES UNDER CONSTRUCTION
WIND-INDUCED LOADS OVER DOUBLE CANTILEVER BRIDGES UNDER CONSTRUCTION S. Pindado, J. Meseguer, J. M. Perales, A. Sanz-Andres and A. Martinez Key words: Wind loads, bridge construction, yawing moment. Abstract.
More informationEntrainment and Mixing in Cumulus Clouds
Introduction NCAS, University of Leeds 25 March 2009 Outline Introduction Outline Cloud microphysics Cumulus clouds Start Movie Liquid water content Introduction Outline Cloud microphysics From: Warner
More informationImprovement of the Volumetric and Isentropic Efficiency Due to Modifications of the Suction Chamber
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2008 Improvement of the Volumetric and Isentropic Efficiency Due to Modifications of the
More informationThe Influence of Ocean Surface Waves on Offshore Wind Turbine Aerodynamics. Ali Al Sam
The Influence of Ocean Surface Waves on Offshore Wind Turbine Aerodynamics Ali Al Sam What I m going to wear today? Do I need to leave early to get to work? Taking buss or riding bike? Where will we drink
More informationLecture 7. More on BL wind profiles and turbulent eddy structures. In this lecture
Lecture 7. More on BL wind profiles and turbulent eddy structures In this lecture Stability and baroclinicity effects on PBL wind and temperature profiles Large-eddy structures and entrainment in shear-driven
More informationSingle Phase Pressure Drop and Flow Distribution in Brazed Plate Heat Exchangers
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 2016 Single Phase Pressure Drop and Flow Distribution in Brazed Plate Heat Exchangers
More informationSOLAR 93 THE 1993 AMERICAN SOLAR ENERGY SOCIETY ANNUAL CONFERENCE. Washington, DC April 22028,1993. Editors: S. M. Burley M. E.
SOLAR 93 THE 1993 AMERICAN SOLAR ENERGY SOCIETY ANNUAL CONFERENCE Washington, DC April 22028,1993 Editors: S. M. Burley M. E. Arden American Solar Energy Society U.S. Section of the International Solar
More informationPHSC 3033: Meteorology Stability
PHSC 3033: Meteorology Stability Equilibrium and Stability Equilibrium s 2 States: Stable Unstable Perturbed from its initial state, an object can either tend to return to equilibrium (A. stable) or deviate
More informationMeteorology. Circle the letter that corresponds to the correct answer
Chapter 4 Worksheet 3 Meteorology Name: Circle the letter that corresponds to the correct answer 1) Natural convection and turbulence are most likely to occur when: a) temperature decreases rapidly with
More informationGravity waves in stable atmospheric boundary layers
Gravity waves in stable atmospheric boundary layers Carmen J. Nappo CJN Research Meteorology Knoxville, Tennessee 37919, USA Abstract Gravity waves permeate the stable atmospheric planetary boundary layer,
More informationA Computational Assessment of Gas Jets in a Bubbly Co-Flow 1
A Computational Assessment of Gas Jets in a Bubbly Co-Flow 1 Melissa Fronzeo*, 1 Michael Kinzel 1 The Pennsylvania State University, University Park, PA, USA Abstract In this effort, Computational Fluid
More informationLecture 7. The Indian monsoon: is it a gigantic land-sea breeze?
Lecture 7 The Indian monsoon: is it a gigantic land-sea breeze? In the next set of lectures I am going to discuss the different hypotheses put forth for the basic system responsible for the monsoon. I
More informationWind tunnel effects on wingtip vortices
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 2010, Orlando, Florida AIAA 2010-325 Wind tunnel effects on wingtip vortices Xin Huang 1, Hirofumi
More informationFigure 1. Outer dimensions of the model trucks.
1 Standardisation of test method for salt spreader: Air flow experiments Report 3: Simulations of airflow patterns by Jan S. Strøm, Consultant Aarhus University, Engineering Centre Bygholm, Test and Development
More information2.4. Applications of Boundary Layer Meteorology
2.4. Applications of Boundary Layer Meteorology 2.4.1. Temporal Evolution & Prediction of the PBL Earlier, we saw the following figure showing the diurnal evolution of PBL. With a typical diurnal cycle,
More informationAlongshore wind stress (out of the page) Kawase/Ocean 420/Winter 2006 Upwelling 1. Coastal upwelling circulation
Kawase/Ocean 420/Winter 2006 Upwelling 1 Coastal upwelling circulation We found that in the northern hemisphere, the transport in the surface Ekman layer is to the right of the wind. At the bottom, there
More informationBlast Damage Consideratons for Horizontal Pressure Vessel and Potential for Domino Effects
A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 26, 2012 Guest Editors: Valerio Cozzani, Eddy De Rademaeker Copyright 2012, AIDIC Servizi S.r.l., ISBN 978-88-95608-17-4; ISSN 1974-9791 The Italian
More informationPURE SUBSTANCE. Nitrogen and gaseous air are pure substances.
CLASS Third Units PURE SUBSTANCE Pure substance: A substance that has a fixed chemical composition throughout. Air is a mixture of several gases, but it is considered to be a pure substance. Nitrogen and
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I Fluid statics. Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET
COURSE NUMBER: ME 321 Fluid Mechanics I Fluid statics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Fluid statics Fluid statics is the study of fluids in
More informationATS 351 Lecture 6. Air Parcel. Air Parcel Movement: Why does rising air expand and cool? Stability & Skew-T Diagrams
ATS 351 Lecture 6 Stability & Skew-T Diagrams To demonstrate stability, a parcel of air is used Expands and contracts freely Always has uniform properties throughout Air Parcel Air Parcel Movement: Why
More informationBioreactor System ERT 314. Sidang /2011
Bioreactor System ERT 314 Sidang 1 2010/2011 Chapter 2:Types of Bioreactors Week 4 Flow Patterns in Agitated Tanks The flow pattern in an agitated tank depends on the impeller design, the properties of
More informationBubble-bubble interactions and wall pressures/temperatures produced by the collapse of a bubble pair near a rigid surface
Bubble-bubble interactions and wall pressures/temperatures produced by the collapse of a bubble pair near a rigid surface Abstract Shahaboddin A. Beig*; Eric Johnsen; Mechanical Engineering Department,
More informationAir Flow Exchange Velocity of Urban Canyon Cavities due to Thermal Spatial Differences
Air Flow Exchange Velocity of Urban Canyon Cavities due to Thermal Spatial Differences Marta Oliveira Panão 1, Helder Gonçalves 1 and Paulo Ferrão 1 Renewable Energy Department, INETI, Lisbon, Portugal
More information9 Mixing. I Fundamental relations and definitions. Milan Jahoda revision Radim Petříček, Lukáš Valenz
9 ixing ilan Jahoda revision 14-7-017 Radim Petříček, Lukáš Valenz I Fundamental relations and definitions ixing is a hydrodynamic process, in which different methods are used to bring about motion of
More informationA Simple Horizontal Velocity Model of a Rising Thermal Instability in the Atmospheric Boundary Layer
A Simple Horizontal Velocity Model of a Rising Thermal Instability in the Atmospheric Boundary Layer (In other words, how to take a closer look at downwind thermal drift) Abstract This pilot, as well as
More informationThe Coriolis force, geostrophy, Rossby waves and the westward intensification
Chapter 3 The Coriolis force, geostrophy, Rossby waves and the westward intensification The oceanic circulation is the result of a certain balance of forces. Geophysical Fluid Dynamics shows that a very
More informationDetermination of the wind pressure distribution on the facade of the triangularly shaped high-rise building structure
Determination of the wind pressure distribution on the facade of the triangularly shaped high-rise building structure Norbert Jendzelovsky 1,*, Roland Antal 1 and Lenka Konecna 1 1 STU in Bratislava, Faculty
More informationExperimental Characterization and Modeling of Helium Dispersion in a ¼-Scale Two-Car Residential Garage
Experimental Characterization and Modeling of Helium Dispersion in a ¼-Scale Two-Car Residential Garage William M. Pitts, Jiann C. Yang, Kuldeep Prasad, and Marco Fernandez National Institute of Standards
More informationMicro Channel Recuperator for a Reverse Brayton Cycle Cryocooler
Micro Channel Recuperator for a Reverse Brayton Cycle Cryocooler C. Becnel, J. Lagrone, and K. Kelly Mezzo Technologies Baton Rouge, LA USA 70806 ABSTRACT The Missile Defense Agency has supported a research
More informationA Theory for Strong Long-Lived Squall Lines Revisited
A Theory for Strong Long-Lived Squall Lines Revisited MORRIS L. WEISMAN AND RICHARD ROTUNNO National Center for Atmospheric Research,* Boulder, Colorado (Manuscript received 29 April 2002, in final form
More informationApplication of Simulation Technology to Mitsubishi Air Lubrication System
50 Application of Simulation Technology to Mitsubishi Air Lubrication System CHIHARU KAWAKITA *1 SHINSUKE SATO *2 TAKAHIRO OKIMOTO *2 For the development and design of the Mitsubishi Air Lubrication System
More information3 1 PRESSURE. This is illustrated in Fig. 3 3.
P = 3 psi 66 FLUID MECHANICS 150 pounds A feet = 50 in P = 6 psi P = s W 150 lbf n = = 50 in = 3 psi A feet FIGURE 3 1 The normal stress (or pressure ) on the feet of a chubby person is much greater than
More informationOcean circulation and surface buoyancy fluxes: dynamics and energetics!
Ocean circulation and surface buoyancy fluxes: dynamics and energetics Ross Griffiths with credits: Graham Hughes, Andy Hogg, Kial Stewart, Julia Mullarney, J.Tan Research School of Earth Sciences The
More informationDevelopment of High-speed Gas Dissolution Device
Development of High-speed Gas Dissolution Device Yoichi Nakano*, Atsushi Suehiro**, Tetsuhiko Fujisato***, Jun Ma**** Kesayoshi Hadano****, Masayuki Fukagawa***** *Ube National College of Technology, Tokiwadai
More informationIMAGE-BASED STUDY OF BREAKING AND BROKEN WAVE CHARACTERISTICS IN FRONT OF THE SEAWALL
IMAGE-BASED STUDY OF BREAKING AND BROKEN WAVE CHARACTERISTICS IN FRONT OF THE SEAWALL Weijie Liu 1 and Yoshimitsu Tajima 1 This study aims to study the breaking and broken wave characteristics in front
More informationCFD AND EXPERIMENTAL STUDY OF AERODYNAMIC DEGRADATION OF ICED AIRFOILS
Colloquium FLUID DYNAMICS 2008 Institute of Thermomechanics AS CR, v.v.i., Prague, October 22-24, 2008 p.1 CFD AND EXPERIMENTAL STUDY OF AERODYNAMIC DEGRADATION OF ICED AIRFOILS Vladimír Horák 1, Dalibor
More informationSound scattering by hydrodynamic wakes of sea animals
ICES Journal of Marine Science, 53: 377 381. 1996 Sound scattering by hydrodynamic wakes of sea animals Dmitry A. Selivanovsky and Alexander B. Ezersky Selivanovsky, D. A. and Ezersky, A. B. 1996. Sound
More informationThorsten Mauritsen *, Gunilla Svensson Stockholm University, Stockholm, Sweden
J.1 WAVE FLOW SIMULATIONS OVER ARCTIC LEADS Thorsten Mauritsen *, Gunilla Svensson Stockholm University, Stockholm, Sweden Branko Grisogono Department of Geophysics, Faculty of Science, Zagreb, Croatia
More informationCRITERIA OF BOW-DIVING PHENOMENA FOR PLANING CRAFT
531 CRITERIA OF BOW-DIVING PHENOMENA FOR PLANING CRAFT Toru KATAYAMA, Graduate School of Engineering, Osaka Prefecture University (Japan) Kentarou TAMURA, Universal Shipbuilding Corporation (Japan) Yoshiho
More informationMultiple sources of buoyancy in a naturally ventilated enclosure
J. Fluid Mech. (1996), vol. 311, pp. 177-192 Copyright 0 1996 Cambridge University Press 177 Multiple sources of buoyancy in a naturally ventilated enclosure By P. F. LINDEN' AND P. COOPER' Department
More informationWind Flow Validation Summary
IBHS Research Center Validation of Wind Capabilities The Insurance Institute for Business & Home Safety (IBHS) Research Center full-scale test facility provides opportunities to simulate natural wind conditions
More informationVortical motions under short wind waves
Vortical motions under short wind waves X. Zhang and C. S. Cox Scripps Institution of Oceanography, UCSD, CA 9293-23, USA Abstract The flow structures under ocean wind waves are of great significance in
More informationCOMPUTATIONAL FLOW MODEL OF WESTFALL'S LEADING TAB FLOW CONDITIONER AGM-09-R-08 Rev. B. By Kimbal A. Hall, PE
COMPUTATIONAL FLOW MODEL OF WESTFALL'S LEADING TAB FLOW CONDITIONER AGM-09-R-08 Rev. B By Kimbal A. Hall, PE Submitted to: WESTFALL MANUFACTURING COMPANY September 2009 ALDEN RESEARCH LABORATORY, INC.
More informationAtmospheric Pressure. Conceptual Physics 11 th Edition. Atmospheric Pressure. Atmospheric Pressure. The Atmosphere
Atmospheric Pressure Conceptual Physics 11 th Edition Chapter 14: GASES Atmospheric pressure Caused by weight of air Varies from one locality to another Not uniform Measurements are used to predict weather
More informationVisual Observation of Nucleate Boiling and Sliding Phenomena of Boiling Bubbles on a Horizontal Tube Heater
Proceedings of the 2 nd World Congress on Mechanical, Chemical, and Material Engineering (MCM'16) Budapest, Hungary August 22 23, 216 Paper No. HTFF 146 DOI:.11159/htff16.146 Visual Observation of Nucleate
More informationConditions for Offshore Wind Energy Use
Carl von Ossietzky Universität Oldenburg Institute of Physics Energy Meteorology Group Detlev Heinemann Conditions for Offshore Wind Energy Use Detlev Heinemann ForWind Carl von Ossietzky Universität Oldenburg
More informationPressure coefficient on flat roofs of rectangular buildings
Pressure coefficient on flat roofs of rectangular buildings T. Lipecki 1 1 Faculty of Civil Engineering and Architecture, Lublin University of Technology, Poland. t.lipecki@pollub.pl Abstract The paper
More informationEXPERIMENTAL RESEARCH ON COEFFICIENT OF WAVE TRANSMISSION THROUGH IMMERSED VERTICAL BARRIER OF OPEN-TYPE BREAKWATER
EXPERIMENTAL RESEARCH ON COEFFICIENT OF WAVE TRANSMISSION THROUGH IMMERSED VERTICAL BARRIER OF OPEN-TYPE BREAKWATER Liehong Ju 1, Peng Li,Ji hua Yang 3 Extensive researches have been done for the interaction
More informationMET Lecture 8 Atmospheric Stability
MET 4300 Lecture 8 Atmospheric Stability Stability Concept Stable: Ball returns to original position Neutral: Ball stays wherever it is placed Unstable: Displacement grows with time. Atmospheric Stability
More informationWhat is Air Temperature? Temperature, Buoyancy, and Vertical Motion. How Atmospehric Temperature is Measured. Temperature Scales
Temperature, Buoyancy, and Vertical Motion Temperature, Pressure, and Density Buoyancy and Static Stability Temperature Lapse Rates Rising & Falling Motions in the Air What is Air Temperature? Temperature
More informationGoal: Develop quantitative understanding of ENSO genesis, evolution, and impacts
The Delayed Oscillator Zebiak and Cane (1987) Model Other Theories Theory of ENSO teleconnections Goal: Develop quantitative understanding of ENSO genesis, evolution, and impacts The delayed oscillator
More informationInteracting Turbulent Plumes in a Naturally Ventilated Enclosure
International Journal of Ventilation ISSN 473-33 Volume 4 No 4 Interacting Turbulent Plumes in a Naturally Ventilated Enclosure P. F. Linden and N. B. Kaye Department of Mechanical and erospace Engineering,
More informationHOW FAST/FAR DOES FLY LINE FALL? N. Perkins of the University of Michigan, March 2003
HOW FAST/FAR DOES FLY LINE FALL? N. Perkins of the University of Michigan, March 003 This report summarizes a simple model for the free fall dynamics of a length of fly line. The line is assumed to remain
More informationAbrupt marine boundary layer changes revealed by airborne in situ and lidar measurements
Abrupt marine boundary layer changes revealed by airborne in situ and lidar measurements David A. Rahn 1, Thomas R. Parish 2, and David Leon 2 1 Univeristy of Kansas 2 Univeristy of Wyoming Precision Atmospheric
More informationPlane Turbulent Wall Jets in Limited Tailwater Depth
International Journal of Engineering & Technology IJET-IJENS Vol: 11 No: 6 192 Plane Turbulent Wall Jets in Limited Tailwater Depth Shazy A. Shabayek 1 Abstract This paper presents laboratory study of
More informationMODELING INDIAN OCEAN CIRCULATION: BAY OF BENGAL FRESH PLUME AND ARABIAN SEA MINI WARM POOL
MODELING INDIAN OCEAN CIRCULATION: BAY OF BENGAL FRESH PLUME AND ARABIAN SEA MINI WARM POOL P. N. Vinayachandran* 1 1, *2 and J. Kurian* * 1 Centre for Atmospheric and Oceanic Sciences, Indian Institute
More informationSection 1. Global Wind Patterns and Weather. What Do You See? Think About It. Investigate. Learning Outcomes
Chapter 5 Winds, Oceans, Weather, and Climate Section 1 Global Wind Patterns and Weather What Do You See? Learning Outcomes In this section, you will Determine the effects of Earth s rotation and the uneven
More informationANALYSIS OF HEAT TRANSFER THROUGH EXTERNAL FINS USING CFD TOOL
ANALYSIS OF HEAT TRANSFER THROUGH EXTERNAL FINS USING CFD TOOL B. Usha Rani 1 and M.E Thermal 2 1,2 Asst.Professor, Dadi Institute of Engineering and Technology, India Abstract-The rate of heat transfer
More informationMethods The experiments were carried out in the wind tunnel in the Air physics laboratory at the Engineering Centre Bygholm as shown in figure 1.
1 Standardisation of test method for salt spreader: Air flow experiments Report 2: Visualization of airflow patterns by Jan S. Strøm, Consultant Aarhus University, Engineering Centre Bygholm, Test and
More information