Turbulence in Fluids Plumes and Thermals enoi Cushman-Roisin Thayer School of Engineering Darmouh College Why do hese srucures behave he way hey do? How much mixing do hey accomplish? 1
Plumes Plumes are common feaures in environmenal fluids. They occur whenever a persisen source of buoyancy creaes a rising moion of he buoyan fluid upward and away from he source. The cleares example is ha of hydrohermal vens a he boom of he ocean. hp://www.livescience.com/ www.panoramio.com/phoo/8607 Phoo by by Dalibor Ribičić Anoher occurrence is he rising of freshwaer from he boom of he sea a submarine springs in karsic regions such as along he Dalmaian Coas of Croaia, where such feaures are called vrulje. The common urban smokesack plume is, however, somewha differen because he warm gas rises no only under is own buoyancy bu also under he propulsion of momenum (ineria). Such plume is more properly caegoried as a buoyan je or forced plume. uoyancy imporan a laer sages Je momenum imporan a he beginning Wha drives a plume is is hea flux, defined as he amoun of hea (expressed in joules) being discharged hrough he exi hole per uni ime. ecause i is more pracical in laer mahemaical developmens, his quaniy is divided by 0 C p (he fluid's reference densiy imes is hea capaciy a consan pressure) and hen muliplied by (he fluid's hermal expansion coefficien) and he graviaional acceleraion g, giving rise o he buoyancy flux: g hea F 0C p ime Noe ha because hea per ime is expressed in J/s, he buoyancy flux is measured in m 4 /s.
Le us consider a hree-dimensional radially symmeric plume progressing verically from he boom hrough a homogeneous and resing fluid. If we denoe by T 0 he emperaure of he ambien fluid, hen he emperaure inside he plume has he value T 0 + T', in which T' denoes he emperaure anomaly (posiive in a rising plume, negaive in a sinking plume). To his emperaure anomaly corresponds a densiy anomaly ' = 0 T'. I is convenien o define he local buoyancy, or reduced graviy, g' as: g g gt 0 The buoyancy flux F can be expressed as he inegral across he secion of he plume of he produc of he verical velociy w wih he buoyancy g'. In erms of he mean values across he plume, a reasonably good approximaion is: F w g R gt w R 1) Hea budge along he plume reduces o F = consan = upsream value. We need now o esablish he mass and momenum budges. For his, consider an infiniesimal slice along he plume. ) Mass budge: Mass exiing a op = Mass enering a boom + Mass enrained hrough he side d ( R w) R u d ) Verical-momenum budge: Verical momenum exiing a op = verical momenum enering a boom + Verical momenum enrained hrough side + Upward buoyancy force = 0 d ( R w ) R g d 4) Closure assumpion: Enrainmen velociy proporional o verical velociy inside plume u a w wih coefficien a o be deermined empirically
The problem consiss of 4 equaions for 4 unknowns, each a funcion of he verical coordinae : R, w, u and g' The soluion is self-similar: R 0.157 F w.14 F u 0.8 F g 6.08 1/ 1/ 1/ 1/ / 5/ Coefficien was fied o observaions (cone of x 8.9 o ) Comparison wih he urbulen je: angle = 11.8 o R = x/5 = 0.0 x A plume is narrower. The fi of he cone angle ses he value of he coefficien a in he enrainmen assumpion: a 0.10 u 0. 10 w Comparing his coefficien o ha for a urbulen je (= 0.10), we deduce ha mixing is slighly more vigorous in a plume han in a je. Pu anoher way, buoyancy is a more effecive mixer han momenum. Plumes in a Sraified Environmen When a plume rises (or sinks) in a sraified environmen, i encouners a emperaure becoming closer o is own and progressively loses buoyancy. A some level, i will have los all buoyancy and will begin o spread horionally. Such is he case of a smokesack plume in a calm (no wind) and sraified amosphere ypical of he early morning. The obvious quesion o ask is how high does he plume reach? 4
Thermals A hermal is a finie parcel of fluid consising of he same fluid as is surroundings bu a a differen emperaure. ecause of is buoyancy, a cold hermal sinks (negaive buoyancy), while a warm hermal rises (posiive buoyancy). The name was given by glider pilos o wha hey perceived as regions of warm air rising above a heaed ground in which hey could soar. Convecion in he amosphere does indeed proceed by means of rising hermals. The siuaion, however, can be quie chaoic, wih a collecion of hermals rising here and here a various imes, some of hem smaller and slower, and ohers larger and faser. Here, for he sake of undersanding he basic mechanism, we are be concerned wih a single hermal immersed in an infinie homogeneous fluid a res. Experimens have been conduced in he laboraory, and i has been found ha all hermals roughly behave in similar ways. As hey rise (or sink), hey enrain surrounding fluid and become more dilue, hereby slowing down in heir ascen (or descen). The acual shape of a hermal, however, can vary considerably from one se of observaions o anoher. Here, basic dynamics supplemened by a few dimensionless numbers gleaned from experimens will be used o esablish a simple heory for he predicion of a hermal's behavior over ime. 5
Descending hermals in a laboraory experimen. These hermals in waer are made visible by barium sulfae. The sem lef behind by each hermal is due o he manner a spherical cap was roaed o provoke he release. The second hermal (boom row) has a larger negaive buoyancy han he firs (op row). [From Scorer, 1997] Define he hermal s average radius R and volume V, hen relae he wo by: V m R in which m is a shape facor (would be = 4/ = 4.19 if he hermal were a sphere). 1) Mass budge: dv dv Au a R w d d in which A is he enclosing surface area, u he enrainmen velociy, and a a dimensionless consan ) Momenum budge: d d hermal V w hermal added mass effec gv d d V w gv in which g is he reduced graviy experienced by he hermal ) Hea budge: gv consan in which is he oal buoyancy of he hermal 6
Soluion: 4a R 9m 64a V 79m 1/ 9m w 4a 1/ 79m g 64a / 4 / / Then, obain elevaion by inegraing velociy over ime: d w d 6m a 1/ and find proporionaliy beween radius R and heigh : R a m Experimenal daa reveal an angle of 14 o. Thus, R o an14 0.5 a 0.75 m Oher experimenal evidence: 5.80 1/ 6m 5.80 a m a 0.94 Togeher wih a = 0.75 m, we obain: a = 1.90 and m =.54. R 0.60 V 0.55.41 w 1.0 g 1.81 1/ / 4 / 1/ 1/ / Plo of he quaniy / (a ime consan during he life of hermal) versus he square roo of he hermal's buoyancy. Each numbered do refers o a differen laboraory experimen, and he solid line shows he bes linear fi. [From Scorer, 1997] 7
Conclusions Jes, plumes and hermals are complex urbulen srucures, bu heir dynamical behavior can be modeled by a few relaively simple equaions, in which a dimensionless facor or wo can laer be adjused by laboraory experimenaion or field daa. Jes, plumes and hermals accomplish mixing, which is easily deermined from he soluion o he dynamical equaions. Mixing is accomplished by enrainmen of ambien fluid ino he coheren srucure. The (laeral) enrainmen velociy is on he order of 10% of he average speed of he fluid wihin he srucure. No one has ye proposed heories ha predic he numerical values of - he angles relaing sie o disance, - he raio of enrainmen velociy o mean velociy. 8