CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   Artificial Compressibility (https://www.cfd-online.com/Forums/main/9996-artificial-compressibility.html)

Ruben Paredes October 2, 2005 23:04

Artificial Compressibility
 
Hi all

I am trying to implement this method with FVM, in an unstructured cell centered grid, for a 2D viscous fluid. The inviscid flux is calculated with the Roe Method.

The problem is that my code doesn't converge, even when in the first iterations(5000) the residue seems to get smaller. Any idea?

Thanks in advance,

Ruben

A.S. October 3, 2005 00:42

Re: Artificial Compressibility
 
Hi,

Have a look at NSCKE it might help.

Regards

A.S.

Ruben October 3, 2005 11:55

Re: Artificial Compressibility
 
Sorry but where can I find this?

buch October 3, 2005 12:09

Re: Artificial Compressibility
 
Hi,

1) Can you please describe your resolution method ?

2) what is the value of the artificial compressibility parameter you choose ?

Regards

Buch

Ruben October 3, 2005 13:00

Re: Artificial Compressibility
 
Hi,

1) Can you please describe your resolution method ?

I am trying to solve the dimensionless N-S Eq. for an incompressible flow using the Artificial Compressibility concept with an triangular unstructured grid. In order to discretizate this Eqs. I used the Finite Volume method, and I get the following system:

dq/dt + Suma(E-Ev)

where E is the inviscid flux normal to the face Ev is the viscous flux normal to the face q the vector with (p,u,v) variables

The inviscid flux is calculated using the Roe method The viscous flux is computed by applying Gauss Theorem For the time integration it is used four-stage Runge-Kutta.

Boundary Conditions: I use ghost cells inflow : u=1.0 v=0.0 Pghost= P neig outflow: Ughost= Uneig , Vghost= V neig, P = 0.0 wall : Ughost= -Uneig , Vghost= -V neig, Pghost= Pneig

2) what is the value of the artificial compressibility parameter you choose ?

First, I tried to implement a formula proposed by Hino, Betha = (u2+v2), but now I am trying to keep this parameter constant.

Thanks in advance,

Ruben


versi October 3, 2005 21:51

Re: Artificial Compressibility
 
in Roe scheme, Q_L and Q_R at the cell face have to be reconstructed. It may be the ways of reconstruction that affect the convergence. How do you obtain Q_L and Q_R?

Ruben October 4, 2005 10:23

Re: Artificial Compressibility
 
QL is just the value of the neighbour cell and QR is the cell analized. is that right?


All times are GMT -4. The time now is 05:24.