Turbomachines

CFD Simulation approaches for Turbomachines

Multiple Reference Frame, Sliding Mesh Motion

Table of Contents: MRF Approach - Limitations | Roll, Pitch and Yaw | Flow inside a centrifugal blower for HVAC applications | Dynamic Mesh Motion | Cavitation in Pumps | Mixing Tanks or Stirred Tanks | Propellers

CFD simulation approach for turbomachines such as centrifugal pump and blowers, appropriateness of various modeling approaches namely Single Reference Frame (SRF), Multiple Reference Frame (MRF) or Frozen Rotor Method, Sliding Mesh Motion (SMM) along with applications to industrial problems are described in this page. Before moving to the CFD simulation aspects, following two graphics summarize the types of turbomachines available in industries. The simulation can be just one phase (liquid or gas) or multi-phase applications such as investigation into mixing effects in Gas-Liquid-Solid Stirred Reactor. The explanation on multi-phase flow simulations can be found here.

Pumps need to transfer energy to the liquids and their suction (inlet) lines are at lower pressures than the discharge lines. In general, it is assumed that the inlet line is always filled with water without any traces of air or gas. However, this is not always the case. If the centreline of the pump is above the water supply reservoir and there is no check valve or the pump is started first time, there will be air present in the suction line. Centrifugal pumps with basic design features cannot remove the air from the suction line and can only start to pump the fluid after it has initially been primed with the fluid.

The term 'priming' refers to the phenomena of removing the air from the system. Types of self-priming pumps are centrifugal pumps with "separation chamber", "side channel" or "water ring pumps" and "two casing chambers and an open impeller pumps". Some self-priming pumps come with an integrated vacuum pump that ensure the pump primes at unfavorable suction pipe layouts. The self-priming feature can be imparted to a centrifugal pump by ensuring impeller is submerged in water or retains enough water when it stops. It can be achieve simply by designing a suction and discharge cavity above the centerline of the impeller or installing a check valve near the suction eye. Note that even self-priming pumps will need initial priming after commissioning.

CFD simulations of a pump is carried with known mass flow rate and known pressure at one of the boundaries. However, any simulation of a self-priming process needs to be carried out with atmospheric boundary conditions at both inlet and outlet as well as keeping the gravity on - to separate gas from liquid. Alternatively, pressure inlet and the opening at the outlet can be used as inlet and outlet boundaries in order to approximate the actual self-priming operation This makes the simulation process special and a transient simulation is necessary.

Pumps are characterized by a single value known as specific speed. There are some dimensonal variants of specific speeds and some non-dimensional. Hence, utmost care should be taken to check and use appropriate units used to derive the specific speeds. ISO 5801:2017 specifies procedures for the determination of the performance of fans of all types except those designed solely for air circulation, e.g. ceiling fans and table fans. Testing of jet fans is described in ISO 13350. The aeroacoustic measurements are performed using a test rig (the so-called In-Duct method) according to the industry norm ISO 5136 where farfield noise levels are recorded inside a circular duct using three slit-tube microphones.

CFD Simulation Methods for Turbomachines

SRF: This method is used when the computational domain is axi-symmetric. This is called 'single' reference frame because only one reference frame (which is rotating) needs to be defined. This method can be used when whole geometry (the computational domain) can be assumed rotating.

Lumped Fan (FP) or Body Force Model: This method does not model the blades and hub, instead uses an interface at the location where the fan blades would be located. An experimentally obtained fan curve (pressure drop vs. flow rate) is applies as pressure rise (pressure jump or porous jump) over the interface according to the mass flow rate. The drawbacks of this model are that the velocity vectors exiting the interface of the standard LF model only have an axial component, no swirl (tangential component) and no reduced velocity where the hub would be located. These drawbacks can be reduced by adding a swirl to the outlet flow or geometrically leaving out the hub region from the interface.

This method uses more that one reference frames - at least one stationary (outer) and 1 rotating (inner). This is also known as Frozen Rotor Method (FRM) as the rotating parts are kept frozen in position and rotation is accounted for by the additional source terms through inclusion of centrifugal and Coriolis forces. Instead of the blades moving physically through the air, the air moves around the rotor blades with a corresponding angular velocity. Even for the cases where transient simulation is required, MRF method is useful for attaining initial values for time-dependent simulations because the pseudo-steady state can be reached within a few revolutions starting from zero initial velocity. MFR approach is appropriate if the flow seen from the rotating/moving frame of reference is steady along the interface. In other words, when flow is relatively uniform at the interface between the moving and stationary zones. E.g. in mixing tanks the impeller-baffle interactions are relatively weak, large-scale transient effects are not present and the MRF model can be used.

This paper gives a comparison of MRF and SMM for a centrifugal fan including description of the boundary conditions used for simulations.

