首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 219 毫秒
1.
赵滨  张博  钟青 《气象科学》2010,30(1):8-14
蛙跳(leapfrog)时间差分格式采用Asselin-Robert时间滤波方案去除计算解能够降低原始方程组的时间差分格式的计算精度,采用二阶Adams-Bashforth格式构造的欧拉前差方案可弥补蛙跳格式的不足,即:在不存在计算解的条件下去除滤波的影响,更大程度的保持方程组的计算准确性。本文基于NCAR CAM3.0(Community Atmosphere Model 3.0)完善的软件平台,将原模式的三时间层蛙跳时间差分方案修改为两时间层二阶Adams-Bashforth时间差分格式,对与重力波有关项使用中央差隐式处理,以此构建半隐式大气环流谱模式。利用动力检验的方法探讨模式对垂直分辨率的敏感性,从而寻找模式在较小计算代价下提高计算效果的可能。通过斜压波实验发现,提高垂直分辨率使模式具有更强的斜压波模拟能力,其模拟效果甚至已经与提高水平分辨率的效果相当,可以作为一种弥补模式运算效率不足的可行方案。  相似文献   

2.
陈嘉滨  江野 《大气科学》1993,17(3):274-282
本文研究了具有参考大气的谱模式半隐式时间差分格式的计算稳定性问题.通过试验与分析,找到了影响这种半隐式格式计算稳定性的关键因子.另外本文还进一步试验了第一作者过去提出的修改的半隐式时间差分格式,这种格式采用了修改半隐式项而不改变重力波性状的做法.取得了稳定的计算结果.  相似文献   

3.
陈雄山  邢如楠 《气象学报》1981,39(4):495-502
用低分辨率的初始方程二层模式在给定下垫面温度下模拟了一月及七月全球海平面气压场,其中加热项主要是感热垂直输送。计算采用守恒的Lilly空间差分格式及Matsuno的时间差分格式。海平面气压的初值给成常数值,高空的风速及气温的初值给成纬向平均的气侯值。在数值模拟中,下垫面温度用一月及七月的气候实际值。计算的冬半球海平面气压与观测值相当一致,而夏半球的结果稍差些。本模式所求得的全球海平面气压场能与实况相比,结果并不差于高分辨率的复杂模式。  相似文献   

4.
斜压原始方程半隐式全能量守恒格式的构造问题长期没有解决。本研究在成功地构造实现其全能量完全守恒的半隐式方案基础上,进行了此守恒方案与欧洲中期天气预报中心(ECMWF)的σ-坐标原始方程全球谱模式半隐式方案间的实际资料对比实验。实验表明,850hPa平均预报高度场RMS误差在积分一周以后得到明显改进,到第30天其预报误差降低达到了50%,进一步的对比实验表明,对流层中部和下部的月预报平均高度场RMS误差也显降低,而且一些明显的系统性误差也得到大幅度改进。更加详细的分析显示,这些收益的很大一部分是从超长波成分的改进中得到的。这说明,通过构造守恒性时间差分方案消除了响应的计算性系统误差源汇,进而能够使模式气候漂移得到显改进,而这种误差源汇存在于传统的,现仍被普遍采用的斜压原始方程天气气候模式中。  相似文献   

5.
显式自由表面模型时间分裂格式的一个积分方案   总被引:5,自引:0,他引:5  
讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中,由于采用蛙跳格式,斜压动量方程具有三层的时间离散结构,但是,正压动量方程在相应的时间层上只有二层的离散结构。因此,正压方程与斜压方程的时间离散结构是不相容的。分析结果表明,这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。  相似文献   

6.
钟青 《气象学报》1997,55(6):641-661
文中构造并证明了一般二次和三次物理守恒律时间差分保真格式两个构造定理,以往一些主要时间离散守恒格式构造方案可作为两个定理特例给出。它们不仅可为解决更加广泛类别的时间离散保真格式构造基本问题提供适用数学基础,而且也为结合已有瞬时空间离散守恒格式,解决更加广泛类别的时-空离散意义下保真格式构造基本问题提供适用的数学基础。此外,文中两个定理还可解决两大类问题的线性和非线性计算不稳定性问题。斜压原始方程传统半隐式全球谱-垂直有限差分模式目前是世界上许多国家的业务预报和大气环流模式。本工作利用文中新构定理,构造并且实现了斜压原始方程全球谱-垂直有限差分模式半隐式高阶全能量守恒方案。以往该项基本问题无论在理论还是实践上长期以来一直都未能得到解决。该项全能量守恒半隐式全球谱模式方案适用于实测资料的长时间数值预报积分。使用FGGE夏季资料进行的13个个例30d数值积分实验表明:新型全能量半隐式保真方案可以有效地改进传统预报方案中关于能量质量守恒性质的系统性偏差。值得注意的是,实验统计分析还显示:在本文实验条件下,传统方案中由于时间离散过程中原物理守恒律性质破坏导致的系统误差(简称Z类误差),对于实验总体均方根系统误差的贡献  相似文献   

