[ufocfd]

I am developing a uniform compressible finite volume solver,

Runs ok when using a cell size of 2.0 but somethings
going wrong when I reduce the cell size to 0.1 units.

I think the CFL & timestep specification is ok but
I think its related to the update of the Q vector:

RHS = dE/dx + dF/dy + dG/dz

Q(n+1) = Q(n) + RHS*dt/vol, where dt = timestep and vol = cell volume

Is this correct ? Have I missed something obvious ? Thanks.

[FMDenaro]

Originally Posted by ufocfd

I am developing a uniform compressible finite volume solver,

Runs ok when using a cell size of 2.0 but somethings
going wrong when I reduce the cell size to 0.1 units.

I think the CFL & timestep specification is ok but
I think its related to the update of the Q vector:

RHS = dE/dx + dF/dy + dG/dz

Q(n+1) = Q(n) + RHS*dt/vol, where dt = timestep and vol = cell volume

Is this correct ? Have I missed something obvious ? Thanks.

I see a discrepancy, dimensionally your equation is not correct ...

If V is your Volume, then the equation writes:

d/dt [Integral_V Q dV] + Integral_V Div Flux dV = 0

Therefore:

d/dt [Q_average] + (1/vol) Integral_V Div Flux dV = 0

Hence, your RHS must be vol*(dE/dx + dF/dy + dG/dz)

[ufocfd]

thanks very much FMD !