首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.  相似文献   

2.
The Weather Research and Forecasting (WRF) model can be used to simulate atmospheric processes ranging from quasi-global to tens of m in scale. Here we employ large-eddy simulation (LES) using the WRF model, with the LES-domain nested within a mesoscale WRF model domain with grid spacing decreasing from 12.15 km (mesoscale) to 0.03 km (LES). We simulate real-world conditions in the convective planetary boundary layer over an area of complex terrain. The WRF-LES model results are evaluated against observations collected during the US Department of Energy-supported Columbia Basin Wind Energy Study. Comparison of the first- and second-order moments, turbulence spectrum, and probability density function of wind speed shows good agreement between the simulations and observations. One key result is to demonstrate that a systematic methodology needs to be applied to select the grid spacing and refinement ratio used between domains, to avoid having a grid resolution that falls in the grey zone and to minimize artefacts in the WRF-LES model solutions. Furthermore, the WRF-LES model variables show large variability in space and time caused by the complex topography in the LES domain. Analyses of WRF-LES model results show that the flow structures, such as roll vortices and convective cells, vary depending on both the location and time of day as well as the distance from the inflow boundaries.  相似文献   

3.
Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.  相似文献   

4.
The influence of mesoscale circulations induced by urban-rural differential surface sensible heat flux and roughness on convective boundary-layer (CBL) flow statistics over an isolated urban area has been examined using large-eddy simulation (LES). Results are analyzed when the circulations influence the entire urban area under a zero background wind. For comparison, the CBL flow over an infinite urban area with identical urban surface characteristics under the same background meteorological conditions is generated as a control case (without circulations). The turbulent flow over the isolated urban area exhibits a mix of streaky structure and cellular pattern, while the cellular pattern dominates in the control case. The mixed-layer height varies significantly over the isolated urban area, and can be lower near the edge of the urban area than over the rural area. The vertical profiles of turbulence statistics over the isolated urban area vary horizontally and are dramatically different from the control case. The turbulent kinetic energy (TKE) sources include wind shear, convergence, and buoyancy productions, compared to only buoyancy production in the control case. The normalized vertical velocity variance is reduced compared to the control case except in the central urban area where it is little affected. The low-level flow convergence is mainly responsible for the enhanced horizontal velocity variance in the central urban area, while wind shear is responsible for the additional local maximum of the horizontal velocity variance near the middle of the CBL outside the central area. Parameterizations in the prognostic equation for TKE used in mesoscale models are evaluated against the LES results over the isolated urban area. We also discuss conditions under which the urban-induced circulations occur and when they may affect the entire urban area. Given that urban-induced circulations can influence the entire urban area within hours for an urban area of a realistic size, it is inappropriate to directly apply empirical relations of turbulence statistics derived under horizontally-homogenous flow conditions to an urban area.  相似文献   

5.
We examine the dynamics of turbulence subgrid (or sub-filter) scales over a lake surface and the implications for large-eddy simulations (LES) of the atmospheric boundary layer. The analysis is based on measurements obtained during the Lake-Atmosphere Turbulent EXchange (LATEX) field campaign (August–October, 2006) over Lake Geneva, Switzerland. Wind velocity, temperature and humidity profiles were measured at 20 Hz using a vertical array of four sonic anemometers and open-path gas analyzers. The results indicate that the observed subgrid-scale statistics are very similar to those observed over land surfaces, suggesting that the effect of the lake waves on surface-layer turbulence during LATEX is small. The measurements allowed, for the first time, the study of subgrid-scale turbulent transport of water vapour, which is found to be well correlated with the transport of heat, suggesting that the subgrid-scale modelling of the two scalars may be coupled to save computational resources during LES.  相似文献   

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

7.
A theoretical approach suggests that the surface heterogeneity on a scale of tens of kilometres can generate mesoscale motions that are not in a quasi-stationary state. The starting point of the theoretical approach is the equations of horizontal velocity and potential temperature that are low-pass filtered with a mesoscale cut-off wavelength. The transition of the generated mesoscale motions from a quasi-stationary state to a non-stationary state occurs when horizontal advection is strong enough to level out the potential temperature gradient on the surface heterogeneity scale. Large-eddy simulations (LES) suggest that the convective boundary layer (CBL) changes to a non-stationary state when forced by a surface heat-flux variation of amplitude of 100W m−2 or higher and a wavelength of the order of 10 km. Spectral analysis of the LES reveals that when the mesoscale motions are in a quasi-stationary state, the energy provided by the surface heat-flux variation remains in organized mesoscale motions on the scale of the surface variation itself. However, in a non-stationary state, the energy cascades to smaller scales, with the cascade extending down into the turbulence scale when the wavelength of the surface heat-flux variation is on a scale smaller than 100 times the CBL height. The energy transfer from the generated mesoscale motions to the CBL turbulence results in the absence of a spectral gap between the two scales. The absence of an obvious spectral gap between the generated mesoscale motions and the turbulence raises questions about the applicability of mesoscale models for studies on the effect of high-amplitude surface heterogeneity on a scale of tens of kilometres. The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

8.
Measurements and Computations of Flow in an Urban Street System   总被引：1，自引：1，他引：0  
We present results from laboratory and computational experiments on the turbulent flow over an array of rectangular blocks modelling a typical, asymmetric urban canopy at various orientations to the approach flow. The work forms part of a larger study on dispersion within such arrays (project DIPLOS) and concentrates on the nature of the mean flow and turbulence fields within the canopy region, recognising that unless the flow field is adequately represented in computational models there is no reason to expect realistic simulations of the nature of the dispersion of pollutants emitted within the canopy. Comparisons between the experimental data and those obtained from both large-eddy simulation (LES) and direct numerical simulation (DNS) are shown and it is concluded that careful use of LES can produce generally excellent agreement with laboratory and DNS results, lending further confidence in the use of LES for such situations. Various crucial issues are discussed and advice offered to both experimentalists and those seeking to compute canopy flows with turbulence resolving models.  相似文献   

9.
Large-eddy simulations were performed of a neutrally-stratified turbulent flow within and above an ideal, horizontally- and vertically-homogeneous plant canopy. Three simulations were performed for shear-driven flows in small and large computational domains, and a pressure-driven flow in a small domain, to enable the nature of canopy turbulence unaffected by external conditions to be captured. The simulations reproduced quite realistic canopy turbulence characteristics, including typical ramp structures appearing in time traces of the scalar concentration near the canopy top. Then, the spatial structure of the organised turbulence that caused the scalar ramps was examined using conditional sampling of three-dimensional instantaneous fields, triggered by the occurrence of ramp structures. A wavelet transform was used for the detection of ramp structures in the time traces. The ensemble-averaged results illustrate that the scalar ramps are associated with the microfrontal structure in the scalar, the ejection-sweep structure in the streamwise and vertical velocities, a laterally divergent flow just around the ramp-detection point, and a positive, vertically-coherent pressure perturbation. These vertical structures were consistent with previous measurements made in fields or wind tunnels. However, the most striking feature is that the horizontal slice of the same structure revealed a streamwise-elongated region of high-speed streamwise velocity impacting on another elongated region of low-speed velocity. These elongated structures resemble the so-called streak structures that are commonly observed in near-wall shear layers. Since elongated structures of essentially similar spatial scales were observed in all of the runs, these streak structures appear to be inherent in near-canopy turbulence. Presumably, strong wind shear formed just above the canopy is involved in their formation. By synthesis of the ensemble-averaged and instantaneous results, the following processes were inferred for the development of scalar microfronts and their associated flow structures: (1) a distinct scalar microfront develops where a coherent downdraft associated with a high-speed streak penetrates into the region of a low-speed streak; (2) a stagnation in flow between two streaks of different velocities builds up a vertically-coherent high-pressure region there; (3) the pressure gradients around the high-pressure region work to reduce the longitudinal variations in streamwise velocity and to enhance the laterally-divergent flow and lifted updrafts downstream of the microfront; (4) as the coherent mother downdraft impinges on the canopy, canopy-scale eddies are formed near the canopy top in a similar manner as observed in conventional mixing-layer turbulence.  相似文献   

10.
In nested nonhydrostatic mesoscale model MM5,the characteristic quantities of atmosphericturbulence,i.e.,the standard deviations of the turbulent fluctuated speeds for three directions inPBL are computed by Mellor-Yamada's level 2.5 closure scheme.The magnitudes and the verticalprofiles of these quantities computed from the model are closely connected with temperature andwind speed profiles as well as the type of the ground with a significant diurnal variation,and are inagreement with known magnitudes and regularities in different stratification conditions.Hence themethod in this paper is reasonable and convincible.Their horizontal distribution depends on thehorizontal distribution of the stratification.The method of predicted characteristic quantities ofturbulence from mesoscale model in this paper can be used in the problem of atmospheric diffusionand atmospheric environment.  相似文献   

11.
  相似文献   

12.
In nested nonhydrostatic mesoscale model MM5,the characteristic quantities of atmospheric turbulence,i.e.,the standard deviations of the turbulent fluctuated speeds for three directions in PBL are computed by Mellor-Yamada's level 2.5 closure scheme.The magnitudes and the vertical profiles of these quantities computed from the model are closely connected with temperature and wind speed profiles as well as the type of the ground with a significant diurnal variation,and are in agreement with known magnitudes and regularities in different stratification conditions.Hence the method in this paper is reasonable and convincible.Their horizontal distribution depends on the horizontal distribution of the stratification.The method of predicted characteristic quantities of turbulence from mesoscale model in this paper can be used in the problem of atmospheric diffusion and atmospheric environment.  相似文献   

13.
The structure of turbulent flows along a transition between tall-forested canopies and forest clearings continues to be an active research topic in canopy turbulence. The difficulties in describing the turbulent flow along these transitions stem from the fact that the vertical structure of the canopy and its leaf area distribution cannot be ignored or represented by an effective roughness length. Large-eddy simulation (LES) runs were performed to explore the effect of a homogeneous variation in the forest leaf area index (LAI) on the turbulent flow across forest edges. A nested grid numerical method was used to ensure the development of a deep boundary layer above the forest while maintaining a sufficiently high resolution in the region close to the ground. It was demonstrated that the LES here predicted first-order and second-order mean velocity statistics within the canopy that agree with reported Reynolds-Averaged Navier–Stokes (RANS) model results, field and laboratory experiments. In the simulations reported here, the LAI was varied between 2 and 8 spanning a broad range of observed LAI in terrestrial ecosystems. By increasing the forest LAI, the mean flow properties both within the forest and in the clearing near the forest edge were altered in two fundamental ways: near the forest edge and into the clearing, the flow statistical properties resembled the so-called back-facing step (BFS) flow with a mean recirculation zone near the edge. Another recirculation zone sets up downstream of the clearing as the flow enters the tall forest canopy. The genesis of this within-forest recirculation zone can be primarily described using the interplay between the mean pressure gradients (forcing the flow) and the drag force (opposing the flow). Using the LES results, a simplified analytical model was also proposed to explain the location of the recirculation zone inside the canopy and its dependence on the forest LAI. Furthermore, a simplified scaling argument that decomposes the mean velocity at the outflow edge into a superposition of ‘exit flow’ and BFS-like flow with their relative importance determined by LAI was explored.  相似文献   