• The equations for the inner rotating region are solved in a rotating reference frame (when a rotating reference frame is used the centrifugal and coriolis forces are included in the momentum equation). The equations for the outer stationary region are solved using a stationary reference frame. The solutions from both regions are matched at the interface between the rotating and stationary regions via velocity transformation from one frame to the other. This velocity-matching step involves the implicit assumption of a steady flow condition at the interface. The local velocity value may depend on the position of the impeller blade.
• Domain and Mesh: The accuracy of simulation result is sensitive to the selection of MRF region. It is defined so that it includes rotating parts and extends to upstream and downstream sections where a mixed-out (relatively uniform) flow field exists. MRF region can also include stationary parts if they are rotationally symmetric. The mesh (elements and nodes) inside the rotating domain including external or internal boundaries rotate as a solid body with the rotation angular speed specified about the axis of rotation.
• Walls: In order to account for wall shear, solver needs to know the speed of the wall. Any external or internal wall inside the rotating domain is assumed to rotate at the speed of the domain. If a wall defining the boundary of rotating domain is required to be defined stationary, it either needs to be specified as "counter-rotating wall" or its angular velocity needs to be defined as "0 [rad/s]". In FLUENT, this can be achieved using the feature motion "Relative to Adjacet Cell Zone" or 'Absolute'.
• Governing Equations: In both SRF and MRF methods, the momentum conservation is governed by the Navier-Stokes equations and the mass conservation is governed by the continuity equation. Both SRF and MRF methods are called steady-state approach (or pseudo-transient due to usage of psuedo time stepping) where the solution is time independent and rotation is achieved through mesh fixed in space and time. As a result, it does not model unsteady (transient) effects at the frame change interface.
• Result is influenced by 'frozen' position of rotating components (blade) with respect to the stationary parts (baffles, involute). Hence, the validity of result should be thoroughly checked when either the number of blade count is low or the speed of rotation is relatively low.
• MRF result will be inaccurate when flow crosses the interface from both directions - that is when flow enters and leaves the outer boundary of rotating domain.
• As per ANSYS user manaul, the use of the realizable k-ε turbulence model with multiple reference frames is not recommended. The PRESTO! scheme for the discretization of pressure equation is recommended for solving rotary machines.

MRF Approach - Limitations:

The fan performance curve, describing the pressure rise over volume flow rate through the fan, is commonly under-predicted when using the MRF model. This under-prediction mainly occurs at low to medium volume flow rates (radial and transitional regime), therefore it is concluded that the MRF method works best when the fan is operating in axial conditions. Liu et al. (2016) found that the thrust of a tidal current turbine was equally under-predicted with the MRF approach as the pressure rise. Furthermore, they found that the frozen rotor position has a significant influence on the flow field in the near wake. Kobayashi and Kohri (2011) showed that uniform inflow conditions are necessary in order to facilitate the transition to the moving reference frame. Apart from the operating conditions, it is also important that flow structures originating from the blades (e.g. tip vortices) are completely encompassed by the MRF domain, in order not to be split or in other ways being hindered from developing. As was pointed out by Gullberg and Löfdahl (2011) and Kobayashi et al. (2014) the performance of the MRF approach is highly dependent on the users choice of the size of the MRF. A large MRF domain has therefore been shown to give better agreement with experimentally obtained fan curves. For open fans (i.e. no ring connecting the blade tips), the radial extent is more important, due to the occurring tip vortices, while for both the open and closed fan type a sufficiently large axial extent can lead to more uniform inflow conditions and hence a better performance prediction. Therefore it can be concluded that the users choice of the MRF domain has a substantial impact on the accuracy of the results.

• Note that the mesh motion can be constant speed or accelerating – the solver accommodate both situations.
• In the SMM formulation, the motion(s) of moving zone(s) is tracked relative to the stationary reference frame where the motion of any point or node in the domain is given by a time rate of change of the position vector - known as grid speed.
• SMM and DMM uses same equations where in case of DMM, has additional feature for nodes to move relative to each other. Hence, SMM can be assumed to be a subset of more general DMM method.
• Since this approach involves mesh motion it is necessary that the interface between the zones be two overlapping faces, one of which is attached to the rotating region while the other is attached to the stationary region. The interaction between the two regions is modeled by interpolating the information across these faces. No additional source terms needs to be added into the momentum equations. For each time step, nodes are rotated and the fluxes at the sliding interfaces (interface at the stationary and rotating) are recomputed.
• Time steps for transient simulation is a function of element size (Δs) at the sliding interface. The time steps should always be less that it requires a moving / rotating cells to cross past a stationary point at the interface that is Δt ≤ Δs/ω/r where ω is the rotational speed of moving domain and r is the radius of sliding interface.
• The mesh interface between rotating and stationary domains must be situated such that there is no motion normal to it.
• The common boundary or mesh interface can be any shape (including non-planar surfaces) provided the two interface boundaries are derived from a set of common geometrical entities (lines in 2D and surfaces in 3D).
• By default, the velocity of wall is zero relative to the motion of the mesh or cell zone it is attached to. For walls bounding a moving mesh this imposes the "no-slip" condition in the reference frame of the mesh. Hence, one need to modify the wall velocity boundary conditions only if it is stationary in the absolute frame and therefore moving in the relative frame. For example, shroud of a fan is a stationary wall bounding the rotating reference frame defined for the blades. Hence, its velocity should be explicitly set to ZERO.

• An interface is a MUST between rotating and stationary regions (domains or zones). Sometimes, such surface zones can be set as 'interior' instead of an interface.
• However, such 'interfaces' must be a surface of revolution having axis of revolution coinciding with the axis of revolution of the rotating zones and walls.
• The recommended practice is to chose the location of such interface(s) at the mid-way between the tip of rotating walls (or the blade tip) and the nearest stationary housing walls. Sometimes, the location of interfaces are also governed by meshing considerations (such as number of boundary layers on the walls of the blades) and free mesh size beyond this region.
• For rotating domains embedded in relative large domains (such as a fan relative very small as compared to room), the recommendation is to have the interface at a location where flow is likely to be uniform.
• Note that there is a significant difference between "Rotating Frame" and "Rotating Mesh".

DMM:

