首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The Regional Atmospheric Modeling System （RAMS）, which is a non-hydrostatic numerical model, has been used to investigate the impact of terrain shape and large-scale forcing on the Antarctic surface-wind regime, focusing on their roles in establishing favorable flow conditions for the formation of katabatic flow jumps. A series of quasi-2D numerical simulations were conducted over idealized slopes representing the slopes of Antarctica during austral winter conditions. Results indicate that the steepness and variations of the underlying slope play a role in the evolution of near-surface flows and thus the formation of katabatic flow jumps. However, large-scale forcing has a more noticeable effect on the occurrence of this small-scale phenomenon by establishing essential upstream and downstream flow conditions, including the upstream supercritical flow, the less stably stratified or unstable layer above the cold katabatic layer, as well as the cold-air pool located near the foot of the slope through an interaction with the underlying topography. Thus, the areas with steep and abrupt change in slopes, e.g. near the coastal areas of the eastern Antarctic, are preferred locations for the occurrence of katabatic flow jumps, especially under supporting synoptic conditions.  相似文献   

2.
We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin–Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.  相似文献   

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

4.
Second-Order Modelling of Turbulence in Katabatic Flows   总被引：1，自引：1，他引：0  
A complete one-dimensional second-order closure model is used to simulate katabatic flows observed on glaciers and ice caps. The model is tested with two different closure assumptions for the viscous dissipation, one based on a prognostic equation for and the other on a diagnostic buoyant length scale. Both formulations give quite similar results. Model simulations are compared to observations made over sloping ice surfaces during periods dominated by katabatic flow. In general, good agreement is found for both mean wind and temperature profiles as well as eddy correlation measurements. It is also found that the turbulent transport terms play an important role in katabatic flows as opposed to the classical stable boundary layer where these terms are usually ignored. Even the turbulent transport of temperature variance, which leads to the well-known countergradient term in unstable boundary layers, is relatively important for modelling the observed temperature profiles. The effect of these terms on the flux-profile relationships, using observed and simulated profiles, is also discussed.  相似文献   

5.
A non-hydrostatic numerical model, the Regional Atmospheric Modeling System (RAMS), has been used to investigate the development of katabatic jumps in Coats Land, Antarctica. In the control run with a 5 m s-1downslope directed initial wind, a katabatic jump develops near the foot of the idealized slope. The jump is manifested as a rapid deceleration of the downslope flow and a change from supercritical to subcritical flow, in a hydraulic sense, i.e., the Froude number (Fr) of the flow changes from Fr > 1 to Fr> 1. Results from sensitivity experiments show that an increase in the upstream flow rate strengthens the jump, while an increase in the downstream inversion-layer depth results in a retreat of the jump. Hydraulic theory and Bernoulli's theorem have been used to explain the surface pressure change across the jump. It is found that hydraulic theory always underestimates the surface pressure change, while Bernoulli's theorem provides a satisfactory estimation. An analysis of the downs balance for the katabatic jump indicates that the important forces are those related to the pressure gradient, advection and, to a lesser extent, the turbulent momentum divergence. The development of katabatic jumps can be divided into two phases. In phase I, the t gradient force is nearly balanced by advection, while in phase II, the pressure gradient force is counterbalanced by turbulent momentum divergence. The upslope pressure gradient force associated with a pool of cold air over the ice shelf facilitates the formation of the katabatic jump.  相似文献   

6.
For the first time, results from a high-resolution numerical simulation (with horizontal grid spacing of 35m) were used to reveal the detailed structure near an atmospheric katabatic jump over an idealized slope. The simulation represents flow over the slopes of Coats Land, Antarctica for austral winter conditions. The katabatic jump is characterised by an updraft with vertical velocities of order 1ms−1 and serves as a possible forcing mechanism for the gravity waves frequently observed over the ice shelves around the Antarctic. Results also indicate that strong turbulence is generally confined within a mixing zone near the top of the katabatic layer upstream of the jump and extends downstream through the top of the strong updraft associated with the jump. Detailed analyses of momentum and heat budgets across the katabatic jump indicate that, upstream of the jump, turbulent mixing is important in decelerating the upper part of the katabatic layer, while within the jump the upslope pressure gradient force associated with the pool of cold air plays a role in decelerating the flow near the surface. The heat budget near the jump reveals a simple two-term balance: the turbulent heat flux divergence is balanced by the advection. A comparison of model results with available theories indicates that mixing between layers of different potential temperature structure indeed plays some role in the development of katabatic flow jumps, especially for strong jumps. Theories used to study katabatic jumps should include this mixing process, of which the amount depends on the intensity of the jump. A conceptual model of a katabatic jump, including the main dynamical processes, is constructed from these detailed analyses.  相似文献   

