Hi Swagat, The MovingReferenceFrame works only with frames moving with constant frame speed. What are the additional terms you need to add? Cheers, Stan On pią, 2019-05-31 at 18:28 +0900, Swagat Kumar Nayak wrote:

Hi all,

I intend to solve the Incompressible Navier Stokes equations in a rotating reference frame. Hence, I am adding the frame velocity of the rotating frame as a Force Type = "MovingReferenceFrame". Will it add the appropriate Coriolis term? Also, does the implementation seem correct?

<NEKTAR> <!-- Setting the Expansion Bases --> <EXPANSIONS> ............................. </EXPANSIONS>

<!-- Setting the Filters to output Energy Details --> <FILTERS> ............................. </FILTERS>

<!-- Setting the Boundary conditions and Solver information for the flow --> <CONDITIONS> <SOLVERINFO> .............................. </SOLVERINFO>

<PARAMETERS> ............................. </PARAMETERS>

<!-- Setting the Variables -->

<VARIABLES> <V ID="0"> u </V> <V ID="1"> v </V> <V ID="2"> w </V> <V ID="3"> p </V> </VARIABLES>

<!-- Setting the Boundary Conditions --> <BOUNDARYREGIONS> ........................................ </BOUNDARYREGIONS>

<BOUNDARYCONDITIONS> ................................................. </BOUNDARYCONDITIONS>

<!-- Setting the Initial conditions and the frame velocity (Precession frame) -->

<FUNCTION NAME="InitialConditions"> <E VAR="u" VALUE="F(y)" /> <E VAR="v" VALUE="G(x)" /> <E VAR="w" VALUE="0" /> <E VAR="p" VALUE="0" /> </FUNCTION>

<FORCE TYPE="MovingReferenceFrame"> <FRAMEVELOCITY> <FUNCTION NAME="MovingReferenceFrame"> <E VAR="u" VALUE="f(y,z,t)" /> <E VAR="v" VALUE="g(y,z,t)" /> <E VAR="w" VALUE="0" /> </FUNCTION> <FRAMEVELOCITY/> </FORCE>

</CONDITIONS> </NEKTAR>

Thank you ! -- Swagat Kumar Nayak Department of Mechanical Engineering IIT Kharagpur

_______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/nektar-users