CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   MacCormick (https://www.cfd-online.com/Forums/main/13784-maccormick.html)

Shuo July 8, 2007 07:51

MacCormick
 
In the MacCormick upwind scheme.

U^(n+1) = U^n - delta_t/delta_x*(f_(i+1) - f_i)

is it f_(i+1) - f_i = A_(i+0.5)*(U_(i+1) - U_i) where A = df/dU

or f = {rho*u, rho*u^2 + P...}^T


Shuo July 9, 2007 02:22

Re: MacCormick
 
it appears that the second order upwind schemes doesn't work for the shock tube problem.

Predictor:

U^(n+1) = U^n - delta_t/delta_x*(f_(i+1)^n - f_i^n)

Corrector:

U^(n+2) = 0.5*(U^n + U^(n+1)) - delta_t/delta_x*(f_(i+1)^(n+1) - f_i^(n+1))

However the flux-splitting version does work. Is there anything i am doing wrong?

Implicit:



All times are GMT -4. The time now is 11:50.