1. Prescribed translation (such as motion of piston and valves in engines), rotation (motion of teeth in gear pumps) or combination of translation and rotation (motion of teeth in helical pumps). In ANSYS FLUENT, "DEFINE_CG_MOTION” is used to describe the motion of a face or cell zone. It provides solver with a velocity at every time step. ANSYS Fluent then translates this to the position of the zone given the current simulation time. DEFINE_CG_MOTION(udf_name, dt, vel, omega, time, dtime), where name defines the name of the UDF. 'dt' is a pointer to a storage that contains the dynamic mesh attributes specified by the user or calculated by FLUENT. 'vel' and 'omega' are the linear and angular velocity arrays respectivly. 'time' defines the current time and 'dtime' defines the time step. These 6 variables are passed to the UDF by FLUENT, the UDF then calculates the linear and angular velocities and returns the values to ANSYS FLUENT.
   • Since movement is known a priori, a sequence of meshes can be generated to accommodate the motion of the boundaries. In between those generated meshes, the mesh can be deformed (for example stretching and compression) according to the prescribed boundary movement.
   • In situation where sliding interfaces need to be created in very tight clearances such as screw pumps, continuous mesh adaptation - even at every time step - might be required.
   • It is when the deformation of the mesh is excessive that the deformed mesh can be replaced with a new generated one (e.g. layer addition and deletion in case of translation of boundary or linear deformation of the fluid zone).
   • Note that the motion is prescribed at the solid boundaries and not at the mesh representing fluid volume which shape and size changes.
   • The motion of each point and cell in the mesh needs to be calculated so that the mesh remains topologically consistent to the changes in shape and size of the fluid domain.
   • During the mesh deformation process, a finite element method is used based to the Laplacian equation (∇.κ∇u = 0) for estimating the velocity and location of the mesh points, either with a constant or variable diffusivity κ. The added diffusivity controls the mesh deformation and quality.
   • If the zone is a rigid body, a 'profile' or 'user-defined function (UDF)' can be used to define the motion of the rigid body or use the 6DOF solver. The motion of a solid-body can be specified either as a profile or as a user-defined function (UDF). A profile may be defined by the following profile fields: time (time), crank angle (angle) (in-cylinder flows only), position (x, y, z), linear velocity (v_x, v_y, v_z), angular velocity (omega_x, omega_y, omega_z), orientation (theta_x, theta_y, theta_z).
2. Unprescribed motion where the subsequent motion is determined based on the solution at each time step(for example, the linear and angular velocities are calculated from the force balance on a swing valve), as is done by the six degree of freedom (6DOF) solver. Motion induced by flow field which is dependent of the solution of flow field itself and found application in phenomena such as aero-elesticity and fluid-structure interactions (FSI). A noticeable difference with prescribed mesh motion case is that the deformation in this category of applications are low and as well as not known a-priori. If there is excessive deformation in solid, a high quality hexa mesh would be required.
3. Overset mesh: this is a dynamic mesh motion option available in Star-CCM+. The Overset mesh consists always of at least two regions partially or fully overlapping each other and is useful when working with moving bodies such as swing valves. It does not require mesh modification after initialization or during the simulation. The solutions in the overlapping mesh zones are interpolated and hence cell sizes in the two meshes should be nearly same.
   • The boundaries and interfaces are defined such that the solver divides the cells into active, inactive and acceptor groups. The inactive cells are deactivated and the flow is not solved for them while the acceptor cells are used to couple the solution between overset and background region.
   • The overset region moves into the background region, and the solver automatically selects appropriate cells (active or acceptor) for the solution.
   • This meshing technique can simulate a moving body with up to six degrees of freedom as well as it also permits the creation of a zero-gap interface required to simulate the full closure of valves.
   • For cases such as swing valves, the closing valve disc is part of the overset region, which shall rotate according to the momentum exerted by all the forces / torques as per Newton’s second law for momentum.
   • Cells in the background overlapping with valve disk becomes inactive because it is occupying that space. Later on the those cells shall get active again when disc changes its position.

Example profile for motion of the solid-body

(
(movement_linear 3 point)
(time  0   1   2 )
(x     2   3   4 )
(v_y   0   -5   0 )
)

(
(movement_angular 3 point)
(time     0   1   2 )
(omega_x  2   3   4 )
) 

Standard Transient Profiles

(
(profile-name transient N periodic?)
(field_name-1 a1 a2 a3 .... aN)
(field_name-2 b1 b2 b3 .... bN)
.
.
.
.
(field_name-r r1 r2 r3 .... rN)
) 

(
(time_vel_curve transient 3 0)
(time 1  2  3)
(u   10 20 30)
)

The 6DOF solver in ANSYS Fluent uses the forces and moments acting on the object in order to compute the translational and angular motion of the center of gravity of an object. The governing equation for the translational motion of the center of gravity is solved in the inertial coordinate system. Once the angular and the translational accelerations are computed from angualr momentum balance and linear momentum balance respectively, the rates (angular velocity / angular displacement and translation velocity / displacement) are derived by numerical integration. The angular and translational velocities are used in the dynamic mesh calculations to update the rigid body position." For the rigid body motion of a body such as valves and diaphragms, the 6DOF solver is which computes external forces and moments on the valve by computing a numerical integration of the pressure and shear stress over the valve’s surface. It can also add additional forces or moments such as e.g. spring and mass inertia forces. When the forces and moments acting on a rigid body is estimated, it calculates the translational and rotational motion of the center of gravity of the body using equation a = 1/m . ΣF, dω/dt = 1/Inertia . Στ

Roll, Pitch and Yaw

Zones and Mesh Updates in DMM

