首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 718 毫秒
1.
Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The turbulent fields are obtained directly by solving the three-dimensional time-dependent Navier-Stokes equations (using the Boussinesq approximation to account for buoyancy effects), and no turbulence model is needed. Two levels of heating are considered, one quite vigorous, the other more moderate. Statistics for the vigorously heated case are found to agree reasonably well with laboratory, field, and large-eddy simulation results, when Deardorff's mixed-layer scaling is used. No indication of large-scale longitudinal roll cells is found in this convection-dominated flow, for which the inversion height to Obukhov length scale ratio –z i /L *=26. However, when heating is more moderate (so that –z i /L *=2), evidence of coherent rolls is present. About 10% of the total turbulent kinetic energy and turbulent heat flux, and 20% of the Reynolds shear stress, are estimated to be a direct consequence of the observed cells.  相似文献   

2.
A general parameterization for solid and liquid hydrometeors is presented. hydrometeors basically are viewed as porous spheroids with the following variable parameters: diameter, axial ratio, mass, and porosity. Based on this parameterization a functional dependence on the Reynolds number of the drag of hydrometeors is derived, which is based on boundary layer theory. The basic form of this functional dependence is consistent with viscous theory and the inertial drag at low Reynolds numbers is predicted with good accuracy by matching the results from the boundary layer theory with Oseen's theory of creeping motion. Based on this solution a general semi-empirical expression for the Reynolds number and fall speed of particles is found. The results from the present theory are in remarkable agreement with experiments: The errors generally are < 5–10% for a wide variety of hydrometeors in the range of Reynolds numbers 0<NRe<5×105, including columnar and variously branched planar ice crystals, rimed and unrimed aggregates, lump, conical, and hexagonal graupel, hail, and rain drops. The present parameterisation aims far beyond the limits of the conventional methods since it is suitable for mixed-phase models of the microphysics of precipitation with continuously varying particle mass and shape characteristics and including processes such as depositional growth of ice crystals under varying environmental conditions, collisional growth of particles, and melting.  相似文献   

3.
The spray content in the surface boundary layer above an air—water interface was determined by a series of measurements at various feteches and wind speeds in a laboratory facility. The droplet flux density N(z) can be described in terms of the scaling flux density N* and von Karman constant K throguh the equation, N(z)/N* = −(1/K) ln(z/z0d) where z is height above the mean water level and z0d is the droplet boundary layer thickness. N* is given by a unique relationship in terms of the roughness Reynolds number u*σ/ν where σ is the root-mean-square surface displacement. Spray inception occurred for u* 0.3. The dominant mode of spray generation in the present and most other laboratory tests, as well as in available field data, appears to be bubble bursting.  相似文献   

4.
In order to investigate flows over topography in an atmospheric context, we have studied experimentally the wake structure of axi-symmetric Gaussian obstacles towed through a linearly stratified fluid. Three dimensionless parameters govern the flow dynamics: F, the Froude number based on the topography height h; Re, the Reynolds number and the aspect ratio r = h/L, where L is the topography horizontal scale. Two-dimensional (2-D), saturated lee wave (SLW) and three-dimensional (3-D) regimes, as defined in Chomaz et al. (1993), are found to be functions of F and r only (Fig. 1) as soon as Re is larger than Rec ≈ 2000. For F < 0.7 the flow goes around the obstacle and the motion in the wake is quasi-two-dimensional. This 2-D layer is topped by a region affected by lee wave motions with amplitude increasing with r and F. For 0.7 < F < 1/r, the flow is entirely dominated by a lee wave of saturated amplitude which suppresses the separation of the boundary layer from the obstacle. Above the critical value 1/r, the lee wave amplitude decreases with F and a recirculating zone appears behind the obstacle. Simultaneously, coherent large-scale vortices start to be shed periodically from the wake at a Strouhal number which decreases as 1/F until it reaches its neutral asymptotic value.  相似文献   

5.
The inertial subrange Kolmogorov constant C 0, which determines the effective turbulent diffusion in velocity space, plays an important role in the Lagrangian modelling of pollutants. A wide range of values of the constant are found in the literature, most of them determined at low Reynolds number and/or under different assumptions. Here we estimate the constant C 0 by tracking an ensemble of Lagrangian particles in a planetary boundary layer simulated with a large-eddy simulation model and analysing the Lagrangian velocity structure function in the inertial subrange. The advantage of this technique is that it easily allows Reynolds numbers to be achieved typical of convective turbulent flows. Our estimates of C 0 is C 0=4.3±0.3 consistent with values found in the literature  相似文献   