7.
关于半拉格朗日半隐式大气模式的时步问题   总被引:1,自引:0,他引:1  
讨论半拉格朗日半隐式大气模式中时步的限制,展示时步过大时用一般的首尾两点平均方案计算非线性源汇项的严重误差,提出精确的源汇项计算格式并作了特例计算和比较,讨论了被模拟的大气过程和波动的特征对时步的要求,指出柯朗数NC=C·△t/△X是大气模式时空步长匹配的重要参数.指出云降水和大气化学过程特征时间对相应的正定变量的计算时步的限制.  相似文献   

8.
台湾海峡海陆风数值模式与数值模拟试验   总被引:5,自引:1,他引:5       下载免费PDF全文
研制了一个包括水平及垂直扩散、牛顿冷却的二维46层非弹性运动方程组的台湾海峡海陆风数值模式,并用此模式来模拟及研究台湾海峡两岸海陆风的生成与变化特征.模式中考虑了太阳辐射、长波辐射及其日变化、地表向大气的感热与潜热输送以及向土壤层的热传导等.数值计算中采用了分解算法及隐式时间差分方案.用此模式得出的模拟结果与闽东南及台湾海陆风的观测事实比较吻合,表明了此模式能够较好地描述海峡两岸的海陆风变化规律.  相似文献   

9.
基于浅水波方程组建立一维溃坝流模型,并给出数值模拟结果.其中,空间离散采用HWENO(Hermit Weighted Essentially Non-Oscillatory) 格式,时间离散采用四步TVD(Total Variation Diminishing)Runge-Kutta模拟堤坝溃决时洪水演进过程。模拟结果表明:较采用WENO格式所得数值解更精确-同时,相比WENO格式的相应算法,该算法解决一维溃坝流问题能更有效地减弱振荡,对间断具有更高的分辨率  相似文献   

10.
设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法(GCR)求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。   相似文献   

11.
The Asselin-Robert time filter used in the leapfrog scheme can degrade the accuracy of calculations. The second-order Adams-Bashforth method with the same accuracy as the leapfrog scheme is not subject to time splitting instability. A new semi-implicit atmospheric general circulation spectral model is developed on the basis of NCAR (National Center for Atmospheric Research) CAM3.0 (Community Atmosphere Model3.0). In this new model, the second-order Adams-Bashforth method is used as an alternative to the leapfrog scheme, and a Crank-Nicholson scheme is incorporated for the treatment of fast gravity modes. In this paper, the new model is tested by the Held-Suarez test and an idealized baroclinic wave test. Results of the Held-Suarez test show that the second-order Adams-Bashforth model has similar climate states to those of many other global models and it converges with resolutions. Based on the idealized baroclinic wave test, the capability of di?erent time di?erencing methods for keeping the initial steady-state are compared.This convinces a better ability of the second-order Adams-Bashforth method in maintaining the stability of the initial state. Furthermore, after the baroclinic wave is triggered through overlaying the steady-state initial conditions with the zonal perturbation, the second-order Adams-Bashforth method has an excellent property of convergence, and can represent the process of the baroclinic wave development much better than the original scheme in CAM3.0. A long-term integration of the new model during the period of 1980-1999 is also carried out and compared with that of CAM3.0. It is found that due to the reduction of simulation errors of prognostic variables, the second-order Adams-Bashforth method also has a better simulation ability for the diagnostic variables, such as precipitation.  相似文献   

12.
The Asselin-Robert time filter used in the leapfrog scheme does degrade the accuracy of calculations. As an attractive alternative to leapfrog time differencing, the second-order Adams-Bashforth method is not subject to time splitting instability and keeps excellent calculation accuracy. A second-order Adams-Bashforth model has been developed, which represents better stability, excellent convergence and improved simulation of prognostic variables. Based on these results, the higher-order Adams-Bashforth methods are developed on the basis of NCAR (National Center for Atmospheric Research) CAM 3.1 (Community Atmosphere Model 3.1) and the characteristics of dynamical cores are analyzed in this paper. By using Lorenz nonlinear convective equations, the filtered leapfrog scheme shows an excellent pattern for eliminating 2Δt wave solutions after 20 steps but represents less computational solution accuracy. The fourth-order Adams-Bashforth method is closely converged to the exact solution and provides a reference against which other methods may be compared. Thus, the Adams-Bashforth methods produce more accurate and convergent solution with differencing order increasing. The Held-Suarez idealized test is carried out to demonstrate that all methods have similar climate states to the results of many other global models for long-term integration. Besides, higher-order methods perform better in mass conservation and exhibit improvement in simulating tropospheric westerly jets, which is likely equivalent to the advantages of increasing horizontal resolutions. Based on the idealized baroclinic wave test, a better capability of the higher-order method in maintaining simulation stability is convinced. Furthermore, after the baroclinic wave is triggered through overlaying the steady-state initial conditions with the zonal perturbation, the higher-order method has a better ability in the simulation of baroclinic wave perturbation.  相似文献   

