共查询到20条相似文献,搜索用时 46 毫秒
1.
数值模拟圆柱绕流问题是计算流体力学和计算气动声学中的一个经典问题,要准确模拟这一类问题,远场无反射边界条件是一个关键环节。本文采用适用于主控方程为Navier-Stokes方程的完全耦合层无反射边界条件对低雷诺数条件下圆柱绕流进行直接数值模拟,研究了计算域大小的选取对数值结果的影响。数值模拟研究表明:(1)由完全耦合层边界条件产生的相对反射误差很小;(2)如果仅是为了获得圆柱附近的流场和声场信息,采用完全耦合层边界条件可以选用较小的计算域;(3)适当增加计算域的上下边界宽度,可以很大程度地减小右边界的宽度,从而在保证精度的情况下减小计算量。 相似文献
2.
研究了数值求解双曲方程外问题的人工边界方法.在圆形的人工边界上,得到了三类等价的完全无反射的人工边界条件.给出的数值例子验证了方法的有效性. 相似文献
3.
基于光强传播方程的相位恢复 总被引:3,自引:3,他引:3
对于圆形孔径的光学系统,利用光强传播方程进行相位恢复时很难得到相位的边界径向斜率值,另外,要获得精确的圆域边界采样值也并非易事。为了克服上述困难,提出了一种相位恢复的改进方法,即改变了方程的表示形式、计算区域和边界条件,并用多重网格法进行求解获得重构相位,最后再将其修正。为了对重构相位进行修正,还得到了重构相位与原始相位之间泽尼克系数的传递矩阵。对均匀照明的情形进行了仿真,发现该方法不仅可以避免复杂的边界条件,减少运算时间,而且还能够较好地恢复出原始波前的泽尼克系数,即便在加噪的情况下,修正相位与原始相位的均方根误差也在可以接受的范围内。 相似文献
4.
在傍轴近似条件下, 可以利用光强传播方程(ITE)对畸变波前进行相位恢复。对于衍射受限的光学系统, 很难获得相位的边界径向斜率值作为边界条件, 此外, 要获得精确的圆域边界采样值也并非易事。为了克服上述困难, 进一步研究了相位恢复改进方法, 即改变了方程的表示形式、计算区域和边界条件, 并用多重网格法进行求解获得重构相位, 最后再将其修正。为了验证算法的精确性, 搭建了实验系统对算法进行实验验证, 将由CCD探测的光强分布进行计算得到的相位分布与相位恢复算法(PR算法)的结果进行比较, 证明在光强分布非均匀的情况下, 该方法恢复的相位均方根误差能够达到017 λ, 可以适用于波前传感精度要求不是很高的相位恢复。 相似文献
5.
6.
7.
提出了一种便于分析光纤布拉格光栅(FBG)、基于强制边界条件的简化TD-FD-BPM模型.采用交替隐式法(ADIM)和综合道格拉斯(GD)法,离散化时域束传输方程,获得比普通差分格式低二阶的误差。提出并采用强制边界条件,简化了光波导分析中边界处理问题.重点讨论和确定了在此条件下激励源位置、脉冲形状和径向分布函数、连续波激励条件等.在修正的正弦脉冲激励情况下,通过离散傅利叶变换,直接得到较为理想的短FBG反射谱.此法可减少数值计算时间,直观反映光脉冲传输. 相似文献
8.
9.
对于非线性Euler方程,提出一类基于完美匹配层(PML)技术的吸收边界条件。首先对线性化的Euler方程设计出PML公式,然后将线性化Euler方程中的通量函数替换成相对应的非线性通量函数,得到非线性的PML方程。考虑到PML方程中包含有一个刚性的源项,文中采用一种隐显Runge-Kutta方法来求解空间半离散后得到的ODE系统。数值实验表明设计的非线性PML吸收边界条件优于传统的特征边界条件。 相似文献
10.
11.
12.
本文采用声类比法探讨了射流中小尺度湍流远场噪声特性.此声类比理论基于线化欧拉方程及声源互易原理,建立了关联声源和远场声压的伴随格林函数,适于声波折射现象及小尺度湍流噪声计算.声场的边界采用完全匹配层边界条件,有效地避免了声波产生数值反射.为减小频散及耗散误差,时空导数均采用高阶预测-校正格式离散控制方程.计算结果与前人的实验相一致. 相似文献
13.
14.
介绍了有限区域水平风场分解的调和-余弦计算方法,该方法把函数分成两部分之和.第一部分是Laplace方程在给定边界条件下的解,由于Laplace方程的解是调和函数,这个部分可称为调和部分,又因为其与区域内部值无关,也称外部部分.第二部分是原始函数与调和部分之差,这个函数是齐次边条件下Poisson方程的解,只与区域内部的涡度或散度有关,故称为内部部分,可以展开成双傅氏的余旋函数系列.调和-余弦计算方法的求导都是用谱系数进行,计算精度比常用的差分方法高两阶以上.而且,由于外部部分对应的边界条件物理意义清楚,边界光滑,成功克服了有限区域流函数和速度势迭代求解出现的计算不稳定、原始风场无法还原、边界上的系统缺失等问题,可以准确分解和重建有限区域的风场.利用NCEP/NCAR 1°×1°的实时分析资料和日本气象厅区域谱模式(RSM)20km分辨率的再分析资料,利用调和-余弦算法得到的无辐散风分量和无旋风分量,对2006年的8号超强台风“桑美"(SAOMEI)进行风场结构的比较分析.结果发现,低层无辐散风比原始风场与台风中心的对应关系更好;同时,无旋风分量能更好地显示原始风场上并不明显的低层辐合高层辐散的特征,大尺度无辐散风分量可以更清晰地显示出台风的水汽输送通道.从与台风中心的对应关系看,台风在海上发展阶段,SAOMEI台风的旋转中心与辐合中心并不是时时重合,这个特点只能通过风场分解才能得到.此外,SAOMEI登陆以后,南部洋面上发展起来的对流活动从水汽和能量补充方面都不利于SAOMEI的维持.可见,分解后的无辐散风场和无旋风场能更清楚地体现出SAOMEI的风场结构,在台风结构分析中有重要的推广应用价值.
关键词:
台风
水平风场分解
调和-余弦算法 相似文献
15.
16.
17.
18.
运用Landau-Lifshitz方程,边界处动态磁化强度由有效偶极边界条件约束,研究了无限长金属磁条中自旋波传播的特征,得到了该体系抽运微波磁场的阈值曲线以及色散曲线的解析式,揭示出自旋波激发谱与磁条宽度的具体关系.结果表明,平面约束下磁条的阈值曲线中出现了明显的扭结现象.并且随着磁条宽度的增加,其扭结数近似地呈指数形式增加,而最低扭结处两阈值的跳跃差值几乎呈反比形式变化.
关键词:
Landau-Lifshitz方程
有效偶极边界条件
阈值曲线
色散曲线 相似文献
19.
S_1流面跨声流场流函数矩阵解 总被引:1,自引:0,他引:1
跨声速叶栅流的计算,可采用时间相关法求解Euler方程,或用松弛方法求解势函数方程和流函数方程。一般说来,时间相关法耗费机时较多,势函数方法仅对无旋流适用。流函数方法适用于二元有旋流的计算,并且边界条件也较为简单,可方便地进行S_1和S_2两类流面迭代得到三元解。流函数方法的跨声计算最大的困难是密度双值问题 相似文献
20.
耦合随机湍流速度生成模型与线化欧拉方程技术,形成了一套具备模拟噪声在非均匀流场中传播能力的气动噪声混合预测方法。该混合方法的随机湍流速度生成模型采用了快速随机粒子网格法,为声传播模拟提供了可靠的源项。而噪声的传播计算选用线化欧拉方程,其空间离散采用9点5阶的色散保持关系格式,时间推进选用了高精度大时间步长的6级4阶龙格库塔格式,远场边界应用了无分裂形式的理想匹配层边界条件。首先,选用高斯脉冲传播算例对线化欧拉方程的时空离散格式、远场无反射边界条件进行了验证分析。然后,计算分析各向同性湍流的空间相关性验证湍流速度生成模型的可靠性。最后,基于已搭建的气动噪声混合预测方法进行了30P30N三段翼缝翼噪声的计算分析。计算分析可知:监测点处功率谱密度曲线、噪声指向性等计算结果与参考文献结果取得了较好的一致性。数值计算结果表明所建立的气动噪声混合预测方法能有效预测二维复杂构型的气动噪声问题。 相似文献