Rotating reference frame
I have a question about the rotating reference frame source code for the space centered schemes in numerics_convective.cpp.
For instance, in numerics_convective.cpp: L3198-3208.
Thanks for the question. In SU2, we are using a version of the equations expressed in a non-interial frame while still solving for the absolute velocities as part of the conserved variables (an 'absolute velocity' formulation). The transformation of the equations to this form results in not only an adjustment of the convective fluxes to take into account the rotational velocity of the mesh faces, but also changes in the momentum and energy equations. More specifically, the momentum equations gain a source term, and the energy flux is slightly modified to become density*enthalpy*(velocity - rotational_velocity) + rotational_velocity*pressure, which could also be written as density*enthalpy*velocity - rotational_velocity*density*energy after substituting the definition of enthalpy. The latter expression is implemented in SU2.
If you are interested in more details on the formulation, there is a nice description in "A Three-Dimensional Euler Solver for Turbomachinery Blade Rows" by D. G. Holmes & S. S. Tong, or please see the following paper using the formulation within SU2: http://adl.stanford.edu/papers/AIAA-2012-3018.pdf.
Hope this helps!
|All times are GMT -4. The time now is 13:16.|