Mixer Module Users Guide PDF
Mixer Module Users Guide PDF
User’s Guide
Mixer Module User’s Guide
© 1998–2018 COMSOL
Protected by patents listed on www.comsol.com/patents, and U.S. Patents 7,519,518; 7,596,474;
7,623,991; 8,457,932; 8,954,302; 9,098,106; 9,146,652; 9,323,503; 9,372,673; and 9,454,625. Patents
pending.
This Documentation and the Programs described herein are furnished under the COMSOL Software License
Agreement (www.comsol.com/comsol-license-agreement) and may be used or copied only under the terms
of the license agreement.
COMSOL, the COMSOL logo, COMSOL Multiphysics, COMSOL Desktop, COMSOL Server, and
LiveLink are either registered trademarks or trademarks of COMSOL AB. All other trademarks are the
property of their respective owners, and COMSOL AB and its subsidiaries and products are not affiliated
with, endorsed by, sponsored by, or supported by those trademark owners. For a list of such trademark
owners, see www.comsol.com/trademarks.
Version: COMSOL 5.4
Contact Information
Visit the Contact COMSOL page at www.comsol.com/contact to submit general
inquiries, contact Technical Support, or search for an address and phone number. You can
also visit the Worldwide Sales Offices page at www.comsol.com/contact/offices for
address and contact information.
If you need to contact Support, an online request form is located at the COMSOL Access
page at www.comsol.com/support/case. Other useful links include:
Chapter 1: Introduction
Tank Parts 41
Impeller Parts 42
Impeller Parameters. . . . . . . . . . . . . . . . . . . . . . 58
CONTENTS |3
Chapter 4: Rotating Machinery, Fluid Flow
4 | CONTENTS
The Rotating Machinery, Mixture Model, v2-f Interface . . . . . . . 102
Domain, Boundary, Pair, and Point Nodes for the Rotating
Machinery, Mixture Model Interfaces . . . . . . . . . . . . . 103
Free Surface. . . . . . . . . . . . . . . . . . . . . . . . 104
Contact Angle . . . . . . . . . . . . . . . . . . . . . . . 106
Rotating Shaft . . . . . . . . . . . . . . . . . . . . . . . 106
CONTENTS |5
Reacting Flow Interfaces
6 | CONTENTS
2
Introduction
This guide describes the Mixer Module, an optional add-on package for
COMSOL Multiphysics® designed to assist you in setting up and solving transport
problems in mixers and stirred vessels. The module is an add-on to the CFD
Module and provides additional support for modeling fluid flow in rotating
machinery.
This chapter introduces you to the capabilities of the module. A summary of the
physics interfaces and information about where to find documentation and model
examples is also included. This is followed by a brief overview with links to each
chapter in the guide. The last section in this introduction contains a tutorial
example from the Mixer Module application library.
7
About the Mixer Module
The development of new products and processing lines continuously places high
demands on mixer design. The Mixer Module provides engineers and scientists with
the necessary simulation tools for predicting and understanding the interactions
between fluid flow and chemical processes in mixers and stirred vessels. The detailed
information gained from computational fluid dynamics (CFD) simulations
substantially reduces the number of empirical studies needed in new product and
process designs, and has eminent impact on both costs and responses to environmental
considerations.
Figure 2-1: Streamlines and velocity magnitude in a cross section of a dished bottom
baffled mixer with a four-bladed impeller. The color scale on the baffles represents the
pressure.
The Mixer Module complements the CFD Module with additional functionality for
the Rotating Machinery, Fluid Flow branch. The added functionality includes
extended capability for modeling turbulence in the Rotating Machinery interfaces.
In order to facilitate fast and efficient setup of mixer geometries, the Mixer Module
Part Library includes predefined geometry components typical of mixer equipment.
8 | CHAPTER 2: INTRODUCTION
The part library includes impeller parts for axial impellers, radial impellers, and
impellers designed for highly viscous fluids. In addition to impellers, three types of
different tank geometries and a cylindrical impeller shaft geometry are available in the
part library. All mixer parts are modularized through a number of input parameters
corresponding to important geometrical properties of each part. These can be adjusted
in order to fit the mixer system under investigation.
A high rotation rate or a strong acceleration of the rotation rate may induce a
substantial deformation of the free surface in an open vessel. This topology change in
turn influences the flow pattern inside the vessel. The Mixer Module includes
free-surface features to capture the displacement of the liquid-air interface induced by
the bulk motion in the domain, by the walls, and by the rotating shaft.
The physics interfaces define a fluid-flow problem using physical quantities such as
pressure, flow rate, temperature, and species composition, as well as physical
properties, such as viscosity, thermal diffusivity, and density. The different physics
interfaces cover a wide range of laminar and turbulent mixer flows. The conservation
laws formulated by the physics interfaces are expressed in terms of partial differential
equations along with corresponding initial and boundary conditions. The equations
are solved by the module using stabilized finite element formulations for fluid flow in
combination with damped Newton methods and, for time-dependent problems, in
combination with various time-dependent solver algorithms. The Mixer Module’s
general capabilities include frozen-rotor and time-dependent flows in two- and
three-dimensional spaces. For a so-called frozen-rotor flow, the topology relative to
the rotating reference frame is fixed (“frozen”). When the flow field is, or can be
approximated to be, of this type the computational time (CPU time) can be
substantially reduced using the Frozen Rotor (see the CFD Module User’s Guide)
study type.
The workflow in the Mixer Module is quite straightforward. Set up a simulation using
one of the Rotating Machinery interfaces, described by the following steps: define the
geometry, select the fluid to be modeled, select the type of flow, define boundary and
initial conditions, define the finite element mesh, select a solver, compute the solution,
and visualize the results. All these steps are performed from the COMSOL Desktop.
The mesh and solver steps are usually carried out automatically using default settings,
which are tuned specifically for each Fluid Flow interface.
The models available in the Mixer Module application library describe the physics
interfaces and their features through examples for different types of mixer flows. Here
you find examples of industrial equipment and devices, tutorial models for practice,
To help you get started, this introduction contains a list of the physics interfaces and
an example, Tutorial Example — Nonisothermal Mixer, to introduce you to the
workflow.
10 | CHAPTER 2: INTRODUCTION
The physics interfaces in the Mixer Module are able to perform all steps in mixer
analyses, from the initial idea and qualitative simulations to the final optimization of
the product or process.
Figure 2-2: The Rotating Machinery interfaces for laminar and turbulent flow.
The Rotating Machinery, Turbulent Flow interfaces ( ) are used to model flow at high
Reynolds numbers. These physics interfaces solve the Reynolds-averaged
Navier-Stokes (RANS) equations for the averaged velocity field and averaged pressure.
The different physics interfaces in this branch have different models for the turbulent
viscosity. There are several turbulence models available — two algebraic turbulence
models, the Algebraic yPlus and L-VEL models, and seven transport-equation models,
including a standard k-ε model, the Realizable k-ε model, a k-ω model, a SST model,
a low Reynolds number k-ε model, the Spalart-Allmaras model, and the v2-f model.
Similarly to the Rotating Machinery, Laminar Flow interface, compressibility is set to
incompressible by default.
The Algebraic yPlus and L-VEL turbulence models are so-called enhanced viscosity
models. A turbulent viscosity is computed from the local distance to the nearest wall.
For this reason, the algebraic turbulence models are best suited for internal flows, such
as in electronic cooling applications. Algebraic turbulence models are computationally
economical, and more robust but, in general, less accurate than transport-equation
models. Among the transport-equation turbulence models, the standard k-ε model is
the most widely used turbulence model because it is often a good compromise between
accuracy and computational cost (memory and CPU time). The Realizable k-ε model
is similar to the standard k-e model but has built-in realizability constraints, resulting
in improved performance for certain flows, such as turbulent jets. The k-ω model is an
alternative to the standard k-ε model and often gives more accurate results, especially
in recirculation regions and close to solid walls. However, the k-ω model is also less
robust than the standard k-ε model. The SST model combines the robustness of the
k-ε model with the accuracy of the k-ω model, making it applicable to a wide variety
of turbulent flows. The low Reynolds number k-ε model is more accurate than the
standard k-ε model, especially close to walls, but requires higher resolution in the
near-wall region. The Spalart-Allmaras model is specifically designed for aerodynamic
applications, such as flow around wing profiles, but is also widely used in other
applications due to its high robustness and decent accuracy. In the v2-f model, the
turbulent viscosity is based on the wall-normal velocity fluctuations, whereby wall
12 | CHAPTER 2: INTRODUCTION
blockage effects and low Reynolds number effects are captured separately. The v2-f
model also includes non-local effects of the fluctuating pressure on the turbulent fields.
The Rotating Machinery, Nonisothermal Flow, Turbulent Flow interfaces ( ) solve the
Reynolds-averaged Navier-Stokes (RANS) equations together with the equations for
heat transfer in fluids and in solids. There is support for all the fluid-flow turbulence
models - the Algebraic yPlus model, the L-VEL model, the standard k-ε model, the
Realizable k-ε model, the k-ω model, the SST model, a low Reynolds number k-ε
model, the Spalart-Allmaras model, and the v2-f model.
Single-Phase Flow
14 | CHAPTER 2: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
16 | CHAPTER 2: INTRODUCTION
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
18 | CHAPTER 2: INTRODUCTION
To open the Help window:
• In the Model Builder or Physics Builder click a node or window and then
press F1.
• On the main toolbar, click the Help ( ) button.
• From the main menu, select Help>Help.
• Press Ctrl+F1.
• From the File menu select Help>Documentation ( ).
• Press Ctrl+F1.
• On the main toolbar, click the Documentation ( ) button.
• From the main menu, select Help>Documentation.
Once the Application Libraries window is opened, you can search by name or browse
under a module folder name. Click to view a summary of the model or application and
its properties, including options to open it or its associated PDF document.
To include the latest versions of model examples, from the Help menu
select ( ) Update COMSOL Application Library.
20 | CHAPTER 2: INTRODUCTION
COMSOL Video Gallery www.comsol.com/video
Support Knowledge Base www.comsol.com/support/knowledgebase
For general information about setting up and solving CFD applications, see the CFD
Module User’s Guide.
The last section in this chapter features a model you can access from The Application
Libraries Window. The Tutorial Example — Nonisothermal Mixer solves a mixer-flow
problem using the Laminar Flow interface in the Rotating Machinery, Nonisothermal
Flow branch.
The following chapter documents the Mixer Module Part Library. It includes an
overview of the all included impellers and tanks.
The third chapter introduces you to the Theory for the Free Surface Features. It
includes descriptions of the mesh deformation within the free surface domain and the
conditions that need to be satisfied at a free surface, at a three-phase boundary, and on
a rotating shaft within the free surface domain.
The remaining chapters describe the physics interfaces and features under the Rotating
Machinery, Fluid Flow branch, and the additional branches exclusive to the Mixer
Module: Rotating Machinery, Multiphase Flow Interfaces, Rotating Machinery,
Nonisothermal Flow Interfaces, and Rotating Machinery, Reacting Flow Interfaces.
22 | CHAPTER 2: INTRODUCTION
Tutorial Example — Nonisothermal
Mixer
This model demonstrates how to obtain the temperature distribution in a simplified
tabletop lab mixer using the Rotating Machinery, Nonisothermal Flow branch in the
Mixer Module. The key instructive element is a demonstration of the Frozen Rotor
method, which substantially reduces the computational time for a mixing study.
Immersion heater
Rotational
direction
Figure 2-3: Tabletop lab mixer geometry, showing the rotational direction of the impeller
and the location of the immersion heater.
The tank is made of steel and is subjected to cooling by natural convection on the
outside of the mixer vessel. The surrounding conditions corresponds to a temperature
of 20°C and a pressure of 1 atmosphere. The total height of the reactor is 0.2 m. These
conditions are needed as input for the natural convection correlations, which are used
to calculate the heat transfer coefficient from the tank wall to the surroundings.
MODEL SETUP
The Reynolds number for a mixer is commonly calculated as
ND a2
Re = ------------- (2-1)
ν
where N is the number of rotations per second, Da the impeller diameter, and ν the
kinematic viscosity. A high Reynolds number means that the flow has a tendency to
become turbulent. Evaluating Equation 2-1 using ν at 60 °C gives Re = 6944. This
Reynolds number indicates that the flow is at least partly turbulent. For simplicity in
this instructional model, the flow is assumed to be two-dimensional and no turbulence
model is used. Possible extensions of the model includes to resolve it using the full
three dimensional geometry, and also to apply a turbulence model to investigate the
effect of turbulent structures occurring in the flow.
A computationally more efficient method is to first simulate the flow using the
frozen-rotor approach. The frozen-rotor approach is a modeling concept that treats
the rotor as fixed, or frozen in space. The flow in the rotating domain is assumed to be
stationary in terms of a rotating coordinate system. The effect of the rotation is then
accounted for by the Coriolis and centrifugal forces. The flow in the non-rotating parts
is also assumed to be stationary, but in a non-rotating coordinate system. (See Frozen
Rotor in the CFD Module User’s Guide for more information.) The result of a
frozen-rotor simulation is an approximation to the flow induced by the impeller. The
result depends on the angular position of the impeller and cannot represent transient
effects. However, it is still a very good starting point to reach operating conditions.
24 | CHAPTER 2: INTRODUCTION
Here the frozen-rotor result is used as input to a time-dependent simulation and the
progress towards the operating conditions is monitored using probe plots.
Figure 2-4: Velocity field obtained from the frozen rotor simulation.
Figure 2-5 shows the temperature distribution obtained from the frozen-rotor
simulation. Streamlines are also included to visualize the flow field. The temperature
is relatively homogeneous throughout the mixer. There are some cold spots in
connection to the recirculation zones adjacent to the outer wall. This is expected
because the fluid there has a longer residence time close to the solid wall, and therefore
has less contact with the heated fluid closer to the center of the mixer.
The progress of a solution can be monitored using probes (see Probes in the COMSOL
Multiphysics Reference Manual). The velocity magnitude and temperature are
probed at (x,y) = (−0.05,0.065). The location is indicated in Figure 2-5, just outside
the recirculation zone along the top wall.
The probe plots produced during the time-dependent simulation are shown in
Figure 2-6. The velocity probe plot shows that the flow pattern, after an initial
transient, oscillates around the frozen-rotor result with an amplitude of about 10%.
The deviations in temperature are much smaller.
The velocity probe plot clearly shows the frequency that corresponds to the passing of
the blades (0.75 Hz) and an additional frequency (about 0.25 Hz) related to the
intermittent boundary-layer separation. Evaluating the frequency from the Strouhal
number for the immersion heater gives a value of about 0.25 Hz, suggesting that the
boundary-layer separation may be induced by the Kármán vortex street behind the
heater. The temperature value does not display the same frequencies and instead
changes more smoothly.
26 | CHAPTER 2: INTRODUCTION
Figure 2-6: Probe plots of velocity and temperature in the mixer from the time-dependent
simulation.
Looking at Figure 2-5, it can be seen that the recirculation zone along the top wall
contains a single, large vortex. As the simulation progresses (t = 20 s to t = 40 s), the
size and strength and position of the vortices along the top wall varies as a result of the
interaction between the disturbance, produced by the immersion heater, and the outer
wall.
Figure 2-7: Evolution of the temperature field from frozen rotor solution toward operating
conditions.
The results obtained in this model are typical for rotating-machinery models: The
frozen-rotor approach can, with minimal computational effort, deliver a decent
approximation of the flow and temperature fields. But transient effects can only be
28 | CHAPTER 2: INTRODUCTION
captured with a time-dependent simulation, and these effects can change local
temperature and velocity values significantly.
MODEL WIZARD
The first step to build a model is to open COMSOL Multiphysics, then select the
physics interface and specify the type of analysis you want to do — in this case, a
frozen-rotor analysis using a Rotating Machinery, Nonisothermal Laminar Flow
interface.
1 Open COMSOL Multiphysics. On the New page, click Model Wizard . Then click
the 2D button .
2 In the Select Physics tree, under Fluid Flow>Nonisothermal Flow>Rotating Machinery,
Nonisothermal Flow click Laminar Flow .
3 Click the Add button.
4 Click the Study button .
5 In the tree, select Preset Studies for Selected Physics Interfaces>Frozen Rotor .
6 Click the Done button .
MATERIALS
On the Home toolbar, click Add Material .
There are many ways to select geometric entities. When you know the domain to add,
such as in this exercise, you can click the Paste Selection button and enter the
information in the Selection text field. In this example, enter 1 in the Paste Selection
window. For more information about selecting geometric entities in the Graphics
window, see the COMSOL Multiphysics Reference Manual.
30 | CHAPTER 2: INTRODUCTION
DEFINITIONS
Rotating Domain 1
1 In the Model Builder window, under Component 1 (comp1)>Definitions>Moving Mesh
click Rotating Domain 1.
1 In the Settings window for Rotating Domain, locate the Domain Selection section.
2 In the list, choose 1 and 2.
3 Click Remove from Selection . Only domain 3, the impeller domain, is selected.
4 In the Settings window for Rotating Domain locate the Rotation section.
5 In the Revolutions per time text field, type -20[rpm].
Flow Continuity 1
1 In the Physics toolbar, in the Boundary section, click Pairs and choose
Flow Continuity .
2 In the Settings window for Flow Continuity, locate the Pair Selection section.
3 In the Pairs list, select Identity Boundary Pair 1 (ap1).
Interior Wall 1
1 In the Physics toolbar, click Boundaries and choose Interior Wall .
Solid 1
1 In the Physics toolbar, click Domains and choose Solid .
2 Select Domain 1, the mixer vessel wall, only.
Continuity 1
1 In the Physics toolbar, in the Boundary section, click Pairs and choose Continuity .
2 In the Settings window for Continuity, locate the Pair Selection section.
3 In the Pairs list, select Identity Boundary Pair 1 (ap1).
Temperature 1
1 In the Physics toolbar, click Boundaries and select Temperature .
2 Select Boundaries 13–16, which make up the boundary of the immersion heater.
3 In the Settings window for Temperature, locate the Temperature section.
4 In the T0 text field, type 60[degC].
Heat Flux 1
1 In the Physics toolbar, click Boundaries and choose Heat Flux .
2 Select Boundaries 1, 2, 7, and 12, which make up the circular outer boundary of the
mixer vessel.
3 In the Settings window for Heat Flux, locate the Heat Flux section.
4 Click the Convective heat flux button.
5 From the Heat transfer coefficient list, choose External natural convection.
6 In the L text field, type 0.2[m].
32 | CHAPTER 2: INTRODUCTION
MESH 1
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1 .
2 In the Settings window for Mesh from the Element size list, select Fine.
The default mesh generates isotropic mesh elements across the thin vessel wall. Edit
the Physics-controlled mesh settings to increase the element size in the circumferential
direction of the solid domain.
Free Triangular 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click
Free Triangular 1 .
Now, add a separate Size node with reduced resolution in narrow regions.
Size 1
1 Right-click Component 1 (comp1)>Mesh 1>Free Triangular 1 and choose Size .
2 In the Settings window for Size 1, locate the Geometric Entity Selection section.
3 From the Geometric entity level list, choose Domain.
4 Select Domain 1, the vessel wall, only.
5 Locate the Element Size section. Click the Custom button.
6 Locate the Element Size Parameters section. Select the Resolution of narrow regions
check box.
7 In the associated text field, type 0.1.
Build the Free Triangular 1 mesh node in order to inspect the mesh at this point.
8 In the Settings window for Size 1, click the Build Selected button
Now add the mesh boundary layers as well. Use boundary layer mesh also in the solid
domain to increase the resolution there.
Boundary Layers 1
1 In the Model Builder window, under Component 1 (comp1)>Mesh 1 click
Boundary Layers 1 .
2 In the Settings window for Boundary Layers, locate the Domain Selection section.
3 From the Selection list, choose All domains.
Now that the physics is defined and the geometry is meshed, solve the model using the
frozen rotor approach.
STUDY 1
In the Home toolbar, click Compute .
RESULTS
Re-create Figure 2-4 using the following steps.
Velocity (spf)
1 In the Model Builder window, under Results , click the Velocity (spf) node.
1 Right-click Velocity (spf) and choose Streamline .
2 In the Settings window for Streamline, locate the Streamline Positioning section.
Temperature (ht)
1 n the Model Builder window, expand the Results>Temperature (ht) node, then click
Surface .
2 In the Settings window for Surface, locate the Coloring and Style section.
34 | CHAPTER 2: INTRODUCTION
3 From the Color table list, choose Wave.
4 Go to the Temperature (ht) toolbar and select Streamline .
5 In the Settings window for Streamline, locate the Streamline Positioning section.
6 From the Positioning list, choose Uniform density.
7 In the Separating distance text field, type 0.02.
8 Locate the Coloring and Style section. Find the Arrow style subsection. From the
Color list, choose Gray.
9 In the Temperature (ht) toolbar, click Plot .
Add a probe to follow the development of the flow during the time-dependent
simulation.
DEFINITIONS
Add a Time Dependent study in order to perform a transient simulation, using the
previous solution as initial condition.
ADD STUDY
1 In the Home toolbar, click Add Study to open the Add Study window.
2 Go to the Add Study window.
STUDY 2
36 | CHAPTER 2: INTRODUCTION
4 Click to expand the Values of Dependent Variables section. Find the
Initial values of variables solved for subsection. From the Settings list, choose
User controlled.
5 From the Method list, choose Solution.
6 From the Study list, choose Study 1, Frozen Rotor.
For the transient simulation, add a restriction on the time step. This will make sure that
the impeller rotation in each time step is bounded, and that a high accuracy is
maintained throughout the simulation. First generate the solver sequence.
Solution 2
1 In the Study toolbar, click Show Default Solver .
2 In the Model Builder window, expand the Solution 2 (sol2) node, then click
Time-Dependent Solver 1 .
Apply a maximum time step of 0.05 s. This is equivalent to an impeller rotation of
6º.
3 In the Settings window for Time-Dependent Solver, click to expand the Time Stepping
section.
4 From the Maximum step constraint list, choose Constant.
5 In the Maximum step text field, type 0.05.
6 In the Study toolbar, click Compute .
Two probe plots are automatically displayed when you start the calculation.
RESULTS
The following steps create an animation that contains the plots in Figure 2-7.
Temperature (ht) 1
1 In the Model Builder window, under Results , click Temperature (ht) 1 .
2 In the Settings window for 2D Plot Group, locate the Plot Settings section.
Plot the data set edges in the spatial frame to make them follow the rotation.
3 From the Frame list, choose Spatial (x, y, z).
Surface
1 In the Model Builder window, expand the Temperature (ht) 1 node, then click
Surface .
2 In the Settings window for Surface, locate the Coloring and Style section.
3 From the Color table list, choose Wave.
Animation 1
1 Go to the Results toolbar, click Animation and choose File .
2 In the Settings window for Animation, locate the Target section.
3 From the Target list, choose Player.
4 Locate the Scene section. From the Subject list, choose Temperature (ht) 1.
5 Locate the Frames section. From the Frame selection list, choose All.
6 Locate the Playing section. In the Display each frame for text field, type 0.25.
7 In the Graphics window, click Play .
38 | CHAPTER 2: INTRODUCTION
3
In this chapter:
39
Part Library Contents
The dedicated part library included in the Mixer Module contains geometrical parts
typical of mixer and stirred vessel equipment. The included parts consists of impellers,
tanks, and a straight cylindrical shaft. When accessing the library, the parts are sorted
under the categories seen in the figure below
The dimensions and configuration of the different parts can be controlled to a great
detail using the input parameters included in each part. For example, adding a pitched
impeller part instance, you can cut different parts of the impeller, round its edges,
define the pitch angle, and define asymmetrical blades, which are wider above the
rotation plane than below. Also several impellers can be added, and also rotate them
slightly around the rotational axis in order to evenly distribute the blades in the
circumferential direction.
Figure 3-2: Tank parts available in the Mixer Module Part Library.
All tank parts can be defined to include an optional number of baffles in order to
improve mixing. Baffles are specially required when radial impellers are used, in order
to avoid a near solid body rotation of the fluid.
By default, the tank’s cylindrical part is positioned with its base at z = 0. This implies
that the dished and cone bottom sections are positioned just below the xy-plane. If you
define a flat bottom tank with a fillet radius on the lower cylinder edge, then also the
fillet section created is placed just below the xy-plane. If you specify a zero fillet radius
for a flat bottom tank, then the base of the tank is positioned in the xy-plane, that is,
at z = 0.
The impeller shaft’s lowest part is also by default positioned at z = 0. You can change
the position, length, and diameter of the impeller shaft. However, the impeller shaft is
always placed in the middle of the tank perpendicular to the xy-plane.
TANK PARTS | 41
Impeller Parts
As seen in Figure 3-1, the impeller parts are sorted into three categories, depending
on the principal direction of the flow generated by the impeller rotation (axial or radial
flow), or if the impeller type is used to mix highly viscous fluids. In total, the Mixer
Module Part Library includes eleven impellers: five axial impellers, four radial
impellers, and two impellers for highly viscous fluids.
Radial impellers force the flow in the radial direction and only give an axial component
once the flow hits the walls of the vessel. In order to achieve good mixing, these
impellers rely on high shear rates and the presence of baffles, that disrupts tangential
flows that would lead to poor mixing.
Figure 3-3: Examples of an axial impeller (pitched impeller) and a radial impeller
(Rushton turbine).
IMPELLER PARTS | 43
be varied between 0 and 90 degrees. The outer vertices of the outer edges can also be
rounded using fillets.
Pitch angle
Figure 3-5: Asymmetrical blade. The upper part of the blade is wider than the lower part.
The figure below contains the notations used for the pitched impeller in the part
library. The full list of parameters together with their description is found in Table 3-1
The vertices of the outer vertical edge can be rounded using fillets. The figure below
shows an impeller with asymmetrical blades and rounded vertices.
Cut upper and Fillets create
inner edges a rounded bend
Figure 3-8: Asymmetrical impeller blades with rounded outer (away from the hub)
vertices.
IMPELLER PARTS | 45
The figure below contains the notations used for the pitched impeller with bent blades
in the part library. The full list of parameters together with their description is found
in Table 3-1
Figure 3-9: Notations used for the pitched impeller with bent blades.
As for the regular pitched impellers, the pitched impeller with constant pitch can be
designed with isosceles trapezoid-shaped blade projections (the blades are not flat).
Figure 3-10: Pitched impellers with constant pitch with different blade shapes.
The inner pitch angle can be varied in order to change the axial and radial components
of the flow induced by the impeller. The pitch angle is then automatically calculated to
give a constant pitch.
IMPELLER PARTS | 47
The outer vertices of the blades can be rounded, which gives additional freedom in the
design of the impeller blades. The figure below shows a propeller created by using
fillets and cuts and applying those on a pitched impeller with constant pitch.
Fillet on upper
and outer vertex
Fillet on lower
and outer vertex
Figure 3-11: Impeller of propeller type created using isosceles trapezoid-shaped blade
projections and large fillet radii for the outer vertices of the blades.
The figure below contains the notations used for the pitched impeller with constant
pitch in the part library. The full list of parameters together with their description is
found in Table 3-1
Figure 3-12: Notations used for the pitched impeller with constant pitch.
HYDROFOIL IMPELLER
Impellers with blades that are curved along the vertical edges are denoted hydrofoil
impellers. These impellers are similar to the pitched blade impellers, but give an
additional pressure difference across the impeller blade with an additionally lowered
pressure on the convex side of the blade.
As in the case for the pitched impeller type, the pitch angle can be varied to control the
relation between the radial and axial flow created as the impeller rotates.
Rectangular curved Asymmetrical blade with wider
blade with attachment at the hub upper part of the blade
Flip
vertical
Figure 3-13: The curved blades can be shaped to create a large variety of impeller designs.
IMPELLER PARTS | 49
The figure below contains the notations used for the hydrofoil impeller in the part
library. The full list of parameters together with their description is found in Table 3-1
As in the case for the regular hydrofoil impellers, the hydrofoil impeller with constant
pitch gives an additional pressure difference across the impeller blade with an
additionally lowered pressure on the convex side of the blade.
Cut wider
outer edge
The figure below contains the notations used for the hydrofoil with constant pitch
impeller in the part library. The full list of parameters together with their description
is found in Table 3-1
Figure 3-16: Notations used for the hydrofoil with constant pitch.
IMPELLER PARTS | 51
the outer radius of the impeller. This impeller can be used to create the so-called
Intermig® impeller.
The impeller can be designed with different pitch angles on the arms of the impeller.
Also the angles of the c-shaped part can be varied. The profile can be changed by
changing the vertical angle of the c-shaped part but keeping the upper and lower
blades parallel. In addition, also the angle of the back of the “c” can be varied, keeping
the lower and upper blades parallel but displacing them in the x-y plane.
C-shaped Adapt for clockwise
double blade and counterclockwise rotation
Impeller arm
Figure 3-18: Notations used for the c-shaped double blade impeller.
RUSHTON TURBINE
Axial impellers, such as the Rushton turbine, work by pumping the fluid toward the
walls of the vessel and then let the collision of the fluid with the wall supply the axial
mixing and the turbulence required for mixing. In order to avoid a tangential flows
that would result in a poor mixing, baffles are often used in combination with Rushton
turbines. Rushton turbines induce a relatively high shear rate and therefore appropriate
in processes where the fluid is not sensitive to shear rates.
The turbine can be designed with different blade lengths and also with different disk
diameters.
Impeller blade
Impeller disk
IMPELLER PARTS | 53
The figure below contains the notations used for the Rushton turbine in the part
library. The full list of parameters together with their description is found in Table 3-1
Direction of
rotation Impeller disk
Figure 3-22: Notations used for the Rushton turbine with backswept blades.
Impeller blade
Direction of
rotation Impeller disk
IMPELLER PARTS | 55
The figure below contains the notations used for the Smith turbine in the part library.
The full list of parameters together with their description is found in Table 3-1
Impeller Curvature
blade radius
Direction of
rotation
ANCHOR IMPELLER
This impeller type is suitable for mixing fluids of very high viscosity. For example, this
impeller is common when mixing Portland concrete and paints in small scales.The
lower part of the impeller is shaped as an ellipse that is cut in the middle. The major
axis of the ellipse equals the impeller diameter while the minor axis equals to the dished
tank where the impeller is placed minus the clearance between the impeller and tank
wall.
Upper hub
Major axis
Minor axis
IMPELLER PARTS | 57
The figure below contains the notations used for the anchor impeller in the part library.
The full list of parameters together with their description is found in Table 3-1
Impeller Parameters
The table below contains the names and descriptions of the input parameter used for
the impellers in the Mixer Module Part Library.
TABLE 3-1: MIXER MODULE PART LIBRARY PARAMETERS
IMPELLER PARTS | 59
60 | CHAPTER 3: MIXER MODULE PART LIBRARY
4
In this chapter:
61
Theory for the Free Surface Features
In this section:
• Deforming Domain
• Free Surface Conditions
• Contact Angle
• Rotating Shaft Conditions
• References for the Free Surface Features
Theory for the Rotating Machinery Interfaces and Theory for the
Nonisothermal Flow and Conjugate Heat Transfer Interfaces in the CFD
Module User’s Guide.
The links to the physics features described in other guides do not work in
the PDF, only from the online help in COMSOL Multiphysics.
Deforming Domain
The mesh within the free surface domain is deformed to account for the movement of
the free surface. This mesh movement is accomplished using a moving mesh approach.
The software perturbs the mesh nodes so that they conform with the free surface and
with other moving or stationary boundaries in the model. The boundary displacement
is propagated throughout the domain to obtain a smooth mesh deformation
everywhere. This is done by solving PDEs for the mesh displacements, using the
movement of the free surface as boundary condition. Taking two dimensions as an
example, a location in the deformed mesh with coordinates (x, y) can be related to its
coordinates in the original undeformed mesh (X,Y) by a function on the form:
x = x ( X, Y, t ), y = y ( X, Y, t )
The original, undeformed, mesh is referred to as the material frame (or reference
frame), while the deformed mesh is called the spatial frame. COMSOL Multiphysics
also defines geometry and mesh frames, which are coincident with the material frame
for this physics interface.
A number of mesh smoothing methods are available: Yeoh (the default), Laplace,
Winslow, and Hyperelastic. Note that the equations used for each smoothing type have
different properties.
Here pext is the pressure outside the free surface domain (SI unit: Pa) and fst denotes
the surface tension forces (SI unit: N/m2). In the surface tension terms, ∇s is the
surface gradient operator (∇s = (I − ni niT)∇ where I is the identity matrix) and σ is
the surface tension coefficient (SI unit: N/m).
The mesh velocity at the free surface is defined as the fluid velocity in the direction
normal to the surface:
( u mesh ⋅ n ) = ( u ⋅ n ) (4-2)
Contact Angle
At a three-phase boundary, it is necessary to add force terms to ensure that the fluid
maintains a consistent contact angle. The forces acting at the contact point are applied
to the fluid by the Contact Angle node (added per default under a Free Surface node).
In equilibrium, the surface tension forces and the normal restoring force from the
surface are in balance at a contact angle (θc), as shown in Figure 4-1. This equilibrium
is expressed by Young's equation, which considers the components of the forces in the
plane of the surface:
σ cos ( θ c ) + γ s1 = γ s2 (4-3)
Figure 4-1: The forces acting at a contact point. In equilibrium, the surface tension forces
and the normal restoring force from the surface are in balance at a contact angle θc.
The normal force balance at the solid surface is handled by the wall boundary
condition, which automatically prevents fluid flow across the solid boundary through
a no-penetration condition. The wall fluid interface feature applies a force, fwf, on the
fluid at the interface:
f wf = σ ( cos ( θ c ) – cos ( θ ) )m s
where θ is the actual contact angle and ms is the binormal to the solid surface.
r ax
u w = w ---------- × ( x – r bp ) (4-4)
r ax
2. W. Ren and E. Weinan, “Boundary Conditions for the Moving Contact Line
Problem,” Physics of Fluids, vol. 19, p. 022101, 2007.
3. W. Ren and D. Hu, “Continuum Models for the Contact Line Problem,” Physics of
Fluids, vol. 22, p. 102103, 2010.
The Rotating Machinery, Laminar Flow (spf) and Rotating Machinery, Turbulent Flow
(spf) interfaces, found under the Single-Phase Flow>Rotating Machinery branch ( )
when adding a physics interface, are used for modeling flow where one or more of
the boundaries rotate in a periodic fashion. This is used for mixers and propellers.
In this chapter:
• Moving Mesh
• The Rotating Machinery, Laminar Flow Interface
• The Rotating Machinery, Turbulent Flow, Algebraic yPlus Interface
• The Rotating Machinery, Turbulent Flow, L-VEL Interface
• The Rotating Machinery, Turbulent Flow, k-ε Interface
• The Rotating Machinery, Turbulent Flow, Realizable k-ε Interface
67
• The Rotating Machinery, Turbulent Flow, k-ω Interface
• The Rotating Machinery, Turbulent Flow, SST Interface
• The Rotating Machinery, Turbulent Flow, Low Re k-ε Interface
• The Rotating Machinery, Turbulent Flow, Spalart-Allmaras Interface
• The Rotating Machinery, Turbulent Flow, v2-f Interface
• Domain, Boundary, Point, and Pair Nodes for the Rotating Machinery Interfaces
ROTATING DOMAINS
When a Rotating Machinery, Fluid Flow interface is added using the Model Wizard, a
Rotating Domain node is added automatically in the Model Builder under the
Definitions>Moving Mesh node. It contains settings for specifying the rotation of one or
several domains. By default, All domains are selected. Under the Rotating Domain, you
can exclude any non-rotating fluid domains. For both 3D and 2D components, the
Rotation type should be set to Specified rotational velocity.
Select a Rotational velocity expression — General revolutions per time (the default),
Constant revolutions per time, General angular velocity, or Constant angular velocity.
• For Constant revolutions per time, or General revolutions per time, enter a value or
expression respectively in the Revolutions per time input field and optionally an Initial
angle. The angular velocity in this case is defined as the input multiplied by 2·π.
• For Constant angular velocity, or General angular velocity enter a value or an
expression respectively in the Angular velocity input field and optionally an Initial
angle.
3D Components
For 3D components, select the Rotation axis base point, rax. The default is the origin.
Select the Rotation axis, urot, the z-axis is the default. If the x-axis is selected, it
corresponds to a rotational axis (1, 0, 0) with the origin as the base point, and
correspondingly for the y-axis and z-axis.
2D Components
For 2D components, enter coordinates for the Rotation axis base point, rax. The
default is the origin (0, 0).
If there is more than one deforming domain, these must not intersect. Also a
deforming domain cannot intersect with a rotating domain.
Free Surface
There are two study types available for this physics interface. Using the Time Dependent
study type, rotation is achieved through moving mesh functionality, also known as
sliding mesh. Using the Frozen Rotor study type, the rotating parts are kept frozen in
position, and rotation is accounted for by the inclusion of centrifugal and Coriolis
forces. In both types, the momentum balance is governed by the Navier-Stokes
equations, and the mass conservation is governed by the continuity equation. See
Theory for the Rotating Machinery Interfaces in the CFD Module User’s Guide.
Pseudo Time Stepping for Laminar Flow Models in the CFD Module
User’s Guide
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is spf.
In addition to the settings described below, see The Creeping Flow, Laminar Flow,
Turbulent Flow, and Large Eddy Simulation Interfaces in the CFD Module User’s
Guide for all the other settings available. See Domain, Boundary, Point, and Pair
Nodes for the Rotating Machinery Interfaces for links to all the physics nodes.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
The equations solved by the Turbulent Flow, Algebraic yPlus interface are the
Reynolds-averaged Navier-Stokes (RANS) equations for conservation of momentum,
the continuity equation for conservation of mass, and an algebraic equation for the
scaled wall distance. Turbulence effects are included using an enhanced viscosity model
based on the scaled wall distance, and the physics interface therefore includes a wall
distance equation.
Two study types are available for this physics interface. Using the Time Dependent with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. Both study types includes an initial Wall
Distance Initialization step in order to compute the wall distance. See Theory for the
Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, Algebraic yPlus — Fluid Properties, Wall and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions and volume forces. You can
also right-click Turbulent Flow, Algebraic yPlus to select physics features from the
context menu. For information about the available physics features, see Domain,
Boundary, Point, and Pair Nodes for the Rotating Machinery Interfaces.
A part for the settings below, see the advanced settings for The Rotating Machinery,
Laminar Flow Interface, and The Creeping Flow, Laminar Flow, Turbulent Flow, and
Large Eddy Simulation Interfaces in the CFD Module User’s Guide.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
• Velocity field u
• Pressure p
• Reciprocal wall distance G
• Wall distance in viscous units yPlus (dimensionless)
The equations solved by the Turbulent Flow, L-VEL interface are the
Reynolds-averaged Navier-Stokes (RANS) equations for conservation of momentum,
the continuity equation for conservation of mass, and an algebraic equation for the
scaled fluid flow velocity tangential to the nearest wall. Turbulence effects are included
using an enhanced viscosity model based on the scaled wall distance, and the physics
interface therefore includes a wall distance equation.
Two study types available for this physics interface. Using the Time Dependent with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. Both study types includes an initial Wall
Distance Initialization step in order to compute the wall distance. See Theory for the
Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, L-VEL — Fluid Properties, Wall, and Initial Values. A
Moving Mesh interface with a Rotating Domain node is added automatically in the Model
Builder under the Definitions node. Then, from the Physics toolbar, add other nodes
A part for the interface settings defined below, see the advanced settings for The
Rotating Machinery, Laminar Flow Interface, and the settings in The Creeping Flow,
Laminar Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the CFD
Module User’s Guide.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
• Velocity field u
• Pressure p
• Reciprocal wall distance G
• Tangential velocity in viscous units uPlus (dimensionless)
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using the standard two-equation k-ε model with realizability constraints. Flow close to
walls is modeled using wall functions.
There are two study types available for this physics interface. Using the Time Dependent
study type, the rotation is achieved through moving mesh functionality, also known as
sliding mesh. Using the Frozen Rotor study type, the rotating parts are kept frozen in
position, and the rotation is accounted for by the inclusion of centrifugal and Coriolis
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, k-ε — Fluid Properties, Wall and Initial Values. A
Moving Mesh interface with a Rotating Domain node is added automatically in the Model
Builder under the Definitions node. Then, from the Physics toolbar, add other nodes
that implement, for example, boundary conditions and volume forces. You can also
right-click Turbulent Flow, k-ε to select physics features from the context menu. For
information about the available physics features, see Domain, Boundary, Point, and
Pair Nodes for the Rotating Machinery Interfaces.
A part for the defined settings below, see the advanced settings for The Rotating
Machinery, Laminar Flow Interface, and the settings in The Creeping Flow, Laminar
Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the CFD Module
User’s Guide.
For this physics interface the parameters are Ce1, Ce2, Cμ, σk, σe, κv, and B.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
There are two study types available for this physics interface. Using the Time Dependent
study type, the rotation is achieved through moving mesh functionality, also known as
sliding mesh. Using the Frozen Rotor study type, the rotating parts are kept frozen in
position, and the rotation is accounted for by the inclusion of centrifugal and Coriolis
forces. See Theory for the Rotating Machinery Interfaces in the CFD Module User’s
Guide.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, Realizable k-ε — Fluid Properties, Wall and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions and volume forces. You can
also right-click Turbulent Flow, Realizable k-ε to select physics features from the context
menu. For information about the available physics features, see Domain, Boundary,
Point, and Pair Nodes for the Rotating Machinery Interfaces.
A part for the defined settings below, see the advanced settings for The Rotating
Machinery, Laminar Flow Interface, and the settings in The Creeping Flow, Laminar
Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the CFD Module
User’s Guide.
For this physics interface the parameters are C2, A0, σk, σe, κv, and B.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using Wilcox revised two-equation k-ω model with realizability constraints. Flow close
to walls is modeled using wall functions.
There are two study types available for this physics interface. Using the Time Dependent
study type, the rotation is achieved through moving mesh functionality, also known as
sliding mesh. Using the Frozen Rotor study type, the rotating parts are kept frozen in
position, and the rotation is accounted for by the inclusion of centrifugal and Coriolis
forces. See Theory for the Rotating Machinery Interfaces in the CFD Module User’s
Guide.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, k-ω — Fluid Properties, Wall and Initial Values. A
Moving Mesh interface with a Rotating Domain node is added automatically in the Model
Builder under the Definitions node. Then, from the Physics toolbar, add other nodes
that implement, for example, boundary conditions and volume forces. You can also
right-click Turbulent Flow, k-ω to select physics features from the context menu. For
information about the available physics features, see Domain, Boundary, Point, and
Pair Nodes for the Rotating Machinery Interfaces.
A part for the settings defined below, see the advanced settings for The Rotating
Machinery, Laminar Flow Interface, and the settings in The Creeping Flow, Laminar
Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the CFD Module
User’s Guide.
For this physics interface the parameters are α, σk∗, σw, β0, β0∗, κv, and B.
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using the SST two-equation model with realizability constraints. The SST model is a
so-called low-Reynolds number model, which means that it resolves the flow all the
way down to the wall. The SST model depends on the distance to the closest wall. The
physics interface therefore includes a wall distance equation.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. See Theory for the Rotating
Machinery Interfaces in the CFD Module User’s Guide. In both study types
initialization is performed by first solving for the distance to the closest wall.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, SST — Fluid Properties, Wall and Initial Values. A
Moving Mesh interface with a Rotating Domain node is added automatically in the Model
Builder under the Definitions node. Then, from the Physics toolbar, add other nodes
that implement, for example, boundary conditions and volume forces. You can also
right-click Turbulent Flow, SST to select physics features from the context menu. For
information about the available physics features, see Domain, Boundary, Point, and
Pair Nodes for the Rotating Machinery Interfaces.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using the AKN two-equation k-ε model with realizability constraints. The AKN model
is a so-called low-Reynolds number model, which means that it resolves the flow all the
way down to the wall. The AKN model depends on the distance to the closest wall.
The physics interface therefore includes a wall distance equation.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, Low Re k-ε— Fluid Properties, Wall and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions and volume forces. You can
also right-click Turbulent Flow, Low Re k-ε to select physics features from the context
menu. For information about the available physics features, see Domain, Boundary,
Point, and Pair Nodes for the Rotating Machinery Interfaces.
A part for the interface settings defined below, see the advanced settings defined for
The Rotating Machinery, Laminar Flow Interface, and the settings in The Creeping
Flow, Laminar Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the
CFD Module User’s Guide.
For this physics interface the parameters are Ce1, Ce2, Cμ, σk, σe, and κv.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using the Spalart-Allmaras one-equation model. The Spalart-Allmaras model is a
so-called low-Reynolds number model, which means that it resolves the flow all the
way down to the wall. The Spalart-Allmaras model depends on the distance to the
closest wall. The physics interface therefore includes a wall distance equation.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. See Theory for the Rotating
Machinery Interfaces in the CFD Module User’s Guide. In both study types
initialization is performed by first solving for the distance to the closest wall.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, Spalart-Allmaras — Fluid Properties, Wall, and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions and volume forces. You can
also right-click Turbulent Flow, Spalart-Allmaras to select physics features from the
context menu. For information about the available physics features, see Domain,
Boundary, Point, and Pair Nodes for the Rotating Machinery Interfaces.
A part for the interface settings defined below, see the advanced settings defined for
The Rotating Machinery, Laminar Flow Interface, and the settings in The Creeping
Flow, Laminar Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the
CFD Module User’s Guide.
For this physics interface the parameters are Cb1, Cb2, Cv1, σv, Cw2, Cw3, κv, and Crot.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
The momentum balance is governed by the Navier-Stokes equations, and the mass
conservation is governed by the continuity equation. Turbulence effects are modeled
using the v2-f transport-equation model with realizability constraints. The v2-f model
is a so-called low-Reynolds number model, which means that it can resolve the flow all
the way down to the wall. The physics interface includes a wall distance equation.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. See Theory for the Rotating
Machinery Interfaces in the CFD Module User’s Guide. In both study types
initialization is performed by first solving for the distance to the closest wall.
When this physics interface is added, the following physics nodes are also added in the
Model Builder under Turbulent Flow, v2-f — Fluid Properties, Wall and Initial Values. A
Moving Mesh interface with a Rotating Domain node is added automatically in the Model
Builder under the Definitions node. Then, from the Physics toolbar, add other nodes
A part for the interface settings defined below, see the advanced settings defined for
The Rotating Machinery, Laminar Flow Interface, and the settings in The Creeping
Flow, Laminar Flow, Turbulent Flow, and Large Eddy Simulation Interfaces in the
CFD Module User’s Guide.
For this physics interface the parameters are Cμv2f, Cτ, Cη, CL, C1, C2, Cε1v2f, Cε2v2f,
A1, σkv2f, σεv2f, σζv2f, κv, and B.
DEPENDENT VARIABLES
The following dependent variables (fields) are defined for this physics interface:
• Contact Angle
• Free Surface
• Rotating Shaft
• Stationary Free Surface
The nodes are available under the option Rotating Machinery from the Physics ribbon
toolbar (Windows users), Physics context menu (Mac or Linux users), or right-click to
access the context menu (all users).
The following nodes, available in the Fluid Flow interfaces, are described for the
Laminar Flow interface in the CFD Module User’s Guide:
ˆ ⋅ ∇p
p ( x 0 ) – p ext + n η FS = – σ ∇ S2 η FS (5-1)
x = x0
The Stationary Free Surface feature is supported by laminar and turbulent flow, and is
applicable for small surface deformations (ηFS and its gradient should be small
compared to the dimensions of the computational domain). If the surface deformation
is large, a time-dependent study with a Deforming Domain and a Free Surface feature
should be used instead.
Iterative solvers are necessary to reduce the cost of models with large number of
degrees of freedom. The fluid flow interfaces use a Smoothed aggregation AMG solver
per default. It is recommendable to switch to a Geometric multigrid solver when
working with the Stationary Free Surface feature to ensure robustness of the model.
The Choice of Solver and Solver Settings in the CFD Module User’s
Guide and Studies and Solvers in the COMSOL Multiphysics Reference
Manual
• For Library coefficient, liquid/gas interface select an option from the list — Water/Air
(the default), Acetone/Air, Acetic acid/Air, Ethanol/Air, Ethylene glycol/Ethylene glycol
vapor, Diethyl ether/Air, Glycerol/Air, Heptane/Nitrogen, Mercury/Mercury vapor, or
Toluene/Air.
• For Library coefficient, liquid/liquid interface select an option from the list —
Benzene/Water, 20°C, Corn oil/Water, 20°C, Ether/Water, 20°C, Hexane/Water, 20°C,
Mercury/Water, 20°C, or Olive oil/Water, 20°C.
Only one Stationary Free Surface feature can be applied in a Rotating Machinery, Fluid
Flow interface, and it requires the study to be stationary (Frozen Rotor).
Ribbon
Physics tab with the physics interface selected:
Free Surface
This feature is used to define the two-phase interface separating the modeled fluid in
a free surface domain from the outer fluid. By applying this boundary condition the
motion of the two-phase interface can be tracked. The free surface follows the fluid
motion in the normal direction. In the tangential direction the free surface motion is
governed by the fluid motion in combination with the surface tension of the fluid-fluid
interface.
The Free Surface node can only be applied to outer boundaries of a Deforming Domain
that do not have any boundary condition assigned under Moving Mesh.
To control the attachment angle between the free surface and adjacent walls, a Contact
Angle node is by default added under the Free Surface node.
The Free Surface feature is supported by laminar flow and by turbulence models using
wall functions or automatic wall treatment but not low-Reynolds number wall
FREE SURFACE
Enter an External pressure pext (SI unit: Pa) corresponding to the pressure level in the
external fluid phase.
SURFACE TENSION
The Surface tension coefficient σ (SI unit: N/m) can be specified from predefined
libraries, by selecting Library coefficient, liquid/gas interface (the default) or Library
coefficient, liquid/liquid interface, or set to User defined.
• For Library coefficient, liquid/gas interface select an option from the list — Water/Air
(the default), Acetone/Air, Acetic acid/Air, Ethanol/Air, Ethylene glycol/Ethylene glycol
vapor, Diethyl ether/Air, Glycerol/Air, Heptane/Nitrogen, Mercury/Mercury vapor, or
Toluene/Air.
• For Library coefficient, liquid/liquid interface select an option from the list —
Benzene/Water, 20°C, Corn oil/Water, 20°C, Ether/Water, 20°C, Hexane/Water, 20°C,
Mercury/Water, 20°C, or Olive oil/Water, 20°C.
Context Menus
Right click the physics interface and select Rotating Machinery>Free Surface
Ribbon
Physics tab with the physics interface selected:
Contact Angle
This feature specifies the contact angle between the free surface and a solid wall. The
Contact Angle feature is a subfeature to, and added under a Stationary Free Surface or
a Free Surface feature.
CONTACT ANGLE
Select an option from the Specify contact angle list — Directly (the default) or Through
Young’s equation.
• For Directly enter a Contact angle θw (SI unit: rad). The default is π/2.
• For Through Young’s equation enter values or expressions for Phase 1-Solid surface
energy density γs1 (SI unit: J/m2) and Phase 2-Solid surface energy density γs2 (SI
unit: J/m2).
The remainder of this section applies when the Contact Angle feature is a subfeature
to the Free Surface feature.
The condition applies forces on the wall that balance the surface tension for the
prescribed contact angle. For laminar flow, the Contact Angle node should be used on
edges or points where a Wall feature with the Navier Slip option is prescribed adjacent
to the free surface. In this case, the Navier Slip option is required for the contact line
to move along the wall.
The representation of the contact angle is dependent on the mesh resolution in the
region where the fluid-fluid interface attaches to the wall. If, during the simulation, the
contact angle is found to fluctuate around a given value, this indicates that the
resolution needs to be improved. This can for example be achieved by adding more
mesh boundary layers on the wall.
ROTATING SHAFT
Axis of Rotation
Select the Axis of rotation. If x-axis is selected, this corresponds to a rotational axis
(1, 0, 0) with the origin as the base point. Correspondingly, if y-axis or z-axis is
selected, the rotational axis is defined by a unit vector in the selected direction together
with a base point at the origin. For User defined enter values into these fields: Rotation
axis base point and Rotation axis direction. The z-axis is the default choice and the
default values in the fields for User defined correspond to those for the z-axis.
Angular Velocity
Select the rotating shaft Rotational frequency from the list:
• Select the Revolution per time option to manually define the rotational velocity (SI
unit: rad/s) and select the Rotational direction as Positive angular velocity or Negative
angular velocity. This is the default setting.
• Select the Angular velocity option to manually define the angular velocity (SI unit:
rad/s).
Context Menus
Right click the physics interface and select Rotating Machinery>Rotating Shaft
Ribbon
Physics tab with the physics interface selected:
Rotating Machinery,
Multiphase Flow Interfaces
This chapter describes the physics interfaces found under the Fluid Flow>Multiphase
Flow>Rotating Machinery, Multiphase Flow branch ( ). In this chapter:
91
The Rotating Machinery, Multiphase
Flow, Mixture Model, Laminar Flow
and Turbulent Flow Interfaces
The Rotating Machinery, Mixture Model, Laminar Flow Interface
The Laminar Flow version of the Rotating Machinery, Mixture Model (mm) interface ( )
is used to model low to moderate Reynolds number flow of a multiphase mixture in
equipment containing one or more rotating parts. The mixture consists of a liquid with
dispersed particles, bubbles, liquid droplets, or solid particles, which are assumed to
always travel with terminal velocity.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. The physics interface can also model the distribution of the
number density, which in turn can be used to calculate the interfacial area, useful when
simulating chemical reactions in the mixture.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
When this physics interface is added, the following default nodes are also added in the
Model Builder under Mixture Model, Laminar Flow — Mixture Properties, Wall, and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
SETTINGS
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics
interface. Refer to such physics interface variables in expressions using the pattern
<name>.<variable_name>. In order to distinguish between variables belonging to
different physics interfaces, the name string must be unique. Only letters, numbers, and
underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is mm.
Most of the physics interface settings are the same as described for The Mixture Model,
Laminar Flow Interface in the CFD Module User’s Guide, and for Moving Mesh.
Only unique settings are discussed below.
ADVANCED SETTINGS
To display this section, click the Show button ( ) and select Advanced Physics Options.
Normally these settings do not need to be changed.
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
CFLCMP which in turn triggers a PID regulator for the CFL number. For Manual enter
a Local CFL number CFLloc (dimensionless).
Pseudo Time Stepping for Laminar Flow Models in the CFD Module
User’s Guide
DEPENDENT VARIABLES
Enter values for the dependent variables (field variables):
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are included using and enhanced viscosity
model based on the scaled wall distance, and the physics interface therefore includes a
wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, Algebraic yPlus — Mixture Properties, Wall, and
Initial Values. A Moving Mesh interface with a Rotating Domain node is added
automatically in the Model Builder under the Definitions node. Then, from the Physics
toolbar, add other nodes that implement, for example, boundary conditions. You can
also right-click Mixture Model, Algebraic yPlus to select physics features from the context
menu.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, Algebraic yPlus Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are included using and enhanced viscosity
model based on the local wall distance, and the physics interface therefore includes a
wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. All study types include an initial
Wall Distance Initialization step in order to compute the wall distance. See Theory for
the Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, L-VEL — Mixture Properties, Wall, and Initial Values.
A Moving Mesh interface with a Rotating Domain node is added automatically in the
Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions. You can also right-click
Mixture Model, L-VEL to select physics features from the context menu.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, L-VEL Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the standard
two-equation k-ε model with realizability constraints. Flow close to walls is modeled
using wall functions.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, k-ε — Mixture Properties, Wall, and Initial Values.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, k-ε Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the realizable
two-equation k-ε model. Flow close to walls is modeled using wall functions.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, Realizable k-ε — Mixture Properties, Wall, and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
nodes that implement, for example, boundary conditions. You can also right-click
Mixture Model, Realizable k-ε to select physics features from the context menu.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, Realizable k-ε Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the Wilcox revised
two-equation k-ω model with realizability constraints. The k-ω model is a so-called
low-Reynolds-number model, which means that it can resolve the flow all the way
down to the wall.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, k-ω — Mixture Properties, Wall, and Initial Values.
A Moving Mesh interface with a Rotating Domain node is added automatically in the
Model Builder under the Definitions node. Then, from the Physics toolbar, add other
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, k-ω Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the Menter shear-stress
transport (SST) two-equation model from 2003 with realizability constraints. The SST
model is a so-called low-Reynolds-number model, which means that it can resolve the
flow all the way down to the wall. The SST model depends on the distance to the
closest wall. The physics interface therefore includes a wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. All study types include an initial
Wall Distance Initialization step in order to compute the wall distance. See Theory for
the Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, SST — Mixture Properties, Wall, and Initial Values.
A Moving Mesh interface with a Rotating Domain node is added automatically in the
Model Builder under the Definitions node. Then, from the Physics toolbar, add other
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
nodes that implement, for example, boundary conditions. You can also right-click
Mixture Model, SST to select physics features from the context menu.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, SST Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the AKN two-equation
k-ε model with realizability constraints. The AKN model is a so-called
low-Reynolds-number model, which means that it can resolve the flow all the way
down to the wall. The AKN model depends on the distance to the closest wall. The
physics interface therefore includes a wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. All study types include an initial
Wall Distance Initialization step in order to compute the wall distance. See Theory for
the Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, Low Re k-ε — Mixture Properties, Wall, and Initial
Values. A Moving Mesh interface with a Rotating Domain node is added automatically in
the Model Builder under the Definitions node. Then, from the Physics toolbar, add other
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, Low Re k-ε Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the Spalart-Allmaras
one-equation model. The Spalart-Allmaras model is a so-called low-Reynolds-number
model, which means that it can resolve the flow all the way down to the wall. The
Spalart-Allmaras model depends on the distance to the closest wall. The physics
interface therefore includes a wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. All study types include an initial
Wall Distance Initialization step in order to compute the wall distance. See Theory for
the Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, Spalart-Allmaras — Mixture Properties, Wall, and
Initial Values. A Moving Mesh interface with a Rotating Domain node is added
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
automatically in the Model Builder under the Definitions node. Then, from the Physics
toolbar, add other nodes that implement, for example, boundary conditions. You can
also right-click Mixture Model, Spalart-Allmaras to select physics features from the
context menu.
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, Spalart-Allmaras Interface in the CFD Module User’s Guide.
The interface solves one set of Navier-Stokes equations for the momentum of the
mixture. The pressure distribution is calculated from a mixture-averaged continuity
equation and the velocity of the dispersed phase is described by a slip model. The
dispersed particles are tracked by solving a transport equation for the volume fraction
of the dispersed phase. Turbulence effects are modeled using the v2-f
transport-equation model with realizability constraints. The v2-f model is a so-called
low-Reynolds-number model, which means that it can resolve the flow all the way
down to the wall. The v2-f model depends on the distance to the closest wall. The
physics interface therefore includes a wall distance equation.
The physics interface can also model the distribution of the number density, which in
turn can be used to calculate the interfacial area, useful when simulating chemical
reactions in the mixture.
There are two study types available for this physics interface. Using the Time Dependent
with Initialization study type, the rotation is achieved through moving mesh
functionality, also known as sliding mesh. Using the Frozen Rotor with Initialization
study type, the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces. All study types include an initial
Wall Distance Initialization step in order to compute the wall distance. See Theory for
the Rotating Machinery Interfaces in the CFD Module User’s Guide.
When this physics interface is added, the following default nodes are also added in the
Model Builder, under Mixture Model, v2-f — Mixture Properties, Wall, and Initial Values.
A Moving Mesh interface with a Rotating Domain node is added automatically in the
The physics interface settings are the same as described for The Rotating Machinery,
Mixture Model, Laminar Flow Interface, for Moving Mesh, and for the The Mixture
Model, v2-f Interface in the CFD Module User’s Guide.
• Contact Angle
• Free Surface
• Rotating Shaft
The nodes are available from the Physics ribbon toolbar (Windows users), Physics
context menu (Mac or Linux users), or right-click to access the context menu (all
users).
The following nodes (listed in alphabetical order) are described for the Mixture Model
interface in the CFD Module User’s Guide:
• Gravity • Outlet
• Initial Values • Periodic Flow Condition
• Inlet • Pressure Point Constraint1
• Interior Wall • Symmetry
• Mass Transfer • Volume Force1
• Mixture Properties • Wall
• Mixture Continuity
1
Described for the Laminar Flow interface in the CFD Module User’s Guide.
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
Free Surface
The Free Surface node is used to define the free surface separating an interior modeled
fluid from the outer fluid. The free surface follows the fluid motion, in the normal
direction, and is governed by the fluid motion in the free surface domain and the
surface tension of the fluid-fluid interface. The Free Surface node can only be applied
to outer boundaries of a Deforming Domain that do not have any boundary condition
assigned under Moving Mesh.
To control the attachment angle between the free surface and adjacent walls, a Contact
Angle node is by default added under the Free Surface node.
The Free Surface feature is supported by laminar flow and by turbulence models using
wall functions or automatic wall treatment, but not low-Reynolds number wall
treatment. It also requires the study to be time dependent. The node is inactive for
Frozen Rotor studies.
FREE SURFACE
Enter an External pressure pext (SI unit: Pa) corresponding to the pressure level in the
external fluid phase.
SURFACE TENSION
The Surface tension coefficient σ (SI unit: N/m) can be specified from predefined
libraries, by selecting Library coefficient, liquid/gas interface (the default) or Library
coefficient, liquid/liquid interface, or set to User defined.
• For Library coefficient, liquid/gas interface select an option from the list — Water/Air
(the default), Acetone/Air, Acetic acid/Air, Ethanol/Air, Ethylene glycol/Ethylene glycol
vapor, Diethyl ether/Air, Glycerol/Air, Heptane/Nitrogen, Mercury/Mercury vapor, or
Toluene/Air.
• For Library coefficient, liquid/liquid interface select an option from the list —
Benzene/Water, 20°C, Corn oil/Water, 20°C, Ether/Water, 20°C, Hexane/Water, 20°C,
Mercury/Water, 20°C, or Olive oil/Water, 20°C.
Context Menus
Right click the physics interface and select Rotating Machinery>Free Surface
Ribbon
Physics tab with the physics interface selected:
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
Contact Angle
This feature specifies the contact angle between the free surface and a solid wall. The
Contact Angle feature is a subfeature to, and added under a Free Surface feature.
CONTACT ANGLE
Select an option from the Specify contact angle list—Directly (the default) or Through
Young’s equation.
• For Directly enter a Contact angle θw (SI unit: rad). The default is π/2.
• For Through Young’s equation enter values or expressions for Phase 1-Solid surface
energy density γs1 (SI unit: J/m2) and Phase 2-Solid surface energy density γs2 (SI
unit: J/m2).
This condition applies forces on the wall that balance the surface tension for the
prescribed contact angle.
For laminar flow, the Contact Angle node should be used on edges or points where a
Navier Slip wall boundary condition is applied adjacent to the free surface. In this case,
the Navier Slip condition is required for the contact line to move along the wall.
The representation of the contact angle is dependent on the mesh resolution in the
region where the fluid-fluid interface attaches to the wall. If, during the simulation, the
contact angle is found to fluctuate around a given value, this indicates that the
resolution needs to be improved. This can for example be achieved by adding more
mesh boundary layers on the wall.
Rotating Shaft
Use the Rotating Shaft node to define boundaries where the fluid is subjected to a solid
shaft rotation. This feature should be applied to boundaries outside of, but connecting
to, a Rotating Domain. This feature can be applied to rotating shaft boundaries in a
free-surface domain modeled by a Deforming Domain.
Axis of Rotation
Select the Axis of rotation. If x-axis is selected, this corresponds to a rotational axis
(1, 0, 0) with the origin as the base point. Correspondingly, if y-axis or z-axis is
selected, the rotational axis is defined by a unit vector in the selected direction together
with a base point at the origin. For User defined enter values into these fields: Rotation
axis base point and Rotation axis direction. The z-axis is the default choice and the
default values in the fields for User defined correspond to those for the z-axis.
Angular Velocity
Select the rotating shaft Rotational frequency from the list:
• Select the Revolution per time option to manually define the rotational velocity (SI
unit: rad/s) and select the Rotational direction as Positive angular velocity or Negative
angular velocity. This is the default setting.
• Select the Angular velocity option to manually define the angular velocity (SI unit:
rad/s).
THE ROTATING MACHINERY, MULTIPHASE FLOW, MIXTURE MODEL, LAMINAR FLOW AND TURBULENT FLOW INTERFACES
LOCATION IN USER INTERFACE
Context Menus
Right click the physics interface and select Rotating Machinery>Rotating Shaft
Ribbon
Physics tab with the physics interface selected:
Rotating Machinery,
109
• Rotating Machinery, Nonisothermal Flow, Turbulent Flow, v2-f Interface
• Domain, Boundary, Edge, Point, and Pair Nodes for the Rotating Machinery,
Nonisothermal Flow Interfaces
This physics interface combines the capabilities of the Rotating Machinery, Laminar
Flow interface and the Heat Transfer in Fluids interface, and can be used to simulate
fluid flows where the fluid properties depend on the temperature. Furthermore it is
also possible to include heat transfer in solids, stationary and rotating, as well as
surface-to-surface radiation and radiation in participating media. The physics interface
is available in 2D and 3D and supports low Mach number (typically less than 0.3)
flows, as well as non-Newtonian fluids.
Four study types are available for this physics interface. Using the Time Dependent study
type, the rotation is achieved through moving mesh functionality, also known as sliding
mesh. Using the Frozen Rotor study, the rotating parts are kept frozen in position, and
the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.The
Time Dependent, One-Way Coupled, NITF and Frozen Rotor, One-Way Coupled, NITF study
are used for frozen rotor and time dependent nonisothermal flow computations where
the model is one-way coupled in the sense that the temperature variations are so small
that they do not affect the flow.
When adding the multiphysics interface, the Laminar Flow and Heat Transfer in Fluids
interfaces are added to the Model Builder. A Moving Mesh interface with a Rotating
Domain node is added automatically in the Model Builder under the Definitions node. In
addition, the Multiphysics node is added, which automatically includes the multiphysics
coupling feature Nonisothermal Flow. The following default nodes are also added in the
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 111
Model Builder under Laminar Flow — Fluid Properties, Wall, and Initial Values; and under
Heat Transfer in Fluids — Fluid, Thermal Insulation, and Initial Values.
• Domain, Boundary, Edge, Point, and Pair Nodes for the Rotating
Machinery, Nonisothermal Flow Interfaces
• Moving Mesh
• Theory for the Free Surface Features
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, Algebraic yPlus interface and the Heat Transfer in Fluids interface, and can be
used to simulate high Reynolds number flows where the fluid properties depend on
the temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
The physics interface solves for conservation of energy, mass, and momentum in fluids,
and for conservation of energy in solids. The interface also solves an algebraic equation
for the scaled wall distance. Turbulence effects are included using an enhanced
viscosity model based on the scaled wall distance, and the physics interface therefore
includes a wall distance equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. The Time Dependent, One-Way Coupled,
NITF and Frozen Rotor, One-Way Coupled, NITF study are used for frozen rotor and time
dependent nonisothermal flow computations where the model is one-way coupled in
the sense that the temperature variations are so small that they do not affect the flow.
All study types include an initial Wall Distance Initialization step in order to compute
the wall distance. See Theory for the Rotating Machinery Interfaces in the CFD
Module User’s Guide.
When adding the multiphysics interface, the Turbulent Flow, Algebraic yPlus and Heat
Transfer in Fluids interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Nonisothermal Flow. The following default
nodes are also added in the Model Builder under Turbulent Flow, Algebraic yPlus — Fluid
Properties, Wall, and Initial Values; and under Heat Transfer in Fluids— Fluid, Thermal
Insulation, and Initial Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 113
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, L-VEL interface and the Heat Transfer in Fluids interface, and can be used to
simulate high Reynolds number flows where the fluid properties depend on the
temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids,
and for conservation of energy in solids. The interface also solves an algebraic equation
for the scaled fluid flow velocity tangential to the nearest wall. Turbulence effects are
included using an enhanced viscosity model based on the scaled wall distance, and the
physics interface therefore includes a wall distance equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. The Time Dependent, One-Way Coupled,
NITF and Frozen Rotor, One-Way Coupled, NITF study are used for frozen rotor and time
dependent nonisothermal flow computations where the model is one-way coupled in
the sense that the temperature variations are so small that they do not affect the flow.
All study types include an initial Wall Distance Initialization step in order to compute
the wall distance. See Theory for the Rotating Machinery Interfaces in the CFD
Module User’s Guide.
When adding the multiphysics interface, the Turbulent Flow, L-VEL and Heat Transfer in
Fluids interfaces are added to the Model Builder. A Moving Mesh interface with a
Rotating Domain node is added automatically in the Model Builder under the Definitions
node. In addition, the Multiphysics node is added, which automatically includes the
multiphysics coupling feature Nonisothermal Flow. The following default nodes are also
added in the Model Builder under Turbulent Flow, L-VEL — Fluid Properties, Wall, and
Initial Values; and under Heat Transfer in Fluids— Fluid, Thermal Insulation, and Initial
Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, k-ε interface and the Heat Transfer in Fluids interface, and can be used to
simulate high Reynolds number fluid flows where the fluid properties depend on the
temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids,
and for conservation of energy in solids. Turbulence effects are modeled using the
standard two-equation k-ε model with realizability constraints. Flow and heat transfer
close to walls are modeled using wall functions.
Four study types are available for this physics interface. Using the Time Dependent study
type, the rotation is achieved through moving mesh functionality, also known as sliding
mesh. Using the Frozen Rotor study, the rotating parts are kept frozen in position, and
the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.The
Time Dependent, One-Way Coupled, NITF and Frozen Rotor, One-Way Coupled, NITF study
are used for frozen rotor and time dependent nonisothermal flow computations where
the model is one-way coupled in the sense that the temperature variations are so small
that they do not affect the flow.
When adding the multiphysics interface, the Turbulent Flow, k-ε and Heat Transfer in
Fluids interfaces are added to the Model Builder. A Moving Mesh interface with a
Rotating Domain node is added automatically in the Model Builder under the Definitions
node. In addition, the Multiphysics node is added, which automatically includes the
multiphysics coupling feature Nonisothermal Flow. The following default nodes are also
added in the Model Builder under Turbulent Flow, k-ε — Fluid Properties, Wall, and Initial
Values; and under Heat Transfer in Fluids — Fluid, Thermal Insulation, and Initial Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 115
Rotating Machinery, Nonisothermal Flow, Turbulent Flow,
Realizable k-ε Interface
The Turbulent Flow, Realizable k-ε version of the Rotating Machinery, Nonisothermal Flow
(rmnitf) interface ( ), found under the Fluid Flow>Nonisothermal Flow>Rotating
Machinery, Nonisothermal Flow branch ( ) when adding a physics interface, is used to
simulate turbulent flow and heat transfer in equipment containing one or more
rotating parts.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, Realizable k-ε interface and the Heat Transfer in Fluids interface, and can be
used to simulate high Reynolds number fluid flows where the fluid properties depend
on the temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids,
and for conservation of energy in solids. Turbulence effects are modeled using the
Realizable two-equation k-ε model. Flow and heat transfer close to walls are modeled
using wall functions.
Four study types are available for this physics interface. Using the Time Dependent study
type, the rotation is achieved through moving mesh functionality, also known as sliding
mesh. Using the Frozen Rotor study, the rotating parts are kept frozen in position, and
the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.The
Time Dependent, One-Way Coupled, NITF and Frozen Rotor, One-Way Coupled, NITF study
are used for frozen rotor and time dependent nonisothermal flow computations where
the model is one-way coupled in the sense that the temperature variations are so small
that they do not affect the flow.
When adding the multiphysics interface, the Turbulent Flow, Realizable k-ε and Heat
Transfer in Fluids interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Nonisothermal Flow. The following default
nodes are also added in the Model Builder under Turbulent Flow, Realizable k-ε — Fluid
Properties, Wall, and Initial Values; and under Heat Transfer in Fluids — Fluid, Thermal
Insulation, and Initial Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, k-ω interface and the Heat Transfer in Fluids interface, and can be used to
simulate high Reynolds number fluid flows where the fluid properties depend on the
temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids
and for conservation of energy in solids. Turbulence effects are modeled using the
Wilcox revised two-equation k-ω model with realizability constraints. The k-ω model
is a so-called low-Reynolds number model, which means that it can resolve the flow
and heat transfer all the way down to the wall.
Four study types are available for this physics interface. Using the Time Dependent study
type, the rotation is achieved through moving mesh functionality, also known as sliding
mesh. Using the Frozen Rotor study, the rotating parts are kept frozen in position, and
the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.The
Time Dependent, One-Way Coupled, NITF and Frozen Rotor, One-Way Coupled, NITF study
are used for frozen rotor and time dependent nonisothermal flow computations where
the model is one-way coupled in the sense that the temperature variations are so small
that they do not affect the flow.
When adding the multiphysics interface, the Turbulent Flow, k-ω and Heat Transfer in
Fluids interfaces are added to the Model Builder. A Moving Mesh interface with a
Rotating Domain node is added automatically in the Model Builder under the Definitions
node. In addition, the Multiphysics node is added, which automatically includes the
multiphysics coupling feature Nonisothermal Flow. The following default nodes are also
added in the Model Builder under Turbulent Flow, k-ω — Fluid Properties, Wall, and
Initial Values; and under Heat Transfer in Fluids — Fluid, Thermal Insulation, and Initial
Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 117
Rotating Machinery, Nonisothermal Flow, Turbulent Flow, SST
Interface
The Turbulent Flow, SST version of the Rotating Machinery, Nonisothermal Flow
interface ( ), found under the Fluid Flow>Nonisothermal Flow>Rotating Machinery,
Nonisothermal Flow branch ( ) when adding a physics interface, is used to simulate
turbulent flow and heat transfer in equipment containing one or more rotating parts.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, SST interface and the Heat Transfer in Fluids interface, and can be used to
simulate high Reynolds number fluid flows where the fluid properties depend on the
temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids
and for conservation of energy in solids. Turbulence effects are modeled using the SST
two-equation model with realizability constraints. The SST model is a so-called
low-Reynolds number model, which means that it can resolve the flow and heat
transfer all the way down to the wall. The SST model depends on the distance to the
closest wall. The physics interface therefore includes a wall distance equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. The Time Dependent, One-Way Coupled,
NITF and Frozen Rotor, One-Way Coupled, NITF study are used for frozen rotor and time
dependent nonisothermal flow computations where the model is one-way coupled in
the sense that the temperature variations are so small that they do not affect the flow.
All study types include an initial Wall Distance Initialization step in order to compute
the wall distance. See Theory for the Rotating Machinery Interfaces in the CFD
Module User’s Guide.
When adding the multiphysics interface, the Turbulent Flow, SST and Heat Transfer in
Fluids interfaces are added to the Model Builder. A Moving Mesh interface with a
Rotating Domain node is added automatically in the Model Builder under the Definitions
node. In addition, the Multiphysics node is added, which automatically includes the
multiphysics coupling feature Nonisothermal Flow. The following default nodes are also
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, Low Re k-ε interface and the Heat Transfer in Fluids interface, and can be used
to simulate high Reynolds number fluid flows where the fluid properties depend on
the temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids
and for conservation of energy in solids. Turbulence effects are modeled using the
AKN two-equation k-ε model with realizability constraints. The AKN model is a
so-called low-Reynolds number model, which means that it can resolve the flow and
heat transfer all the way down to the wall. The AKN model depends on the distance
to the closest wall. The physics interface therefore includes a wall distance equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. The Time Dependent, One-Way Coupled,
NITF and Frozen Rotor, One-Way Coupled, NITF study are used for frozen rotor and time
dependent nonisothermal flow computations where the model is one-way coupled in
the sense that the temperature variations are so small that they do not affect the flow.
All study types include an initial Wall Distance Initialization step in order to compute
the wall distance. See Theory for the Rotating Machinery Interfaces in the CFD
Module User’s Guide.
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 119
When adding the multiphysics interface, the Turbulent Flow, low Re k-ε and Heat
Transfer in Fluids interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Nonisothermal Flow. The following default
nodes are also added in the Model Builder under Turbulent Flow, low Re k-ε — Fluid
Properties, Wall, and Initial Values; and under Heat Transfer in Fluids— Fluid, Thermal
Insulation, and Initial Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, Spalart-Allmaras interface and the Heat Transfer in Fluids interface, and can be
used to simulate high Reynolds number flows where the fluid properties depend on
the temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids
and for conservation of energy in solids. Turbulence effects are modeled using the
Spalart-Allmaras one-equation model. The Spalart-Allmaras model is a so-called
low-Reynolds number model, which means that it can resolve the flow and heat
transfer all the way down to the wall. The Spalart-Allmaras model depends on the
distance to the closest wall, and the physics interface therefore includes a wall distance
equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
When adding the multiphysics interface, the Turbulent Flow, Spalart-Allmaras and Heat
Transfer in Fluids interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Nonisothermal Flow. The following default
nodes are also added in the Model Builder under Turbulent Flow, Spalart-Allmaras — Fluid
Properties, Wall, and Initial Values; and under Heat Transfer in Fluids— Fluid, Thermal
Insulation, and Initial Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
This physics interface combines the capabilities of the Rotating Machinery, Turbulent
Flow, v2-f interface and the Heat Transfer in Fluids interface, and can be used to
simulate high Reynolds number flows where the fluid properties depend on the
temperature. Furthermore it is also possible to include heat transfer in solids,
stationary and rotating, as well as surface-to-surface radiation and radiation in
participating media. The physics interface is available in 2D and 3D and supports low
Mach number (typically less than 0.3) flows, as well as non-Newtonian fluids.
The physics interface solves for conservation of energy, mass, and momentum in fluids
and for conservation of energy in solids. Turbulence effects are modeled using the v2-f
transport-equation model with realizability constraints. The v2-f model is a so-called
low-Reynolds number model, which means that it can resolve the flow and heat
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 121
transfer all the way down to the wall. The physics interface includes a wall distance
equation.
Four study types available for this physics interface. Using the Transient with
Initialization study type, the rotation is achieved through moving mesh functionality,
also known as sliding mesh. Using the Frozen Rotor with Initialization study type, the
rotating parts are kept frozen in position, and the rotation is accounted for by the
inclusion of centrifugal and Coriolis forces. The Time Dependent, One-Way Coupled,
NITF and Frozen Rotor, One-Way Coupled, NITF study are used for frozen rotor and time
dependent nonisothermal flow computations where the model is one-way coupled in
the sense that the temperature variations are so small that they do not affect the flow.
All study types include an initial Wall Distance Initialization step in order to compute
the wall distance. See Theory for the Rotating Machinery Interfaces in the CFD
Module User’s Guide.
When adding the multiphysics interface, the Turbulent Flow, v2-f and Heat Transfer in
Fluids interfaces are added to the Model Builder. A Moving Mesh interface with a
Rotating Domain node is added automatically in the Model Builder under the Definitions
node. In addition, the Multiphysics node is added, which automatically includes the
multiphysics coupling feature Nonisothermal Flow. The following default nodes are also
added in the Model Builder under Turbulent Flow, v2-f — Fluid Properties, Wall, and
Initial Values; and under Heat Transfer in Fluids— Fluid, Thermal Insulation, and Initial
Values.
See Rotating Machinery, Nonisothermal Flow, Laminar Flow Interface for details.
Domain, Boundary, Edge, Point, and Pair Nodes for the Rotating
Machinery, Nonisothermal Flow Interfaces
Physics nodes are available from the Physics ribbon toolbar (Windows users), Physics
context menu (Mac or Linux users), or right-click to access the context menu (all
users).
LAMINAR FLOW
The available physics features for The Rotating Machinery, Laminar Flow Interface are
listed in the section Domain, Boundary, Point, and Pair Nodes for the Rotating
Machinery Interfaces.
T U R B U L E N T F L O W , k-ε
The available physics features for The Rotating Machinery, Turbulent Flow, k-ε
Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
T U R B U L E N T F L O W , R E A L I Z A B L E k-ε
The available physics features for The Rotating Machinery, Turbulent Flow, Realizable
k-ε Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
T U R B U L E N T F L O W , k-ω
The available physics features for The Rotating Machinery, Turbulent Flow, k-ω
Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
THE ROTATING MACHINERY, NONISOTHERMAL FLOW, LAMINAR AND TURBULENT FLOW INTERFACES | 123
T U R B U L E N T F L O W , L O W R E k-ε
The available physics features for The Rotating Machinery, Turbulent Flow, Low Re
k-ε Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
Rotating Machinery,
Reacting Flow Interfaces
This chapter describes the physics interfaces found under the Chemical Species
Transport>Rotating Machinery, Reacting Flow branch ( ).
In this chapter:
125
The Rotating Machinery, Reacting
Flow, Laminar Flow and Turbulent
Flow Interfaces
The Rotating Machinery, Reacting Flow, Laminar Flow Interface
The Laminar Flow version of the Rotating Machinery, Reacting Flow interface ( ) is
used to simulate the laminar flow and transport of chemical species in equipment
containing one or more rotating parts. Reactions between the participating species can
also be studied.
The physics interface combines the functionality of the Rotating Machinery, Laminar
Flow and the Transport of Concentrated Species interfaces. The Reacting Flow
multiphysics coupling, which is added automatically, couples fluid flow and mass
transport. The instantaneous velocity and pressure fields are solved together with an
arbitrary number of mass fractions. The momentum balance is governed by the
Navier-Stokes equations, and the total mass conservation is governed by the continuity
equation. For the chemical species, the governing transport equations include
convection, diffusion and, optionally, migration in an electric field.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type (see Frozen Rotor in the CFD Module User’s
Guide), the rotating parts are kept frozen in position, and the rotation is accounted
for by the inclusion of centrifugal and Coriolis forces.
When adding the multiphysics interface, the Laminar Flow and Transport of
Concentrated Species interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Reacting Flow. The following default nodes
are also added in the Model Builder under Laminar Flow — Fluid Properties, Wall, and
This physics interface combines the functionality of the Rotating Machinery, Turbulent
Flow, k-e and the Transport of Concentrated Species interfaces. The Reacting Flow
multiphysics coupling, which is added automatically, couples fluid flow and mass
transport. The averaged velocity and pressure fields are solved together with an
arbitrary number of averaged mass fractions. The momentum balance is governed by
THE ROTATING MACHINERY, REACTING FLOW, LAMINAR FLOW AND TURBULENT FLOW INTERFACES | 127
the averaged Navier-Stokes equations, and the total mass conservation is governed by
the continuity equation. The fluid-flow turbulence is modeled using the standard
two-equation k-ε model with realizability constraints. For the chemical species, the
governing transport equations include convection, diffusion and, optionally, migration
in an electric field. Flow close to solid walls is modeled using wall functions for the flow
and species transport. Turbulent reactions are modeled using the eddy dissipation
model.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
When adding the multiphysics interface, the Turbulent Flow, k-ε and Transport of
Concentrated Species interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Reacting Flow. The following default nodes
are also added in the Model Builder under Turbulent Flow, k-ε — Fluid Properties, Wall,
and Initial Values; and under Transport of Concentrated Species — Transport Properties,
No Flux, and Initial Values.
See The Rotating Machinery, Reacting Flow, Laminar Flow Interface for details.
The physics interface combines the functionality of the Rotating Machinery, Turbulent
Flow k-w and the Transport of Concentrated Species interfaces. The Reacting Flow
multiphysics coupling, which is added automatically, couples fluid flow and mass
transport. The averaged velocity and pressure fields are solved together with an
arbitrary number of averaged mass fractions. The momentum balance is governed by
the averaged Navier-Stokes equations, and the total mass conservation is governed by
the continuity equation. The fluid-flow turbulence is modeled using the Wilcox
revised k-ω model with realizability constraints. The k-ω model is a so-called
low-Reynolds number model, which means that it can resolve the velocity, pressure
and mass fractions all the way down to the wall. For the chemical species, the
governing transport equations include convection, diffusion and, optionally, migration
in an electric field. Turbulent reactions are modeled using the eddy dissipation model.
There are two study types available for this physics interface. For the Time Dependent
study type the rotation is achieved by moving mesh functionality, also known as sliding
mesh. For the Frozen Rotor study type the rotating parts are kept frozen in position,
and the rotation is accounted for by the inclusion of centrifugal and Coriolis forces.
When adding the multiphysics interface, the Turbulent Flow, k-ω and Transport of
Concentrated Species interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Reacting Flow. The following default nodes
are also added in the Model Builder under Turbulent Flow, k-ω — Fluid Properties, Wall,
and Initial Values; and under Transport of Concentrated Species — Transport Properties,
No Flux, and Initial Values.
See The Rotating Machinery, Reacting Flow, Laminar Flow Interface for details.
THE ROTATING MACHINERY, REACTING FLOW, LAMINAR FLOW AND TURBULENT FLOW INTERFACES | 129
The Rotating Machinery, Reacting Flow, Turbulent Flow, SST
Interface
The Turbulent Flow, SST version of the Rotating Machinery, Reacting Flow interface ( ),
found under the Rotating Machinery, Reacting Flow>Turbulent Flow branch ( ) when
adding a physics interface, is used to simulate the turbulent flow and mixing of
chemical species in equipment containing one or more rotating parts. Reactions
between the participating species can also be studied.
The physics interface combines the functionality of the Rotating Machinery, Turbulent
Flow, SST and the Transport of Concentrated Species interface. The Reacting Flow
multiphysics coupling, which is added automatically, couples fluid flow and mass
transport. The averaged velocity and pressure fields are solved together with an
arbitrary number of averaged mass fractions. The momentum balance is governed by
the averaged Navier-Stokes equations, and the total mass conservation is governed by
the continuity equation. The fluid-flow turbulence is modeled using the SST
two-equation model with realizability constraints. The SST model is a so-called
low-Reynolds number model, which means that it can resolve the velocity, pressure
and mass fractions all the way down to the wall. For this reason the physics interface
can be used to study mass transfer at high Schmidt numbers. The SST model depends
on the distance to the closest wall, and the physics interface therefore includes a wall
distance equation. For the chemical species, the governing transport equations include
convection, diffusion and, optionally, migration in an electric field. Turbulent reactions
are modeled using the eddy dissipation model.
There are two study types available for this physics interface. For the Time Dependent
with Initialization study type the rotation is achieved by moving mesh functionality, also
known as sliding mesh. For the Frozen Rotor with Initialization study type the rotating
parts are kept frozen in position, and the rotation is accounted for by the inclusion of
centrifugal and Coriolis forces. In both study types initialization is performed by first
solving for the distance to the closest wall.
When adding the multiphysics interface, the Turbulent Flow, SST and Transport of
Concentrated Species interfaces are added to the Model Builder. A Moving Mesh interface
with a Rotating Domain node is added automatically in the Model Builder under the
Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Reacting Flow. The following default nodes
are also added in the Model Builder under Turbulent Flow, SST — Fluid Properties, Wall,
and Initial Values; and under Transport of Concentrated Species — Transport Properties,
No Flux, and Initial Values.
The physics interface combines the functionality of the Rotating Machinery, Turbulent
Flow, Low Re k-e and the Transport of Concentrated Species interfaces. It solves for the
averaged velocity and pressure fields, together with an arbitrary number of averaged
mass fractions. The momentum balance is governed by the averaged Navier-Stokes
equations, and the total mass conservation is governed by the continuity equation. The
fluid-flow turbulence is modeled using the AKN low-Reynolds number k-ε model.
The low-Reynolds number model can resolve the velocity, pressure and mass fractions
all the way down to the wall. For this reason the physics interface is suited for studying
mass transfer at high Schmidt numbers. The AKN model depends on the distance to
the closest wall, and the physics interface therefore includes a wall distance equation.
For the chemical species, the governing transport equations include convection,
diffusion and, optionally, migration in an electric field. Turbulent reactions are
modeled using the eddy dissipation model.
There are two study types available for this physics interface. For the Time Dependent
with Initialization study type the rotation is achieved by moving mesh functionality, also
known as sliding mesh. For the Frozen Rotor with Initialization study type the rotating
parts are kept frozen in position, and the rotation is accounted for by the inclusion of
centrifugal and Coriolis forces. In both study types initialization is performed by first
solving for the distance to the closest wall.
THE ROTATING MACHINERY, REACTING FLOW, LAMINAR FLOW AND TURBULENT FLOW INTERFACES | 131
When adding the multiphysics interface, the Turbulent Flow, Low Re k-ε and Transport
of Concentrated Species interfaces are added to the Model Builder. A Moving Mesh
interface with a Rotating Domain node is added automatically in the Model Builder under
the Definitions node. In addition, the Multiphysics node is added, which automatically
includes the multiphysics coupling feature Reacting Flow. The following default nodes
are also added in the Model Builder under Turbulent Flow, Low Re k-ε — Fluid Properties,
Wall, and Initial Values; and under Transport of Concentrated Species — Transport
Properties, No Flux, and Initial Values.
See The Rotating Machinery, Reacting Flow, Laminar Flow Interface for details.
LAMINAR FLOW
The available physics features for The Rotating Machinery, Laminar Flow Interface are
listed in the section Domain, Boundary, Point, and Pair Nodes for the Rotating
Machinery Interfaces.
T U R B U L E N T F L O W , k-ω
The available physics features for The Rotating Machinery, Turbulent Flow, k-ω
Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
T U R B U L E N T F L O W , L O W R E k-ε
The available physics features for The Rotating Machinery, Turbulent Flow, Low Re
k-ε Interface are listed in the section Domain, Boundary, Point, and Pair Nodes for the
Rotating Machinery Interfaces.
THE ROTATING MACHINERY, REACTING FLOW, LAMINAR FLOW AND TURBULENT FLOW INTERFACES | 133
134 | CHAPTER 8: ROTATING MACHINERY, REACTING FLOW INTERFACES
I n d e x
A Application Libraries window 19 R rotating machinery, laminar flow (rmspf)
application library examples interface 70
rmspf interfaces 71 rotating machinery, mixture model, alge-
braic yplus interface 94
B boundary nodes
rotating machinery, mixture model,
rmmf interfaces 103
k-eps interface 96
rmspf interfaces 84
rotating machinery, mixture model,
C CFL number
k-omega interface 98
pseudo time stepping, and 72, 93
rotating machinery, mixture model, lam-
contact angle (node) 88, 106
inar flow interface 92
D documentation 18 rotating machinery, mixture model, Low
domain nodes Re k-eps interface 100
rmmf interfaces 103 rotating machinery, mixture model, l-vel
rmspf interfaces 84 interface 95
rotating machinery, mixture model, real-
E emailing COMSOL 20
izable k-eps interface 97
F free surface (node) 86, 104
rotating machinery, mixture model,
free surface domain 62
Spalart-Allmaras interface 101
G geometry frames 62 rotating machinery, mixture model, SST
interface 99
I internet resources 17
rotating machinery, mixture model, tur-
K knowledge base, COMSOL 21
bulent flow interface 96
L laminar flow (rmnitf) interface 111 rotating machinery, mixture model, v2-f
laminar flow (rmrf) interface 126 interface 102
local CFL number 72, 93 rotating machinery, turbulent flow (rm-
spf) interface 72–75, 77–80, 82
M material frames 62
rotating shaft (node) 89, 106
mesh frames 62
moving mesh 62 S spatial frames 62
MPH-files 19
T technical support, COMSOL 20
P pair nodes turbulent flow k-omega (rmrf) interface
rmmf interfaces 103 128
rmspf interfaces 84 turbulent flow low re k-epsilon (rmnitf)
point nodes interface 119–121
rmmf interfaces 103 turbulent flow low re k-epsilon (rmrf) in-
rmspf interfaces 84 terface 130–131
INDEX| 135
turbulent flow sst interface 118
turbulent flow, algebraic yplus interface
112
turbulent flow, k-epsilon (rmnitf) inter-
face 115
turbulent flow, k-epsilon (rmrf) interface
127
turbulent flow, k-omega (rmnitf) inter-
face 117
turbulent flow, l-vel interface 113
turbulent flow, realizable k-epsilon inter-
face 116
W websites, COMSOL 20
136 | I N D E X