[syoo]

Hello,

I am currently trying to code an 2D Euler-eq solver using the Beam-Warming scheme for a given deforming mesh. So far, I can't preserve freestream using my solver (embarrassing ...) for deforming grids. The freestream is conserved for a non-deforming grid.

My current understanding is that the mesh deformation (mesh velocity) is accounted for in the solution through the Jacobian which transforms the grid from the cartesian coordinate to generalized body-fitted coordinate. I also satisfy the Geometric Conservation Law by calculating the grid metrics written the conservative form, suggested by Thomas and Lombard.

Should I be including the grid velocity in the flux terms? I have went to many literatures that use Beam-Warming scheme and none of them mention this.

Any help, suggestion, or comment will be very much appreciated.

thank you.

-Paul