資源簡介
利用5階WENO有限差分法求解一維Euler方程組,求解Sod型激波管問題的詳細程序。
代碼片段和文件信息
function?[LRLambdaPLambdaNA]=Count1DLR(Ur)
for?i=1:size(U2)
????rho=U(1i);?E=U(3i);?u=U(2i)/U(1i);
????p=(r-1)*(E-rho.*u^2/2);?c=sqrt(r*p/rho);
????L{i}=(r-1)/(2*c^2)*[u^2/2+u*c/(r-1)?-u-c/(r-1)?1;?-u^2+2*c^2/(r-1)?2*u?-2;?u^2/2-u*c/(r-1)?-u+c/(r-1)?1];
????R{i}=[1?1?1;?u-c?u?u+c;?c^2/(r-1)+u^2/2-u*c?u^2/2?c^2/(r-1)+u^2/2+u*c];
????Lambda01=u-c;?Lambda02=u;?Lambda03=u+c;?Lambda=1.1*max(abs([Lambda01?Lambda02?Lambda03]));
????LambdaP(:i)=[Lambda01+Lambda;?Lambda02+Lambda;?Lambda03+Lambda]/2;
????LambdaN(:i)=[Lambda01-Lambda;?Lambda02-Lambda;?Lambda03-Lambda]/2;
????A{i}=[0?1?0;?(r-3)*u^2/2?(3-r)*u?r-1;?-u*c^2/(r-1)+(r-2)*u^3/2?c^2/(r-1)-(r-3/2)*u^2?r*u];
end
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????689??2014-02-06?09:57??Sod\Count1DLR.m
?????文件???????3199??2014-02-06?13:48??Sod\SodProgram.m
?????文件???????1498??2014-02-06?09:59??Sod\WENO51DEulerLu.m
?????目錄??????????0??2014-02-07?10:59??Sod
-----------?---------??----------?-----??----
?????????????????5386????????????????????4
- 上一篇:foa 果蠅優化算法matlab程序
- 下一篇:clippingFBMC
評論
共有 條評論