• Deforming zone: This refers to the mesh which will deform with the simulation. For example, zone or volume in 3D consisting of volume bounded by piston, liner, head and valves. The zone deforms by the valve and piston motion.
• Layering Zone: This zone represents the valve motion and it located above the valve, between valve curtain and valve stem. As the name suggest, mesh layers are added / deleted during the simulation. Whenever the mesh quality in the deforming zone or layering zone deteriorates beyond a pre-defined value, mesh is regenerated by adding or removing layers in the layering zone and remeshing the deforming zone by adding (refining) or removing (coarsening) cells (depending on if the cells are compressed or stretched during piston motion).
• Rigid Motion Zone: This is located between the layering zone and the valve. This zone keeps its geometrical shape throughout the whole engine cycle, but follows the valve motion. The function of the rigid motion zone is to capture the concave valve shape without affecting the layering zone. In this way it is easier to make new meshes at different valve lifts by simply adding or removing layers of cells in the layering zone. It also helps to preserve the mesh quality in the area of the valve exit.
• Static Zone: This zone consists of all volumes that are unaffected by boundary motion and as the name suggests they keep their geometry and shape during the whole engine cycle. This zone includes all remaining volumes in the intake and exhaust ports.
• The mesh deformation is handled in 3 ways:
  • Smoothing: This is method to apply a boundary movement which retains original cells of the mesh and move the nodes to ensure certain mesh quality. This method is appropriate only when the boundaries undergo small deformation. In case of larger deformation, the layer adjacent to boundary may get so skewed that any mesh smoothing operation will not improve the quality to permissible level. The smoothing method includes three different options:
    • Spring-Based Smoothing: use the edges between two mesh nodes as interconnected springs. Equilibrium state is before any boundary motion have occurred. When a displacement happen at a node, it generates a force equal to the displacement between the two nodes. Hooke’s law is used to calculate the force between the nodes. This method is used if the movement of the boundary is primarily only in one direction, resulting in no excessive stretching or compression of the cell zone
    • Laplacian Smoothing: This method adjusts the mesh vertex to the center of the adjacent vertices, does not guarantee improvements on mesh quality and often lead to poor mesh quality.
    • Boundary Layer Smoothing: this is often combined with another mesh motion UDF where the smoothing method is used to deform the boundary layer during the mesh motion.
  • Layering: In this method, the boundary selected as a layering face has certain prescribed motion profile. All the cells adjacent to that face (in the domain marked layring) are allowed to compress and expand. All other cells of the domain perform rigid body motion - compressing or expansion as the same rate as the first layer of cells attached the layering face. This method required the layering zone to consists of either hexahedrons or prism elements. Note that this method is viable only if the moving object (boundary) is not surrounded by cells.
  • Remeshing:The Remeshing method create new cells due to skewness, that is a new mesh is regenerated in the specified zone when the mesh quality deteriorates beyond a pre-defined level. An informative video on YouTube covering 6-DOF (6 Degrees-of-Freedom) solver in ANSYS FLUENT is here.

• Skewness > the value specified by user
• Cell size < specified minimum length scale
• Cell size > specified maximum length scale
• Cell height does not meet length scale specified for adjacent boundary, e.g. a moving wall (valve or piston)

Sample UDF for 6DOF case. DEFINE_SDOF_PROPERTIES (name, properties, dt, time, dtime) specifies custom properties of moving objects for the six degrees of freedom (SDOF) solver which includes mass, moment and products of inertia, external forces and external moments. real *properties - pointer to the array that stores the SDOF properties. The properties of an object which can consist of multiple zones can change in time, if desired. External load forces and moments can either be specified as global coordinates or body coordinates. In addition, you can specify custom transformation matrices using DEFINE_SDOF_PROPERTIES. The boolean properties[SDOF_LOAD_LOCAL] can be used to determine whether the forces and moments are expressed in terms of global coordinates (FALSE) or body coordinates (TRUE). The default value for properties[SDOF_LOAD_LOCAL] is FALSE.

#include "udf.h"
#define PI 3.141592654
DEFINE_TRANSIENT_PROFILE(speed, time) {
  real A = PI/12;   /*Amplitude in [rad]   */
  real f = 10.0;    /*Frequency in [rad/s] */
  real w;           /*Angular displacement */
  w = 2.0*PI*A*cos(f * time);
  return w;
} 

Sign convention of rotor power in FLUENT: 'positive' value implies fluid is supplying energy to the rotor (e.g. turbine), 'negative' value implies rotor is supplying energy to the fluid (e.g. compressor or pump).

DFBI - DMM is known as DFBI in STAR-CCM+ where the DFBI (Dynamic Fluid Body Interaction) module is used to simulate the motion of a rigid body in response to pressure and shear forces the fluid exerts, and to additional user-defined forces such as weight and inertia. STAR-CCM+ calculates the resultant force and moment acting on the body due to all influences, and solves the governing equations of rigid body motion to find the new position of the rigid body. There are multiple type of DFBI features. One such feature is DFBI Superposed Rotation which superimposes an additional fixed body rotation in addition to the DFBI motion. For example, this option can be used to model rotating propellers attached to rotating and/or translating marine boats.

The 6-DOF Solver in STAR-CCM+ computes fluid forces, moments, and gravitational forces on a 6-DOF object where pressure and shear forces are integrated over the surfaces of the 6-DOF bodies. These forces and moments are used to compute the translational motion of the center of mass of the body and the angular motion of the orientation of the body. The 6-DOF body object defines the surface of the floating body for the calculations of the 6-DOF solver.

Mesh Motion in OpenFOAM: The solid body motion is defined by classes derived from common base solidBodyMotionFunction. The motion function returns a septernian which describes the motion of the body. Septernion class used to perform translations and rotations in 3D space. It is composed of a translation vector and rotation quaternion and as such has seven components hence the name 'septernion' from the Latin to be consistent with quaternion rather than 'hepternion' derived from the Greek. Quaternion are the extension of complex number system in 2D geometry to that in 3D geometry (4-dimensional division algebra discovered by Irish mathematician and physicist William Rowan Hamilton).

