implementation of upwind scheme in MATLAB

17 views (last 30 days)
nathan blanc
nathan blanc on 21 Dec 2021
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

Answers (0)

Products


Release

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by