共查询到19条相似文献,搜索用时 140 毫秒
1.
在求解Euler方程时,带TVB限制器的Runge-Kutta间断有限元(RKDG)方法是一种高精度、高并行效率的方法,而浸入边界方法是一种较新颖且对网格要求较低的方法,适用于处理复杂几何外形的边界.尝试了将上述2种方法结合起来求解Euler方程在笛卡尔网格上具有复杂几何外形的物体绕流问题,数个经典算例的数值结果验证了该方法的有效性. 相似文献
2.
3.
基于浅水波方程组建立一维溃坝流模型,并给出数值模拟结果.其中,空间离散采用HWENO(Hermit Weighted Essentially Non-Oscillatory) 格式,时间离散采用四步TVD(Total Variation Diminishing)Runge-Kutta模拟堤坝溃决时洪水演进过程。模拟结果表明:较采用WENO格式所得数值解更精确-同时,相比WENO格式的相应算法,该算法解决一维溃坝流问题能更有效地减弱振荡,对间断具有更高的分辨率 相似文献
4.
5.
用有限区域风速场准确求解流函数和速度势场的方法 总被引:4,自引:1,他引:3
流函数和速度势是气象业务和研究中常用于表述风速的一组变量。用有限区域风速场, 使用有限差分方法求解得到的流函数和速度势场重建初始风速场, 由于受区域边界的限制往往有明显的偏差。虽然有许多求解方法的研究, 但是, 至今仍尚未见到一种真正准确的求解计算方案。首先, 介绍用Arakawa A网格和D网格分布的有限区域风速场求解流函数和速度势场的一般有限差分计算方法, 探讨用它们的解重建风速场产生误差的原因。然后, 针对这些原因, 对给定的有限区域, 通过线性外推初始风速场, 扩展求解计算区域, 使用协调、一致的有限差分格式方案, 准确计算求解区域的边界有旋风速、散度风速和速度势的定解边界条件, 以及恰当选择流函数、速度势、涡度和散度等变量的分布网格, 设计了用上述两种网格分布的风速场准确求解流函数、速度势场的方案, 并对其正确性加以证明, 它们可以推广应用于其他Arakawa网格。用实际资料试验同样显示, 方案避免了重建风速场误差的出现, 与初始风速场相比, 全场风速最大偏差精度达到10-12m/s或以上, 在计算机精度造成的计算误差影响范围内。本文的研究很好解决了长期以来用有限区域风速场、 使用有限差分方法无法准确求解流函数和速度势场的问题。 相似文献
6.
球面上斜压原始方程组保持总有效能量守恒的差分格式 总被引:8,自引:9,他引:8
本文从标准层结近似下球面上的斜压原始方程组出发,针对两类常用的网格系统——C-网格和B-网格,分别设计出可以保持总动能、总有效位能和总有效表面位能之和守恒的差分格式.同时,讨论了定义在交错网格上的差分和平均算子的一些很有用的性质. 相似文献
7.
基于动态自适应网格的开源软件Gerris受到越来越多海洋和水文研究者的关注.概述了Gerris开发背景、研究现状和特点,详细阐述了Gerris数值方案,包括动态自适应网格、动态负载平衡技术原理、广义正交曲线坐标系、内嵌复杂固体边界和地形数据的处理方法,并探讨了Gerris在海洋数值模拟中的初步应用.结果表明,Gerris动态自适应网格在多尺度问题模拟中的优势独特,在海洋数值模拟应用中可通过自适应网格提高地理特征的精度,通过GTS(或KDT)格式的数据来处理地形和网格,达到同时兼顾精确性和易用性的目的,使得Gerris与其他海洋模式进行有机结合成为重要发展方向. 相似文献
8.
由于中尺度模式中采用不同的三维变量配置,必将产生不同的预报和模拟效果.怎样的三维变量分布,才能产生最佳的模拟效果呢?一直没有人从事这方面的研究.为此,在线性非静力滞弹性方程组的基础上,从频率、水平群速和垂直群速方面将中尺度模式中常见的四种三维网格的计算频散性与解析解的情况进行了对比,主要采用图示比较的方法.结果表明:总的来看三维网格C/CP计算频散性能最好,Z/LZ和Z/LY网格次之,C/L网格最差.因此设计非静力斜压模式时,应尽量采用C/CP网格.但C/CP网格对水平波长较短的波,误差相对要大些.如果模拟的波动水平尺度较小时,为减少误差,水平格距要减小.另外C/CP网格主要用在有限差分模式中;要考虑谱模式中的变量配置时,应把C/CP网格和Z/LZ网格或Z/LY网格结合起来使用. 相似文献
9.
城市边界层数值模式研究以及在香港地区复杂地形下的应用 总被引:10,自引:2,他引:8
由于城市表面复杂的下垫面的影响,使得城市边界层风温场结构较其他下垫面有很大不同.作者通过将城市中500 m网格内的各种地表类型,按照各自在网格中所占的百分比及各自的地表参数加权平均,得到此网格的平均的地表参数,以此准确反映下垫面的情况,建立分辨率为500 m的城市边界层能量平衡模式,将此城市边界层能量平衡模式嵌入动力学框架,并用中尺度模式MM5作为初始条件和边界条件,建立一个既考虑中尺度背景场又详细考虑城市下垫面复杂性与多样性的城市边界层模式系统.将模式系统运用于香港复杂地形下的边界层特点的模拟研究.通过与观测值的比较,模式能够较准确的模拟出海陆风、城市热岛等热力过程,及气流过山引起的绕流等动力过程,并且通过对边界层高度的模拟预测污染扩散的条件等.说明模式系统具有模拟在中尺度的背景场的控制下海陆风环流、过山堆积和绕流及城市热力影响的能力. 相似文献
10.
为了研究垂直变量配置对静力适应过程的影响,本文从描写静力适应过程的方程组出发,分别在将所有变量置于整层上的非跳点N网格;将垂直速度和温度放置在整层,水平速度、气压和密度等变量放置在半层的Charney Phillips跳点网格(CP网格);将水平速度、气压和温度放置在整层,将垂直速度和密度放置在半层的Lorenz跳点网格(L网格);将密度变量放置在整层的Charney Phillips跳点网格(CP_N网格);将密度放置在整层的Lorenz跳点网格(L_N网格)上进行离散,垂直格距分1 km、0.5 km、0.2 km和0.01 km,研究了在这5种网格上产生的频率和垂直群速的相对误差。结果表明:(1)L_N网格和CP网格是完全等效的两种网格。(2)不论垂直格距为多少,CP网格和L网格的误差都是最小,N网格次之,CP_N网格的误差最大。(3)随着垂直格距的减少,在这几种网格上产生的误差都在减小。对于CP网格、L网格和N网格,在水平长波和垂直短波处产生的误差较大。而CP_N网格对水平波长变化不敏感,垂直波长越短,误差越大。(4)当垂直格距为0.01 km时,这几种网格都对水平波长的变化不敏感了,仅对垂直波长敏感。(5)CP网格、L_N网格和L网格在描写静力适应过程和斜压地转适应过程都是误差最小的垂直变量配置方案,因此在非静力完全可压缩深层大气数值预报模式中应优先选择这3种方案。 相似文献
11.
自适应网格在大气海洋问题中的初步应用 总被引:15,自引:4,他引:15
自适应网格法是80年代兴起的通过求解椭圆型方程的边值问题来数值生成网格的一种新方法。它是在任意形状的区域上求偏微分方程的数值解的一种非常有效的工具。该方法抛弃了等距均匀的差分网格,代之以能够自动地适应所研究问题中解的特征的疏密程度不均的曲线网格。如在边界上计算网格与实际边界相重合,在区域内部可任意调节网格点的疏密程度等。本文扼要地介绍了自适应网格的原理及其构造方法。并将其应用于生成南海区域的计算网格以及数值预报台风路径的自适应网格。 相似文献
12.
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. 相似文献
13.
针对处理运动界面问题的流体体积函数(VOF)法,给出了一种高分辨率的运动界面捕捉方法.该方法采用五阶高精度和高分辨率的加权本质无振荡(WENO)有限差分格式离散VOF函数的空间导数;采用四阶Runge-Kutta方法离散时间导数;采用LocalLax-Friedrich通量作为数值流通量.用该方法对旋转流场和剪切流场中的运动界面追踪,结果表明该方法有较好的适用性和精确性. 相似文献
14.
1.IntroductionUnsaturatedsoilwaterflowisaflowwherewaterisnotfullofsoilhole,whichisanimportantformofflowinporousmedia.Predictionofanunsaturatedflowisprovidedwithsignificanceinmanybranchesofscienceandengineering.Theseincludeatmosphericscience,soilscien... 相似文献
15.
16.
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid. Two kinds of moments(i.e., point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell. The PV moment is updated by a conventional semi-Lagrangian method, while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation. Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme, a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment. To achieve exactly positive-definite preserving, two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme. The resulting scheme is inherently conservative, non-negative, and allows a Courant number larger than one.Moreover, the spatial reconstruction can be performed within a single cell, which is very efficient and economical for practical implementation. In addition, a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry, which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry. Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity, and it has the potential to transport the tracers accurately in a real atmospheric model. 相似文献
17.
D. Pavlidis G. J. Gorman J. L. M. A. Gomes C. C. Pain H. ApSimon 《Boundary-Layer Meteorology》2010,136(2):285-299
The computational fluid dynamics code Fluidity, with anisotropic mesh adaptivity, is used as a multi-scale obstacle-accommodating
meteorological model. A novel method for generating realistic inlet boundary conditions based on the view of turbulence as
a superposition of synthetic eddies is adopted. It is able to reproduce prescribed first-order and second-order one-point
statistics and turbulence length scales. The aim is to simulate an urban boundary layer. The model is validated against two
standard benchmark tests: a plane channel flow numerical simulation and a flow past a cube physical simulation. The performed
large-eddy simulations are in good agreement with both reference models giving confidence that the model can be used to successfully
simulate urban atmospheric flows. 相似文献
18.
Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial guess or the subsequent updated solution.In the second type,the NBE is rearranged into a quadratic form of the absolute vorticity with the positive root of this quadratic form used in the form of a Poisson equation to solve NBE iteratively.The two methods are rederived by expanding the solution asymptotically upon a small Rossby number,and a criterion for optimally truncating the asymptotic expansion is proposed to obtain the super-asymptotic approximation of the solution.For each rederived method,two iterative procedures are designed using the integral-form Poisson solver versus the over-relaxation scheme to solve the boundary value problem in each iteration.Upon testing with analytically formulated wavering jet flows on the synoptic,sub-synoptic and meso-αscales,the iterative procedure designed for the first method with the Poisson solver,named M1a,is found to be the most accurate and efficient.For the synoptic wavering jet flow in which the NBE is entirely elliptic,M1a is extremely accurate.For the sub-synoptic wavering jet flow in which the NBE is mostly elliptic,M1a is sufficiently accurate.For the meso-αwavering jet flow in which the NBE is partially hyperbolic so its boundary value problem becomes seriously ill-posed,M1a can effectively reduce the solution error for the cyclonically curved part of the wavering jet flow,but not for the anti-cyclonically curved part. 相似文献
19.
Turbulent Transfer Between Street Canyons and the Overlying Atmospheric Boundary Layer 总被引:1,自引:1,他引:0
Pietro Salizzoni Massimo Marro Lionel Soulhac Nathalie Grosjean Richard J. Perkins 《Boundary-Layer Meteorology》2011,141(3):393-414
The turbulent exchange of momentum between a two-dimensional cavity and the overlying boundary layer has been studied experimentally,
using hot-wire anemometry and particle image velocimetry (PIV). Conditions within the boundary layer were varied by changing
the width of the canyons upstream of the test canyon, whilst maintaining the square geometry of the test canyon. The results
show that turbulent transfer is due to the coupling between the instabilities generated in the shear layer above the canyons
and the turbulent structures in the oncoming boundary layer. As a result, there is no single, unique velocity scale that correctly
characterizes all the processes involved in the turbulent exchange of momentum across the boundary layer. Similarly, there
is no single velocity scale that can characterize the different properties of the turbulent flow within the canyon, which
depends strongly on the way in which turbulence from the outer flow is entrained into the cavity and carried round by the
mean flow. The results from this study will be useful in developing simple parametrizations for momentum exchange in the urban
canopy, in situations where the street geometry consists principally of relatively long, uniform streets arranged in grid-like
patterns; they are unlikely to be applicable to sparse geometries composed of isolated three-dimensional obstacles. 相似文献