Immersed Solid Approach: ANSYS CFX uses an immersed solid approach to model to model steady-state or transient simulations involving rigid solid objects that can move through fluid domains such as lobe pump, gear pumps and axial flow fans. The immersed solid is represented as a moving walls (or two counter-rotating moving walls in case of gear pumps) and a source term in the fluid equations that drives the fluid velocity to match the solid velocity. Some of the limitations of this approach are: the immersed solid domain cannot undergo mesh deformation and the surfaces of the immersed solid body are not explicitly resolved by the mesh. In addition, a wall function cannot be applied to the boundary of an immersed solid and hence the accuracy of simulation results may be lower than can be obtained using mesh deformation methods or other techniques that support the use of wall boundaries to directly resolve solid surfaces.

• Like MRF, MPM is a steady-state approach, that is the mixing plane model is useful for predicting steady-state flow in a turbomachine stage when local interaction effects (such as wake and shock wave interaction) are weak.
• After a prescribed number of iterations, the flow data at the mixing plane interface are (by default area-)averaged in the tangential (circumferential) direction at the interface: on both the rotor outlet and the stator inlet faces.
• If rotor-stator interaction effects are important, then a transient sliding mesh calculation is required.
• Note that the mass flux for the rotor portion and stator volume would be different as periodic angles are also different. However, for a complete 360° domain, the mass fluxes for rotor and stator should be equal.

Flow inside Automotive Transmission Gearboxes

The lubrication system used in an automotive gearboxes are of dip-lubricated type or splash-lubricated gearboxes where the gears are partly immersed in oil and oil is transported to the tooth-meshing region while the teeth come out of the oil well. The rotation of gear teeth in highly viscous oil lead to drag known as churning power loss (category of load-independent power losses). Note that in real life, the two faces of gear teeth are in solid-to-solid contact which changes with time.

Thus, in order to apply the CFD method, some arbitrary gaps need to be assumed so that the two interacting teeth never touch each other (sometimes achieved by scaling down the driver and driven gears to 99.5% of their original sizes). It is recommened to start with 2D geometry where VOF multiphase model with dynamic meshing (smoothing, local- or global-remeshing and/or layering) to simulate the fluid flow for a pair of mating gears. Continuous variation of the geometry of the fluid volume during the operation, determined by the mating cycles of the gears, leads to high complexity in the simulation of fluid dynamics inside gearboxes, because it necessarily requires an update of the (distorted) mesh after a few time steps or fraction of degree rotation.

Another source of load-independent losses in gears is squeezing or pocketing generated due to variation of the volume between mating teeth thereby producing axial flows (squeeze effect) of the trapped lubricant. This loss is a lower order of magnitude when compared to churning losses.

Flow inside a centrifugal blower for HVAC applications

The purpose of this demonstration is:

• to gain insight into the operation of a centrifugal blower, effect of casing on pressure recovery
• develop an optimized throat shape, size and location
• execute basic optimization using 2D simulation
• use the 2D mesh to generate the 3D mesh, a novel use of mesh extrusion
• assess the improvement in results using Sliding Mesh Model (SMM) over Multiple Reference Frame (MRF) model
• study the effect of clearances between blade and casing on overall performance - 3D simulation
• check for the issues observed when solution for SMM is initialized with a converged MRF solution vs. full transient start where flow field is uniform.

