共查询到20条相似文献，搜索用时 421 毫秒
1.
《Dynamics of Atmospheres and Oceans》1997,26(1):1-25
Principles for incorporating the upstream effects of deep sills into numerical ocean circulation models using nonlinear analytical hydraulic models are discussed within the context of reduced gravity flow. A method is developed allowing the upstream influence of a numerically unresolvable deep sill or width contraction to be reproduced. The method consists of placing an artificial boundary in the numerical model's overflowing layer at some distance upstream of the actual sill or width contraction of the deep strait. Given the model state at time t, the dependent flow variables are then predicted at the artificial boundary at time t + Δt by using the method of characteristics in combination with quasi-steady hydraulic laws. The calculation requires the use of Riemann invariants and examples are given for a simple nonrotating flow and for rotating channel flow with uniform potential vorticity. The computation is considerably simplified by linearizing the relevant equations in the vicinity of the artificial boundary, resulting in a linear wave reflection problem. The reflection coefficients for the two cases are calculated and these can be used directly to numerically satisfy the boundary condition in a straightforward way. 相似文献
2.
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. 相似文献
3.
Modeling nonhydrostatic atmospheric flow requires the solution of the vertical equation of motion and a prognostic or diagnostic equation for pressure. If the nonhydrostatic components of the flow are relatively small, they can be approximated and incorporated into a purely hydrostatic model, which usually is conceptually simpler and computationally more efficient. A method to do this for a linear model of local thermally-induced circulations is further developed and adapted to a non-linear numerical model of the neutral atmospheric boundary layer. A hydrostatic model and the quasi-nonhydrostatic version were used to simulate neutral flow over simple terrain features. One set of observations taken over a simple change in roughness and another set taken over a change in both roughness and terrain were simulated by both models to assess the capabilities of the quasi-nonhydrostatic technique.It is found that (as expected) the pressure deviation from the hydrostatic state is negligible for the roughness change, but it is an important aspect of neutral flow over terrain. Thus, for flow encountering a simple roughness change, the hydrostatic approximation is good, even for small horizontal scales. However, the quasi-nonhydrostatic model qualitatively produces the features in the observations for flow over a terrain change that the hydrostatic model cannot produce.Journal Paper No. J-12737 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa. Project No. 2779. 相似文献
4.
5.
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. 相似文献
6.
A two-dimensional turbulent diffusion equation is derived in a streamline coordinate system, defined for rotational flow over complex terrain and limited aloft by an elevated, impenetrable inversion. In the first instance, the steady-state equation is solved for an inner region of the boundary layer, in which the effect of curvature is negligible and, for simplicity, it is assumed that vorticity has a power-law dependence upon stream function. A variational method of solution is also discussed, in which vorticity may have a more general representation. A numerical calculation is performed for a special case of symmetrical flow over an isolated hill. The dependence of pollutant concentration upon the flow field, downwind distance and source is examined and the effect of wind acceleration in the neighbourhood of the top of the hill is discussed. It is pointed out that the diffusion model can be applied to any realistic flow field, provided that the streamlines are specified. 相似文献
7.
1.IntroductionUnsaturatedsoilwaterflowisaflowwherewaterisnotfullofsoilhole,whichisanimportantformofflowinporousmedia.Predictionofanunsaturatedflowisprovidedwithsignificanceinmanybranchesofscienceandengineering.Theseincludeatmosphericscience,soilscien... 相似文献
8.
复杂地形风场的精细数值模拟 总被引：1，自引：0，他引：1
风能是一种重要气候资源,随着我国风电规模的迅速增大,发展风能资源评估系统和风功率预测系统已成为一项重要的研究内容。国内外对复杂地形风场结构的数值模拟有大量研究,随着计算机能力增强,以往用于空气动力学精细流场计算的计算流体力学(Computational Fluid Dynamics,CFD)模式越来越多地在气象领域得到应用,人们开始研究用中尺度预报模式和CFD模式结合进行复杂地形风场的数值模拟。本文的耦合模式系统采用中尺度气象模式(WRF),通过嵌套网格到内层尺度(一般是几公里),然后通过耦合CFD模式Fluent软件获得高分辨率(水平30~100 m,垂直150 m高度以下10 m)的风速分布资料,得到精细化的风场信息。通过对鄱阳湖北部区域和云南杨梅山复杂地形的风场模拟,提供了风能评估和预报的一种可行的方法。 相似文献
9.
基于CFD技术的陡峭山体风场模拟方法研究 总被引：1，自引：0，他引：1
为更好将计算流体力学(CFD)技术应用于陡峭山体的风场模拟研究,文中对其中两个关键技术问题进行了研究,一是适宜的CFD解域顶高问题,二是入流边界条件设置方法问题。研究采用了香港国际机场的多普勒雷达径向速度观测资料,通过CFD模拟结果与观测结果对比验证展开研究。对于解域顶高,一般环境流体力学或建筑风工程领域要求顶高越高越好,例如5倍于地面障碍物的垂直高度,究其原因,主要是为了获得地面障碍物周边的整个流场结构,避免顶边界过低影响障碍物顶部流场的准确性。通过数值试验发现,若模拟关心的区域为较低海拔高度时,不必完全拘泥于该要求,对一座高度为1000m量级的山体进行模拟,顶高为3000和6000m的两组试验在边界层中低部给出的模拟结果没有显著差别。进行了3组边界条件设置方法的对比试验,结果表明单纯采用廓线法不足以充分描述来流信息,二维插值法尽管提供了更高分辨率的边界数据,但其结果总体上逊于0维插值法。研究表明,CFD可以很好地描述山地激发的涡旋/波动脱体运动,比单纯采用中尺度模式效果要好得多。在利用CFD工具研究复杂山地风场时,应本着"实用主义"的精神设置物理模型、参数和选取适当的方法、边界条件,以达到模拟精度和计算量的优化平衡,并足以用于研究所关心的具体问题。 相似文献
10.
利用Ｇａｌｅｒｋｉｎ有限元集中质量法,讨论了非饱和土壤水流动的数值计算模型,用以模拟均质土壤,地下水埋藏很深且在不同的初始和边界条件下的水分运动。通过将第一类边界条件齐次化,而对于第二类非线性边界条件采用变分的办法将其直接转化为已知边界通量的计算,从而容易且较好地处理了带通量约束边界条件。对于入渗或蒸发问题,用通常的离散方法容易产生振荡非物理解。我们采用有限元集中质量法,将质量矩阵的非对角元集中到主对角元,避免了这种非物理解的产生。利用该方案统一地模拟如下情况：稳定入渗和蒸发、蒸发一段时间后达到风干率,以及波动入渗蒸发时土壤水势的变化状况。模拟表明,在已知有关土壤参数的情况下,该模型能较好地模拟非饱和流水势及含水量分布。 相似文献
11.
The adjoint sensitivity related to explosive cyclogenesis in a conditionally unstable atmosphere is investigated in this study.The PSU/NCAR limited-area,nonhydrostatic primitive equation numerical model MM5 and its adjoint system are employed for numerical simulation and adjoint computation,respectively.To ensure the explosive development of a baroclinic wave,the forecast model is initialized with an idealized condition including an idealized two-dimensional baroclinic jet with a balanced three-dimensional moderateamplitude disturbance,derived from a potential vorticity inversion technique.Firstly,the validity period of the tangent linear model for this idealized baroclinic wave case is discussed,considering different initial moisture distributions and a dry condition.Secondly,the 48-h forecast surface pressure center and the vertical component of the relative vorticity of the cyclone are selected as the response functions for adjoint computation in a dry and moist environment,respectively.The preliminary results show that the validity of the tangent linear assumption for this idealized baroclinic wave case can extend to 48 h with intense moist convection,and the validity period can last even longer in the dry adjoint integration.Adjoint sensitivity analysis indicates that the rapid development of the idealized baroclinic wave is sensitive to the initial wind and temperature perturbations around the steering level in the upstream.Moreover,the moist adjoint sensitivity can capture a secondary high sensitivity center in the upper troposphere,which cannot be depicted in the dry adjoint run. 相似文献
12.
《Dynamics of Atmospheres and Oceans》1988,12(1):47-70
Two local implementations of no-slip boundary conditions are investigated for both the vorticity—streamfunction and momentum—pressure formulations of the time-dependent planar incompressible Navier-Stokes equations, as applied to barotropic ocean circulation modelling. The objective is to determine the extent to which the local accuracy and numerical consistency of these conditions affects the global solution. The effects of a non-local implementation of no-slip conditions for the vorticity—streamfunction equations are also studied. In all cases, boundary condition effects are measured by comparing time-averaged dynamics of turbulent solutions of numerical models based on the two formulations.In the model interior, the energy and enstrophy conserving Arakawa Jacobian is used for the vorticity—streamfunction equations while an extension of the energy and potential enstrophy conserving Arakawa and Lamb finite difference scheme is used for the momentum-pressure equations.Numerical experiments performed with a non-linear model similar to Bryan's barotropic ocean reveal no significant differences between the time-averaged solutions obtained with either of the two formulations, with each using either of the two local boundary conditions. A simple one-dimensional analogue of the vorticity—streamfunction equations is solved algebraically to explain the experimental results. A similar analogue suggests that an apparent inconsistency in the no-slip boundary conditions within the Cox stratified, primitive equation, ocean circulation model should not affect the accuracy or convergence of the global solution. 相似文献
13.
空气污染预报属于正问题，而从污染物浓度来求解扩散系数则属于反问题。正问题和反问题有着本质的不同，在解的定义和求解方法上也有很大的区别。从最优控制的角度定义了大气边界层中垂直扩散系数反演问题的解，用伴随模式方法得到目标函数的梯度并求解反问题。研究中发现，反演的结果与模式差分格式的选取有关，与测试源的设置也有直接的关系。在经过多次数值试验后，对于误差的来源进行了理论分析，发现了反演结果与差分格式及测试源之间的联系，得到了满意的反演结果，并为实验测定扩散系数提供了依据。 相似文献
14.
Large-Eddy Simulation of the Daytime Boundary Layer in an Idealized Valley Using the Weather Research and Forecasting Numerical Model 总被引：1，自引：1，他引：0
A three-dimensional numerical meteorological model is used to perform large-eddy simulations of the upslope flow circulation
over a periodic ridge-valley terrain. The subgrid-scale quantities are modelled using a prognostic turbulence kinetic energy
(TKE) scheme, with a grid that has a constant horizontal resolution of 50 m and is stretched along the vertical direction.
To account for the grid anisotropy, a modified subgrid length scale is used. To allow for the response of the surface fluxes
to the valley-flow circulation, the soil surface temperature is imposed and the surface heat and momentum fluxes are computed
based on Monin–Obukhov similarity theory. The model is designed with a symmetrical geometry using periodic boundary conditions
in both the x and y directions. Two cases are simulated to study the influence of along-valley geostrophic wind forcing with different intensities.
The presence of the orography introduces numerous complexities both in the mean properties of the flow and in the turbulent
features, even for the idealized symmetric geometry. Classical definitions for the height of the planetary boundary layer
(PBL) are revisited and redefined to capture the complex structure of the boundary layer. Analysis of first- and second-moment
statistics, along with TKE budget, highlights the different structure of the PBL at different regions of the domain. 相似文献
15.
Massimiliano Burlando Luigi Carassale Emilia Georgieva Corrado F. Ratto Giovanni Solari 《Boundary-Layer Meteorology》2007,125(3):417-439
In spite of recent progress in the prognostic numerical simulation of the atmospheric boundary layer, the explicit simulation
of turbulent flows in actual complex terrain is generally still very complicated and time consuming for many environmental
applications. In an attempt to develop simpler and more efficient application oriented techniques, although less refined,
we propose a multi-step procedure for simulating wind fields. Once obtained the necessary meteorological input, the mass-consistent
modelling technique is used to perform high-resolution mean wind flow simulations taking into account recent developments
in the atmospheric boundary-layer theory. Besides, a procedure based on a generalisation of the local logarithmic law-of-the-wall
over complex terrain is used to estimate the effective parameters characterising the simulated wind profiles. Turbulence intensities
and spectral properties are then calculated through the estimated effective parameters, in particular through the effective
friction velocity parameter. Finally, time series of the instantaneous velocity field are simulated by the Monte Carlo technique.
Two applications of the proposed approach are discussed briefly: the first one is related to a coastal area in southern Italy
(the Messina Straits), where the construction of the world’s longest central span bridge is being planned; the second one
corresponds to the flow in a mountainous area in northern Italy (the Albenga Airport). 相似文献
16.
J. Berg J. Mann A. Bechmann M. S. Courtney H. E. Jørgensen 《Boundary-Layer Meteorology》2011,141(2):219-243
17.
Large sudden wind-direction shifts and submeso variability under nocturnal conditions are examined using a micrometeorological network of stations in north-western Victoria, Australia. The network was located in an area with mostly homogeneous and flat terrain. We have investigated the main characteristics of the horizontal propagation of events causing the wind-direction shift and not addressed in previous studies. The submeso motions at the study site exhibit behaviour typical of flat terrain, such as the lower relative mesovelocity scale and smaller cross-wind variances than that for complex terrain. The distribution of wind-direction shifts shows that there is a small but persistent preference for counter-clockwise rotation, occurring for 55% of the time. Large wind-direction shifts tend to be associated with a sharp decrease in air temperature (74% of the time), which is associated with rising motion of cold air, followed by an increase in turbulent mixing. The horizontal propagation of events was analyzed using the cross-correlation function method. There is no preferred mean wind direction associated with the events nor is there any relationship between the mean wind and propagation directions. The latter indicates that the events are most likely not local flow perturbations advected by the mean flow but are rather features of generally unknown origin. This needs to be taken into account when developing parametrizations of the stable boundary layer in numerical models. 相似文献
18.
A method for solving the planetary boundary-layer equations 总被引：1，自引：0，他引：1
A method for solving the nonlinear three-dimensional steady-state equations which govern planetary boundary layer flow is described. The method is applicable to air motions over terrain with horizontally varying surface roughness, temperature and moisture. It can also be applied to a physical system consisting of the air and the sea (or earth) boundary layers taken together. Examples of calculations for selected cases of terrain variations are presented. Convergence of the method has been assessed by determining the degree to which the solutions satisfy the set of equations.The results of the calculations show that the method produces qualitatively realistic distributions of the different meteorological variables.Contribution No. 1251 from the University of Miami, Rosenstiel School of Marine and Atmospheric Sciences. 相似文献
19.
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. 相似文献
20.
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. 相似文献