首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 765 毫秒
1.
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.  相似文献   

2.
A comparison is made of numerical and experimental results for flow over two-dimensional hills in both neutral and stably stratified flow. The numerical simulations are carried out using a range of one-and-a-half order and second-order closure schemes. The performance of the various turbulence schemes in predicting both the mean and turbulent quantities over the hill is assessed by comparing the results with new wind-tunnel measurements. The wind-tunnel experiments include both neutral and stably stratified flow over two different hills with different slopes, one of which is steep enough to induce flow separation. The dataset includes measurements of the mean and turbulent parts of the flow using laser Doppler anemometry. Pressure measurements are also made across the surface of the hill. These features make the dataset an excellent test of the model performance. In general second-order turbulence schemes provide the best agreement with the experimental data, however, they can be numerically unstable for steep hills. Some modifications can be made to the standard one-and-a-half order closure scheme, which results in improved performance at a fraction of the computation cost of the second-order schemes.  相似文献   

3.
Within the framework of the semiempirical theory of turbulence for stratified fluids some aspects of the problem of internal wave-turbulence interaction in the upper layer of the ocean are discussed. The conditions of amplification and sustaining of turbulence by internal waves are investigated. Stationary distributions of turbulent energy are found for a stratified fluid with a shear flow produced, for example, by a low-frequency internal wave. The internal wave damping due to both turbulent viscosity and turbulent diffusion in the thermocline is studied. For a two-layer model damping constant is determined as a function of the wave number. The variation of surface turbulence by internal waves is estimated and the role of this process in slick formation is considered.  相似文献   

4.
Laboratory and numerical model simulations of turbulent circulations within the wake regions of thunderstorm outflows have been done with the assumption that there is no turbulence within the ambient airmass. Furthermore, many observational studies have used Doppler radar data that have been filtered so that turbulent structures are reduced in amplitude or eliminated altogether. This study presents unique Doppler radar observations of the collision of a roll-like boundary-layer circulation with a gust flow. The boundary-layer circulation is seen to interact with the circulation within the gust flow head and to reappear within the wake region. It is suggested that the ambient boundary layer may be an energy source for the generation and/or maintenance of turbulence in the wake region.  相似文献   

5.
The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, having a significant impact on medium-range weather forecasts and climate integrations. To pursue this further, a moderately stratified Arctic case is simulated by nineteen single-column turbulence schemes. Statistics from a large-eddy simulation intercomparison made for the same case by eleven different models are used as a guiding reference. The single-column parameterizations include research and operational schemes from major forecast and climate research centres. Results from first-order schemes, a large number of turbulence kinetic energy closures, and other models were used. There is a large spread in the results; in general, the operational schemes mix over a deeper layer than the research schemes, and the turbulence kinetic energy and other higher-order closures give results closer to the statistics obtained from the large-eddy simulations. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.  相似文献   

6.
We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.  相似文献   

7.
A new approach to investigations of the structure of the boundary layer above waves is discussed. The approach is based on direct numerical simulation of wave motions in the boundary layer produced by a moving curved surface. Model equations are derived, which are the Reynolds equations in a curvilinear nonstationary system of co-ordinates, evolution equations for turbulent kinetic energy, and Kolmogorov's approximate similarity formulae relating the coefficient of turbulent viscosity to the dissipation of turbulent energy; the length scale is assumed to grow linearly with increasing distance from the surface. Principles of constructing the model numerical scheme are described. Results are given of modelling the structure of the boundary layer above a nonsteady surface, which, in a general case, is a superposition of progressive waves with assigned dispersion relations and amplitudes. Mechanisms of energy and momentum transfer to the surface, effects of density stratification and energy structure in the boundary layer are studied. Merits and demerits of the approach are discussed.  相似文献   

8.
Along-Coast Features of Bora-Related Turbulence   总被引：1，自引：1，他引：0  
The along-coast, offshore turbulence structure of the Bora flow that occurred on 7 November 1999 during the Mesoscale Alpine Programme (MAP) Intensive Observation Period 15 is examined. In this analysis we employ the aircraft and dropsonde data obtained over the Adriatic Sea, where the turbulence structure is determined by estimating turbulent kinetic energy (TKE) and its dissipation rate along the flight legs. The turbulence characteristics of Bora in the lee of the Dinaric Alps is greatly influenced by the mesoscale Bora flow structure over the Adriatic Sea, which in the cross-wind direction features an interchange of jets and wakes related to mountain gaps and peaks. In order to establish the origin of turbulence, the Weather Research and Forecasting—Advanced Research WRF (WRF-ARW) numerical model is used and its results are compared to the measurements. All five TKE-prediction parametrization schemes available in the model show reasonable agreement with the measured values. Since these parametrization schemes do not have horizontal advection included, they suggest that the along-flight structure of the Bora turbulence is principally generated by the local vertical wind shear. Further evidence is needed to support this hypothesis.  相似文献   

9.
The air flow above breaking monochromatic Stokes waves is studied using a numerical nonlinear model of the turbulent air flow above waves of finite amplitude. The breaking event (spilling breaker) is parameterized by increasing the local roughness at the downwind slope of the wave, just beyond the crest. Both moderate slope waves and steep waves are considered. Above steep breaking waves, a large increase (typically 100%) in the total wind stress — averaged over the wave profile — is found compared to nonbreaking moderate slope waves. This is due to the drastic increase of the form drag, which arises from the asymmetrical surface pressure pattern above breaking waves. Both increase of wave slope (sharpening of the crest) and increase of local roughness in the spilling breaker area cause this asymmetrical surface pressure pattern. A comparison of the numerical results with the recent experimental measurements of Banner (1990) is carried out and a good agreement is found for the structure of the pressure pattern above breaking waves and for the magnitude of enhanced momentum transfer. Also: Dept. of Applied Physics, Techn. Univ. Delft, Netherlands.  相似文献   

10.
In this work, three turbulence closure models, Mellor andYamada level 2.5, E - l and E - implemented in a circulation model, are compared in neutral condition over complex terrain. They are firstly applied to a one-dimensional case on flat terrain and then to a schematic two-dimensional valley. The simulation results, in terms of wind field and turbulent kinetic energy, are tested against measurements from a wind-tunnel experiment. The empirical constants defining the characteristic length scales of the closures are modified based on turbulence parameters estimated in the experiment. The formulation of the diffusion coefficients is analysed to explain the differences among the various closures in the simulation results. Regarding the mean flow, both on flat and complex terrain, all the closures yield satisfactory results. Concerning the turbulent kinetic energy, the best results are obtained by E - l and E - closures.  相似文献   

11.
A mesoscale planetary boundary layer (PBL) numerical model has been developed to study airflow over complex topography. Turbulence closures using the turbulent kinetic energy (TKE) and dissipation () equations are investigated in combination with the level 2.5 scheme of Mellor and Yamada (1982) to determine eddy diffusivities for momentum and heat. This modified E- closure is simpler than the level 3 one which requires more prognostic equations for moist turbulent transport.One-dimensional (1-D) model results show that the PBL mean flows under various stability conditions are not significantly sensitive to the modified Blackadar and Kolmogorov eddy mixing-length formulations used in this E- model, although the latter yields excessively large mixing lengths in the entrainment region of the upper PBL. Eddy mixing lengths in the Kolmogorov-type formulation can be better defined by introducing background dissipation. Using the same prognostic TKE equation, the 1-D model results are not significantly affected by different diagnostic formulations in the closures. The simulated results compare well with large-eddy simulations and those obtained using higher-order closure schemes including the level 3 one. The results are found to be insensitive to eddy Prandtl number, in contrast to the 2-D model results (see Part II).  相似文献   

12.
We investigate a Cartesian-mesh immersed-boundary formulation within an incompressible flow solver to simulate laminar and turbulent katabatic slope flows. As a proof-of-concept study, we consider four different immersed-boundary reconstruction schemes for imposing a Neumann-type boundary condition on the buoyancy field. Prandtl’s laminar solution is used to demonstrate the second-order accuracy of the numerical solutions globally. Direct numerical simulation of a turbulent katabatic flow is then performed to investigate the applicability of the proposed schemes in the turbulent regime by analyzing both first- and second-order statistics of turbulence. First-order statistics show that turbulent katabatic flow simulations are noticeably sensitive to the specifics of the immersed-boundary formulation. We find that reconstruction schemes that work well in the laminar regime may not perform as well when applied to a turbulent regime. Our proposed immersed-boundary reconstruction scheme agrees closely with the terrain-fitted reference solutions in both flow regimes.  相似文献   

13.
Further development of the non-linear mixed spectral finite difference (NLMSFD) model of turbulent boundary-layer flow over topography is documented. This includes modifications and refinements to the solution procedure, the incorporation of second-order turbulence closures to the model and the three-dimensional extension of the model. Based on these higher order closures, linear limitations, boundary-layer approximation and non-linear effects are discussed. The impact of different turbulence closures on the prediction of the NLMSFD model is also demonstrated. Furthermore, sample results for 3D idealized topography (sinusoidal) are presented. The parameterization of drag over small-scale topography is also addressed.  相似文献   

