共查询到20条相似文献，搜索用时 359 毫秒
1.
Given incident logarithmic profiles of wind and pollutant concentration above a rough, absorbing surface, the three-dimensional distribution of pollutant concentration over a hill of gentle slope is determined from a linearized model. The model is applied in neutrally stratified flow, without chemistry, and is integrated using spectral methods in the horizontal and a finite-difference scheme in the vertical. This approach allows for flexibility in choosing a closure scheme and a variety of surface boundary conditions. This was not possible in the analytic approach of Padro (1987) who added pollutant concentration and flux to the MS3DJH/1 model of Walmsley et al. (1980). The present model requires as input the turbulent kinetic energy, E, dissipation, , and the perturbation vertical velocity, w, from the three-dimensional boundary-layer flow model of Beljaars et al. (1987), hereinafter referred to as MSFD, The latter model also supplies wind velocity perturbations at the upper boundary, as input to upper boundary conditions on the pollutant flux perturbations.The present study describes applications of the model to idealized terrain features: isolated two- and three-dimensional hills and ridges and an infinite series of ridges. (Application to real terrain, however, presents no difficulties.) Comparisons were made with different (though uniform) surface roughnesses. Tests were performed to examine the effect of upstream terrain features in the periodic domain and to illustrate the importance of the vertical resolution of the output for interpreting results from the sinusoidal terrain case. 相似文献
2.
Keith W. Ayotte 《Boundary-Layer Meteorology》1997,83(2):285-309
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. 相似文献
3.
4.
5.
Prof. Dr. Z. Petkovšek 《Meteorology and Atmospheric Physics》1990,43(1-4):197-202
Summary For the inclined stationary internal atmospheric boundary as is the upper boundary of the bora flow down the lee-side slope of an orographic obstacle, the dynamic forces must be appropriately balanced, similar as at the atmospheric frontal surfaces. Starting from this idea, a model for the bora flow conditions was developed. Introducing some special suppositions, the model finally consists of four partial differential equations, which are solved numerically.Numerical experiments with the model by changing parameters such as friction and warm air stratification, as well as initial values of the flow: thickness, velocity, inclination and air mass temperature difference, show interesting results and finings. The shape of upper bora boundary and the distribution of the bora velocity above and downward the slope are graphycally presented, as they develop under the influence of different values of the mentioned parameters.With 5 Figures 相似文献
6.
S. R. Karpik 《Boundary-Layer Meteorology》1988,43(3):273-286
In their Mixed Spectral Finite Difference (MSFD) model for flow over complex terrain, Beljaars et al. (1987) solve a set of coupled, second-order ordinary differential equations (ODEs) for the first-order perturbations to the logarithmic velocity profile caused by nonuniform surface roughness and topography. To solve this set of ODEs, they employ a Forward Euler Shooting Method. It is demonstrated here that the shooting method is computationally unstable for this problem. An absolutely stable finite-difference method based on a block tridiagonal LU factorization of the finite-difference matrix is presented. The advantages of the present algorithm over the method used by Beljaars et al. are demonstrated both by theoretical argument and numerical experiment. 相似文献
7.
Roy M. Endlich 《Boundary-Layer Meteorology》1984,30(1-4):375-386
A diagnostic model is a relatively simple and practical tool for modeling the wind flow of the boundary layer in complex terrain. The model begins with a wind analysis based on available surface wind reports and geostrophic winds (computed from pressure data). The height of the boundary layer top (upper surface of the computational domain) is prescribed to fit local conditions. Using the continuity equation in terrain-following coordinates, the winds at mesh points are adjusted to produce nondivergence while maintaining the original vertical component of vorticity. The method of computing the nondivergent winds uses direct alterations. This method may be useful for other modeling purposes and will be described. Data for a long period (usually a year) are analyzed to obtain eigenvectors and the associated time series of their coefficients at each observation time. The model is run only for the five or six eigenvectors that explain most of the variance. The wind field at any particular time is reconstructed from the eigenvector solutions and their appropriate coefficients. Comparisons of model results with measured winds at sites representing different types of terrain will be shown. The accuracy and economy of the model make it a useful tool for estimating wind energy and also for giving wind fields for low-level diffusion models. 相似文献
8.
Daniel Reinert Volkmar Wirth Joachim Eichhorn Walter-Georg Panhans 《Boundary-Layer Meteorology》2007,125(1):109-132
This paper presents the dry version of a new large-eddy simulation (LES) model, which is designed to simulate air flow and
clouds above highly complex terrain. The model is three-dimensional and nonhydrostatic, and the governing equations are sound
filtered by use of the anelastic approximation. A fractional step method is applied to solve the equations on a staggered
Cartesian grid. Arbitrarily steep and complex orography can be accounted for through the method of viscous topography. The
dynamical model core is validated by comparing the results for a spreading density current against a benchmark solution. The
model accuracy is further assessed through the simulation of turbulent flow across a quasi two-dimensional ridge. The results
are compared with wind-tunnel data. The method of viscous topography is not restricted to moderately sloped terrain. Compared
to models using curvilinear grids, it allows this model to be applied to a much wider range of flows. This is illustrated
through the simulation of an atmospheric boundary-layer flow over a surface mounted cube. The results show that the dry model
version is able to accurately represent the complex flow in the vicinity of three-dimensional obstacles. It is concluded that
the method of viscous topography was successfully implemented into a micrometeorological LES model. As will be shown in Part
II, this allows the detailed study of clouds in highly complex terrain. 相似文献
9.
Stability Dependence of Height Scales and Effective Roughness Lengths of Momentum and Heat Transfer Over Roughness Changes 总被引：1，自引：0，他引：1
It is widely accepted that the correct formulation of an effective roughness length, defined as the area average of the roughness length in heterogeneous terrain, relies upon the appropriate de-termination of a height scale. At this height a meteorological quantity is approximately in equilibrium with local surface conditions and independent of horizontal position. This research note determines explicitly the different height scales from the perturbation solutions of flow velocity and temperature, as well as the fluxes of momentum and heat, in a stratified boundary layer. These solutions are derived from the asymptotic approximation theory and shown to capture major characteristics of momentum and heat transfer over heterogeneous terrain with changes of the underlying roughness lengths. The effective roughness lengths can then be computed by use of these height scales. The dependence of height scales and effective roughness lengths upon stratification is also discussed briefly. 相似文献
10.
Alessandro Parente Catherine Gorlé Jeroen van Beeck Carlo Benocci 《Boundary-Layer Meteorology》2011,140(3):411-428
We report on a novel approach for the Reynolds-averaged Navier-Stokes (RANS) modelling of the neutral atmospheric boundary
layer (ABL), using the standard k-ek-{\varepsilon} turbulence model. A new inlet condition for turbulent kinetic energy is analytically derived from the solution of the k-ek-{\varepsilon} model transport equations, resulting in a consistent set of fully developed inlet conditions for the neutral ABL. A modification
of the standard k-ek-{\varepsilon} model is also employed to ensure consistency between the inlet conditions and the turbulence model. In particular, the turbulence
model constant C
μ
is generalized as a location-dependent parameter, and a source term is introduced in the transport equation for the turbulent
dissipation rate. The application of the proposed methodology to cases involving obstacles in the flow is made possible through
the implementation of an algorithm, which automatically switches the turbulence model formulation when going from the region
where the ABL is undisturbed to the region directly affected by the building. Finally, the model is completed with a slightly
modified version of the Richards and Hoxey rough-wall boundary condition. The methodology is implemented and tested in the
commercial code Ansys Fluent 12.1. Results are presented for a neutral boundary layer over flat terrain and for the flow around
a single building immersed in an ABL. 相似文献
11.
进一步讨论了有关非线性不稳定的一些问题,其主要内容有： １．考察了有代表性的三类发展方程,指出其对应的差分格式是否出现非线性计算不稳定,与原微分方程解的性质密切相关。 ２．进一步讨论了带周期边条件的守恒型差分格式的非线性计算稳定性问题,总结了克服非线性不稳定的有效措施。 ３．以非线性平流方程为例,着重分析了带非周期边条件的非守恒差分格式的非线性计算稳定性问题,给出了判别其计算稳定性的“综合分析判别法”。 相似文献
12.
Jacob Padro 《Boundary-Layer Meteorology》1987,38(1-2):17-28
Given the incident profiles of wind velocity and pollutant concentration, we seek to determine the 3-dimensional concentration field of a pollutant upon a region with complex terrain. The analytic solution of the wind velocity in a 3-dimensional boundary-layer model by Walmsley et al. (1980) is utilized as a forcing function in the simplified concentration perturbation equation for a pollutant. The resulting solution applies to an isolated cosine-squared hill in a neutrally stratified boundary-layer flow with a surface type which absorbs the pollutant totally. The solution shows that the concentration perturbation field is organized in accordance with the wind field. In particular, the east-west cross-section is 180° out of phase with the velocity perturbation field. The vertical profiles of the concentration perturbations for selected grid points approach the value of the upper boundary condition very rapidly. 相似文献
13.
14.
地形对边界层影响的有限元数值模式 总被引：2，自引：0，他引：2
本文在定常和中性层结条件下应用有限元法初步建立了一个能够处理二线复杂地形上边界层问题的数值模式。模式使八节点抛物型等参元,有限元方程组用“波前法”求解。应用模式计算了几种典型地形上的近地层流动。模拟出了风速、气压、湍流动量通量以及湍流交换系数等气象要素的主要分布特征,并将计算结果与有关的观测事实做了比较。 相似文献
15.
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. 相似文献
16.
André Tremblay 《大气与海洋》2013,51(3):567-603
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. 相似文献
17.
小尺度地形引起的切变重力波 总被引：2，自引：0，他引：2
在大气边界层中,特别是山区的边界层中,经常可以观测到重力波的活动。例如Bull和Neisser分析了1974年4月至7月的2300时的微气压计记录。他们发现约38％的观测时间里都存在振幅至少为10—20 dPa的重力波活动,其中大多数重力波的振幅为25—70 dPa。由于风速切变和温度层结的不同,重力波引起的垂直位移的大小可能为百米,也可能达到整个夜间边界层厚度的量级。 很多大气活动都可能引起重力波,例如锋面、低空急流、局地风切变、逆温及强对流等等。地形特别是在夜间也是产生重力波的重要来源。当稳定度随高度急剧减少或风速随高度急剧增加时,在山脊的背风坡经常可以发现一系列的波动活动。 相似文献
18.
复杂地形下北京雷暴新生地点变化的加密观测研究 总被引：3，自引：1，他引：2
2008 年8 月14 日北京发生了雷暴群形式的局地暴雨,雷暴新生地点复杂多变,形成了多个γ 中尺度的强降水中心。本文利用近年来北京气象现代化建设取得的加密地面自动站、多普勒雷达、风廓线仪、微波辐射计等多种新型高时空分辨率观测资料及雷达四维变分同化系统（VDRAS）反演资料,通过精细分析地面（边界层）风场、温度场等的演变特征,讨论了雷暴新生地点变化的机制。结果表明:复杂地形与雷暴冷池出流作用相结合,主导了雷暴新生地点的变化,进而影响γ 中尺度强降水中心的变化;天气尺度高低空涡、槽的配合不一致,并且系统移动缓慢,以及对流层低层的弱的环境垂直风切变,是雷暴冷池结合复杂地形发挥雷暴新生地点主导作用的重要前提;复杂地形使得冷空气在一定范围内流动,在边界层产生碰撞和辐合,起到触发和增强对流作用,并使得对流风暴的形态和走向与地形呈现出紧密相关性;一定强度的冷池出流、边界层前期的暖湿空气和对流不稳定能量的积累,是冷池出流触发雷暴新生和演变的必要条件;北京周边地区的雷暴,通过其雷暴冷池出流沿着沟谷地形或向平原地区流动,与北京山谷或城区的边界层暖湿空气形成辐合抬升机制,触发雷暴新生。 相似文献
19.
1.IntroductionUnsaturatedsoilwaterflowisaflowwherewaterisnotfullofsoilhole,whichisanimportantformofflowinporousmedia.Predictionofanunsaturatedflowisprovidedwithsignificanceinmanybranchesofscienceandengineering.Theseincludeatmosphericscience,soilscien... 相似文献
20.
This paper analyses a western boundary current striking a solid boundary. Interest is concentrated on the case where inertial effects are sufficient to modify the flow from its ‘Stommel-layer’ form with Ekman friction relatively unimportant in the interior. It is shown that, beyond a critical inflow speed, a complicated system of four asymptotic regions forms near the corner, turning some of the flow along the blocking boundary and then returning it westwards to rejoin the western boundary current. A comparison with the results of a simple ocean model shows that many of the features in that flow can be explained through the asymptotic theory. 相似文献