Notes on the Turbulence Closure Model for Atmospheric Boundary Layers
|
|
- Julian Mosley
- 6 years ago
- Views:
Transcription
1 466 Journal of the Meteorological Society of Japan Vol. 56, No. 5 Notes on the Turbulence Closure Model for Atmospheric Boundary Layers By Kanzaburo Gambo Geophysical Institute, Tokyo University (Manuscript received 24 March 1978, in revised form 12 July 1978) Abstract A turbulence closure model for atmospheric boundary layers is examined, under the assumption that the turbulent flow is steady in its ensemble average and the advection and diffusion terms in the turbulent Reynolds stress and heat flux equations are neglected. The constants which are introduced in order to obtain a closure system are determined referring the experimental results obtained in the wind tunnel. The validity of the closure model thus obtained is checked referring the observational results obtained in the constant-flux layer. In order to apply our model to the planetary boundary layer, simple forms for the eddy transport coefficients of momentum and heat are formulated. In this case the parametarization of the scale-length l=l(z) (z: heihgt) which is introduced in estimating the dissipation rate of turbulent kinetic energy with height is discussed in order to satisfy the observational result. 1. Introduction Recently the so-called closure models which use exact equations for the mean field and approximate ones for the turbulence have been successfully applied to predict the detailed structure of the planetary boundary layer (for example, Mellor and Yamada (1974), Yamada and Mellor (1975), Wyngaad, Cote and Rao (1974), Deardorff (1974a, 1974b), Wyngaad (1975), Schemm and Lipps (1976), Sommeria (1976)). In the closure model, however, there is the choice how to select the empirical constants which are used to represent the approximate equations for the Reynolds stress and heat flux due to the turbulence. We examine the choice of the empirical constants in the closure model in this paper. In this note we proceed our discussion refering Mellor's (1973) paper, because he is the first to apply the closure models which were developed in the field of engineering to the prediction of the properties of the constant flux layer in the atmosphere. Mellor determined the empirical constants from the data under the neutral stability condition. Generally speaking, the estimation of the pressure-strain correlation in the Reynolds stress (uti uj) equation and the pressure-temperature correlation in the heat-flux (ui*) equation is essential in the success of the closure model. Therefore, the pressure-strain correlation in the Reynolds stress equation is re-examined in the first part of this paper, comparing the treatment proposed by Mellor. After that the pressuretemperature correlation is also re-examined. In Mellor's case, it is assumed that there is no direct influence of gravity on the pressure-containing correlation in the equations for the Reynolds stress and heat-flux. The present note argues that such an influence should be included, refering Launder's (1975) proposal. In the latter part of this paper, our model is extended to predict the properties of the planetary boundary layer and we present the parametarized scheme for the eddy coefficients of momentum and heat, KM and KH, introducing the scale-length l. In the constant flux layer, the scale-length l is l*kz (k: Karman's constant) outside the roughness sublayer. In the planetary boundary layer, however, the scale-length l is l(z). Therefore the stipulation of the functional form of l(z) should be considered in order to satisfy the observational data such as w2 (vertical velocity variance) and *2 (potential temperature variance).
2 October 1978 Kanzaburo Gambo The basic equations The equations of motion for the mean velocity and mean potential temperature may be written i n the form they are small in the treatment of planetary boundary layer. On the other hand, the equation of *2 may be written in the form U* and u* are mean and fluctuating velocity component (i, j, k run from 1 to 3, *1*, 2*y and *3*z in (x, y, z)-coordinate * system), and * are mean and fluctuating potential * temperature and P is the mean pressure. The overbars represent ensemble averages. The Boussinesq approximation and the adiabatic motion have been made here, and *0 and. *0 are constant representative values of density and potential temperature respectively. *ij* is the alternating unit tensor, fi= (0, 0, f) is the Coriolis parameter, g is the gravitational acceleration, and * is the Kronecker delta. Molecular diffusive terms have been neglected. The equations govering the transport of the Reynolds stress u* uj and heat flux ui * may be written in the form 3. Approximation of the pressure-strain correlations. The equations govering the Reynolds stress and heat flux contain unknown correlations of fluctuating quantities such as shown in (2.3), (2.4) and (2.5). These must be approximated in terms of the main dependent variables if we wish to obtain the closed system for the govering equations (2.1)*(2.5). In this section we consider the case there is no buoyant force for the sake of simplicity and discuss the parameterization of the pressure-strain correlation in (2.3), Concerning this problem many proposals have been discussed and each proposal seems to be based on the different assumption. Since the detailed review of individual proposal is not main problem of this paper, we omitt the reference of these proposals. However, it might be generally accepted that the principle of all proposals is based on the parameterization of the pressure-strain correlation proposed by Rotta (1951). Therefore we pick up here Rotta's original proposal in order to discuss the closed system for (2.3) systematically. For the estimation of pressure-strain correlation, Rotta proposed the following form for (p/*0)(*ui/*j) in the homogeneous turbulence: p is the pressure fluctuation about mean value, is the kinematic viscosity and * is the kinematic * thermal diffusivity. The terms which contain the Coriolis term have been neglected here because Here C1 is a constant and E and * are the timeaveraged turbulent kinetic energy and energy dissipation rate respectively, i.e.,
3 468 Journal of the Meteorological Society of Japan Vol. 56, No. 5 in a case of isotropic turbulence in (2.3). The quotient E/* thus represents a characteristic decay time of the turbulence. The second term on the right hand side of (3.1), *ij2 arises from the presence of mean rate of strain and the coefficient aljmi is approximated by determined to satisfy the experimental results of the nearly homogeneous shear flow obtained by Champagne et al. (1970), Ui= {U, 0, 0} and U varies linearly with *3 (*3*z). And the most reasonable constants are given by the *'s are the Cartesian components of the position vector X-Y and term with and without a prime relate to values at Y and X respectively (the integration being carried out over Y space). Rotta pointed out that the fourthorder tensor {aljmi} should satisfy the following kinematic constrains: In a special case Ui*{U, 0, 0}, U= U(*3), and the homogeneous turbulence is assumed, Rotta computed approximately the values of a1131, a1232, * as follows: provided that C1 takes a value of about 1.5*. Considering that C2*0.4 in (3.9), it is easily understood that the first term on the right-hand side of (3.9) predominates compared with other terms. Then Launder et al. proposed the approximate expression for (3.9) as follows: * is a constant. It is expected that the constant * will differ somewhat in magnitude from the coefficient of the first term in (3.9) (*8.4/11 for C2*0.4) to compensate in part for the neglected terms. In a case Ui= {U, 0, 0} and U=U(*3), we have P=-u1u3(*U/*3) in (3.12) and then Later Hanjalic and Launder (1972), and Launder, Reece and Rodi (1975) extended Rotta's treatment in a general form. Since the forms of (3.5)*(3.7) suggest that aljmi might be satisfactorily approximated by a linear combination of Reynolds stress, Launder et al. proposed the following form for *ij2 in the homogeneous shear flow: Comparing these values with those in (3.8) obtained by Rotta, the appropriate value of * will be given by 0.6 * In this case (*11)2, (*22)2 and (*33)2 give the same result as Rotta's estimation. However there is a little difference in the estimation of [(*13)2+ (*31)2], because Rotta's estimation is given by and P denotes the rate of production of turbulence energy. C2 in (3.9) is considered to be a constant. The constants C1 and C2 in (3.2) and (3.9) were Considering the discussion mentioned above, we use the equation (3.12) in the present paper for estimation of the pressure-strain correlations. Now we consider to determine * from the another point of view. If we neglect the transport and buoyancy terms in (2.3), incorporate the pressure-strain approximation such as shown in (3.12) and the dissipation term is expressed by 2/3(*ij), the equation (2.3) reduces to a set of following equations for the normal stresses: * Rotta originally proposed that C1 should be about 1.4.
4 October 1978 Kanzaburo Gambo 469 may be written in the form Here we use the relation of *=-u1u3(*u1/*3) and the values denoted by parentheses of the last terms in the above equations give the experimental result of Champagne et al. Under the assumption of C1=1.5 given by Launder et al., we have from the above relations between ui2 and * Thus the pressure fluctuation due to the effect of temperature fluctuations is expected in the case of buoyant flow and Launder proposed the pressure-strain interaction as follows: Considering the above estimation of *, we use * =0.65 and C1=1.5 tentatively in this paper. As pointed by Launder et al., the equation (3.9) is reduced, irrespective of the value of C2, to the following equation: provided that uiuj=2/3 (E*ij). This equation is the same as that obtained by Crow (1968) for the case of isotropic turbulence subjected to sudden distortion. If *=0.6 and uiuj=2/3(e*ij), the equation (3.12) is also reduced to the same equation as (3.15). In this sense, the expression such as (3.15) seems to be a special form for the estimation of (*ij2+*ji2). Here we note that the equation similar to (3.15) is usually used in the treatment of turbulence in the planetary boundary layer (Mellor, Mellor and Yamada, Deardorff (1975), Som- (1974a), Wyngaad et al., Yamada meria, Schemm and Lipps). 4. Analysis in buoyant shear flows In the preceeding section we discussed the approximate expression of the pressure-strain interaction for the case of non-buoyant flows. In this section we consider the extention of the above discussion to the case of buoyant flow. The main part of this section is due to Launder's proposal and some part of his paper is referred here, because we are particularly concerned with the process how the empirical constants are reasonably introduced in the parameterization of the pressure-strain and pressure-temperature correlations in (2.3) and (2.4). The divergence of the momentum equation for ui under the assumption of incompressible fluid Here Pij and P are considered to stand for the total generation of uiuj and E due to the combined effects of shear and buoyancy. The value of C2 is put as C2=*. Consistently, the same consideration is adopted in approximating the pressure scrambling term, (p/*0)(*/ *i) in (2.4). And we put C1T and C2T are constants, and In the usual closure mode, the second term of the right-hand side in (4.5) is neglected, i.e., C2T =0. However, we assume that C2T*0 in parallel with (4.2). Next we consider the dissipation term in (2.3), (2.4) and (2.5). Assuming isotropy of the fine scale motion, the dissipation term in each equation may be written CT' is a proportional constant. In esti-
5 470 Journal of the Meteorological Society of Japan Vol. 56, No. 5 mating CT', Launder used the following relation: which is derived from (2.5) and (4.9) in equilibrium flows, neglecting the diffusion term in (2.5). From the decay characteristics of ui2(=2e) and * 2 behind the grid, CT' is estimated by Launder as and CT'*1.6 (4.11) Here we note some comment about CT'. We have from (4.7) and (4.9) The hydrostatic relation is used in (5.3). Now we consider to rewrite the equations (2.3), (2.4) and (2.5), referring the analysis in section 4. In section 4, we use E/ * as the parameter which represents a characteristic decay time of Hinze (1975) estimated as CT'=4/3*1.33 under turbulence. In the so-called closure model, the the assumption of the "linear" velocity decay. prognostic equations for * is usually used to close On the other hand, Mellor assumed that CT'= the system which we use. In this paper, however, 16/15*1.07*, referring the boundary layer measurements of Johnson (1956). In this sense, there we assume the following simple relation between and E, * remains uncertain factor in determing the value of CT'.** In this paper, however, we use tentatively the value of CT' such as CT'=1.6. l is the characteristic scale length and C Concerning the estimation of C1T and C2T in is a constant. Then we have (4.5) Launder used the following relations Comparing these relations with Webster's (1964) experiment (for example, 1/CT'C1T*0.2), Launder estimated as C1T*3.2, C2T*0.5 (4.15) 5. Govering equations in the boundary layer If we apply the discussion mentioned in sections 2, 3 and 4 to a boundary layer the vertical scale is much less than the horizontal scale, equations (2.1) and (2.2) may be written as follows: and the equations (2.3), (2.4) and (2.5) may be written in the form * Following Mellor's notation we have CT'/2= B2/B1=8/15 in (4.12). ** The value of CT' is sensitive for the estimation
6 October 1978 Kanzaburo Gambo The constant flux layer If the mean velocity field is assumed as Ui= {U, 0, 0) and U=U(z), the Reynolds stresses and heat fluxes in the constant flux layer may be put in the following way: As mentioned in sections 3 and 4, the parameters C1, C2, CT', C1T and C2T are given in this paper as follows: In the constant flux layer, we assume further that the diffusion terms in (5.8), (5.9) and (5.10) are neglected. Then the equations for u2=u12, *2= u22, w2=u32 and uw=u1u3 in (5.8), u*=u1* and w*=u3* in (5.9) and *2 in (5.10) may be rewritten in the nondimensional form as follows: In order to fix the value of C in (5.8), (5.9) and (5.10), we use the equation (5.8), assuming that /*t=0, Ui= {U, 0, 0}, C1T=C2T=0 (no bouyant * flow) and the diffusion term is neglected. Since we assume no buoyance, we have From the equation for u1u3(=uw) in (5.8) we have On the other hand we have from (3.14) Elimination of w2 and E from (5.14)*(5.16) gives the following relation between uw and U/*z: * If l=kz (k: Karman's constant) and G(C)=1, the equation (5.17) corresponds to the formula, given by Prandtl's mixing length theory. Thus the value of C is given by G(C)=1, i.e.,
7 472 Journal of the Meteorological Society of Japan Vol. 56, No. 5 As mentioned in section 5 (cf. (5.17) and (5.19)), the value of C is determined under the assumption that l is proportional to z outside the roughness sublayer and we put = kz l (6.18) k is the Karman's constant. Now we consider to obtain the value of *M and *H as the function of *, eliminating u*2, *2 *, w*2, u* **, **2 from (6.4)*(6.10). By making the addition of (6.4), (6.5) and (6.6), we have the energy equation such as M =*+CE*3/2 * (6.19) Therefore we have in a strongly unstable case (-*>>1) provided that C=0.203 as mentioned in section 5. If we eliminate **, w* and *M from (6.6), (6.9), (6.10) and (6.19), we have Fig. 1 Calculated values of *M (solid curve). Dotted curve is obtained by Mellor. The equation for *M is obtained by eliminating u*** and w* from (6.6), (6.7) and (6.8) as follows: Eliminating E* from (6.19), (6.21) and (6.23), we obtain *M(*) and *H(*) such as shown in * The value of CT' is very sensitive for the estimation of *H on the stable side. A variation, CT'/CT'=0.2 products a variation *H/*H*0.3. * Fig. 2 Calculated values of *H (solid curve). Dotted curve is obtained by Mellor. Figs. 1 and 2.* In order to compare our results with those obtained by Mellor we show his results by dotted line in Figs. 1 and 2. In Mellor's case the type of (3.15) is used to estimate the pressurestrain correlation and the assumption that C2T= 0 in (4.5) is used. In Figs. 3 and 4, the values of *M and *H on the unstable side (*<0) are shown by the contineous line, while Mellor's case corresponds to the dotted line. For the sake of comparison, the empirical formulas such as *M4-9*M3=1 (KEYPS) and *M=(1-15*)-1/4 (Bussinger et al. (1971)) are shown by dot-chain lines in Fig. 3. The observed data plotted in Fig. 4 are quoted from the paper of Bussinger et al.
8 October 1978 Kanzaburo Gambo 473 Fig. 3 Comparison of the calculated values of *M (solid curve) with interpolation formulas (dot-chain curve) under unstable conditions (*<0). Dotted curve is obtained by Mellor. Fig. 5 Comparison of vertical velocity variance *W*2 with data assembled] by Wyngaad et al. (1971). Dotted line is obtained by Mellor. Fig. 4 Comparison of the calculated values of *H (solid curve) with interpolation formula (dot-chain curve) and observed data (from Businger et al., 1971) under unstable conditions (*<0). Dotted curve is obtained by Mellor. (1971a) and the dot-chain line in Fig. 4 corresponds to the empirical formula *H=0.74 (1-9*)-1/2 (Businger et al.). The values of *w*z and **2 are also shown in Figs. 5 and 6, refering the Mellor's case by dotted line. The observed data in these figures are quoted from the paper of Wyngaad et al. (1971). As may be seen in Figs. 3*6, our results show the good agreement with the observational data in unstable conditions (*<0). It may be suggested that the treatment of the pressure-strain correlation such as (4.2) plays an important role to obtain the reasonal values of *H(*) and *M(*) compared with those obtained by Mellor. In stable conditions (*>0), there is discrepancy about the values of *H and *M between Mellor's result and ours. Concerning the value of *H we Fig. 6 Comparison of temperature variance with data assembled by Wyngaad **2 et al. (1971). Dotted line is obtained by Mellor. notice that *H is sensitive with respect to CT' in stable conditions and *M is not sensitive with the value of CT'. As mentioned in section 4 there is uncertain factor in determining the value of CT'. If we put CT'=1.2 instead of CT'=1.6, the value of *H on the stable side in Fig. 2 gives the almost same value as that shown by dotted line in the figure. Concerning the discrepancy of *M in Mellor's case and ours, we must mention that there is no unique empirical formula of *M in stable conditions (*>0) to compare the theoretical result. In the discussion of *M and *H mentioned above, we have only refered the empirical formulas proposed by Businger et al, as one
9 474 Journal of the Meteorological Society of Japan Vol. 56, No. 5 example. Concerning the empirical formulas of M and *H, however, Pruit et al. (1973) and * Kondo (1962, 1975) proposed the following ones: (b) stable case: Generally speaking, these empirical formulas including that of Businger et al. show the difference in considerable magnitude in stable conditions (*>0), though the difference is small in the unstable conditions (*<0). For example, we have for *=1 M*5.7, * *H*5.44 Businger et al. M*2.2, *HN*3.0 Pruit et al. * Yamamoto and Shimanuki (1966) and Yamamoto (1975) proposed the generalization of the KEYPS formula in order to explain the discrepancy mentioned above. In these circumstances, we do not discuss here in detail about the difference between the empirical formula and the theoretical result in stable conditions. Concerning the value of **2 on the stable side in Fig. 6, both Mellor's case and ours show that **2 increases with *, while the observation shows the slowly decrease of **2 with *. This comes from the fact that **2 is proportional to H/E*1/2 such as shown in (6.10) and the increase of *H is larger than that of E*1/2 with * respect to * (>0). Anyway the problem to be revised still remains in stable conditions. In order to understand the characteristic feature of *M and *H we compute the values of *M and *H when -*>>1 and *=0. In case -*>>1, we have from the energy equation such as shown in (6.19) CE*3/2*-* (6.24) Replacement of (6.24) to (6.21) and (6.23) gives the following relations: The coefficient * in (6.26) is about two times large compared with that given by Mellor, while the coefficient * in (6.26) is almost the same as Mellor's value.* If *=0, we have from (6.19) and (6.23) the following two relations: Thus we have Replacement of C1=1.5 and C2=0.65 reduces to E**2.90, *M*1.00 (6.30) The value of *H at *=0 is estimated by making use of the following turbulent Prandtl number Pr, which is derived from (6.21) and (6.23). In our case we have Pr=0.75, *H=0.75 (6.31) If we use the values of E*, *M and *H at *=0 mentioned above, we have at *=0 As mentioned above, the values of *M, *H, u*2, w*2, ** at *=0 (or z*0) are determined based on parameters of C1, C2, C1T, CT' * If C1=1.5, C2=0.65, C1T=3.2, C2T=0.5, CT' =1.6, we have *0.22 and *0.39. **In the neutral case the strong shear of mean flow exists we have u*2>**2>w*2. In this sense the application of our discussion is limited in the case the shear of mean flow is weak.
10 October 1978 Kanzaburo Gambo which are determined from the experiment of the homogeneous shear flow or the buoyant flow in the wind tunnel. On the other hand Wyngaad et al. (1974, 1975) proposed the similar equations for uiuj, uj* and *2 with ours, though the expression of *ij2 is the type of (3.15). However they used the observed boundary layer values of *M, *H, u*2, w*2, * at *=0 in order to determine the values of C1, C2, C1T, CT' In this sense, there is no guarantee that the govering equations for uiuj, Ui* and *2 proposed by Wyngaad et al. (1974) could be safely applied to the constant flux layer. According to numerical experiment of the planetary boundary layer, however, the result seems to be reasonable compared qualitatively with Deardorff's (1972) experiment. The comparison of their result with ours will be briefly explained in section Eddy coefficients far momentum and heat, KM and KH. The prognostic equations in the planetary boundary layer are given by (5.1)*(5.10) as mentioned in section 5. However, the prognostic equations for uiuj, ui* and *2 may be simplified if we assume that locally steady state for the quantities of uiuj, ui* and *2 is satisfied and their diffusive processes are neglected. Mellor and Yamada (1974) classified such a system as the level II in their terminology. In this case the equations for w2, uw, u*, w* and *2 may be written, refering (5.8), (5.9) and (5.10), as follows: Here we assume that Ui= {U, 0, 0} and U=U(z). The energy equation for E=(u2 +v2 +w2)/ 2 is given by After some algebraic manipulation the solutions of the set of equations (7.1) to (7.6) may be expressed as follows: KM and KH are the eddy coefficients for momentum and heat and Rf is the flux Richardson number, i.e., In our case we have Rf1*0.29, Rf2*0.33, Rf 3*0.45 Thus the critical flux Richardson number Rfc is
11 476 Journal of the Meteorological Society of Japan Vol. 56, No. 5 given by Rfc=0.29. Concerning the value of Rfc, Yamamoto pointed out that it depends on the functional form of *M and *H with respect to and showed that we may be able to expect * the non-existence of Rfc in a special case. In our case the value of Rfc=0.29 is close to those suggested by Businger et al. and Kamimal and Izumi (1965). Since the functional forms of KM and KH in (7.18) and (7.19) are complicated ones, Takano (1977) rewrite (7.18) and (7.19) in the following way Here we use the relation such as And a, b2(ri) and b2(ri) are computed refering the results obtained by Mellor and Yamada. These formulas for KM and KH are originally derived by Yamamoto and he finds the empirical constant values of b1=25, b2=40 and a=1.25, under the assumption that b1(ri)=const., b2(ri)= const.. In order to compare our result with those obtained by Takano and Yamamoto, we discuss the case -*>>1. When -*>>1, we have from (7.20) and (7.21) because b1(ri) and b2(ri)>>1, On the other hand Mellor and Yamada obtain the following approximate expression (refer Fig. 1 in their paper): If we compare (7.24) and (7.26) with (7.28) and (7.30) the values of KM and KH estimated by Mellor and Yamada are about two times large compared with ours for given value of Rf. If we use the empirical constants of a=1.25, b1=25 and b2=40 in (7.22) and (7.23) suggested by Yamamoto, we have These formulas correspond to (7.29) and (7.31). Our results of KM and KH are larger than the values estimated by Yamamoto about 10% and 25% respectively when -Ri>>1. 8. Scale-length l The eddy coefficients for momentum and heat, KM and KH are expressed with the parameter l as shown in (7.8) and (7.9). In the constant flux layer we have l*kz outside the roughness sublayer. If we consider to apply the eddy coefficients, KM and KH to the planetary boundary layer, we must stipulate the function l(z). In the present paper, we assume Blackadar's (1962) interpolation formula l=kz/1+kz/l0 (8.1) Here we use Rf*1.446 Ri. In our case we have from (7.16) and (7.17) Therefore the approximate expressions for KN and KH in (7.18) and (7.19) are given by which interpolate between two limits l*kz as z* 0 and l*l0 as z*. Various expressions for l0 are proposed. For example, l0=ug/ U9 is the geostrophic velocity or l0=u*/f in the neutral boundary layer. However it is shown that these representations are not satisfactory in the convective turbulence field (Deardorff, 1974a), In the present paper we propose the following form for l0:
12 October 1978 Kanzaburo Gambo 477 l0=h/n (8.2) Here we note that h is the height of convective layer and N is a constant to be determined. Roughly speaking, N corresponds to the value of 1/* which is defined by Mellor and Yamada (cf. (7.2) in their paper.). Now we consider to determine the appropriate value of N. Introducing the expression of (8.1) into (5.8), (5.9) and (5.10), refering the discussion of the constant flux layer in section 6, we have the following equations corresponding to the equations (6.4)*(6.10): The energy equation is written from (8.3)*(8.5) as follows: In a strongly unstable case (-*>>1), the equations (8.3)*(8.11) are simplified just as we computed the values of *M and *H in section 6. After same algebraic manipulation we have * is given by (6.27). On the other hand we have from (8.11) Then we have from (8.9) *= *CT'/C2/3.* Eliminating *M from (8.11) and (8.5), we have In order to compare with the observational data in the planetary boundary layer, we compute the following values: * In our case *l.02 and *2.83
13 478 Journal of the Meteorological Society of Japan Vol. 56, No. 5 w* is the convective velocity scale defined by Then we have k is the Karman's constant. Assuming that we compute W2 and *2 for N=5, 10 and 15 and the results are shown in Figs. 7 and 8. For the sake of comparison, the observational data in the figures are quoted from Deardorff's (1972) paper. In the figures the results of numerical experiment in a case of -h/l=45 computed by Deardorff (1972) are also shown by the dotted line, although his example corresponds to the case w*/w*0=(1-z/h)2. Comparing the values of W2 with the observational data in Fig. 7, the value of N seems to be around N=5*10 i.e., h/5*10. the differences In the case of *2 in Fig. 8, however, l0* of *2 due to the value of N is rela- Fig. 7 Mean profiles of vertical velocity variance, W2=w*2u*02/w*2 for N. Ordinate is z/h (h: height of convective layer). Observed data are quoted from Deardorff's (1972) paper. Dotted curve is obtained by Deardorff (1972) for h/l= Fig. 8 Mean profiles of potential temperature variance, *2=**2(W*2/u*02) for N. Ordinate is z/h (h: height of convective layer). Observed data are quoted from Deardorff's (1972) paper. Dotted curve is obtained by Deardorff (1972) for -h/l=4.5 tively small and it is difficult to determine the reasonable value of N from the observational data. Here we note that Wyngaad et al. (1974) obtained the similar results with those in Figs. 7 and 8. In their case the prognostic equation for is used instead of defining the scale-length and * they determined the constants for the parametarization of closure model from the observational data at *=0 as mentioned in section 6. Their results w*/w*0*(1-z/h) and -z/l>>1* are shown in Figs. 9 and 10 by dotted line, while our results in Figs. 7 and 8 are shown by contineous lines. As may be seen in Figs. 9 and 10, the case of N=5 seems to correspond to their result in case of W2. In case of *2, however, N= 15 seems to correspond to their result. If we assume that N=10 and h*1,000m, we have l=0.4z/( z), i.e., l*67m at z= 500m. For the sake of comparison, we refer the observational result of the Eulerian scale length e obtained by Gamo et al. (1975). They showed l that the vertical distribution of le is quite similar to (8.1) and le is about 100m at z=500m* 1,000m. Though there is some problem about the direct comparison of l with 1e, our assumption that N*10 may be reasonable. Since l is proportional factor in estimating the eddy coefficients of momentum and heat, the numerical simulation of the planetary boundary layer should be per- * They computed the cases of -h/l=10, 20 and 50. The results show that there is no difference in each case.
14 October 1978 Kanzaburo Gambo 479 *w*31.7*1/3*** (8.22) Gamo et al. (1976) showed that *w=c0*1/3 and the mean proportional factor C0 which was determined from the observations at the heights of 200m, 300m, 460m, 700m and 1,000m is about 28, i.e., *w*100cm/sec, when *=50cm2/ sec3. Strictly speaking, C0 is regarded as C0= 1/2[l(z)]1/3 such as shown in (8.21). However, * Gamo et al. (1976) estimated the mean value of C0 in the planetary boundary layer. Therefore we omit the further discussion about the quantitative difference between the coefficient in (8.22) and the value of C0. Fig. 9 Comparison of W2 for various values of N (solid line) in case -*>1 with that obtained by Wyngaad et al. (1974) in case -h/l>>1. Ordinate is z/h. Acknowledgement The author wishes to express his thanks to Prof. T. Matsuno and Dr. K. Takano for their kind discussion. Thanks are due to Mr. Y. Fujiki and Mrs. K. Kudo for drawing the figures and typing the draft. The author is grateful to the reviewers for their valuable comments on the original version of this paper. The research reported here was supported by Funds for Scientific Research from the Ministry of Education. Fig. 10 Comparison N of *2 for various values of N (solid line) in case -*>>1 with that obtained by Wyngaad et al. (1974) in case -h/l>>1. Ordinate is z/h. formed in the future in order to check the vertica distribution of l comparing the observationa data. As a short comment to (8.15), here we refer the observations of w2 and * by Gamo et al (1976). The equation (8.15) is rewritten as fol lows: Here we use the relation of *(g/*0)w* uncler the assumption of -*>>1. Putting *2.82* amd l=67m**, we have in unit of cm/sec * cf. footnote of (8.15) N=10 and h=1,000m ** are used in estimating the value of l. l*67m is the value at z=504m. References Blackadar, A. K., 1962: The vertical distribution of wind and turbulent exchange in neutral atmosphere. J. Geophy. Res., 67, Businger, J. A., J. C. Wyngaad, Y. Izumi and E. F. Bradley, 1971: Flux-profile relationships in the atmospheric surface layer. J. Atmos. Sci., 28, Champane, F. H., V. G. Harris and S. Corrsin, 1970: Experiments on nearly homogeneous turbulent shear flow. J. Fluid Mech., 41, Crow, S. C., 1968: Viscoelastic properties of finegained incompressible tu*bulence. J. Fluid Mech. 33, Deardorff, J. W., 1972: Numerical integration of neutral and unstable planetary boundary layers. J. Atmos. Sci., 29, a: Three-dimensional *, numerical study of the height and mean structure of a heated planetary boundary layer. Boundary- Layer Meteor., 7, b: *, Three-dimensional numerical study of turbulence in an entraining mixed layer. Boundary-Layer Meteor. 7, Gamo, M. T. and O. Yokoyama, 1975: Observed characteristics of the standard deviation of the vertical wind velocity in upper part of the atmospheric boundary layer. J. Meteor. Soc. Japan, 53, *** *w123cm/sec when =50cm2/sec3.
15 *,O 480 Journal of the Meteorological Society of Japan Vol. 56, No. 5 from a simplified three-dimensional numerical model of atmospheric turbulence. J. Atmos. Sci., 33, Sommeria, G., 176: Three-dimensional simulation of turbulent processes in an undisturbed tradewind boundary layer. J. Atmos. Sci., 33, Takano, K., 1977: Three-dimensional numerical modeling of the land and sea breezes and the urban heat island in the Kanto Plain. Sc. D. Thesis, Tokyo University (to be published). Webster, C. A. G., 1964: An experimental study of turbulence in a density stratified shear flow. J. Fluid Mech., 19, Wyngaad, J. C., O. R. Cote and Y. Izumi, 1971: Local free convection, similarity, and the budgets of shear stress and heat flux. J. Atmos., Sci., 28, *,* and K. S. Rao, 1974: Modeling the atmospheric boundary layer. Advances in Geophysics, 18A, : Modeling *, the planetary boundary layer extension to the stable case. Boundary- Layer Meteor. 9, Yamada, T., 1975: The critical Richardson number and the ratio of the eddy transport coefficient cbtained from a turbulent closure model. J. Atmos. Sci., 32, and G. *, L. Mellor, 1975: A numerical model prediction of the Wyangara atmospheric boundary layer data. J. Atmos. Sci., 32, Yamamoto, G., 1976: Generalization of the KEYPS formula in diabatic condition and related discussion on the critical Richardson number, J. Meteor. Soc. Japan, 53, A. Shimanuki, *, 1966: Turbulent transfer in diabatic conditions. J. Meteor. Soc. Japan,
IMPROVEMENT OF THE MELLOR YAMADA TURBULENCE CLOSURE MODEL BASED ON LARGE-EDDY SIMULATION DATA. 1. Introduction
IMPROVEMENT OF THE MELLOR YAMADA TURBULENCE CLOSURE MODEL BASED ON LARGE-EDDY SIMULATION DATA MIKIO NAKANISHI Japan Weather Association, Toshima, Tokyo 70-6055, Japan (Received in final form 3 October
More informationMellor-Yamada Level 2.5 Turbulence Closure in RAMS. Nick Parazoo AT 730 April 26, 2006
Mellor-Yamada Level 2.5 Turbulence Closure in RAMS Nick Parazoo AT 730 April 26, 2006 Overview Derive level 2.5 model from basic equations Review modifications of model for RAMS Assess sensitivity of vertical
More informationA Simple Turbulence Closure Model
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =
More informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
More informationA Simple Turbulence Closure Model. Atmospheric Sciences 6150
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationVALIDATION OF LES FOR LOCAL HEAT ENVIRONMENT IN TOKYO -COMPARISON WITH FIELD MEASUREMENT DATA-
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 29, Taipei, Taiwan VALIDATION OF LES FOR LOCAL HEAT ENVIRONMENT IN TOKYO -COMPARISON WITH FIELD MEASUREMENT DATA- Tsuyoshi Nozu 1,
More informationConvective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very
More informationLES of wind turbulence and heat environment around dense tall buildings
EACWE 5 Florence, Italy 19 th 23 rd July 2009 LES of wind turbulence and heat environment around dense tall buildings Flying Sphere image Museo Ideale L. Da Vinci Tsuyoshi Nozu 1, Takeshi Kishida 2, Tetsuro
More informationADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS
ADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS Radi Sadek 1, Lionel Soulhac 1, Fabien Brocheton 2 and Emmanuel Buisson 2 1 Laboratoire de Mécanique des Fluides et d Acoustique,
More informationThe applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization
Julia Palamarchuk Odessa State Environmental University, Ukraine The applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization The low-level katabatic
More informationTURBULENT KINETIC ENERGY
TURBULENT KINETIC ENERGY THE CLOSURE PROBLEM Prognostic Moment Equation Number Number of Ea. fg[i Q! Ilial.!.IokoQlI!!ol Ui au. First = at au.'u.' '_J_ ax j 3 6 ui'u/ au.'u.' a u.'u.'u k ' Second ' J =
More informationOn stably stratified homogeneous shear flows subjected to rotation
Center for Turbulence Research Proceedings of the Summer Program 2000 241 On stably stratified homogeneous shear flows subjected to rotation By B. A. Pettersson Reif, A.Ooi AND P. A. Durbin Theoretical
More information6A.3 Stably stratified boundary layer simulations with a non-local closure model
6A.3 Stably stratified boundary layer simulations with a non-local closure model N. M. Colonna, E. Ferrero*, Dipartimento di Scienze e Tecnologie Avanzate, University of Piemonte Orientale, Alessandria,
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationESS Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1
ESS5203.03 - Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1 Text: J.R.Garratt, The Atmospheric Boundary Layer, 1994. Cambridge Also some material from J.C. Kaimal and
More information(Wind profile) Chapter five. 5.1 The Nature of Airflow over the surface:
Chapter five (Wind profile) 5.1 The Nature of Airflow over the surface: The fluid moving over a level surface exerts a horizontal force on the surface in the direction of motion of the fluid, such a drag
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationτ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization
τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization Note that w *3 /z i is used to normalized the TKE equation in case of free
More informationRoughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationThe atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:
The atmospheric boundary layer: Utrecht Summer School on Physics of the Climate System Carleen Tijm-Reijmer IMAU The atmospheric boundary layer: Where the atmosphere meets the surface Photo: Mark Wolvenne:
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationOn the Velocity Gradient in Stably Stratified Sheared Flows. Part 2: Observations and Models
Boundary-Layer Meteorol (2010) 135:513 517 DOI 10.1007/s10546-010-9487-y RESEARCH NOTE On the Velocity Gradient in Stably Stratified Sheared Flows. Part 2: Observations and Models Rostislav D. Kouznetsov
More informationBuoyancy Fluxes in a Stratified Fluid
27 Buoyancy Fluxes in a Stratified Fluid G. N. Ivey, J. Imberger and J. R. Koseff Abstract Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed
More informationA Further Study of Dynamics of the Four-Day Circulation. By Yoshihisa Matsuda
February 1982 Y. Matsuda 245 A Further Study of Dynamics of the Four-Day Circulation in the Venus Atmosphere By Yoshihisa Matsuda Department of Astronomy and Earth Sciences, Tokyo Gakugei University, Koganei
More informationCOMMENTS ON "FLUX-GRADIENT RELATIONSHIP, SELF-CORRELATION AND INTERMITTENCY IN THE STABLE BOUNDARY LAYER" Zbigniew Sorbjan
COMMENTS ON "FLUX-GRADIENT RELATIONSHIP, SELF-CORRELATION AND INTERMITTENCY IN THE STABLE BOUNDARY LAYER" Zbigniew Sorbjan Department of Physics, Marquette University, Milwaukee, WI 5301, U.S.A. A comment
More informationA Note on the Barotropic Instability of the Tropical Easterly Current
April 1969 Tsuyoshi Nitta and M. Yanai 127 A Note on the Barotropic Instability of the Tropical Easterly Current By Tsuyoshi Nitta and M. Yanai Geophysical Institute, Tokyo University, Tokyo (Manuscript
More informationNumerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model
Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model Hitoshi Sugiyama 1),Naoto Kato 1), Masahiro Ouchi 1) 1) Graduate School of Engineering,Utsunomiya University
More informationA Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds over a Tropical Urban Terrain
Pure appl. geophys. 160 (2003) 395 404 0033 4553/03/020395 10 Ó Birkhäuser Verlag, Basel, 2003 Pure and Applied Geophysics A Note on the Estimation of Eddy Diffusivity and Dissipation Length in Low Winds
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationA Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation by Sullivan
耶鲁 - 南京信息工程大学大气环境中心 Yale-NUIST Center on Atmospheric Environment A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation
More informationImplementation of the Quasi-Normal Scale Elimination (QNSE) Model of Stably Stratified Turbulence in WRF
Implementation of the Quasi-ormal Scale Elimination (QSE) odel of Stably Stratified Turbulence in WRF Semion Sukoriansky (Ben-Gurion University of the egev Beer-Sheva, Israel) Implementation of the Quasi-ormal
More informationRica Mae Enriquez*, Robert L. Street, Francis L. Ludwig Stanford University, Stanford, CA. 0 = u x A u i. ij,lass. c 2 ( P ij. = A k. P = A ij.
P1.45 ASSESSMENT OF A COUPLED MOMENTUM AND PASSIVE SCALAR FLUX SUBGRID- SCALE TURBULENCE MODEL FOR LARGE-EDDY SIMULATION OF FLOW IN THE PLANETARY BOUNDARY LAYER Rica Mae Enriquez*, Robert L. Street, Francis
More informationAtm S 547 Boundary Layer Meteorology
Lecture 8. Parameterization of BL Turbulence I In this lecture Fundamental challenges and grid resolution constraints for BL parameterization Turbulence closure (e. g. first-order closure and TKE) parameterizations
More informationAtm S 547 Boundary Layer Meteorology
Lecture 5. The logarithmic sublayer and surface roughness In this lecture Similarity theory for the logarithmic sublayer. Characterization of different land and water surfaces for surface flux parameterization
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationPart III: Modeling atmospheric convective boundary layer (CBL) Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, USA
Physical modeling of atmospheric boundary layer flows Part III: Modeling atmospheric convective boundary layer (CBL) Outline Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, USA
More informationReynolds Averaging. Let u and v be two flow variables (which might or might not be velocity components), and suppose that. u t + x uv ( ) = S u,
! Revised January 23, 208 7:7 PM! Reynolds Averaging David Randall Introduction It is neither feasible nor desirable to consider in detail all of the small-scale fluctuations that occur in the atmosphere.
More informationStatistical turbulence modelling in the oceanic mixed layer
Statistical turbulence modelling in the oceanic mixed layer Hans Burchard hans.burchard@io-warnemuende.de Baltic Sea Research Institute Warnemünde, Germany Seminar at Bjerknes Centre, October 16, 23, Bergen,
More informationFrictional boundary layers
Chapter 12 Frictional boundary layers Turbulent convective boundary layers were an example of boundary layers generated by buoyancy forcing at a boundary. Another example of anisotropic, inhomogeneous
More informationOn the turbulent Prandtl number in homogeneous stably stratified turbulence
J. Fluid Mech. (2010), vol. 644, pp. 359 369. c Cambridge University Press 2010 doi:10.1017/s002211200999293x 359 On the turbulent Prandtl number in homogeneous stably stratified turbulence SUBHAS K. VENAYAGAMOORTHY
More informationModel description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club
Model description of AGCM5 of GFD-Dennou-Club edition SWAMP project, GFD-Dennou-Club Mar 01, 2006 AGCM5 of the GFD-DENNOU CLUB edition is a three-dimensional primitive system on a sphere (Swamp Project,
More informationSurface layer parameterization in WRF
Surface layer parameteriation in WRF Laura Bianco ATOC 7500: Mesoscale Meteorological Modeling Spring 008 Surface Boundary Layer: The atmospheric surface layer is the lowest part of the atmospheric boundary
More informationJournal of Fluid Science and Technology
Science and Technology Effect of Molecular Diffusivities on Countergradient Scalar Transfer in a Strong Stable Stratified Flow (Study on the Linear and Nonlinear Processes by using RDT) Kouji NAGATA, Takashi
More informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
More informationPhysical Processes in Acoustic Wave Heating of the Thermosphere
Publications 4-9-2005 Physical Processes in Acoustic Wave Heating of the Thermosphere G. Schubert Institute of Geophysics and Planetary Physics, University of California Michael P. Hickey Ph.D. Embry-Riddle
More informationLecture 12. The diurnal cycle and the nocturnal BL
Lecture 12. The diurnal cycle and the nocturnal BL Over flat land, under clear skies and with weak thermal advection, the atmospheric boundary layer undergoes a pronounced diurnal cycle. A schematic and
More informationA NONSTATIONARY NOCTURNAL DRAINAGE FLOW MODEL. 1. Introduction
A NONSTATIONARY NOCTURNAL DRAINAGE FLOW MODEL K. S. RAO and H. F. SNODGRASS* Atmospheric Turbulence and Diffusion Laboratory, National Oceanic and Atmospheric Administration, Oak Ridge, Tennessee 37830,
More informationLogarithmic velocity profile in the atmospheric (rough wall) boundary layer
Logarithmic velocity profile in the atmospheric (rough wall) boundary layer P =< u w > U z = u 2 U z ~ ε = u 3 /kz Mean velocity profile in the Atmospheric Boundary layer Experimentally it was found that
More informationLarge eddy simulation studies on convective atmospheric boundary layer
Large eddy simulation studies on convective atmospheric boundary layer Antti Hellsten & Sergej Zilitinkevich Finnish Meteorological Institute Outline Short introduction to atmospheric boundary layer (ABL)
More informationThe Johns Hopkins Turbulence Databases (JHTDB)
The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES
More informationThe Stable Boundary layer
The Stable Boundary layer the statistically stable or stratified regime occurs when surface is cooler than the air The stable BL forms at night over land (Nocturnal Boundary Layer) or when warm air travels
More informationPart I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows
Physical modeling of atmospheric boundary layer flows Part I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows Outline Evgeni Fedorovich School of
More informationTesting Turbulence Closure Models Against Oceanic Turbulence Measurements
Testing Turbulence Closure Models Against Oceanic Turbulence Measurements J. H. Trowbridge Woods Hole Oceanographic Institution Woods Hole, MA 02543 phone: 508-289-2296 fax: 508-457-2194 e-mail: jtrowbridge@whoi.edu
More informationSimulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
More informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More informationPassive Scalars in Stratified Turbulence
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Passive Scalars in Stratified Turbulence G. Brethouwer Linné Flow Centre, KTH Mechanics, SE-100 44 Stockholm, Sweden E. Lindborg Linné Flow Centre,
More informationP2.1 Scalar spectra in the near-dissipation range derived from structure functions
P2. Scalar spectra in the near-dissipation range derived from structure functions A.F. Moene and A. van Dijk,2 Meteorology and Air Quality Group, Wageningen University, The Netherlands 2 Alterra, Wageningen
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationPrototype Instabilities
Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies
More informationTransformed Eulerian Mean
Chapter 15 Transformed Eulerian Mean In the last few lectures we introduced some fundamental ideas on 1) the properties of turbulent flows in rotating stratified environments, like the ocean and the atmosphere,
More informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationLarge eddy simulation of a forced round turbulent buoyant plume in neutral surroundings
Center for Turbulence Research Annual Research Briefs 1999 239 Large eddy simulation of a forced round turbulent buoyant plume in neutral surroundings By A. J. Basu AND N. N. Mansour 1. Motivation and
More informationOn the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationWind driven mixing below the oceanic mixed layer
Wind driven mixing below the oceanic mixed layer Article Published Version Grant, A. L. M. and Belcher, S. (2011) Wind driven mixing below the oceanic mixed layer. Journal of Physical Oceanography, 41
More informationATMOSPHERIC AND OCEANIC FLUID DYNAMICS
ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-scale Circulation G E O F F R E Y K. V A L L I S Princeton University, New Jersey CAMBRIDGE UNIVERSITY PRESS An asterisk indicates more advanced
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
More informationCharacteristics of Langmuir Turbulence in the Ocean Mixed Layer
Characteristics of Langmuir Turbulence in the Ocean Mixed Layer Article Published Version Grant, A. L. M. and Belcher, S. E. (2009) Characteristics of Langmuir Turbulence in the Ocean Mixed Layer. Journal
More informationGrid-generated turbulence, drag, internal waves and mixing in stratified fluids
Grid-generated turbulence, drag, internal waves and mixing in stratified fluids Not all mixing is the same! Stuart Dalziel, Roland Higginson* & Joanne Holford Introduction DAMTP, University of Cambridge
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationPassive scalars in stratified turbulence
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L06809, doi:10.1029/2007gl032906, 2008 Passive scalars in stratified turbulence G. Brethouwer 1 and E. Lindborg 1 Received 5 December 2007; accepted 29 February 2008;
More informationKinematic Effects of Differential Transport on Mixing Efficiency in a Diffusively Stable, Turbulent Flow
Iowa State University From the SelectedWorks of Chris R. Rehmann January, 2003 Kinematic Effects of Differential Transport on Mixing Efficiency in a Diffusively Stable, Turbulent Flow P. Ryan Jackson,
More informationWind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting
Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting Detlev Heinemann ForWind Center for Wind Energy Research Energy Meteorology Unit, Oldenburg University Contents Model Physics
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationDirect Numerical Simulation of the Neutrally Stratified Turbulent Ekman Boundary Layer
Journal of the Earth Simulator, Volume 6, October 26, 3 15 Direct Numerical Simulation of the Neutrally Stratified Turbulent Ekman Boundary Layer Katsuhiro Miyashita 1, Kaoru Iwamoto 1 and Hiroshi Kawamura
More information3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures
3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures Igor ROGACHEVSKII and Nathan KLEEORIN Ben-Gurion University of the Negev, Beer-Sheva, Israel N.I. Lobachevsky State University
More informationSergej S. Zilitinkevich 1,2,3. Division of Atmospheric Sciences, University of Helsinki, Finland
Atmospheric Planetary Boundary Layers (ABLs / PBLs) in stable, neural and unstable stratification: scaling, data, analytical models and surface-flux algorithms Sergej S. Zilitinkevich 1,,3 1 Division of
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationThe parametrization of the planetary boundary layer May 1992
The parametrization of the planetary boundary layer May 99 By Anton Beljaars European Centre for Medium-Range Weather Forecasts Table of contents. Introduction. The planetary boundary layer. Importance
More informationNonsingular Implementation of the Mellor-Yamada. Level 2.5 Scheme in the NCEP Meso model
Nonsingular Implementation of the Mellor-Yamada Level 2.5 Scheme in the NCEP Meso model by Zaviša I. Janjić 1 UCAR Scientific Visitor December, 2001 National Centers for Environmental Prediction Office
More informationStable Boundary Layer Parameterization
Stable Boundary Layer Parameterization Sergej S. Zilitinkevich Meteorological Research, Finnish Meteorological Institute, Helsinki, Finland Atmospheric Sciences and Geophysics,, Finland Nansen Environmental
More informationLECTURE 28. The Planetary Boundary Layer
LECTURE 28 The Planetary Boundary Layer The planetary boundary layer (PBL) [also known as atmospheric boundary layer (ABL)] is the lower part of the atmosphere in which the flow is strongly influenced
More informationObservation of Falling Motion of Columnar Snow Crystals
276 Journal of the Meteorological Society of Japan Vol. 54, No. 5 Observation of Falling Motion of Columnar Snow Crystals By Masahiro Kajikawa Department of Earth Science, Akita University, Akita, Japan
More informationImproving a Turbulence Scheme for the Terra Incognita in a Dry Convective Boundary Layer
December Journal of the 2016 Meteorological Society of Japan, Vol. Y. 94, KITAMURA No. 6, pp. 491 506, 2016 491 DOI:10.2151/jmsj.2016-028 Improving a Turbulence Scheme for the Terra Incognita in a Dry
More informationBulk Boundary-Layer Model
Bulk Boundary-Layer Model David Randall Ball (1960) was the first to propose a model in which the interior of the planetary boundary layer (PBL) is well-mixed in the conservative variables, while the PBL
More informationReynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.
Reynolds Averaging Reynolds Averaging We separate the dynamical fields into sloly varying mean fields and rapidly varying turbulent components. Reynolds Averaging We separate the dynamical fields into
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationDynamic Meteorology (lecture 13, 2016)
Dynamic Meteorology (lecture 13, 2016) Topics Chapter 3, lecture notes: High frequency waves (no rota;on) Boussinesq approxima4on Normal mode analysis of the stability of hydrosta4c balance Buoyancy waves
More informationFall Colloquium on the Physics of Weather and Climate: Regional Weather Predictability and Modelling. 29 September - 10 October, 2008
1966-10 Fall Colloquium on the Physics of Weather and Climate: Regional Weather Predictability and Modelling 9 September - 10 October, 008 Physic of stable ABL and PBL? Possible improvements of their parameterizations
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationMultichoice (22 x 1 2 % = 11%)
EAS572, The Atmospheric Boundary Layer Final Exam 7 Dec. 2010 Professor: J.D. Wilson Time available: 150 mins Value: 35% Notes: Indices j = (1, 2, 3) are to be interpreted as denoting respectively the
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationLECTURE NOTES ON THE Planetary Boundary Layer
LECTURE NOTES ON THE Planetary Boundary Layer Chin-Hoh Moeng prepared for lectures given at the Department of Atmospheric Science, CSU in 1994 & 1998 and at the Department of Atmospheric and Oceanic Sciences,
More information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More information