14.
We perform large-eddy simulations of neutral atmospheric boundary-layer flow over a cluster of buildings surrounded by relatively flat terrain. The first investigated question is the effect of the level of building detail that can be included in the numerical model, a topic not yet addressed by any previous study. The simplest representation is found to give similar results to more refined representations for the mean flow, but not for turbulence. The wind direction on the other hand is found to be important for both mean and turbulent parameters. As many suburban areas are characterised by the clustering of buildings and homes into small areas separated by surfaces of lower roughness, we look at the adjustment of the atmospheric surface layer as it flows from the smoother terrain to the built-up area. This transition has unexpected impacts on the flow; mainly, a zone of global backscatter (energy transfer from the turbulent eddies to the mean flow) is found at the upstream edge of the built-up area.  相似文献   

15.
An analytical model is developed for the initial stage of surface wave generation at an air–water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417–445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate Γ, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer.For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as ∝ 1/Γ), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.  相似文献   

16.
复杂下垫面地域边界层结构的三维细网格数值模拟   总被引：8，自引：0，他引：8  
建立了一个可供复杂下垫面地域使用的非静力的三维细网格边界层模式，就复杂下垫面条件下的边界层结构和湍流特征作了以实例为对照的数值模拟试验，模式采用能量闭合方案，舍弃了静力近似。以实测资料为初台输入，同时还做了一些数值试验，分别获得了采用静力与非静力模式和不同闭合方案对ＰＢＬ模型的结果。  相似文献   

17.
The dependence of the turbulent airflow over water waves on the angle,, between mean wind and wavedirections is investigated. To this end,an existing semi-analytical model is extended. In this model, the main simplification of the problem is obtained by using the well-established divisionof the wave boundary layer into inner and outer regions for modelling turbulence. The effect of waves on turbulence is restricted to the thin inner region. Simulations show that the influence of the wind speed component transverse to the wave direction on the air flow, and hence on the growth rate of the waves, is small. This is confirmed by calculations with a numerical model that solves the full Reynolds equations using a second-order turbulence closure scheme. The growth rate of slowly moving waves (as compared to the wind speed) is then proportional to cos2, whereas, for faster waves, it has a narrower angular distribution.  相似文献   

18.
The dispersion of heavy particles subjected to a turbulent forcing is often simulated with Lagrangian stochastic models. Although these models have been employed successfully over land, the implementation of traditional LS models in the marine boundary layer is significantly more challenging. We present an adaptation of traditional Lagrangian stochastic models to the atmospheric marine boundary layer with a particular focus on the representation of the scalar turbulence for temperature and humidity. In this new model, the atmosphere can be stratified and the bottom boundary is represented by a realistic wavy surface that moves and deforms. Hence, the correlation function for the turbulent flow following a particle is extended to the inhomogenous, anisotropic case. The results reproduce behaviour for scalar Lagrangian turbulence in a stratified airflow that departs only slightly from the expected behaviour in isotropic turbulence. When solving for the surface temperature and the radius of evaporating heavy water droplets in the airflow, the modelled turbulent forcing on the particle also behaves remarkably well. We anticipate that this model will prove especially useful in the context of sea-spray dispersion and its associated sensible heat, latent heat, and gas fluxes between spray droplets and the atmosphere.  相似文献   

19.
Here we advance the physical background of the energy- and flux-budget turbulence closures based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth’s rotation. In accordance with modern experimental evidence, the closure implies the maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: “strong turbulence” at ${Ri \ll 1}$ typical of boundary-layer flows and characterized by the practically constant turbulent Prandtl number Pr T; and “weak turbulence” at Ri > 1 typical of the free atmosphere or deep ocean, where Pr T asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy viscosity and eddy conductivity, and a general non-gradient model based on prognostic equations for all the basic parameters of turbulence including turbulent fluxes.  相似文献   

20.
为了深入理解非静力近似下的波-湍相互作用问题,本研究在σ坐标的海洋环境研究和预报模型(MERF)中引入常用的Mellor-Yamada两方程湍混合参数化方案(MY2.5),评估垂向湍混合对小尺度背风波传播过程的影响.瞬时状态场的模拟结果表明,无论是否为非静力近似条件,上述湍参数化方案的引入都会减弱背风波传播的模拟效果.从时间平均场的试验结果来看,垂向湍混合过程会显著减小非静力近似和静力近似之间的差异.此外,能量收支分析的诊断结果表明,MY2.5方案会显著抑制陆坡地形下的背风波传播过程,进而将更多的潮能转化到不可逆的湍混合过程中.  相似文献   

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

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