python for CFD(前两步问题总结)

在前面的问题中,如果单方面将nx无限变大,(不妨令nx = 100)则结果画出的图形是发散的。

这是因为有限差分中有一个重要的必要条件,CFL数,具体是差分方程的依赖域必须包含相应微分方程的依赖域,最简单可以理解为时间推进求解的速度(dx/dt)必须大于物理扰动传播的速度(u),只有这样才能将物理上所有的扰动俘获到。

每一个时间间隔dt上“波”传播的距离是dx,而dxnx相关,保持数值求解稳定的关系式为:

101058_afqJ_1427937.jpg

import numpy as npimport matplotlib.pyplot as plt
def linearconv(nx):
    dx = 2./(nx-1)
    nt = 20    #nt is the number of timesteps we want to calculate
    c = 1
    sigma = .5   
    
    dt = sigma*dx     #跟步骤2唯一的变化在这

    u = np.ones(nx) 
    u[.5/dx : 1/dx+1]=2

    un = np.ones(nx)

    for n in range(nt):  #iterate through time
        un = u.copy() ##copy the existing values of u into un
        for i in range(1,nx):
            u[i] = un[i]-c*dt/dx*(un[i]-un[i-1])
        
    plt.plot(np.linspace(0,2,nx),u)

linearconv(121)即可,此时图形和步骤2类似


转载于:https://my.oschina.net/sunxichao/blog/374307

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值