finite volume solver problem
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.
If V is your Volume, then the equation writes:
d/dt [Integral_V Q dV] + Integral_V Div Flux dV = 0
d/dt [Q_average] + (1/vol) Integral_V Div Flux dV = 0
Hence, your RHS must be vol*(dE/dx + dF/dy + dG/dz)
thanks very much FMD !
|All times are GMT -4. The time now is 17:26.|