6.
A laboratory study of scalar diffusion in the convective boundary layer has found results that are consistent with a 1999 large-eddy simulation (LES) study by Jonker, Duynkerke and Cuijpers. For bottom-up and top-down scalars (introduced as ‘infinite’ area sources of passive tracer at the surface and inversion, respectively) the dominant length scale was found to be much larger than the length scale for density fluctuations, the latter being equal to the boundary-layer depth h. The variance of the normalized passive scalar grew continuously with time and its magnitude was about 3–5 times larger for the top-down case than for the bottom-up case. The vertical profiles of the normalized passive scalar variance were found to be approximately constant through the convective boundary layer (CBL) with a value of about 3–8c*2 for bottom-up and 10–50c*2 for top-down diffusion. Finally, there was some evidence of a minimum in the variance and dominant length scale for scalar flux ratios (top-down to bottom-up flux) close to −0.5. All these convection tank results confirm the LES results and support the hypothesis that there is a distinct difference in behaviour between the dynamic and passive variables in the CBL.  相似文献   

7.
Mean wind velocity profiles were measured by means of radio-windsondes over the Landes region in southwestern France, which consists primarily of pine forests with scattered villages and clearings with various crops. Analysis of neutral profiles indicated the existence of a logarithmic layer between approximately z – d 0 = 67(±18)z 0 and 128(+-32)z 0 (z is the height above the ground, z 0 the surface roughness and d 0 the displacement height). The upper limit can also be given as z – d 0 = 0.33 (±0.18)h, where h is the height of the bottom of the inversion. The profiles showed that the surface roughness of this terrain is around 1.2 m and the displacement height 6.0 m. Shear stresses derived from the profiles were in good agreement with those obtained just above the forest canopy at a nearby location with the eddy correlation method by a team from the Institute of Hydrology (Wallingford, England).  相似文献   

8.
The influence of an idealized moving wavy surface on the overlying airflow is investigated using direct numerical simulations (DNS). In the present simulations, the bulk Reynolds number is Re = 8000 (; where U0 is the forcing velocity of the flow, h the height of the domain and v the kinematic viscosity) and the phase speed of the imposed waves relative to the friction velocity, i.e., the wave age varies from very slow to fast waves. The wave signal is clearly present in the airflow up to at least 0.15λ (where λ is the wave length) and is present up to higher levels for faster waves. In the kinetic energy budgets, pressure transport is mainly of importance for slow waves. For fast waves, viscous transport and turbulent transport dominate near the surface. Kinetic energy budgets for the wave and turbulent perturbations show a non-negligible transport of turbulent kinetic energy directed from turbulence to the wave perturbation in the airflow. The wave-turbulent energy transport depends on the size, tilt, and phase of the wave-induced part of the turbulent Reynolds stresses.According to the DNS data, slow waves are more efficient in generating isotropic turbulence than fast waves.Despite the differences in wave-shape as well as in Reynolds number between the idealized direct numerical simulations and the atmosphere, there are intriguing similarities in the turbulence structure. Important information about the turbulence above waves in the atmosphere can be obtained from DNS—the data must, however, be interpreted with care.  相似文献   

