首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
Laboratory experiments were carried out to investigate the interaction between turbulent line buoyant plumes and sharp density interfaces, with the aim of using the results to interpret oceanic observations pertinent to crack openings in the polar ice-cap (leads). These openings take the form of long narrow channels, and are often modeled as line bouyant plumes. The plumes descend as in a homogenoous fluid, impinge on the density interface, and then spread horizontally as gravity currents. Depending on the Richardson number , where Δb is the buoyancy jump across the interface, lD is the half-width of the plume before the impingement and q0is the buoyancy flux per unit length of the source, different flow patterns were identified. When Ri < 0.5, the plumes penetrate deep into the bottom layer, deflect horizontally and then spread while showing little vertical rise. When 0.6 < Ri < 5, the penetration is significant, but the fluid bounces back after entraining heavy fluid from the lower layer and then spreads horizontally above the interface as a gravity current. Appreciable mixing between this current and the lower layer was detected when Ri <1. When Ri > 10, the penetration was small and a sharp-nosed gravity current emerged some time after the impact. Measurements were made on the penetration depth, the velocities of the gravity current and the subsurface flow towards the plume, the entrainment rate and other wave parameters. Possible implications of the results for oceanic cases are also discussed.  相似文献   

2.
Turbulent mixing across heat-stratified density interfaces was studied in the laboratory using oscillating-grid generated turbulence. The aim was to study the transition between the entrainment regimes dominated by interfacial wave-breaking and molecular diffusion, and to study the characteristics of the latter. It was observed that, above a critical Richardson number Ric, which depends on the Peclet number Pe, the mixing due to wave breaking disappears and that Ric Pe−n, where the mean value of the exponent n is approximately . Above Ric, the entrainment is molecular-diffusion dominated and takes place through a sequence of events: the buoyancy gradient of the initially sharp density interface is weakened by molecular diffusion until the mixed-layer eddies can engulf a portion of the interfacial layer wherefore the interface sharpens again. Thus, the entrainment events are recurrent with a rate-controlling diffusion stage between them. An entrainment law of the form E Ri−2Pe−2, where E is the entrainment coefficient and Ri is the Richardson number, is suggested for the diffusion-dominated entrainment regime.  相似文献   

3.
Turbulent convection forced by a surface heat flux into a stably stratified region is a feature of both the atmospheric and oceanic planetary boundary layers. Of particular interest is the interface between the convective layer and the stable stratification, where the entrainment of fluid into the convective layer by penetrating plumes may lead to a reverse buoyancy flux, and an enhancement of the stable stratification. Whereas in the atmosphere the influence of rotation on this penetrative convection is negligible, oceanic convection may be subjected to lower Rossby numbers and hence greater rotational influence. To isolate the effects of rotation, we present three numerical solutions for turbulent penetrative convection, characterised by different rotation rates, with all other parameters being held constant. Our results indicate that at lower Rossby numbers the lateral scale of the plumes is reduced, whereas the vertical vorticity of the plumes is much enhanced. Vertical transports of buoyancy and kinetic energy across the convective layer are reduced, leading to less efficient penetration at the interface with the stratified layer, and hence less reverse buoyancy flux in this region.  相似文献   

4.
For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w *. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).  相似文献   

5.
The transfer processes within and above a simulated urban street canyon were investigated in a generic manner. Computational fluid dynamics (CFD) was used to aid understanding and to produce some simple operational parameterisations. In this study we addressed specifically the commonly met situation where buoyancy effects arising from elevated surface temperatures are not important, i.e. when mechanical forces outweigh buoyancy forces. In a geophysical context this requires that some suitably defined Richardson number is small. From an engineering perspective this is interpreted as the important case when heat transfer within and above urban street canyons is by forced convection. Surprisingly, this particular scenario (for which the heat transfer coefficient between buildings and the flow is largest), has been less well studied than the situation where buoyancy effects are important. The CFD technique was compared against wind-tunnel experiments to provide model evaluation. The height-to-width ratio of the canyon was varied through the range 0.5–5 and the flow was normal to the canyon axis. By setting the canyon’s facets to have the same or different temperatures or to have a partial temperature distribution, simulations were carried out to investigate: (a) the influence of geometry on the flow and mixing within the canyon and (b) the exchange processes within the canyon and across the canyon top interface. Results showed that the vortex-type circulation and turbulence developed within the canyon produced a temperature distribution that was, essentially, spatially uniform (apart from a relatively thin near-wall thermal boundary layer) This allowed the temperatures within the street canyon to be specified by just one value T can , the canyon temperature. The variation of T can with wind speed, surface temperatures and geometry was extensively studied. Finally, the exchange velocity u E across the interface between the canyon and the flow above was calculated based on a heat flux balance within the canyon and between the canyon and the flow above. Results showed that u E was approximately 1% of a characteristic wind velocity above the street canyon. The problem of radiative exchange is not addressed but it can, of course, be introduced analytically, or computationally, when necessary.  相似文献   

