- **Main CFD Forum**
(*http://www.cfd-online.com/Forums/main/*)

- - **4th order Runge-Kutte & uncoupled method Navier Stokes equations**
(*http://www.cfd-online.com/Forums/main/72564-4th-order-runge-kutte-uncoupled-method-navier-stokes-equations.html*)

4th order Runge-Kutte & uncoupled method Navier Stokes equationsHello CFD's users!
I have programmed a code in order to solve Navier-Stokes equations with Finite Element Method. In order to solve it, I use an uncoupled method (first I calculate the velocity, secondly, the PPE to obtain the pressure), to make time discretization I have applied 4th explicit Runge Kutta method. But I don't know if I am applying rightly Runge-Kutta with the uncoupled method. I do: 1) Calculate velocity u_n+1 doing: 1.a) k1=incrT*(-Advection_Term(n,u_n)-Difusion_Term(n,u_n)-Gradient*p_n) 1.b) k2=incrT*(-Advection_Term(n+1/2,u_n+k1/2 )-Difusion_Term(n+1/2,u_n+k1/2 )-Gradient*p_n) 1.c) k3=incrT*(-Advection_Term(n+1/2,u_n+k2/2 )-Difusion_Term(n+1/2,u_n+k2/2 )-Gradient*p_n) 1.c) k4=incrT*(-Advection_Term(n+1,u_n+k3 )-Difusion_Term(n+1,u_n+k3 )-Gradient*p_n) In that way, we have Mu_n+1=Mu_n+k1/6+k2/3+k3/3+k4/6, where M is the mass matrix. 2) Calculate pressure p_n+1 with velocity u_n+1, throughout the Pressure Poisson Equations (PPE). I would like to know if I am doing it in the right way, does it make sense? Thanks in advance! Best regards! Isi |

All times are GMT -4. The time now is 00:41. |