9.
A model is developed to simulate the potential temperature and the height of the mixed layer under advection conditions. It includes analytic expressions for the effects of mixed-layer conditions upwind of the interface between two different surfaces on the development of the mixed layer downwind from the interface. Model performance is evaluated against tethersonde data obtained on two summer days during sea breeze flow in Vancouver, Canada. It is found that the mixed-layer height and temperature over the ocean has a small but noticeable effect on the development of the mixed layer observed 10 km inland from the coast. For these two clear days, the subsidence velocity at the inversion base capping the mixed layer is estimated to be about 30 mm s–1 from late morning to late afternoon. When the effects of subsidence are included in the model, the mixed-layer height is considerably underpredicted, while the prediction for the mean potential temperature in the mixed layer is considerably improved. Good predictions for both height and temperature can be obtained when values for the heat entrainment ratio,c, 0.44 and 0.68 for these two days respectively for the period from 1000 to 1300 LAT, were used. These values are estimated using an equation including the additional effects on heat entrainment due to the mechanical mixing caused by wind shear at the top of the mixed layer and surface friction. The contribution of wind shear to entrainment was equal to, or greater than, that from buoyant convection resulting from the surface heat flux. Strong wind shear occurred near the top of the mixed layer between the lower level inland flow and the return flow aloft in the sea breeze circulation.Symbols c entrainment parameter for sensible heat - c p specific heat of air at constant pressure, 1010 J kg–1 K–1 - d 1 the thickness of velocity shear at the mixed-layer top, m - Q H surface sensible heat flux, W m–2 - u m mean mixed-layer wind speed, m s–1 - u * friction velocity at the surface, m s–1 - w subsidence velocity, m s–1 - W subsidence warming,oC s–1 - w e entrainment velocity, m s–1 - w * convection velocity in the mixed layer, m s–1 - x downwind horizontal distance from the water-land interface, m - y dummy variable forx, m - Z height above the surface, m - Z i height of capping inversion, m - Z m mixed-layer depth, i.e.,Z i–Zs, m - Z s height of the surface layer, m - lapse rate of potential temperature aboveZ i, K m–1 - potential temperature step atZ i, K - u h velocity step change at the mixed-layer top - m mean mixed-layer potential temperature, K  相似文献   

10.
Measurements of the temperature and zonal velocity fields which develop in a rotating annulus of fluid with an upper surface, differentially heated from the inner to outer cylinder, are described for the lower symmetric regime (small radial temperature differences). The temperature field is essentially conductive for moderate to large rotation rates, Ω (>1.0 sec−1). The zonal velocity field is poorly approximated by the thermal wind equation.Measurements of the transition to waves from the lower symmetric regime at very large rotation rates are presented for positive and negative radial temperature differences. They suggest that the centrifugal buoyancy force and the free surface curvature may be important factors for the lower symmetric-wave transition at large Ω. By varying the stratification of the fluid over a range of 103 independently of the radial temperature difference, ΔrwT, it is conclusively shown that several theories are correct in predicting that the lower symmetric transition is independent of the stratification at small ΔrwT > 0 for large enough Ω.  相似文献   

11.
A new method for obtaining instantaneous vertical profiles of two components of velocity and temperature in thermally stratified turbulent shear flows is presented. In this report, the design and construction of the traversing system will be discussed and results to date will be presented. The method is based on rapid vertical sampling whereby probe sensors are moved vertically at a high speed such that the measurement is approximately instantaneous. The system is designed to collect many measurements for the calculation of statistics such as vertical wave number spectra, mean square vertical gradients, and Thorpe scales. Results are presented for vertical profiles of temperature and compared to vertical profiles measured by single-point Eulerian time averages. The quality of the vertical profiles is found to be good over many profiles. Some comparisons are made between vertical measurements and standard single-point Eulerian measurements for three cases of stably stratified turbulent shear flow in which the initial microscale Reynolds number, Reλ≈30. In case 1, the mean conditions are characterized by a gradient Richardson number, Rig=0.015, for which the flow is “unstable”, meaning the spatially evolving turbulent kinetic energy (Ek) grows. In case 2, Rig=0.095, for which the evolving turbulent kinetic energy is almost constant. In case 3, the flow is highly stable, where Rig=0.25 and Ek decays with spatial evolution. The measurements indicate anisotropy in the small scales for all cases. In particular, it is found that the ratio grows initially to a maximum and then decays with further evolution. Maximum Thorpe displacements are measured and compared to single-point measures of the vertical scales. It is found that vertical length scales derived from single-point measurements, such as the Ozmidov scale, LO=(ε/N3)1/2 and the overturn scale, Lt=θ′/(dT/dz), do not represent well the wide range of overturning scales which are actually present in the turbulence.  相似文献   

12.
  相似文献   

13.
Spectra and co-spectra of the streamwise (u) and normal or vertical (w) velocity fluctuations have been measured in the inner region of a large Reynolds number laboratory boundary layer over a rough wall. There is reasonable evidence of ak 1 –1 range in theu spectrum (wherek 1 is the streamwise wavenumber). Such a range results from an overlap between a spectral region dominated by largescale, inactive motion, which scales on the boundary-layer thickness, and a region dominated by smaller-scale, active motion which scales on the distance from the wall. Spectra ofw, anduw cospectra, scale in a manner consistent with the dominance by active motion. The present spectral data do not support local isotropy over the inertial subrange. A comparison between measuredw spectra and calculations based on isotropy indicates that the inertial subrange anisotropy is only slightly affected by the magnitude of the non-dimensional mean shear.  相似文献   

