共查询到20条相似文献,搜索用时 31 毫秒
1.
The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry. 相似文献
2.
针对线性椭圆型问题的,发展了一种无重叠区域分裂方法,首先引入了Lagrange乘子,以使子域交界处解的连续性约束条件获得弱满足。接着重点讨论了Lagrange乘子的近似空间,即粘接元(mortar elements)空间的建立,以及所引起的离散鞍点问题的共轭梯度迭代解法。这种区域分裂方法非常适用于网格不匹配情形。然后,对基于二阶鼓包(bump)函数近似的当地事后误差估算方法提出改进,以使之适应于区分裂情形,我们还将上述方法应用于广义Stokes问题。最后,给出了以误差估算值为准确则的自适应网格上的数值结果。 相似文献
3.
4.
一种包含运动边界的高精度流场数值计算方法 总被引:1,自引:0,他引:1
为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。 相似文献
5.
6.
对Goldstein提出的反馈力浸入边界法进行了新的思考,改进了其对力源项的计算,拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时,含有对速度误差的时间积分项,只能用于含时间项的Navier-Stokes (N-S)方程的求解,且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项,避免了时间相关的参数,使其不仅能适合非定常隐式时间推进,还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性,对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算,计算结果均与文献结果符合较好,表明本文改进的方法是有效的。得出的结论为:可以直接基于迭代次数进行反馈力源项的计算,改进的反馈力浸入边界法不仅可与非定常N-S方程结合,进行隐式求解,还可以与定常N-S方程结合用于定常流动的模拟,可将改进的方法运用到更多的流动问题当中。 相似文献
7.
适用于非周期流固耦合问题的时间谱方法 总被引:2,自引:1,他引:1
针对非定常流固耦合设计、分析问题所面临的计算效率与精度、鲁棒性之间的矛盾,采取将Chebyshev谱方法与基于非定常面元和几何非线性梁有限元模型的流固耦合分析方法相结合的方式,建立了针对大展弦比机翼非定常流固耦合优化设计问题的,可与伴随方法相结合的时间谱方法。Chebyshev谱方法直接对流固耦合的控制方程进行处理,利用Chebyshev算子替换系统状态变量,将非定常问题转化为Chebyshev控制点处耦合的定常问题。通过这种方式建立的流固耦合时间谱方法具有较高的计算精度、效率和足够的鲁棒性。验证算例及Goland机翼颤振速度计算实例表明,Chebyshev谱方法的计算精度随着Chebyshev控制点个数的增加而不断增大。只需选取较少的控制点,Chebyshev谱方法便可以达到满足精度要求的计算结果。与此同时,建立的流固耦合时间谱方法不仅适用于周期性非定常问题还适用于非周期性非定常问题。 相似文献
8.
9.
将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。 相似文献
10.
11.
研究了时间谱方法求解周期性非定常流场的计算效率,并对时间谱方法应用于周期性非定常流动的隐式求解方法进行探讨。当采样点数增加或减缩频率增大时,时间谱方法对应的雅可比矩阵对角占优性质迅速恶化,导致很多传统的迭代方法失效。为了解决上述问题,论文采用带预处理的广义极小残差(GMRES)算法来提高雅可比系数矩阵的计算收敛性。使用时间谱方法对NACA0012翼型强迫振荡算例进行计算,并与时域差分方法的计算效率和精度进行对比。研究表明在保证计算精度的同时,时间谱方法普遍可将计算效率提高一个量级左右。对于跨声速周期性流动,广义极小残差算法不论是稳定性还是收敛性都优于对称SGS迭代算法。 相似文献
12.
基于BVD原理的高保真空间重构方法 总被引:1,自引:0,他引:1
简要综述了一类基于单元边界变差最小化(Boundary Variation Diminishing,BVD)原理,设计双曲守恒律高保真数值格式的空间重构方法。BVD原理要求尽量减少通过重构得到的网格边界两侧物理量之间的差,从而能够有效地控制黎曼求解器中的数值黏性。BVD方法针对数值解的空间分布特征,选择多个函数作为空间重构的候补函数,并根据BVD判定准则从候补函数中选取最合适的函数进行空间重构。BVD判据不需要根据求解对象进行经验参数(阈值)的调整。选用适当的候补函数和BVD准则,可以完全避免现有算法中为抑制数值振荡而必须采用的非线性限制。BVD格式能在抑制数值振荡的同时,有效地控制数值耗散,可以对光滑解与间断解都获得高保真的计算结果。本文概述了BVD方法的基本思想、设计相关格式的基本思路,以及一些具有很强实用价值的BVD格式。并通过单相和两相可压缩流动的一些典型算例验证BVD格式的特点和优势。 相似文献
13.
14.
15.
16.
17.
18.
19.
基于自适应小波配点法和人工黏性,构造了二维激波问题的小波数值计算格式。利用最细尺度的小波系数构造两组激波定位函数,分别用以控制两个正交方向人工黏性的大小和分布。对三类二维激波问题进行计算和验证,结果表明:自适应小波配点法利用小波阈值滤波删除变化平缓区域大量网格点,而保留变化急剧区域的网格点,比传统方法计算效率高,且分辨率越高,计算效率更高;利用最细层小波系数构建的幂函数形式的激波定位函数能准确判断激波位置和控制X、Y方向人工黏性的大小和分布,从而捕捉不同方向和不同强度的激波。 相似文献
20.
基于谱方法的管路充填过程仿真 总被引:4,自引:1,他引:3
基于一维瞬变流理论,建立了推进剂供应系统管路充填过程数学模型.提出采用Fourier谱方法求解充填管道中液体瞬变流控制方程的新方法.利用坐标变换方法,很好地处理了解算管内气液移动边界时,非线性迭代收敛速度较慢这一难题.对由贮箱、隔离阀、管道组成的系统充填过程进行了仿真计算,并将计算结果与已发表的采用特征线方法获得的结果进行对比分析,表明该方法和仿真结果可信.进一步研究了预存气体压力、贮箱压力、多变指数等影响因素对最大充填压力峰的影响规律. 相似文献