14.
The phase-averaged characteristics of the turbulent velocity fields beneath steep short wind waves are investigated. A scheme was developed to compute the phase of individual wind waves using spatial surface displacement data. This information was used to analyze the two-dimensional velocity data acquired using particle image velocimetry (PIV) in a wind-wave tank. The experiments were conducted at a fetch of 5.5m and at wind speeds that ranged from 4 to 10ms−1. Under these conditions previous studies have shown that a significant percentage of the waves are microscale breaking waves. An analysis of the phase-averaged results suggests under these conditions (short fetches and moderate wind speeds) a wind-driven water surface can be divided into three regions based on the intensity of the turbulence. These are the crests of microscale breaking waves, the crests of non-breaking waves and the troughs of all waves. The turbulence is most intense beneath the crests of microscale breaking waves. In the crest region of microscale breaking waves coherent structures were observed that were stronger and occurred more frequently than beneath the crests of non-breaking waves. Beneath the crests of non-breaking waves the turbulence is a factor of two to three times weaker and beneath the wave troughs it is a factor of six weaker. These findings provide additional support for the hypothesis that approximately two-thirds of the gas and heat fluxes occur across the turbulent wakes produced by microscale breaking waves.  相似文献   

15.
A series of idealized model simulations are analyzed to determine the sensitivity of model results to different configurations of the lateral boundary conditions (LBCs) in simulating mesoscale shallow convection over hilly terrain, In the simulations with steady thermal forcing at the model surface, a radiation condition at both boundaries is the best choice under high wind conditions, and the best results are produced when both the normal velocities and the temperature are treated with the radiation scheme in which the phase speed is the same for different variables, When the background wind speed is reasonably small, the LBC configuration with either the radiation or the zero gradient condition at both boundaries tends to make the numerical solution unstable. The choice of a constant condition at the inflow boundary and a radiation outflow boundary condition is appropriate in most cases, In the simulations with diurnal thermal forcing at the model surface, different LBC schemes are combined together to reduce spurious signals induced by the outflow boundary, A specification inflow boundary condition, in which the velocity fields at the inflow boundary are provided using the time-dependent results of a simulation with periodic LBCs over a flat domain, is tested and the results indicate that the specification condition at the inflow boundary makes it possible to use a smaller model domain to obtain reasonable results. The model horizontal domain length should be greater than a critical length, which depends on the domain depth H and the angle between gravity wave phase lines and the vertical, An estimate of minimum domain length is given by [(H-zi)/πU]√N^2L2x-4π^2u^2, where N and U are the background stability and wind speed,respectively, Lx is the typical gravity wavelength scale, and zi is the convective boundary layer (CBL)depth.  相似文献   

16.
In weak wind stable conditions, eddy-correlation fluxes calculated using conventional averaging times of 5 min or longer to define the perturbations are severely contaminated by poorly sampled mesoscale motions. A method is developed to identify the averaging time for each individual data record that captures the turbulence while excluding most of the mesoscale motions. The method is based on multiresolution decomposition of the heat flux, and provides an objective procedure for selecting the averaging time for calculating eddy-correlation fluxes. Eddy-correlation data collected in weak turbulence conditions over grass, snow, a pine forest and the ocean are used to demonstrate the approach.When the small-scale turbulence and mesoscale motions are clearly separated by a gap region in the heat flux cospectra, the variable window width reduces the influence of nonstationarity by more effectively filtering out mesoscale motions compared to traditional methods using constant averaging time. For records where turbulence and mesoscale motions overlap in scale, the method is not well posed, although such records occur infrequently for our datasets. These ambiguous cases correspond to significant nonstationarity at scales that overlap with turbulence scales. The improved turbulence fluxes calculated with the proposed method are the appropriate fluxes for evaluating flux-gradient relationships and Monin–Obukov similarity theory for developing improved model parameterizations of turbulence for weakly turbulent flows  相似文献   

17.
Abstract

Numerical simulations of the severe squall line of 14 July 1987 are discussed within the context of semi‐Lagrangian and semi‐implicit integrations. The fully compressible non‐hydrostatic Euler set of equations constitutes the basic dynamical framework of the numerical model. With elementary precipitation physics and with a generalized treatment of lateral boundary conditions, nested integrations simulate the observed structure of the squall line. The numerical solution and observations show a well organized precipitation system including a mesoscale fast‐propagating prefrontal squall line and a slow‐propagating system moving with the synoptic wave. The prefrontal squall line is seen to be a manifestation of the organization of an inertia‐gravity wave and has characteristics of a wave‐CISK mechanism. Owing to the interaction between the upper jet stream dynamics and moisture, the prefrontal perturbation is initiated locally and subsequently irradiates away from its point source. The high computational efficiency and the accuracy of the model emphasize its potential and demonstrate its value as an interesting tool for mesoscale modelling.  相似文献   