14.
To investigate the stability of the bottom boundary layer induced by tidal flow (oscillating flow) in a rotating frame, numerical experiments have been carried out with a two-dimensional non-hydrostatic model. Under homogeneous conditions three types of instability are found depending on the temporal Rossby number Rot, the ratio of the inertial and tidal periods. When Rot < 0.9 (subinertial range), the Ekman type I instability occurs because the effect of rotation is dominant though the flow becomes more stable than the steady Ekman flow with increasing Rot. When Rot > 1.1 (superinertial range), the Stokes layer instability is excited as in the absence of rotation. When 0.9 < Rot < 1.1 (near-inertial range), the Ekman type I or type II instability appears as in the steady Ekman layer. Being much thickened (100 m), the boundary layer becomes unstable even if tidal flow is weak (5 cm/s). The large vertical scale enhances the contribution of the Coriolis effect to destabilization, so that the type II instability tends to appear when Rot > 1.0. However, when Rot < 1.0, the type I instability rather than the type II instability appears because the downward phase change of tidal flow acts to suppress the latter. To evaluate the mixing effect of these instabilities, some experiments have been executed under a weak stratification peculiar to polar oceans (the buoyancy frequency N2  10−6 s−2). Strong mixing occurs in the subinertial and near-inertial ranges such that tracer is well mixed in the boundary layer and an apparent diffusivity there is evaluated at 150–300 cm2/s. This suggests that effective mixing due to these instabilities may play an important role in determining the properties of dense shelf water in the polar regions.  相似文献   

15.
A laboratory study in a rotating stratified basin examines the instability and long time evolution of the geostrophic double gyre introduced by the baroclinic adjustment to an initial basin-scale step height discontinuity in the density interface of a two-layer fluid. The dimensionless parameters that are important in determining the observed response are the Burger number S=R/R0 (where R is the baroclinic Rossby radius of deformation and R0 is the basin radius) and the initial forcing amplitude (H1 is the upper layer depth). Experimental observations and a numerical approach, using contour dynamics, are used to identify the mechanisms that result in the dominance of nonlinear behaviour in the long time evolution, τ>2−1 (where τ is time scaled by the inertial period TI=2π/f). When the influence of rotation is moderate (0.25≤S≤1), the instability mechanism is associated with the finite amplitude potential vorticity (PV) perturbation introduced when the double gyre is established. On the other hand, when the influence of rotation is strong (S≤0.1), baroclinic instability contributes to the nonlinear behaviour. Regardless of the mechanism, nonlinearity acts to transfer energy from the geostrophic double gyre to smaller scales associated with an eddy field. In the lower layer, Ekman damping is pronounced, resulting in the dissipation of the eddy field after only 40TI. In the upper layer, where dissipative effects are weak, the eddy field evolves until it reaches a symmetric distribution of potential vorticity within the domain consisting of cyclonic and anticyclonic eddy pairs, after approximately 100TI. The functional dependence of the characteristic eddy lengthscale LE on S is consistent with previous laboratory studies on continuously forced geostrophic turbulence. The cyclonic and anticyclonic eddy pairs are maintained until viscous effects eventually dissipate all motion in the upper layer after approximately 800TI. The outcomes of this study are considered in terms of their contribution to the understanding of the energy pathways and transport processes associated with basin-scale motions in large stratified lakes.  相似文献   

16.
A study has been made of the effects of varying the (uniform) grid resolution of a one-dimensional finite-difference numerical model of the dry convective boundary layer. The resolution of the inversion at the top of the boundary layer, and representation of the entrainment at the inversion, are found to influence the development of the momentum and buoyancy flux profiles. The modelled change in potential energy in a developing mixed layer is used to define a mixed layer scale, h m, which is found to vary systematically with resolution. The discretization errors (which can be large for resolutions poorer than a few tens of metres, particularly in the early stages of mixed-layer development) are quantified.  相似文献   

