CFD Online Logo CFD Online URL
Home > Forums > SU2

Rotating reference frame

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   February 25, 2013, 21:48
Default Rotating reference frame
New Member
Join Date: May 2012
Posts: 6
Rep Power: 5
amphy404 is on a distinguished road

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.

    /*--- Adjustment for a rotating frame ---*/
    if (rotating_frame) {
        ProjVelocity = Rot_Flux;
        for (iVar = 0; iVar < nVar; iVar++) {
            val_resconv[iVar] -= ProjVelocity * 0.5*(U_i[iVar]+U_j[iVar]);
            if (implicit) {
                val_Jacobian_i[iVar][iVar] -= 0.5*ProjVelocity;
                val_Jacobian_j[iVar][iVar] -= 0.5*ProjVelocity;
From this code, it looks like the rotation is accounted for by subtracting the rotational velocity (Rot_Flux) from the projected flux (val_resconv). This makes sense for the continuity and momentum equations. Perhaps it is just my misunderstanding, but for the energy equation shouldn't it be subtracting the enthalpy, rather than the internal energy?

That is,
val_resconv[nVar-1] -= ProjVelocity * MeanDensity * MeanEnthalpy;
amphy404 is offline   Reply With Quote

Old   February 26, 2013, 07:10
Super Moderator
Thomas D. Economon
Join Date: Jan 2013
Location: Stanford, CA
Posts: 260
Rep Power: 5
economon is on a distinguished road

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:

Hope this helps!
economon is offline   Reply With Quote


moving mesh, rotating frame

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Second Derivative Zero - Boundary Condition fu-ki-pa OpenFOAM 10 May 1, 2014 16:26
OpenFOAM 1.6-ext git installation on Ubuntu 11.10 x64 Attesz OpenFOAM Installation 45 January 13, 2012 13:38
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 00:35
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08
G95 + CGNS Bruno Main CFD Forum 1 January 30, 2007 01:34

All times are GMT -4. The time now is 14:52.