7.
Large-eddy simulation (LES) of a stable atmospheric boundary layer is performed using recently developed dynamic subgrid-scale (SGS) models. These models not only calculate the Smagorinsky coefficient and SGS Prandtl number dynamically based on the smallest resolved motions in the flow, they also allow for scale dependence of those coefficients. This dynamic calculation requires statistical averaging for numerical stability. Here, we evaluate three commonly used averaging schemes in stable atmospheric boundary-layer simulations: averaging over horizontal planes, over adjacent grid points, and following fluid particle trajectories. Particular attention is focused on assessing the effect of the different averaging methods on resolved flow statistics and SGS model coefficients. Our results indicate that averaging schemes that allow the coefficients to fluctuate locally give results that are in better agreement with boundary-layer similarity theory and previous LES studies. Even among models that are local, the averaging method is found to affect model coefficient probability density function distributions and turbulent spectra of the resolved velocity and temperature fields. Overall, averaging along fluid pathlines is found to produce the best combination of self consistent model coefficients, first- and second-order flow statistics and insensitivity to grid resolution.  相似文献   

8.
Impacts of different terrain configurations on the general behaviour of idealised katabatic flows are investigated in a numerical model study. Various simplified terrain models are applied to unveil modifications of the dynamics of nocturnal cold drainage of air as a result of predefined topographical structures. The generated idealised terrain models encompass all major topographical elements of an area in the tropical eastern Andes of southern Ecuador and northern Peru, and the adjacent Amazon. The idealised simulations corroborate that (i) katabatic flows develop over topographical elements (slopes and valleys), that (ii) confluence of katabatic flows in a lowland basin with a concave terrainline occur, and (iii) a complex drainage flow system regime directed into such a basin can sustain the confluence despite varying slope angles and slope distances.  相似文献   

9.
A realistic simulation of katabatic flows is not a straightforward task for numerical models. One complicating factor is that katabatic flows develop within a stably stratified boundary layer, which is poorly resolved and described in many numerical models. To capture the jet-shaped shallow flow a model set-up with high vertical resolution is also required. In this study, ‘a state of the art’ mesoscale numerical model is applied in a simulation of katabatic flow over a melting glacier. A basic agreement between observations and model results is found. From scale analysis, it is concluded that the simulated flow can be classified as katabatic. Although the background flow varies in strength and direction, the simulated katabatic flow over Breidamerkurjökull is persistent. Two factors vital for this persistence are identified. First, the melting snow maintains the surface temperature close to 0 °C while the air temperature warms adiabatically as it descends the slope. This provides a ‘self enhanced’ negative buoyancy that drives the flow to a balance with local friction. Second, the jet-like shape of the resulting flow gives rise to a large ‘curvature term’ in the Scorer parameter, which becomes negative in the upper jet. This prevents vertical wave propagation and isolates the katabatic layer of the influence from the free troposphere aloft. Our results suggest that the formation of local microclimates dominated by katabatic flow is a general feature over melting glaciers. The modelled turbulence structure illustrates the importance of non-local processes. Neglecting the vertical transport of turbulence in katabatic flows is not a valid assumption. It is also found that the local friction velocity remains larger than zero through the katabatic jet, due to directional shear where the scalar wind speed approaches its maximum.  相似文献   

10.
Large-Eddy Simulation of Windbreak Flow   总被引：13，自引：10，他引：3  
A large-eddy simulation has been performed of turbulent flow around multiple windbreaks set within a wheat canopy under neutral stability conditions. The simulation is validated against a wind tunnel data set taken under similar conditions. Velocity profiles and second-order statistics are presented and compared to those found in the wind tunnel. From the numerical simulation, we discuss spatial distributions of instantaneous velocity fields and pressure statistics, which are important and telling features of the flow that are difficult to measure experimentally. We present a discussion of the momentum balance at various locations with respect to the windbreak, and similarly, we introduce the budget of a passive scalar. These discussions show the importance of the terms in each budget equation as they vary upstream and downstream of the windbreak.  相似文献   