17.
We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s–1) - u g, vg components of the geostrophic wind (cm s–1) - V=(u2+v2)1/2 (cm s–1) - (û v)= 1/(h – k) k h(u, v)dz(cm s–1) - V * friction velocity (cm s–1) - z 0 roughness length (cm) - k 1 von Karman constant =0.4 - V d deposition velocity (cm s–1) - V g gravitational settling velocity (cm s–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - gr potential temperature at ground (°K) - K potential temperature at top of surface layer (°K) - P pressure (mb) - P 0 sfc pressure (mb) - C p/Cv - (t)= /z lapse rate of potential temperature (°K cm–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C HO transfer coefficient for sensible heat - C dust concentration (g m–3) - C K dust concentration at top of surface layer (g m–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s–1) - A 1 constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s–2 - Re Reynolds number = - D s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.  相似文献   

18.
Packet Structure of Surface Eddies in the Atmospheric Boundary Layer   总被引：6，自引：6，他引：0  
A smoke visualization experiment has beenperformed in the first 3,m ofneutral and unstable atmospheric boundary layersat very large Reynolds number(Re > 106). Under neutral atmosphericconditions mean wind profiles agreewell with those in the canonical flatplate zero-pressure-gradient turbulentboundary layer. The experiment was designedto minimize the temperaturedifference between the passive marker (smoke)and the air to ensure that anyobserved structures were due to vortical, ratherthan buoyant, motions. Imagesacquired in the streamwise–wall-normal planeusing a planar laser light-sheetare strikingly similar to those observed inlaboratory experiments at low to moderate Reynolds numbers. They reveal large-scaleramp-like structures withdownstream inclination of 3°–35°.This inclination isinterpreted as the hairpin packet growthangle following the hairpin vortexpacket model ofAdrian, Meinhart, and Tomkins.The distribution of this characteristicangle agrees with the results of experiments at far lower Reynolds numbers,suggesting a similarity in structures among low, moderate, and high Reynoldsnumber boundary layers at vastly different scales. These results indicate thatthe hairpin vortex packet model extends over a large range of scales. Theeffect of vertical heat transport in an unstable atmosphere on wall structuresis investigated in terms of the hairpin vortex packet model.  相似文献   

19.
The geostrophic Ekman boundary layer for large Rossby number (Ro) has been investigated by exploring the role played by the mesolayer (intermediate layer) lying between the traditional inner and outer layers. It is shown that the velocity and Reynolds shear stress components in the inner layer (including the overlap region) are universal relations, explicitly independent of surface roughness. This universality of predictions has been supported by observations from experiment, field and direct numerical simulation (DNS) data for fully smooth, transitionally rough and fully rough surfaces. The maxima of Reynolds shear stresses have been shown to be located in the mesolayer of the Ekman boundary layer, whose scale corresponds to the inverse square root of the friction Rossby number. The composite wall-wake universal relations for geostrophic velocity profiles have been proposed, and the two wake functions of the outer layer have been estimated by an eddy viscosity closure model. The geostrophic drag and cross-isobaric angle predictions yield universal relations, which are also supported by extensive field, laboratory and DNS data. The proposed predictions for the geostrophic drag and the cross-isobaric angle compare well with data for Rossby number Ro ≥ 105. The data show low Rossby number effects for Ro < 105 and higher-order effects due to the mesolayer compare well with the data for Ro ≥ 103.  相似文献   

20.
Functional relationships linking at λ0=351 nm aerosol extinction αλ0aer and backscatter coefficient βλ0aer of maritime and desert type aerosols are determined to allow for inversion of the single-wavelength lidar signals. Such relationships are derived as mean behavior of 20,000 extinction versus backscatter computations, performed for aerosol size distributions and compositions whose describing parameters are randomly chosen within the naturally observed variability. For desert-type aerosols, the effect of the particle non-sphericity is considered and it is shown that the extinction to backscatter ratio of non-spherical dust particles can be up to 60% larger than the values obtained for spherical particles. Aerosol extinction and backscatter coefficient profiles obtained inverting the single-wavelength lidar signal with the modeled relationships are then compared to the same profiles measured by a combined elastic-Raman lidar operating at 351 nm. Analytical back trajectories and satellite images are used to characterize advection patterns during lidar measurements and to properly choose the modeled functional relationship. A good accordance between the two techniques is found for advection patterns over the lidar site typical of maritime and dust conditions. Maximum differences between the model-based αλ0aer and βλ0aer vertical profiles and the corresponding ones measured by the combined elastic-Raman lidar technique are of 30% and 40% in maritime and desert dust conditions, respectively. The comparison of elastic-Raman lidar measurements and model-based results also reveals that particle non-sphericity must be taken into account when mineral dust-type aerosols are directly advected over the measurement site.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号