18.
The wake characteristics of a wind turbine in a turbulent atmospheric boundary layer under different thermal stratifications are investigated by means of large-eddy simulation with the geophysical flow solver EULAG. The turbulent inflow is based on a method that imposes the spectral energy distribution of a neutral boundary-layer precursor simulation, the turbulence-preserving method. This method is extended herein to make it applicable for different thermal stratification regimes (convective, stable, neutral) by including suitable turbulence assumptions, which are deduced from velocity fields of a diurnal-cycle precursor simulation. The wind-turbine-wake characteristics derived from simulations that include the parametrization result in good agreement with diurnal-cycle-driven wind-turbine simulations. Furthermore, different levels of accuracy are tested in the parametrization assumptions, representing the thermal stratification. These range from three-dimensional matrices of the precursor-simulation wind field to individual values. The resulting wake characteristics are similar, even for the simplest parametrization set-up, making the diurnal-cycle precursor simulation non-essential for the wind-turbine simulations. Therefore, the proposed parametrization results in a computationally fast, simple, and efficient tool for analyzing the effects of different thermal stratifications on wind-turbine wakes by means of large-eddy simulation.  相似文献   

19.
Both observational and numerical studies of the convective boundary layer (CBL) have demonstrated that when surface heat fluxes are small and mean wind shear is strong, convective updrafts tend to organize into horizontal rolls aligned within 10–20\(^\circ \) of the geostrophic wind direction. However, under large surface heat fluxes and weak to negligible shear, convection tends to organize into open cells, similar to turbulent Rayleigh-Bénard convection. Using a suite of 14 large-eddy simulations (LES) spanning a range of \(-z_i/L\) between zero (neutral) and 1041 (highly convective), where \(z_i\) is the CBL depth and L is the Obukhov length, the transition between roll- and cellular-type convection is investigated systematically for the first time using LES. Mean vertical profiles including velocity variances and turbulent transport efficiencies, as well the “roll factor,” which characterizes the rotational symmetry of the vertical velocity field, indicate the transition occurs gradually over a range of \(-z_i/L\); however, the most significant changes in vertical profiles and CBL organization occur from near-neutral conditions up to about \(-z_i/L \approx \) 15–20. Turbulent transport efficiencies and quadrant analysis are used to characterize the turbulent transport of momentum and heat with increasing \(-z_i/L\). It is found that turbulence transports heat efficiently from weakly to highly convective conditions; however, turbulent momentum transport becomes increasingly inefficient as \(-z_i/L\) increases.  相似文献   

20.
Summary The dynamical effect of land surface heterogeneity on heat fluxes in the atmospheric boundary layer (ABL) is investigated using numerical simulations with a non-hydrostatic model over a wide range of grid resolutions. It is commonly assumed that mesoscale or dynamical fluxes associated with mesoscale and convective circulations simulated by a high-resolution model (subgrid (SG) model) on the subgrid scale of a climate model (large-scale (LS) model) represent additional processes in the ABL, which are not considered by the turbulence scheme of the LS-model, and which can be parameterized using the SG-model. The present study investigates the usefulness of this methodology for small-scale and large-scale idealized heterogeneities using a SG-model resolving mesoscale or even microscale circulations to compute the mesoscale fluxes on the scale of the LS-model. It is shown that the dynamical transports as derived from the SG-model should not be used to correct the parameterized turbulent fluxes of the LS-model. The reason is that the subgrid circulations simulated by the SG-model interact with the fields of wind and scalars in the ABL, which results in reduced turbulent fluxes in the ABL. Thus the methodology of previous studies to use mesoscale/dynamical fluxes for the correction of flux profiles simulated by climate models seems to be questionable.  相似文献   

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

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