Rotation in MAESTROeX
Rotation
To handle rotation in MAESTROeX we move to the co-rotating reference frame. Time derivatives of a vector in the inertial frame are related to those in the co-rotating frame by
where i(rot) refers to the inertial(co-rotating) frame and
Plugging this into the momentum equation and making use of the continuity equation we have a momentum equation in the co-rotating frame:
The Coriolis and Centrifugal terms are force terms that will be added to right hand side of the equations in mk_vel_force. Note that the Centrifugal term can be rewritten as a gradient of a potential and absorbed into either the pressure or gravitational term to create an effective pressure or effective gravitational potential. Furthermore, because it can be written as a gradient, this term would drop out completely in the projection if we were doing incompressible flow. In what follows we will include the Centrifugal term explicitly.
Using Spherical Geometry
In spherical geometry implementing rotation is straightforward as we don’t have
to worry about special boundary conditions. We assume a geometry as shown in
Fig. 27 where the primed coordinate system is
the simulation domain coordinate system, the unprimed system is the primed
system translated by the vector
The direction of
The magnitude of
The Coriolis term is a straightforward cross-product: