“UNSTEADY_NAVIER_STOKES”脚本求解二维稳定Navier-Stokes方程。
空间离散化是通过标准伽辽金方法来实现的。
压力采用双线性四边形单元,速度采用双二次四边形单元。
这种压力速度插值满足所谓的LBB条件,保证了代数系统的可解性。
对于时间积分,采用了theta方法。根据theta的值,得到下列格式:
0 ->正向欧拉
1/2 ->曲克-尼科尔森
3/4 ->伽辽金
1 ->后向欧拉
有限元参数如有限元单元数、高斯积分点数等可以很容易地选择。
这些函数和例子是根据Jean Donea和Antonio Huerta的《流动问题的有限元方法》一书的第6章“粘性不可压缩流动”给出的。
如果你喜欢这个文件,请给我反馈。
引用作为
安德里亚·拉·斯皮纳(2023年)。2D非定常Navier-Stokes(//www.tatmou.com/matlabcentral/fileexchange/60869-2d-unsteady-navier-stokes), MATLAB中央文件交换。检索.