对流方程的有限差分数值解法（步长定律、 固有差分格式、matlab 程序.PDF

对流方程的有限差分数值解法（步长定律、 固有差分格式、matlab 程序和输出图形） mainc@126.com /m/user_content.aspx?id=404234 1 一维对流方程 1.1 一维对流方程的形式： ， 其中，u 代表物质的量或者其他参量，a 代表物质的运动速度。此一维对流方程仅仅表示物 质的运动情况，而与边界条件或是约束条件无关。当a 为常数时，此一维对流方程为一维常 系数对流方程，当a 不为常数时，方程为一维变系数对流方程。在不考虑边界条件或约束条 件的情况下，无论a 是否为常数，此对流方程本身的数值解法存在固有的有限差分格式，条 件是要满足步长定律。 1.2 一维常系数对流方程的步长定律： ，其中 ∆t 为时间步长，∆x 为空间步长。 1.3 一维常系数对流方程的固有差分格式： 或者 1.4 有以下方程： 满足初始条件: x=1:1:113， u0=[0 0 1.8393e-005 3.7259e-005 7.2772e-005 1.3704e-004 2.4885e-004 4.3574e-00 4 7.3575e-004 0.0012 0.0019 0.0028 0.0042 0.0059 0.0080 0.0104 0.0132 0.016 1 0.0189 0.0215 0.0235 0.0248 0.0253 0.0248 0.0235 0.0215 0.0189 0.0161 0.0 132 0.0104 0.0080 0.0059 0.0042 0.0028 0.0019 0.0012 7.3575e-004 4.3574e-00 4 2.4885e-004 1.3704e-004 7.2772e-005 3.7259e-005 1.8393e-005 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]， 和边界条件： u(1,t)=u(113,t)=0. MATLAB 程序： u0 = [0 0 1.8393e-005 3.7259e-005 7.2772e-005 1.3704e-004 2.4885e-004 4.3574e-004 7.3575e-004 0.0012 ... 0.0019 0.0028 0.0042 0.0059 0.0080 0.0104 0.0132 0.0161 0.0189 0.0215 ... 0.0235 0.0248 0.0253 0.0248 0.0235 0.0215 0.0189 0.0161 0.0132 0.0104 ... 0.0080 0.0059 0.0042 0.0028 0.0019 0.0012 7.3575e-004 4.3574e-004 2.4885e-004 1.3704e-004 ... 7.2772e-005 3.7259e-005 1.8393e-005 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0 0 0