The aim of this test case is to validate the following parameters of incompressible laminar couette flow in the annular region of two infinite cylinders (one stationary, one rotating):

Velocity profile

Pressure profile

The presence of both stationary and rotating components in the geometry necessitates the need for distinct representative zones. SimScale provides two options to approach this problem: the Multi-Reference Frame (MRF) method and the Arbitrary Mesh Interface (AMI) method. Briefly stated, MRF uses source terms for coriolis and centrifugal forces to model the rotation, while AMI explicitly ‘rotates’ the appropriate mesh, interpolating values at each iteration over the interface between the fixed and rotating zones. AMI is inherently transient, but MRF allows for cheaper and quicker results. Thus, MRF uses the steady-state simpleFoam solver in contrast to the pimpleDyMFoam solver used by AMI.

The geometry was uploaded on to the SimScale platform, and appropriate meshing techniques used to generate stationary and rotating zones. Simulation results of SimScale were compared to analytical results obtained from methods elucidated in [1].

The infinite cylinders were modeled as a a slice of their annulus (see Fig.1.). The dimensions are shown in Table 1. The outer cylinder remains stationary, while the inner rotates with an angular velocity of 0.001rad.s−1

$0.001\text{}rad.{s}^{-1}$

.

Analysis type and Domain

The geometry was uploaded on the SimScale platform, and appropriate meshing techniques were employed to generate fixed and rotating zones in the mesh.

For both MRF and AMI the rotating mesh zone is assigned an angular velocity of 0.001rad.s−1

$0.001rad.{s}^{-1}$

about the y-axis, although it is interpreted differently by the two cases. The velocity and pressure bounary conditions for MRF and AMI have been stated in Table 3.

Table 3: Boundary Conditions for MRF and AMI Cases¶

Parameter

Inner Wall (A)

Outer Wall (B)

Sides (C)

Velocity

0.0ms−1$0.0\text{}m{s}^{-1}$

(Moving wall velocity)

0.0ms−1$0.0\text{}m{s}^{-1}$

Symmetry

Pressure

Zero Gradient

Zero Gradient

Symmetry

Results

The analytical solution for Taylor-Couette flow is computed from the simplified Navier-Stokes in cylindrical coordinates [1]. The velocity profile is given by:

A comparison of the velocity and pressure variation in the radial direction obtained with SimScale (using MRF and AMI) with analytical results is given in Fig.3A and 3B.

Fig.3. Visualization of velocity and pressure (A, B) in the radial direction for MRF and AMI

This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software and owner of the OPENFOAM® and OpenCFD® trade marks. OPENFOAM® is a registered trade mark of OpenCFD Limited, producer and distributor of the OpenFOAM software.

Contents

Cookie Policy

SimScale uses cookies to improve your user experience. By using this website you consent to our cookie policy.

This website uses cookies so that we can provide you with the best user experience possible. Cookie information is stored in your browser and performs functions such as recognising you when you return to our website and helping our team to understand which sections of the website you find most interesting and useful.

Strictly Necessary Cookies

Strictly Necessary Cookie should be enabled at all times so that we can save your preferences for cookie settings.

If you disable this cookie, we will not be able to save your preferences. This means that every time you visit this website you will need to enable or disable cookies again.