6.
A stainless steel soil corer which was filled with homogenized soil was used to measure the flux (J) of NO between soil and atmosphere and the vertical profile of the NO mixing ratios (m) in the soil atmosphere, both as function of the NO mixing ratio (mm a ) in the atmosphere of the headspace. The NO emission flux decreased linearly with increasing NO mixing ratio and turned into a deposition flux after passage of the compensation point (m c) at about 400 ppbv NO. Almost the same compensation point was obtained when the turnover of NO was measured in flask-incubated soil samples as function of the NO mixing ratio. The flux (J) of NO at the soil-atmosphere interface was calculated from the production rate (P) of NO and the NO uptake rate constant (k) that were measured in these flask-incubated soil samples using the diffusion model of Galbally and Johansson (1989). The calculated fluxes agreed within <15% with those actually measured. The vertical profiles of NO were fitted to an exponential function and analyzed by Fick's first law of diffusion. The shape of the profiles indicated a net production of NO in the upper 10 cm soil layer when the atmospheric NO mixing ratio was below the compensation point and in a net consumption of NO when the atmospheric NO mixing ratio was above the compensation point. In soil layers below 10 cm depth, the turnover of NO resulted in compensation of production and consumption rates. Measurement of the actual diffusion coefficient using SF6 showed that gas transport in the soil core was not only due to molecular diffusion but in addition due to a bidirectional gas flow. The experimentally determined diffusion coefficient was smaller than that computed from soil porosities, but resulted together with the additional transport term in NO fluxes that were close (< ±15%) to those measured. This is the first comprehensive study of NO concentration profiles and turnover rates in soil providing a theoretical basis for modelling NO fluxes at the soil-atmosphere interface.  相似文献   

7.
A series of experiments were performed in a rotating cylindrical tank over a wide range of rotation rates in which convective turbulence was generated by a bottom-mounted heated plate in both homogeneous and stratified fluids. Measurements were made of the turbulent velocities in all three axes over the full depth of the chamber, and of the temperatures at the mid-depth near the centre of the tank. For even small rotation rates, the measurements showed that the turbulent velocities were weakly affected by rotation at all depths, but as the rotation rate increased, the deviation from the non-rotational scaling slowly and progressively increased until eventually the turbulent velocities were fully rotationally controlled. The results indicated that there was no sudden transition of the turbulent field from the non-rotational state (a function only of the surface buoyancy flux B and the depth z) to the rotational state (where the strength of the turbulent field is a function of only B and the Coriolis parameter f). Rather the transition was a smooth asymptotic one from one state to the other. Nevertheless, it was possible to parametrize this transition by a single value of the turbulent or small scale Rossby number, defined by Ro = (B/f3z2)1/3. Our measurements suggested a critical value of Roc ≈ 0.1, below which the turbulence was fully rotationally controlled and which was equivalent to a critical depth zc = (35 ± 15)(B/f3)1/2. Using typical oceanic values for B and f, the oceanic turbulence driven by surface cooling events becomes rotationally controlled only for depths greater than about 10 km, a depth which is greater than that of the bulk of the world's oceans. Thus, convective turbulence actively being generated by cooling of the ocean surface is best described by non-rotating turbulent velocity and length scales and is a function only of the surface buoyancy flux and the depth.  相似文献   

8.
In a series of laboratory experiments, a partially mixed patch was produced in thick linear concentration gradients favorable to salt-finger convection. Salt-fingers, which give rise to an up-gradient flux of buoyancy, can reduce and invert the density gradient in the initial imposed patch. This leads to overturning convection within the patch if (a) the ratio of ambient T and S gradients, Rρ=αTz/βSz, is near one; (b) the initial imposed turbulence results in a nearly well-mixed patch; and (c) the patch thickness is large enough that convective eddies are able to transport T and S faster than salt-fingers. Once overturning occurs, subsequent turbulent entrainment can lead to growth of the patch thickness. Experimental results for one-dimensional patches (layers) agree well with the theoretical prediction. This thickening is in contrast to the collapse that a partially mixed three-dimensional patch would experience due to lateral intrusion in a wide tank.  相似文献   

9.
  相似文献   

10.
Turbulent organized structures (TOS) above building arrays were investigated using a large-eddy simulation (LES) model for a city (LES-CITY). Square and staggered building arrays produced contrasting behaviour in terms of turbulence that roughly corresponded to the conventional classification of ‘D-type’ and ‘K-type’ roughness, respectively: (1) The drag coefficients (referred to the building height) for staggered arrays were sensitive to building area density, but those for square arrays were not. (2) The relative contributions of ejections to sweeps (S2/S4) at the building height for square arrays were sensitive to building area density and nearly equalled or exceeded 1.0 (ejection dominant), but those for staggered arrays were insensitive to building area density and were mostly below 1.0 (sweep dominant). (3) Streaky patterns of longitudinal low speed regions (i.e., low speed streaks) existed in all flows regardless of array type. Height variations of the buildings in the square array drastically increased the drag coefficient and modified the turbulent flow structures. The mechanism of D-type and K-type urban-like roughness flows and the difference from vegetation flows are discussed. Although urban-like roughness flows exhibited mixed properties of mixing layers and flat-wall boundary layers as far as S2/S4 was concerned, the turbulent organized structures of urban-like roughness flows resembled those of flat-wall boundary layers.  相似文献   

11.
A striking feature observed during the turbulent entrainment across a density interface at moderately-large Richardson numbers is the presence of an intermediate layer sandwiched between the well-mixed and interfacial layers. In this note, a possible mechanism for its origin, and its role in buoyancy transfer during the entrainment is discussed. Estimates are made for the thickness of the intermediate layer and are compared with experimental results.Also Environmental Engineering Science, California Institute of Technology, Pasadena, CA 91125, U.S.A.  相似文献   

12.
Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies   总被引：1，自引：8，他引：1  
In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (FT) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating FT to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (ST) and FT have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if ST and FT are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.  相似文献   

13.
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.  相似文献   

14.
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.  相似文献   

15.
Abstract

Estimates of the rate of dissipation of turbulent energy,?v, were made on the continental shelf off the southwest coast of Vancouver Island throughout the diurnal tidal cycle at six stations. If turbulent mixing takes place away from boundaries, and a constant fraction of the energy supplied to turbulence is converted into potential energy, this dissipation rate per unit volume, ?v, is shown to be proportional to the rate of turbulent mass and nutrient flux across the isopycnals.

We compute these fluxes to determine the source of nutrients in the upper mixed layer on the shelf. It is found that tidal mixing on the shelf throughout the water column contributes less than 10% of the flux of nutrients supplied by the estuarine outflow out of Juan de Fuca Strait. Strong winds during upwelling events supply nutrients at rates greater than those due to tidal mixing, but at rates that are likely smaller than the Juan de Fuca source. Therefore, the nutrient‐rich waters observed in the euphotic zone in spring and summer on this shelf derive mainly from Juan de Fuca Strait. Nutrients in bottom water derive from upwelling along shelf‐break canyons.  相似文献   

