共查询到20条相似文献，搜索用时 203 毫秒
1.
针对一类大气、海洋问题，从算子方程出发，在显式完全平方守恒差分格式的基础上提出两种变时间步长的差分算法，使数值计算更加有效。在人工调节时间步长法中，通过人为控制方程中算子能量变化的办法调整时间步长；在自动调节时间步长法中，通过数值试验自动比较出所构造差分格式稳定性的好坏，将它与泰勒展开法及分解算法进行混合应用，验明了自动调节时间步长法具有综合判断差分格式好坏的功能，并得出一些有意义的结论。 相似文献
2.
3.
本文用差分模拟绝对涡度守恒定律,构成正压涡度方程的两类差分格式:积分守恒型格式和相当于准拉格朗日型格式。前一类型格式的显式部分是条件稳定的,隐式部分是无条件稳定的。后者至为无条件稳定的显式格式,可以依据实际需要和可能,来放大无条件稳定格式的时间步长。 本文还根据格式的计算要求,讨论了“风”场的计算和平滑问题,以及差分方程的边、初条件(计算边、初条件)的给法。并通过理想场计算和实例预报,检验了几种显式格式的计算效果。 相似文献
4.
大气海洋环境数值模拟中的若干计算问题 总被引：1，自引：0，他引：1
为了适应发展气候数值模拟和环境数值模拟的需要,综合介绍了与此密切有关的若干计算问题。首先,指出这类问题易于出现计算紊乱、非物理解和非线性计算不稳定,同时,简要介绍造成这些问题出现的三种机理：虚假频散、能量关系破坏和能谱的非线性转移。然后,把这类问题归结为一种“发展方程”,给出了与此有关的几个定理,阐述了计算稳定性、能量守恒性和算子非负性之间的密切关系。接着,分别介绍了多种完全平方守恒的差分格式的构造方法,其中包括隐式的和显式的完全平方守恒的差分格式,高时间精度和高阶紧致的完全平方守恒的差分格式,同时也给出若干具体算例说明这些格式的计算效果。最后,对全文作了小结,并指出今后应进一步研究的若干重要问题。 相似文献
5.
天气数值预报工作中所用的初始方程显式差分格式,时间步长取得很短,实际应用很不方便。比如,空间步长取300公里,时间步长只能用10分钟。本文提出一种新型显式差分格式,其计算稳定性不受时间步长(τ)和空间步长(d)大小的限制,由|λ_e|τ/d(e=1,2,3)的整数部分m_e=[|λ_e|τ/d],控制执行“浮动”计算,保持计算稳定。这种差分格式曾用104电子计算机作检验计算,空间步长用25公里,时间步长取20分钟,计算过程稳定,结果精度较高,可与文献[1]中的结果相比较。这表明,初始方程的差分计算,适当延长时间步长是有可能的。 相似文献
6.
一类计算稳定性好的显式平流差分格式 总被引：1，自引：0，他引：1
通常的显式平流差分格式，如迎风格式，Lax-Wendroff格式等，均是有条件稳定的,其稳定条件与差分网格的时、空步长有关。本文对线性和拟线性平流方程分别构造了一种计算稳定性好的显式差分格式。对前一格式，本文严格证明了它的无条件稳定性及收敛性，并具有一阶精度；对后一格式，由于非线性方程的限制，本文用数值试验研究了它的计算稳定性。 相似文献
7.
本文首先讨论了如何依据空间步长来选取平滑系数。其次考虑在涡度方程中加上粘性项后对微分方程解的影响,讨论如何选取适当的粘性系数。第三讨论涡度方程加粘性项后对差分格式计算稳定性的影响,证明了在常用时间中心差分格式中应把粘性项取在第(n-1)时间层,这样计算才是稳定的。反之,如把粘性项取在第n时间层,则差分格式计算不稳定。并且时间向前差的格式在加上粘性项后也可以稳定,但要求的粘性系数过大。最后讨论了用涡度方程作数值预告时加上平滑过程与加上粘性项之间的关系,从计算上看,二者并不等价。但可以把平滑过程与扩散方程类比,以确定适当的平滑系数。 相似文献
8.
吴诚鸥 《南京气象学院学报》1981,(2)
时间序列分析中的AR、MA、ARMA、ARIMA模型都是一些线性随机差分方程。本文将用推移算子理论,求得线性差分方程的通解(定理4和定理5)。并应用于时间序列分析。 相似文献
9.
发展方程差分格式的构造和应用 总被引：5，自引：5，他引：5
本文把许多拟线性方程归结为算子形式的“发展方程”:F/t+AF=0.证明了带有非负算子A的“强隐式(即1/2≤θ≤1)格式”是绝对隐定的;而带有反对称算子A的“弱隐式和显式(即0≤θ≤1/2)格式”是绝对不稳定的。文中又以一维非线性平流方程为例,具体地构造了具有非负算子A的三种差分格式和相应的带人工耗散项的稳定格式,给出一个计算实例,说明后一类格式对消除虚假的寄生波也是很有效的,适于计算有间断的问题。 相似文献
10.
线性与非线性发展方程计算稳定性的比较分析 总被引：21，自引：0，他引：21
针对性与非线性发展方程的几种差分格式，以一维线性和非线性平流方程为例，对线性与非线笥发展方程差分格式的计算稳定性进行了比较分析，揭示了差分格式结构和初值形式与计算稳定性的关系。理论分析和数值试验证明，线性与非线性发展方程差分格式的计算稳定性在本质上是完全不同的。 相似文献
11.
12.
文中指出了分解算法中的两个问题:(1)为了使算法获得最佳效果,除按快、慢过程分解方程外,进一步的分解应按什么原则?(2)时间离散是否可以任意,如果不然,应采取什么形式?针对这两个问题得到的研究结果是:和一般设计简化模式一样,分解后的方程应保持原方程的重要物理性质,特别是其整体性质(如能量守恒和质量守恒等),而且,其分解后方程的解还应和原方程的解一致,至少对其中重要的解应一致.至于时间离散,在差分的情况则应采用二时间层格式;如采用三时间层格式,则会导致差分方程(或差分微分方程)和偏微分方程的不相容. 相似文献
13.
运用判定非线性发展方程差分格式计算稳定性的Hirt启发性分析方法,对一类非线性Schrdinger方程差分格式的计算稳定性进行分析,得到了保证差分格式计算稳定的必要条件.数值试验结果进一步表明,得到的稳定性判据不仅是保证差分格式计算稳定的必要条件,而且在实际中也是非常有效的. 相似文献
14.
分析了新一代非静力中尺度数值模式中常用的三阶龙格-库塔时间分裂显式算法(RK3)的稳定性和误差性质,特别是分析了空间中央差分和迎风偏斜两种不同情况下该算法不同的稳定性和误差性质。运用数学软件先进的符号计算功能,分析了该算法涉及的复杂高阶、高次幂振幅矩阵的特征值性质;并通过一维线性声波-平流方程组的数值模拟实验,检验了时间分裂算法的模拟效果。对振幅矩阵特征值模的表达式进行高阶的级数展开,得到了该算法的分裂误差项的公式;而且,由于特征值模的公式保留了较高阶项,可以同时分析迎风偏斜和中央差两种空间差分格式的分裂误差性质。根据分裂误差项公式,定量地比较了三阶和二阶龙格-库塔格式(RK2)的分裂误差大小以及误差与小时间步数的关系,发现迎风格式RK3的分裂误差明显小于RK2的误差,并具有更好的稳定性质。空间中央差格式的分裂误差项具有更高阶数,比迎风格式具有更小的时间分裂误差。对于各种不同波长的特征值分析和采用中央差格式的数值模拟,也进一步证实空间差分采用中央差时,RK3时间分裂显式算法在不同方向传播的声波振幅几乎没有差别。另外,误差公式以及数值试验结果说明RK3的分裂误差也略小于Adams-Bashforth-Moulton分裂显式法的分裂误差。 相似文献
15.
强迫耗散非线性发展方程准完全平方守恒格式的构造 总被引：5，自引：0，他引：5
从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散性发展方程的显式准完全平方守恒差分格式，理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的．值得推广应用。 相似文献
16.
变步长显式完全平方守恒差分格式 总被引：8，自引：0，他引：8
综合隐式完全平方守恒差分格式和显式瞬时平方守恒差分格式的优点，针对一类非线性发展方程构造了一种通过自动调节时间步长来保持平方守恒性的显式差分格式。它基于加小耗散的思想，但又与小耗散法有所不同。本文取的耗散项不是一般的人工耗散，而是取能够弥补由于时间离散所产生的截断误差的所谓（时间）协调耗散。因此，该格式具有较高的时间精度。在数值试验中，该类格式可取得满意的效果。 相似文献
17.
在对微分方程初值问题的数值计算中 ,当采用前向或后向差分格式时有两个时间层次 ,而当用中央差分格式时则有 3个时间层次 .由于把预报问题提为初值问题 ,因此只有一个初始场的信息被利用 .本文从扩展信息源的角度出发 ,对微分方程的时间积分提出了一种新格式——回溯时间积分格式 ,并对在大气、海洋、水文等常用的平流方程进行了实例计算 .制约动力系统的微分方程可以写为 : u/ t=F ( u,t) ( 1 )其中 ,u为状态变量 ,F为除局地变化项外的所有项 ,称为源函数 ,设时间间隔为△ t,式 ( 1 )的时间中央差格式为 :u( t △ t) =u( t-△ t) 2△ t… 相似文献
18.
再论发展方程差分格式的构造和应用 总被引：40，自引：8，他引：32
【摘 要】本文把一大类大气、海洋方程归结为一种发展方程,具体构造了若干定时间步长的显式完全平方守恒差分格式。并证明在一定条件下,这类格式也具有能量守恒、“广义能量”守恒和“平均尺度”守恒的特性,它表明这类格式具有较好的计算稳定性和省时性。文中还探讨了显式平方守恒格式与隐式平方守恒格式之间的密切联系。最后给出了令人满意的用四波的R-H波作数值检验的结果。 相似文献
19.
20.
本文用“单向”(one—way)方法进行正压原始方程套网格试验。粗、细网格均采用空间分解准拉格朗日格式和由廖洞贤等提出的与此差分格式一致的内边界条件(即方案二)。试验结果表明,用方案二作出的24小时细网格预报,与按半动量格式和粗网格直接给细网格提供内边界条件(即方案一)所得的结果相比,去除了在流出处出现的不稳定,系统的流入、流出也较好。这两种方案再加上Davies“松弛区域”的方法(分别为方案四,方案三),其结果与上述类似,且方案四比方案二还要好些。 相似文献