浸入边界法（IBM）二维matlab程序

直接运行 处动画，便与学习

function [uuuuu]=fluid（uff）
global a dt rho mu;
w=u-（dt/2）*skew（u）+（dt/（2*rho））*ff;
w=fft（w[]1）;
w=fft（w[]2）;
uu（::1）=a（::11）.*w（::1）+a（::12）.*w（::2）;
uu（::2）=a（::21）.*w（::1）+a（::22）.*w（::2）;
uu=ifft（uu[]2）;
uu=real（ifft（uu[]1））;
w=u-dt*skew（uu）+（dt/rho）*ff+（dt/2）*（mu/rho）*laplacian（u）;
w=fft（w[]1）;
w=fft（w[]2）;
uuu（::1）=a（::11）.*w（::1）+a（::12）.*w（::2）;
uuu（::2）=a（::21）.*w（::1）+a（::22）.*w（::2）;
uuu=ifft（uuu[]2）;
uuu=real（ifft（uuu[]1））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        506  2008-01-28 03:59  2D IBM 程序\fluid.m

     文件         86  2008-01-27 10:14  2D IBM 程序\Force.m

     文件        486  2008-02-18 05:23  2D IBM 程序\ib2D.m

     文件        841  2011-02-09 04:47  2D IBM 程序\initialize.m

     文件        800  2008-01-28 02:05  2D IBM 程序\init_a.m

     文件        276  2011-02-09 02:45  2D IBM 程序\interp.m

     文件        235  2008-01-27 07:50  2D IBM 程序\laplacian.m

     文件        134  2008-01-28 05:47  2D IBM 程序\phi1.m

     文件        134  2008-01-28 05:48  2D IBM 程序\phi2.m

     文件        184  2008-01-28 02:28  2D IBM 程序\sk.m

     文件        133  2008-01-28 02:02  2D IBM 程序\skew.m

     文件        302  2011-02-09 02:44  2D IBM 程序\spread.m

     目录          0  2015-10-05 14:33  2D IBM 程序

----------- ---------  ---------- -----  ----

                 4117                    13