11.
Second-order closure models for the canopy sublayer (CSL) employ aset of closure schemes developed for `free-air' flow equations andthen add extra terms to account for canopy related processes. Muchof the current research thrust in CSL closure has focused on thesecanopy modifications. Instead of offering new closure formulationshere, we propose a new mixing length model that accounts for basicenergetic modes within the CSL. Detailed flume experiments withcylindrical rods in dense arrays to represent a rigid canopy areconducted to test the closure model. We show that when this lengthscale model is combined with standard second-order closureschemes, first and second moments, triple velocity correlations,the mean turbulent kinetic energy dissipation rate, and the wakeproduction are all well reproduced within the CSL provided thedrag coefficient (CD) is well parameterized. The maintheoretical novelty here is the analytical linkage betweengradient-diffusion closure schemes for the triple velocitycorrelation and non-local momentum transfer via cumulant expansionmethods. We showed that second-order closure models reproducereasonably well the relative importance of ejections and sweeps onmomentum transfer despite their local closure approximations.Hence, it is demonstrated that for simple canopy morphology (e.g.,cylindrical rods) with well-defined length scales, standard closureschemes can reproduce key flow statistics without much revision.When all these results are taken together, it appears that thepredictive skills of second-order closure models are not limitedby closure formulations; rather, they are limited by our abilityto independently connect the drag coefficient and the effectivemixing length to the canopy roughness density. With rapidadvancements in laser altimetry, the canopy roughness densitydistribution will become available for many terrestrialecosystems. Quantifying the sheltering effect, the homogeneity andisotropy of the drag coefficient, and more importantly, thecanonical mixing length, for such variable roughness density isstill lacking.  相似文献   

12.
The flow and turbulence quantities governing dispersion in katabatic flows vary with both height and downslope distance. This variation cannot be accounted for in conventional plume dispersion models. In this study, three random-walk models of varying complexity are formulated to simulate dispersion in katabatic flows, and their strengths and weaknesses are discussed. The flow and turbulence parameters required by these models are determined from a high-resolution two-dimensional katabatic flow model based on a turbulent kinetic energy closure. Random-walk model calculations have been performed for several values of source height and slope angle to examine the influence of these parameters on dispersion. Finally, we simulated the perfluorocarbon and heavy methane tracer releases for Night 4 of the 1980 ASCOT field study over a nearly two-dimensional slope in Anderson Creek Valley, California. The observed peak concentrations are generally well-predicted. The effects of the pooling of the drainage air could not be taken into account in our katabatic flow model and, consequently, the predicted concentrations decay much more rapidly with time than the observed values.  相似文献   

13.
The differential equations for first-order (linear) response of the planetary boundary layer are formulated for flow over periodic terrain, for variations in both surface roughness and terrain elevation. A simple second-order closure model of the turbulence is used, and Coriolis forces are neglected. Flow over a periodic terrain produces corresponding periodic structure in all meteorological fields above the surface. The periodic structure consists of two components. The first is very nearly evanescent with height. It corresponds to the motion that would be observed were the atmosphere inviscid. The second component, introduced by turbulent viscosity, exhibits a phase variation with height in addition to a decay in amplitude. W.K.B. [Wentzel-Kramers-Brillouin] approximations for the two components are developed, and the coupling between the components is discussed. The formulation for calculating solutions by numerical integration is developed, including specification of appropriate boundary conditions. Calculations are presented in a companion paper.Wave Propagation Laboratory.Environmental Science Group.  相似文献   

14.
A nonlinear numerical model is developed for turbulent boundary-layer flowover a train of water waves of finite amplitude or slope. The airflow isassumed to be steady, two-dimensional, and neutrally-stratified. The wavesurface is assumed to be aerodynamically rough and flow conditions at thewave surface are prescribed. The numerical model used in this study adoptsthree turbulence closure schemes with different degrees of physicalcompleteness. Two of these are second-order schemes, whichare believed to describe turbulent flow more completely than thesimpler closures used in most previous studies. Although models with all turbulence closures agree qualitatively in the prediction of the amplitude of the surface normal stress perturbation, the lower- and higher-order closures differ significantly in predictions of phase, and hence the form drag and energy transfer rate between wind and waves. Our model results are in reasonable agreement with field and laboratory measurements, although predicted energy transfer rates are generally at the low end of the range of experimental values. Cases with airflow at various angles to the wave direction are also considered.  相似文献   

15.
Spectral Maxima In A Perturbed Stable Boundary Layer   总被引：4，自引：1，他引：4  
Wind velocity data have been collected on Nansen Ice Sheet, Antarctica, close to the base of a steeply sloping glacier along which frequently flow katabatic winds. The aim of this study is to investigate how turbulent energy and momentum flux are perturbed by the flow interaction with topography and by the strong mechanical mixing produced by downslope flows. Spectral and cospectral analyses, performed on the wind velocity components, provide evidence that such a perturbation, at any stability, is restricted to frequencies lower than the inertial subrange. Longitudinal spectra display an energy increment, due to turbulence generated by topography and by mechanical forcing related to the katabatic wind structure. The energy, supplied by the topographic forcing, displaces the turbulent energy maximum toward lower frequencies. In near-neutral stratification the spectral maximum occurs at a reduced frequency, which seems to be consistent with the height of the steepest part of the slope, and seems to shift toward higher frequencies as a linear ,function of the local stability parameter,Ll. The parameterisation of the orographic perturbation by means of a similarity relationship allows us to scale u spectra in the same way as over uniform terrain. The scaled, perturbed spectra collapse onto a unique curve in the mid-frequency as well in the inertial subrange, while maxima are grouped in a cluster. Lateral and vertical velocity spectra exhibit shapes independent of stability, suggesting a topographic perturbation that is predominantly over stability effects.  相似文献   

16.
This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.  相似文献   

17.
Atmospheric flow over complex terrain, particularly recirculation flows, greatly influences wind-turbine siting, forest-fire behaviour, and trace-gas and pollutant dispersion. However, there is a large uncertainty in the simulation of flow over complex topography, which is attributable to the type of turbulence model, the subgrid-scale (SGS) turbulence parametrization, terrain-following coordinates, and numerical errors in finite-difference methods. Here, we upgrade the large-eddy simulation module within the Weather Research and Forecasting model by incorporating the immersed-boundary method into the module to improve simulations of the flow and recirculation over complex terrain. Simulations over the Bolund Hill indicate improved mean absolute speed-up errors with respect to previous studies, as well an improved simulation of the recirculation zone behind the escarpment of the hill. With regard to the SGS parametrization, the Lagrangian-averaged scale-dependent Smagorinsky model performs better than the classic Smagorinsky model in reproducing both velocity and turbulent kinetic energy. A finer grid resolution also improves the strength of the recirculation in flow simulations, with a higher horizontal grid resolution improving simulations just behind the escarpment, and a higher vertical grid resolution improving results on the lee side of the hill. Our modelling approach has broad applications for the simulation of atmospheric flows over complex topography.  相似文献   

18.
We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.  相似文献   

19.
We develop a parameterisation for the effective roughness length of terrain that consists of a repeating sequence of patches, in which each patch is composed of strips of two roughness types. A numerical model with second-order closure in the turbulent stress is developed and used to show that: (i) the normalised Reynolds stress develops as a self-similar profile; (ii) the mixing-length parameterisation is a good first-order approximation to the Reynolds stress. These findings are used to characterise the blending layer, where the stress adjusts smoothly from its local surface value to its effective value aloft. Previous studies have assumed that this adjustment occurs abruptly at a single level, often called the blending height. The blending layer is shown to be characterised by height scales that arise naturally in linear models of surface layer flow over roughness changes, and calculations with the numerical model show that these height scales remain appropriate in the nonlinear regime. This concept of the blending layer allows the development of a new parameterisation of the effective roughness length, which gives values for the effective roughness length that are shown to compare well with both atmospheric measurements and values determined from the second-order model.  相似文献   

20.
Summary The Adélie Land coastal section of East Antarctica is known for strong katabatic winds. Although the primary forcing of these persistent drainage flows has been attributed to the radiative cooling of the sloping ice topography, effects of ambient horizontal pressure gradients can play a central role in shaping the Antarctic surface wind regime as well. Oberrvations of the katabatic wind at the near-coastal Adélie Land station D-10 have been sorted into strong and weak wind classes. Concurrent radiosonde ascents at nearby Dumont D'Urville have been used to depict the timeaveraged large scale conditions accompanying the katabatic wind classes. Results suggest that strong katabatic wind cases are associated with low pressure over the coastal margin and easterly upper level motions. Numerical simulations have been conducted to examine the effect of prescribed large scale forcing on the evolution of the katabatic wind. The model runs indicate that the ambient environment plays a key role in the development and intensity of the katabatic wind regime.With 7 Figures  相似文献   

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

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