16.
Modification of a turbulent flow due to a change from a smooth to a rough surface has been studied by means of a stream function-vorticity model. Results of four models of eddy viscosity (or turbulent exchange coefficient) K mhave been compared. The models are: (1) K m = l2S, where l is the mixing length and S is the deformation of mean flow; (2) K m E/S, which is based on the assumption that turbulent momentum flux is proportional to turbulent kinetic energy E; (3) K m lE1/2, the so called Prandtl-Kolmogoroff approach; and (4) K m E2/, the E — closure, where is the dissipation of turbulent kinetic energy.It is found that net-production, i.e., the difference of production and dissipation of turbulent kinetic energy counteracts the influence of mean shear on turbulent shear stress and diminishes turbulent shear stress. The reduction of mixing-length, being predicted by Model 4 only, adds to this attenuation. As a consequence, in Models 2 and 4, loss of horizontal mean momentum is concentrated close to the ground, which results in an inflexion point in the logarithmic, vertical profile of horizontal mean velocity. By contrast, in Models 1 and 3, modification of turbulent shear stress reaches larger heights causing deeper internal boundary layers. Concerning the existence of an inflexion point in U(lnz), the depth of the internal boundary layer for mean velocity, and the modification of bottom shear stress, Model 4 comes closest to experimental data.A remarkable difference of Models 1, 2, 3 and Model 4 is that only Model 4 predicts a very slow relaxation of eddy viscosity which can be attributed to the reduction of mixing-length.  相似文献   

17.
In spring preceding the record minimum summer ice cover detailed microstructure measurements were made from drifting pack ice in the Arctic Ocean, 110 km from the North Pole. Profiles of hydrography, shear, and temperature microstructure collected in the upper water column covering the core of the Atlantic Water are analyzed to determine the diapycnal eddy diffusivity, the eddy diffusivity for heat, and the turbulent flux of heat. Turbulence in the bulk of the cold halocline layer was not strong enough to generate significant buoyancy flux and mixing. Resulting turbulent heat flux across the upper cold halocline was not significantly different than zero. The results show that the low levels of eddy diffusivity in the upper cold halocline lead to small vertical turbulent transport of heat, thereby allowing the maintenance of the cold halocline in the central Arctic.  相似文献   

18.
Theoretical arguments are developed to describe the effects of a uniform slope on the development of the stably stratified atmospheric boundary layer (SBL). A maximum sustainable surface buoyancy flux exists for the SBL overlying a uniform, non-sloping surface. In this study it is shown that the SBL overlying a uniform shallow slope (with gradient of the order of 1:1000) also supports a maximum sustainable buoyancy flux, B max, but that the value of B max is influenced by the gradient of the slope, γ. It is demonstrated that in the limit γ → 0, results for the SBL over a horizontal surface are recovered.  相似文献   

19.
It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u */w *, where u * and w * are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z i , where z i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.  相似文献   

20.
The El Nin o-Southern Oscillation (ENSO) is modulated by many factors; most previous studies have emphasized the roles of wind stress and heat flux in the tropical Pacific. Freshwater flux (FWF) is another environmental forcing to the ocean; its effect and the related ocean salinity variability in the ENSO region have been of increased interest recently. Currently, accurate quantifications of the FWF roles in the climate remain challenging; the related observations and coupled ocean-atmosphere modeling involve large elements of uncertainty. In this study, we utilized satellite-based data to represent FWF-induced feedback in the tropical Pacific climate system; we then incorporated these data into a hybrid coupled ocean-atmosphere model (HCM) to quantify its effects on ENSO. A new mechanism was revealed by which interannual FWF forcing modulates ENSO in a significant way. As a direct forcing, FWF exerts a significant influence on the ocean through sea surface salinity (SSS) and buoyancy flux (Q B ) in the western-central tropical Pacific. The SSS perturbations directly induced by ENSO-related interannual FWF variability affect the stability and mixing in the upper ocean. At the same time, the ENSO-induced FWF has a compensating effect on heat flux, acting to reduce interannual Q B variability during ENSO cycles. These FWF-induced processes in the ocean tend to modulate the vertical mixing and entrainment in the upper ocean, enhancing cooling during La Nin a and enhancing warming during El Nin o, respectively. The interannual FWF forcing-induced positive feedback acts to enhance ENSO amplitude and lengthen its time scales in the tropical Pacific coupled climate system.  相似文献   

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

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