13.
为了使用以制作数值预报的微分方程中包含更多的信息量、以期提高预报准确度, 该文提出了一种基于记忆动力学的时间积分格式。以天气预报为实例, 计算表明, 平流方程的回溯时间积分格式所得的预报准确度大大高于传统的蛙跃差分格式。此方法在海洋、水文、环境、航空等应用平流方程计算的多种科学中亦是有效的。  相似文献   

14.
AnEconomicalConsistentDisipationOperatorandItsApplicationstotheImprovementofAGCM①WangBin(王斌)andJiZhongzhen(季仲贞)LASG,Institute...  相似文献   

15.
This paper introduces a new consistent dissipation operator. It is based on the explicit square conservation scheme and the theory of consistent dissipation. The operator makes full use of the advantages of the Leap-frog scheme, i.e., its second order time precision and its explicit solution manner. Meanwhile, it overcomes the fatal disad-vantage, the absolute instability in computations, of the scheme. When it is applied to the explicit square conservation scheme, the time precision of the scheme reaches to third order. Especially, the computational stability of this scheme is as good as the third order explicit Runge-Kutta scheme. The CPU time required in computations by the scheme is less than that required by the explicit square conservation scheme with the consistent dissipation operator constructed from the Runge-Kutta method. Therefore, the new operator is an economical one. The application of the operator to the improvement of the dynamical model of the L2IAP AGCM shows its time-saving property and its good effects  相似文献   

16.
In this paper, a special three-step difference scheme is applied to the solution of nonlinear time-evolution equations, whose coefficients are determined according to accuracy constraints, necessary conditions of square conservation, and historical observation information under the linear supposition. As in the linear case, the schemes also have obvious superiority in overall performance in the nonlinear case compared with traditional finite difference schemes, e.g., the leapfrog(LF) scheme and the complete square conservation difference(CSCD) scheme that do not use historical observations in determining their coefficients, and the retrospective time integration(RTI) scheme that does not consider compatibility and square conservation. Ideal numerical experiments using the one-dimensional nonlinear advection equation with an exact solution show that this three-step scheme minimizes its root mean square error(RMSE) during the first 2500 integration steps when no shock waves occur in the exact solution, while the RTI scheme outperforms the LF scheme and CSCD scheme only in the first 1000 steps and then becomes the worst in terms of RMSE up to the 2500th step. It is concluded that reasonable consideration of accuracy, square conservation, and historical observations is also critical for good performance of a finite difference scheme for solving nonlinear equations.  相似文献   

17.
天气雷达测定区域降水量方法的改进与比较   总被引:9,自引:1,他引:9  
讨论变分校准法用于雷达-雨量计系统联合探测降水。由雷达反射率因子Z和地面降水强度I实时地获得最优Z-I关系,在求解欧拉方程时采用多重网格法,不仅可提高计算结果的精度,还可大大提高计算速度。  相似文献   

18.
回溯时间积分格式的有效性研究   总被引:3,自引:0,他引:3       下载免费PDF全文
游性恬  朱禾  曹鸿兴 《大气科学》2002,26(2):249-254
采月正压准地转模式作为动力核,用Rossby-Haurwitz波函数作为理想场,同时使用回溯时间积分格式和普通中央差格式进行数值试验.结果表明,在短时积分的情况下,回溯格式能减小数值解误差1~2个量级,而且它对时间步长的增大不敏感.因此使用回溯格式可以减少积分计算量,延长积分时效.  相似文献   

19.
JFNK方法在求解全隐式一维非线性平流方程中的应用   总被引:1,自引:0,他引:1       下载免费PDF全文
JFNK(Jacobian-free Newton-Krylov)方法是由Newton迭代方法和Krylov子空间迭代方法构成的嵌套迭代方法。作者以全隐式差分的一维非线性平流方程(亦称无粘Burgers方程)探讨采用全隐式格式计算的必要性和JFNK方法的有效性。模拟结果表明, 隐式结果比显式和半隐式结果在稳定度和精度方面较大的优越性, 特别是模拟气流强的系统以及要素空间分布具有较大梯度的系统。  相似文献   

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

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