implementation of upwind scheme in MATLAB

13 次查看(过去 30 天)
More of a CFD or numerical analysis question than a proper MATLAB question but I will try my luck.
I am trying to implement an upwind scheme in MATLAB, in order to solve a 1D compressible navier stokes equation, with a periodic boundary condition:
I am trying to perform the spatial diffrentiation using an upwind scheme, see the attahced file "upWindDiffv2". the scheme should be TVD or numericaly stable but it isn't, I get substantial numerical oscilations and the code diverges. see the attached file "main", it includes a "video" of the solution where the solution clearly diverges.
I am performing the temporal integration using eulers method. I tried other temporal integration methods and they did not work either. note that the CFL criterion is more than satisfied.
I am not sure what I am doing wrong. would be thankfull for any assistance.
Nathan

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Computational Fluid Dynamics (CFD) 的更多信息

产品


版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by