• Incompressible air at 25 [°C] and 1 [atm] resulting in density of 1.185 [kg/m³.
• Realizable k-ε model with enhanced wall treatment
• Coupled solver with 2nd order discretization schemes for mass, momentum and turbulence

The results with Shear Stress Transport (SST) turbulence model is presented in following plots.

• The calculated mass flow rate per unit [that is 1 m] depth of the blade is 4.72 [kg/s].
• Small level of reverse flow observed near the throat area which has been handles by redesign of this section and extending the outlet. Free-slip wall boundaries can be applied to eliminate the effect of extended domain.
• The location of throat is very close to the optimal design and there is no back flow into the blade cascade from the discharge region.

Cyclone Separators

• Cyclone separators are being used in industries for more than a century. This device falls under the category of what is called "Industrial Duct Collectors".
• Dust separation process utilize different methods ranging from fabrics (such as Air Cleaners in Automotive Intake Systems) to Electrostatic Precipitators in Coal-fired power plants
• Cyclone separators fall under the category of inertial separators which uses combination of the 3 most prevalent mechanical forces namely gravitational, centrifugal and inertial.
  • In a cyclone, a high speed flow of fluid is established by tangential entry into a cylindrical geometry followed by a conical section.
  • The tangential entry of fluid stream result in a rotating and translating (swirling) pattern, beginning at the top and ending at the bottom conical end. Heavier particles in the rotating stream gets thrown out due to centrifugal forces and strike the outside wall, falling along the wall then to the bottom exit of the cyclone where they are collected.
  • In the conical section, due to conservation of angular momentum, as the rotational radius of the stream is reduced, the tangential component of velocity is increased throwing out lighter and smaller particles, leading to separation of smaller particles, even up to size of 5 microns.
• [Reference: Particulate formation and control technology -- Bruce Miller, in Fossil Fuel Emissions Control Technologies, 2015] - The pressure loss and collection efficiency are two key performance parameters of this device. Cyclone separators are used for removing particles having diameters ≥ 10 [μm]. However, conventional cyclones seldom remove particles with an efficiency greater than 90% unless the particle size is ≥ 25 [μm]. High-efficiency cyclones can remove particles up to as los as to 5 [μm].

  

  

• A geometry of cyclone separator in STEP and IGES format can be are as available with the respective links. Dimensions in [mm] are Ht = 1250, D = 400, h = 500, S = 200, a = 200, Dx = 100.
• [Reference: Multiphase Flows in Cyclone Separators - Modeling the classification and drying of solid particles using CFD: Master of Science Thesis by Erik A. R. Stendal] Excerpts - When the swirling gas reaches the cone it will be accelerated due to the decreasing cross-sectional area and a vortex going upwards will form. This vortex will move inside the vortex finder and out through the gas outlet. The purpose of the vortex finder is to prevent contact between the inner vortex and the swirling gas in the barrel to prevent large pressure drops. The diameter of the vortex finder is an important parameter which affects the velocities in the barrel and, as a consequence, the total pressure drop of the cyclone. The particles injected into the cyclone separator will experience a centrifugal force towards the outer wall of the cyclone due to the swirling motion of the gas. The centrifugal force will be opposed by the drag force acting towards the core of the cyclone.
• The cut diameter of the cyclone is defined as the size of the particles collected with 50% collection efficiency. Collection efficiency is defined as the ratio of particles of a given size collected in the cyclone to the number of particles of that size entering into the cyclone at inlet.
• Despite such a long history of application in industry, the design principles so far as mostly based on empirical data. Recently, CFD techniques is being used to optimize the designs.
• 4 geometrical parameter which is tightly linked with the performance of cyclone separators are
  • Vortex finder diameter
  • Inlet width
  • Inlet Height
  • Total Height of the Cyclone
• The cone-tip diameter of the separator does not have noticeable effect on its performance.

Design Ratio: Cyclone Separator: The table below summarizes important dimensions of a cyclone separator. The range given are based on information found in reaseach and thesis documents. However, these dimensions are indicative only and are not intended to provide any engineering solutions.

Gear Pumps

• P_Input = input power to drive gear
• ε_α = transverse contact ratio
• ε₁ = addendum contact ratio
• ε₂ = dedendum contact ratio

• grease lubrication - applied in low speed / low load applications with tangential speed < 5 [m/s]
• splash lubrication (oil bath method) - used with an enclosed system. The rotating gears splash lubricant onto the gear system and bearings. It needs at least 3 [m/s] tangential speed to be effective.
• forced oil circulation lubrication: there are many methods, one of them is the oil drop method where an oil pump is used to suck-up the lubricant and then directly drop it on the contact portion of the gears via a delivery pipe.

Centrifugal Pumps

• z [m] = height of the measurement plane with respect to arbitrarily chosen datum plane
• U [m/s] = mean (area-average) velocity at pressure measurement planes
• ρ [kg/m³] = density of fluid
• p [Pa] = static pressure (the value should be consistent at outlet and inlet: either both in gauge pressure or both in absolute pressure
• g [m/s² = 9.806 - acceleration due to gravity
• N [RPM] = rotation speed of impeller in revolutions per minute
• ω [rad/s] = rotation speed of impeller
• Q [m³/s] = volumetric flow rate through the pump (each eye)

Fan Laws: applicable when the efficiency of scaled model and actual model are assumed nearly equal.

Measurement uncertainty for the individual operating parameters: reference - www.ksb.com/centrifugal-pump-lexicon/total-tolerance/191302

Sample Pump Performance Curve

Sample design data of a centrifugal pump [Reference: ¹Numerical 3D RANS simulation of gas-liquid flow in a centrifugal pump with an Euler-Euler two-phase model and a dispersed phase distribution: T. Mueller, P. Limbach, R. Skoda, ²Investigation of the handling ability of centrifugal pumps under air-water two-phase inflow: model and experimental validation: Qiaorui Si et al.]

Losses in Pumps

Slip Losses: The losses due to "imperfect guidance of the flow by the blades" or slip due to fluid not following the solid rotating wall and fluid. As the fluid traverses through the impeller blades the pressure between each adjacent blade (pressure side of one blade and suction side of other blade) will be different due to the adverse tangential or peripheral pressure gradient. This results in a secondary circulation which forces the fluid exiting the blade to flow backward (from tip towards the hub) with respect to the rotational direction of the impeller. To reduce slip losses, it is suggested to increase the inlet blade angle and reduce the outlet angle (Sixsmith). The slip at shut-off head is a measure of the drag (or 'hold') which the blades have on the fluid (Crewdson). The increase in number of blades reduces slip loss but increases frictional loss At the same time, number and spacing of the blades are strongly related to the diameter of the impeller and and the size of the side channel (height of the blades along axial direction).

Shock Losses: Also known as incidence losses, these losses occur at the entry to the blades and predominantly at off-design operating conditions. It is believed that the difference in angle between the blade and the velocity of fluid entering the blade results in difference of angular momenta between the slower moving fluid in the channel and the faster moving fluid in the impeller. The shock effect is quantified as the ratio of the mean peripheral (tangential) fluid velocity at the blade inlets to the velocity of the blade (at the inner diameter of the blades). As the losses from shock and slip are mainly due to misalignment of the blade and fluid angles, these could be minimised by increasing the impeller diameter and decreasing the hub diameter (Raheel and Engeda).

Cavitation in Pumps

• Mass of 1 m³ of water at 20 [°C] = 998.2 [kg]
• Gas constant for water vapour = 8134/18 =451.9 [J/kg-K]
• Volume occupied by water vapour at 20 [°C] and 1 [bar], V = mRT/p = 1360 [m³].
• However, the density ratio of water liquid and vapour is 998.2 [kg/m³] / 0.0173 [kg/m³] = 57655.

Another term associated with cavitation phenomena is NPSH (Net Positive Suction Head). NPSH is the difference "total pressure at pump inlet" - "vapour pressure of liquid at operating temperature" expressed as head of that fluid. That is:

ρ × g × NPSH = P₀ - P_v. Note that total (or stagnation) pressure and not the static pressure is used in calculation.

NPSH is further differentiated in two types: NPSH_REQUIRED or NPSHr and NPSH_AVAILABLE or NPSHa. The former is supplied by pump manufacturers and this refers to the NPSH that must be available at impeller eye. Hence this value must be independent of the the system in which pump is installed and should solely depend on the design (shape and size) of the pump.

• NPSHr indicates the lowest inlet pressure required by the pump at a given flow to avoid cavitation and is independent of the kind of liquid pumped.
• As the liquid flows from the pump inlet port (suction port) to the eye of the impeller, the velocity increases and hence pressure decreases.
• Additional pressure losses occur due to shock and turbulence as the liquid strikes the impeller. The centrifugal force of the impeller vanes further increases the velocity and decreases the pressure of the liquid.
• Cavitation zone predominantly occurs near the suction side of the leading edge of the impeller blades.
• In CFD simulations, the volume fraction plots on the impeller blades can be used to estimate the extent of cavitation.
• CFD simulation for cavitation is primarily a multiphase flow simulation where the volume fraction on non-condensible air can be specified in the range 1 ~ 5 [ppm]. Note that the amount of O₂ dissolved in water at atmospheric pressure at 15 [°C] is approximately 10.0E−6 [g/g]. For N₂, this value is about 15.0E−6 [g/g]. However, the solubility of air in water is not the sum 25.0E−6 [g/g].
• The solubility of O₂ in water is higher than the solubility of N₂. Air dissolved in water contains approximately 36% O₂ compared to 21% in air.
• Lower the NPSHr, higher the suction capacity.
• On the other hand, NPSHa is dependent on the overall system (supply tank or sump) to inlet to the impeller eye.

In manufacturer's cataloques, characteristic curves (Δp-Q curve) of a pump also contain a curve for NPSHr vs. Q. The NPSHr values indicated are based on measurements carried out with cold water as pumping liquid. NPSHr is also referred to as NPSH3 per API 610 and determines the operating point at which a pump will operate at 3% loss of head due to cavitation. In test set-up, the pump is installed with a starving device (flow and pressure regulator) on its suction line so that the test loop can deliver variable NPSHa. Cavitation begins as small bubbles before any indication of loss of head or capacity can be observed. This is called the point of incipient cavitation and corresponding head is denoted by NPSHi. NPSHr ≈ [2 ~ 20] × NPSHi, is sole responsibility of the pump manufacturer.

Reference - www.iso.org/standard/41202.html: ISO 9906:2012 specifies hydraulic performance tests for customers' acceptance of rotodynamic pumps (centrifugal, mixed flow and axial pumps). It is intended to be used for pump acceptance testing at pump test facilities, such as manufacturers' pump test facilities or laboratories. It can be applied to pumps of any size and to any pumped liquids which behave as clean, cold water. It specifies three levels of acceptance:

• grades 1B, 1E and 1U with tighter tolerance
• grades 2B and 2U with broader tolerance
• grade 3B with even broader tolerance

CFD simulations can be used to determine NPSHr by running a series of simulations for a given system and determining when the performance exhibits 3% head loss. For accurate predictions, the effect of vapor formation and collapse (cavitation) and Non-Condensable Gases (NCG) such as dissolved oxygen should also be considered.

The dimensionless parameter that governs the cavitation characteristics of a centrifugal pump is cavitation number described below.

• ρ_L = density of liquid
• p_i = inlet pressure at pump suction port (Pa)
• r_i = impeller tip radius (m)
• Ω = pump rotation speed (rad/s)
• p_V = vapour pressure or saturation pressure at working temperature of the liquid

Refer to the two operating conditions of a pump at same flow rates, pipe diameters, elevation of discharge tanks and water levels in supply tanks. Should the NPSH_REQUIRED be different in the two situations? Would the NPSH_AVAILABLE be same in both of the scenarios?

Cavitation Test

Effect of Non-condensible Gases

• The primary effect of non-condensible gas is due to the expansion of gas at low pressures. This can lead to significant values of local gas volume fraction causing considerable change in density, velocity, viscosity and pressure distributions.
• The secondary effect of a non-condensible gas is on the early inception of cavitation by increasing the number of nucleation sites for bubble formation.

Mixing Tanks or Stirred Tanks

The time of homogenization (mixing time) is defined as the time from the introduction of the tracer to the time when the tracer concentration at the probe position reaches and remains within a certain range of the final value. If the range is set to ± 5% it is designated as t₉₅.

Flash Mixer: An agitator used to mix a small amount of additive into a continuous stream where the Residence Time is extremely short. Residence Time it average time a process component remains in the mixing environment in a continuous process.

Mixing process needs to handle both the immiscible (e.g. water - silicone oil, water - benzene) and miscible liquids (e.g. water - alchohol, water - caustic solution), liquid with large difference in viscosities (e.g. water-0.001 Pa.s and molasses-2 Pa.s) and fluids with large difference in densities. In case of miscible liquids, only transient simulation can be performed and mass or volume fraction of one the phases needs to be monitored at different locations of the tank. If left for a long enough period of time, miscible liquids will any way dissolve in one another and form a homogeneous solution. A sliding mesh model is recommended to transport the velocity of the impeller to the bulk liquid in the stationary domain.

Reference: Hayward Gordon - MASTERING MIXING FUNDAMENTALS: A technical guide from the experts in the industry - Blending:blending of miscible fluids. If left for a long enough period of time, miscible liquids will dissolve in one another and form a homogeneous solution. The majority of liquid/liquid blending applications fall into this category. Applications involving immiscible (insoluble) fluids are classified as dispersion applications which involves very different mixer sizing methods.

There are other types of mixtures known as static mixtures where the mixing elements are inserted in the pipelines and the liquids mix as they flow through them. Kenics static mixer from Chemineer is one such device.

Reference: Mechanical Design of Mixing Equipment, D. S. DICKEY, MixTech, Inc. J. B. FASANO Chemineer, Inc - Dry-solids mixers are normally applied to flowable powdered materials. The action of the mixers can be categorized as summarized below.

1. Tumble mixers
2. Convective mixers which use a ribbon, paddles, or blades to move material
3. High-shear mixers which create a crushing action like a mortar and pestle
4. Fluidized mixers as in fluidized beds and
5. Hopper mixers which use discharge and recirculating flow to cause mixing

The Reynolds number for agitators or mixers are calculated based on blade tip speed. However, the adopted formula has been simplified a bit by droping π and the formula is Re = ρ[kg/m³]×N[rev/s] ×D[m²]/μ [Pa.s]. Flow is assumed turbulent when Re > 10,000 ((McCabe, Unit Operations of Chemical Engineering, 1993). Torque per Equivalent Volume [Torque on impeller / Working Volume of Tank] is an extremely useful ratio which is used as the basis for mixer sizing and describes the level of mixing for any application.

Vapour Pressure of Water

T [°C]	PSAT [Pa]
5	872.60
10	1228.1
15	1705.6
20	2338.8
25	3169.0
30	4245.5
35	5626.7
40	7381.4
45	9589.8
50	12344
55	15752
60	19932
65	25022
70	31176
75	38563
80	47373
85	57815
90	70117
95	84529
100	101320

Surface Water: water on the earth’s surface, including rivers, ponds, water reservoirs (dams), springs, creeks and wetlands/swamps. Most surface water comes from rainfall (precipitation) runoff from the surrounding land area (catchment). Surface water also includes the solid forms of water - snow and ice.

Non-Surface Water: water under the earth’s surface, underground water, aquifer system. It also includes soil water.

Propellers

In other words: The 'face' of a blade is the high-pressure side or pressure face of the blade. This is the side that faces aft (downstream) and pushes the water when the vessel is in forward motion. The 'back' of the blade is the low pressure side or the suction face of the blade. This is the side that faces upstream (incoming water) or towards the front of the vessel.

Propellers are axial flow type. In most of the rotating devices and especially in axial flow machines, blades have two edges: leading edge and trailing (lagging) edge. The blades rotates in the direction formed by rotating trailing edge towards the leading edge. The shape of the blades have speacial twist and the parameters defining the twists are known as rake, skew and pitch.

Propellers operate in non-uniform flow field (wake regions) created by the boat or ship body. However, the theoretical studies of the performance of an impeller is carried out in calm water known as open-water charateristics. Advance ratio, thrust coeffficient, torque coefficient and power coefficient are the key parameters required to define a propeller.

• u_∞ = ship or flight (free-stream) velocity
• u_s = ship velocity
• ρ = density of the liquid
• n = revolution per second (RPS)
• D = propeller tip diamter
• T = Thrust
• Q = Torque
• P = Power
• J = Advance ratio or advance number
• K_T = thrust coefficient
• K_Q = torque coefficient
• η_PROP = propeller efficiency
• η_t = propulsive efficiency
• C_POW = power coefficient
• R_t = total ship resistance
• w = wake fraction
• t = thrust deduction factor

Input

• Thrust, T = 2,000 [N]
• Boat speed, u_∞ = 5[m/s] = 18 [km/hour]
• Fluid Density, ρ = 1000 [kg/m₃]
• Propelle blade tip diameter, D = 0.75 [m]
• RPS, n = 20 [rev/s]
• Number of blades, n_B = 4 [ - ]

Cross-flow Fans

Reciprocating Compressors

Types of Hydraulic pumps

1. Hydrodynamic: a rotating element within a casing imparts radial velocity to the fluid which exits in a discharge nozzle
2. Hydrostatic: also known as positive displacement pumps, displace a fixed volume of fluid per cycle (rotation of shaft)
   • Rotary: a rotating vane, screw or gear collects fluid in its void on the suction side of the pump and forces it to the discharge side
   • Reciprocating: fluid is pulled into a void on the suction side of the pump and then pushed during the compression stroke, through the discharge port of the pump.

Swash-plate axial piston pump

Reference: Axial Piston Pump - Leakage Modelling and Measurement, PhD Thesis by Jonathan Mark Haynes

Reference: Use of CFD Technology in Hydraulics System Design for off-Highway Equipment and Applications by Shivayogi S. Salutagi, Milind S. Kulkarni and Aniruddha Kulkarni - International Journal of Materials, Mechanics and Manufacturing, Vol. 4, No. 1, February 2016

Leakage Paths:

1. The gap between the rotating barrel and fixed port plate introduces a cross-port leakage and a leakage across the barrel edges to the internal casing, and ultimately back to tank.
2. The slipper face consists of a lubrication orifice, through which pressurised fluid can flow, ultimately returning back to tank.

Wind Turbines

Fluid Pressure in Rotating Reference Frame

For a container rotating cylindrical about its axis: the shape of the free surface is a parabola and fluid inside the rotating cylinder forms a paraboloid of revolution, whose volume is one-half of the volume of the "circumscribing cylinder". To calculate angular velocity at which the liquid at the center reaches the bottom of the cylinder just as the liquid at the curved wall reaches the top of the cylinder: ω_spill = (2gH)^0.5/R.

Ball in rotating tube

The centroid of a hemi-sphere is at 3R/8 from the base. Using this value, pressure on the left half of the ball = 1/2.ρ_F ω²(r-3R/8)². The pressure on the right half of the ball = 1/2.ρ_F ω²(r+3R/8)². Net force acting on the ball towards the axis of rotation = 1/2.ρ_F ω²(4.r.3R/8)² × π*R² where projected area = π*R².

Net force due to fluid pressure towards axis of rotation, F_F = 1/2.πR².ρ_F ω²(r.3R/2) = 3/4πR³.ρ_F ω²r

Centrifugal force due to own mass of the ball, F_B = 4/3.πR³.ρ_B ω²r

Root or Lobe Pump and Blowers

References

1. ANSYS FLUENT User Manual
2. Study of mesh deformation features of an open source CFD package and application to a gear pump simulation: Alejandro Roger Ull, ETSEIAT 2012