Pipe Flow Module: Introduction To
Pipe Flow Module: Introduction To
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:
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
The Applications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
The Pipe Flow Module Interfaces . . . . . . . . . . . . . . . . . . . . . . . . . 9
Physics Interface Guide by Space Dimension and Study Type . . . 10
Tutorial Example: Geothermal Heating from a Pond Loop . . . 12
|3
4|
Introduction
The Pipe Flow Module is an optional add-on package for COMSOL Multiphysics
designed to model and simulate fluid flow, heat, and mass transfer in pipes and
channels. Compressible hydraulic transients and acoustic waves can also be
modeled using the Water Hammer interface and Pipe Acoustics interface,
respectively. The Pipe Flow Module can address problems involving flow velocity,
pressure, temperature and sound waves in pipes and channels.
Modeling pipes as curves in 2D or 3D gives a great advantage in computational
efficiency over meshing and computing 3D pipes with finite diameter.
Figure 1: The Pipe Flow interface reduces the 3D flow problem to a 2D or 3D curve.
Pipe systems for which the ratio length/diameter is large enough that the flow
inside of each pipe segment can be considered fully developed are suitable for the
Pipe Flow Module.
The physics interfaces in the module define the conservation of momentum,
energy, and mass of a fluid inside a pipe or channel. The flow rate, pressure,
temperature, and concentration fields are modeled as cross-section averaged
quantities, so that they only vary along the length of the pipes. The pressure losses
along the length of a pipe or in a pipe component are described using friction
factors. A broad range of built-in expressions for Darcy and Fanning friction
factors cover the entire flow regime from laminar to turbulent flow, Newtonian
and non-Newtonian fluids, different cross-sectional geometries, and a wide range
of relative surface roughness values. In addition to the continuous frictional
pressure drop along pipe stretches, pressure drops due to irreversible losses in
components such as bends, contractions, expansions, T-junctions, Y-junctions,
and valves are computed through an extensive library of industry standard loss
coefficients. Pumps are also available as flow inducing devices.
The features in this module are intended for modeling and simulating
incompressible and weakly compressible fluid flow in pipes and channel systems,
as well as compressible hydraulic transients and acoustic waves. Typical simulations
yield the velocity, pressure variation, and temperature in systems of pipes and
channels. Hydraulic transients are also possible to model. These can be the result
of a valve that is closed rapidly in a pipe network, which is referred to as a water
hammer.
|5
The Applications
The module can be used to design and optimize complex cooling systems in
turbines, analyze ventilation systems in buildings, pipe systems in the chemical
process industry, and pipelines in the oil and gas industry, just to mention a few
applications.
Figure 2: A Probe Tube Microphone modeled using the “Pipe Acoustics, Transient” interface in the Pipe
Flow Module.
Any devices in which you find flow, waves, mass or heat transfer phenomena in
narrow pipes or ducts are candidates for simulation with the Pipe Flow Module.
Classical pressure-drop and mass flow calculations through piping with bends,
valves, tanks and so forth are well suited for the Pipe Flow Module.
Figure 3: Tutorial example describing the discharging of a water tank through a simple pipe system with
bends and valves.
For heat transfer studies, the Pipe Flow Module includes several automatic
couplings to the surrounding heat sinks or sources; both shortcut methods with
6|
semi-empirical correlations for forced and convection, and also direct coupling to
a 3D solid, in which the pipe is embedded.
Figure 4: Methods for heat transfer to the surroundings, from left to right: forced convection, convection,
solid conduction.
Figure 5: The heat transfer to the cooling channels embedded in a cooling mold is simulated to
understand the controlled cooling of a polyurethane steering wheel. The Nonisothermal Pipe Flow
interface is used in the model.
|7
chemical reaction modeling. This can include mass transfer, chemical kinetics, heat
transfer, and pressure drop calculations in the same model.
Figure 6: Temperature distribution in an autothermal chemical reactor. The model includes mass
transport, chemical kinetics, heat transfer, and pressure drop and flow calculations.
Figure 7: The Pipeline Insulation model in the application library simulates the effect of viscous heating
on the cooling and transport properties of oil in a pipeline.
8|
The Pipe Flow Module Interfaces
The module includes the following physics interfaces. See also the graphical view
from the Model Wizard, shown in Figure 8.
P IPE F LOW
The Pipe Flow interface ( ) computes the pressure and velocity field in
isothermal pipe systems.
W ATER H AMMER
The Water Hammer interface ( ) solves rapid hydraulic transients in pipe
systems, taking the elastic properties of both the fluid and pipe wall into account.
|9
P IPE A COUSTICS , T RANSIENT
The Pipe Acoustics, Transient interface ( ) models sound waves in flexible pipe
systems, for arbitrary transient variations in pressure.
Figure 8: The physics interfaces for the Pipe Flow Module as shown in the Model Wizard. The interfaces
are located in different groups.
10 |
PHYSICS INTERFACE ICON TAG SPACE AVAILABLE STUDY TYPE
DIMENSION
Acoustics
Acoustic-Structure Interaction
Single-Phase Flow
Nonisothermal Flow
1
Requires both the Pipe Flow Module and the Acoustics Module.
| 11
Tutorial Example: Geothermal Heating from a Pond Loop
Ponds and lakes can serve as thermal reservoirs in geothermal heating applications.
In this example, fluid circulates underwater through polyethylene piping in a
closed system. Inside the pipes, cold water is flowing and picking up low grade
heat from the pond. The return water will have a slightly higher temperature than
the feed. A heat pump (not included in the model) will upgrade the heat for use
in dwellings. Chilled water is then returned back to the pond loop. The pipes are
coiled in a slinky shape and mounted onto sleds. The Nonisothermal Pipe Flow
interface sets up and solves the equations for the temperature and fluid flow in the
pipe system, where the geometry is represented by lines in 3D. A forced
convection correlation is used (see Figure 4 above) for the external cooling effect.
Figure 9: A sled carrying pipe coils shown before the system is submerged.
Problem Specification
High density polyethylene pipe (20 mm diameter) is rolled into sixteen coils.
Groups of eight coils are mounted on two sleds. Each coil has a radius of 1 m and
a length of approximately 75 m. The coil groups are connected to feed and return
piping with a diameter of 50 mm (see Figure 10). The coil groups are 2.4 m in
12 |
height and sit at the bottom of a pond that is 6 m deep. The total length of the
piping is 1446 m.
9m 0m
6m
9m
2.4 m
0m
2.4 m
return
feed
2.4 m
2m
Figure 10: Polyethylene pipe system. Elevation above the pond bottom is indicated. Feed and return
piping (gray) is 50 mm in diameter while coils (black) are 20 mm in diameter. The pipes are insulated
above the pond surface.
The heat exchange between pond water and pipe fluid will depend, among other
things, on the temperature difference between the two. A slow current in the pond
will make the heat transfer more effective than water at rest. The pond is warmer
closer to the surface, as shown by the temperature data in the table below.
POND TEMPERATURE
ELEVATION (M) TEMPERATURE (K)
0 284
2 288
4 291
6 293
| 13
It is easy to set up a function in the software with linear interpolation between
points so that the varying pond temperature can be taken into account in the
simulation.
F LOW E QUATIONS
The stationary momentum and continuity equations below describe the stationary
flow inside the pipe system, and are solved by the Nonisothermal Pipe Flow
interface:
ρ
0 = – ∇p – f D ---------- u u + F (1)
2d h
∇ ⋅ ( Aρu ) = 0 (2)
Above, A (SI unit: m2) is the cross section area of the pipe, ρ (SI unit: kg/m3) is
the density, u (SI unit: m/s) is the fluid velocity in the tangential direction of the
pipe curve segment, and p (SI unit: N/m2) is the pressure. F (SI unit: N/m3) is
a volume force, like gravity.
Gravity can be included explicitly in the model, but since the variation in density
is negligible, and the model is not pressure driven, the only effect of including
gravity is a change in the total pressure level. It is therefore common modeling
practice to exclude gravity from F and interpret the pressure variable as the
reduced pressure p r = p – ρg ( z 0 – z ) , where z0 is the datum level of the free liquid
surface. This reduces the model complexity and yields the same results. For models
where you want to include gravity effects you can add a Volume Force node.
14 |
hydraulic pipe diameter, e/dh. The Nonisothermal Pipe Flow interface provides a
library of built-in expressions for the Darcy friction factor, fD.
Figure 11: Select a shape and different predefined friction models in the Pipe Properties node.
This example uses the Churchill relation (Ref. 1) that is valid for laminar flow,
turbulent flow, and the transitional region in between. The Churchill relation is:
8 12 1 ⁄ 12
f D = 8 ------- + ( A + B )
-1.5
(3)
Re
where
0.9 16
A = -2.457ln ------- + 0.27 ( e ⁄ d )
7
(4)
Re
16
B = ----------------
37530
(5)
Re
| 15
As seen from the equations above, the friction factor is a function of the surface
roughness divided by diameter of the pipe. Surface roughness data can be selected
from a predefined list in the Pipe Properties feature.
Figure 12: Specify the surface roughness in the Pipe Properties node.
The Churchill equation is also a function of the fluid properties, through the
Reynolds number:
ρud
Re = ----------- (6)
μ
where Cp (SI unit: J/(kg·K)) is the heat capacity at constant pressure, T is the
temperature (SI unit: K), and k (SI unit: W/(m·K)) is the thermal conductivity.
The second term on the right-hand side of Equation 7 corresponds to friction heat
dissipated due to viscous shear. Qwall (SI unit: W/m) is a source/sink term due to
heat exchange with the surroundings through the pipe wall:
16 |
Q wall = hZ ( T ext – T ) (8)
where Z (m) is the wetted perimeter of the pipe, h (W/(m2·K)) is an overall heat
transfer coefficient and Text (K) is the external temperature outside of the pipe.
If there are several wall layers, the overall heat transfer coefficient includes
contributions from internal film resistance, layered wall resistance, and external
film resistance. The heat transfer coefficient is automatically calculated taking all
of these resistances into account.
In this model we only have one wall and the properties of the pipe wall are given
in the table below.
Results
Figure 13 shows the pressure (Pa) in the 1446 m pipe system assuming that water
enters the system at a rate of 4 l/s.
Figure 13: Pressure drop over the pipe system due to flow losses.
| 17
The plot below shows the temperature (K) distribution for the pipe fluid. It enters
the pipe system at 5 °C and exits with a temperature of approximately 11 °C.
Turbulent flow conditions in the loop are important for good heat exchange
between the pipes and the surroundings. A plot of the Reynolds number is shown
in Figure 15, confirming that flow is turbulent (Re > 3000) throughout the
system.
Figure 15: The Reynolds number in the pipe loop confirms that the flow conditions are turbulent.
18 |
References
1. S.W. Churchill, “Friction factor equations span all fluid-flow regimes,” Chem. Eng.,
vol. 84, no. 24, p. 91, 1997.
2. F.P. Incropera and D.P. DeWitt, Fundamentals of Heat and Mass Transfer, 4th ed.,
John Wiley & Sons, 1996. Eq 8.62 and Eq 9.34, respectively.
Model Wizard
The first step to build a model is to open COMSOL and then select the physics
interface and specify the type of analysis you want to do—in this case, a stationary,
Nonisothermal Flow analysis.
Note: These instructions are for the user interface on Windows but apply, with
minor differences, also to Linux and Mac.
1 To start the software, double-click the COMSOL icon on the desktop. When
the software opens, you can choose to use the Model Wizard to create a new
COMSOL model or Blank Model to create one manually. For this tutorial, click
the Model Wizard button.
If COMSOL is already open, you can start the Model Wizard by selecting
New from the File menu and then click Model Wizard .
The Model Wizard guides you through the first steps of setting up a model. The
next window lets you select the dimension of the modeling space.
2 In the Select Space Dimension window click 3D .
3 In the Select Physics tree under
Fluid Flow>Nonisothermal Flow, click
Nonisothermal Pipe Flow (nipfl) .
4 Click Add and then click Study .
5 In the Select Study window under
General Studies, click Stationary .
6 Click Done .
| 19
Geometry
Note: The location of the files used in this exercise varies based on your
installation. For example, if the installation is on your hard drive, the file path is
similar to C:\Program Files\COMSOL53a\Multiphysics\applications\.
Start by creating the pipe system geometry. You can simplify this by inserting a
prepared geometry sequence from a file.
1 On the Geometry toolbar click Insert Sequence .
Definitions
Now add some data in the form of interpolation tables and variables.
Interpolation 1
The water in the pond is somewhat stratified and has a temperature gradient along
the depth direction. This has been measured earlier and will be added as an
interpolation table.
1 On the Definitions toolbar click Interpolation .
Note: On Linux and Mac, the Home toolbar refers to the specific set of controls
near the top of the Desktop.
20 |
2 In the Settings window for Interpolation
under Definition, enter the following in
the table:
T F(T)
0 284
2 288
4 291
6 293
Variables 1
1 On the Home toolbar click Variables and choose Local Variables .
2 In the Settings window for Variables locate the Variables section. Click Load
from File .
3 Browse to the application library folder and double-click the file
geothermal_heating_variables.txt, under
Pipe_Flow_Module\Heat_Transfer.
This should be the result of the file import for variable definitions. Note that the expression int1(z) refers
to the interpolation table that was imported in the previous step, and z is the name of the height
coordinate.
| 21
Mate rials
Now is the time to specify the physics feature settings, which include the flow and
temperature conditions, as well as the pressure drop correlations to use in the pipe
segments. Since piping systems typically consist of a large number of different
segments, it can be efficient to define selections to get a manageable model.
Pipe Properties 1
1 In the Model Builder under Nonisothermal Pipe Flow, click Pipe Properties 1.
22 |
Note: In many of the Settings windows in COMSOL, there is an Equation section
that you can expand. This shows you the equations employed by the current
setting. In this case, the correlations for the default friction factor model
Churchill is displayed. Also, if you are interested in the underlying theory of any
Settings window, press F1 or the Help button . This takes you to the
documentation section describing the current Settings window.
2 In the Settings window for Pipe Properties locate the Pipe Shape section.
- From the Shape list, choose Circular.
- In the di text field, type 20[mm].
Temperature 1
1 In the Model Builder under Nonisothermal Pipe Flow, click Temperature 1.
2 In the Settings window for Temperature locate the Temperature section. In the
Tin text field, type 5[degC].
Pipe Properties 2
Now you can assign additional different properties to a new group of pipe
segments.
1 Right-click Nonisothermal Pipe Flow and choose Pipe Properties.
| 23
2 Select the edges 1–12, 15, 17–21, 27–29, 33–35, 39, 41–45, 51–53, 57–59,
63, 65–69, 75–77, 81–83, 87, 89–91, 97, 98, 102, and 103.
This selection represents the feed and return pipes from the coils. It can be
somewhat tedious to select them only by looking at the numbered list above.
However, since we know visually which part this is, you can select it graphically.
- In the Graphics window, click Go to ZX View and then Select Box .
- In the Graphics window, draw a box around the pipes highlighted in blue in
the picture above.
- Click Go to Default 3D View button.
3 In the Settings window for Pipe Properties locate the Pipe Shape section.
- From the list, choose Circular.
- In the di text field, type 50[mm].
24 |
2 Select Edges 7–104 only, using the Select Box tool as before. The parts
corresponding to these edges are shown in blue color below.
3 In the Settings window for Wall Heat Transfer locate the Heat Transfer Model
section. In the Text text field, type T_pond.
Wall Layer 1
Now the polyethene pipe is defined by adding a wall layer to the Wall Heat
Transfer feature.
1 In the Model Builder under Nonisothermal Pipe Flow, right-click Wall Heat
Transfer 1 and choose Wall Layer .
2 In the Settings window for Wall Layer locate the Specification section.
- From the k list, choose User defined.
- In the k text field, type k_wall.
- From the Δw list, choose User defined.
- In the text field, type d_wall.
More than one Wall Layer can be added if you want a multilayered pipe wall.
| 25
External Film Resistance
1 In the Model Builder right-click Wall Heat
Transfer 1 and choose External Film
Resistance .
Inlet 1
1 On the Physics toolbar click Points and choose Inlet .
2 Select Point 1 only.
3 In the Settings window for Inlet locate the Inlet Specification section.
- From the Specification list, choose Volumetric flow rate.
- In the qv,0 text field, type 4[l/s] (as in liters/second).
Heat Outflow 1
1 On the Physics toolbar click Points and choose Heat Outflow .
2 Select Point 2 only.
26 |
Me sh
Solve
Pressure
The default plot groups show the pressure (Figure 13), velocity, and temperature
(Figure 14) in the pipe system. To get a better view, do as follows:
1 Click the Zoom Box button on the Graphics toolbar.
| 27
2 Draw a box in the Graphics window to zoom in on the two coil packages.
Reproduce the Reynolds number plot in Figure 15 with the following steps.
3D Plot Group 4
1 On the Results toolbar click 3D Plot Group .
2 When the 3D Plot Group 4 is selected in the Toolbar, an additional toolbar
containing Plot Tools for the 3D Plot Group displays. On the 3D Plot Group 4
toolbar, click Line .
28 |
Note: In this model the bend features in the pipe corners are deliberately left out
as it is assumed these do not affect the results much. As an exercise, you can add
bend point features and see if this assumption is correct.
| 29
30 |