Crank-Nicolson格式作为隐格式,有无条件稳定的特性,也就是不再受到\Delta t不能过大的影响,比如以下实验,固定住\Delta t, 减小\Delta x, 结果如下: CN格式还是非常稳的 代码在这里: importnumpyasnpimportmatplotlib.pyplotaspltfromscipy.sparseimportdiags,linalga=1u_truth=lambdax,t:x*(x*x+np.exp(t))...
2.2Crank-Nicolson格式 对上式偏导在进行离散, i 表示网格点,共 N+1 个, n 为求解时间步,有: {∂u∂t=uin+1−uinΔt∂2u∂y2=[12(ui+1n+ui+1n+1)−(uin+uin+1)+12(ui−1n+ui−1n+1)]Δy2 从而得到Crank-Nicolson求解格式: −Δt2ReΔy2ui+1n+1+(1+ΔtReΔ...
摘自《FLUENT流体工程仿真计算实例与分析》,程序略有修改两个间距为1cm水平平板,如下图所示:上板匀速平板间流动(Crank-Nicolson格式)【转载】">充满着运动黏度系数υ=1cm2/s的液体。上板做水平运动并在0.1s时间内,速度线性由0线性地增加到10cm/s,如下图
对于扩散方程（包括许多其他方程），可以证明Crank-Nicolson方法无条件稳定。但是，如果时间步长与空间步长平方的比值过大（一般地，大于1/2），近似解中将存在虚假的振荡或衰减。基于这个原因，当要求大时间步或高空间分辨率的时候，往往会采用数值精确较差的后向欧拉方法进行计算，这样即可以保证稳定，又避免...
机器永远是对的,未测试的代码永远是错的 数学- 微分方程数值解 - 第 4 章 抛物型方程的差分解法 - 4.5 Crank-Nicolson 格式 4.5 Crank-Nicolson 格式 本节对于定解问题 (3.1.1)∼(3.1.3)(3.1.1)∼(3.1.3) 建立一个具有 O(τ2+h2)O(τ2+h2) 精度的无条件稳定的差分格式。 注意,对各个符号取...
【内容介绍】本资源主要利用MATLAB的实时脚本编程实现了抛物型偏微分方程数值求解,以图-文-代码三者互相嵌套的形式介绍实现过程,一目了然。包括对迭代的误差分析。 【适用对象】工科生、数学专业等。 【算法涵盖】4点显式差分格式、4点隐式差分格式、Crank-Nicolson格式 感谢支持!
下载得到文件列表 瞬变电磁Crank-Nicolson FDTD三维正演.pdf 相关文档 文档介绍文档介绍:万方数据瞬变电磁狽三维正演孙怀凤,柳尚斌,杨洋地球物理学报狽第卷第年孙怀风,柳尚斌,杨洋..瞬变电磁荩厍蛭锢硌Пǎ:琩:./‘,.篽,—琩:./摘要时域有限差分方法使用衿史值绱懦〉目占洳裳ü奔洳降迪值绱懦∈的D猓...