11
11
11
Laval Nozzle
Introduction
Modeling wet steam is very important in the analysis and design of steam turbines. Wet
steam refers to the two-phase mixture of saturated vapor and fine liquid droplets formed
during rapid condensation of steam.
The purpose of this tutorial is to demonstrate the mechanics of running a non-equilibrium
steam condensation problem using the wet steam model in FLUENT.
This tutorial demonstrates how to do the following:
• Set up a non-equilibrium steam condensation problem using the wet steam model.
• Compare the results obtained for the wet and dry (non-condensing) cases.
Prerequisites
This tutorial assumes that you are familiar with the FLUENT interface and have completed
Tutorial 1 from the FLUENT 6.3 Tutorial Guide. You should also be familiar with the wet
steam model. Refer to Section 23.13: Setting Up the Wet Steam Model in the FLUENT 6.3
User’s Guide for more information.
Problem Description
This problem considers steam condensation in a Laval converging-diverging nozzle. A
schematic of the problem is shown in Figure 1. The nozzle geometry is two-dimensional.
The throat height is 1 cm. Because of the symmetry of the nozzle geometry, only half of
the nozzle is modeled.
c Fluent Inc. September 27, 2006 1
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Axis of symmetry
Inlet Outlet
The fluid entering the nozzle is dry steam. As the flow accelerates and clears the throat
area, rapid condensation occurs. The steam first subcools and then nucleates, forming a
two-phase mixture of steam (gas phase) and fine water droplets (liquid phase). The wet
steam model will solve transport equations for conservation of mass fraction of the condensed
phase and conservation of number density of the droplets per unit volume. These equations
will not be solved for the non-condensing case.
Preparation
2
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Grid
FLUENT 6.3 (2d, dp, pbns, lam)
c Fluent Inc. September 27, 2006 3
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Extra: You can use the right mouse button to check the zone number corresponding
to each boundary. If you click the right mouse button on one of the boundaries
in the graphics window, its zone number, name, and type will be printed in the
FLUENT console. This feature is especially useful when you have several zones
of the same type and you want to distinguish between them.
Step 2: Units
4
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Step 3: Models
The wet steam model can only be used with the density-based solver.
!
(b) Retain the default selection of Implicit in the Formulation list.
(c) Click OK to close the Solver panel.
c Fluent Inc. September 27, 2006 5
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
6
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
c Fluent Inc. September 27, 2006 7
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
(a) Enter 530.5 torr for both, Gauge Total Pressure and Supersonic/Initial Gauge
Pressure.
(b) Retain the selection of Normal to Boundary in the Direction Specification Method
drop-down list.
(c) Retain the default setting of 0 for both, Liquid Mass Fraction and Log 10 (Droplets
Per Unit Volume).
(d) Click the Thermal tab and enter 377 K for Temperature.
(e) Click OK to close the Pressure Inlet panel.
8
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
c Fluent Inc. September 27, 2006 9
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Step 6: Solution
10
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
c Fluent Inc. September 27, 2006 11
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
12
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Residuals
continuity
x-velocity 1e+06
y-velocity
energy
wsb
wsn 1e+04
1e+02
1e+00
1e-02
1e-04
1e-06
0 100 200 300 400 500 600 700 800 900
Iterations
Scaled Residuals
FLUENT 6.3 (2d, dp, dbns imp)
c Fluent Inc. September 27, 2006 13
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Step 7: Postprocessing
14
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
5.30e+02
5.10e+02
4.90e+02
4.70e+02
4.50e+02
4.30e+02
4.10e+02
3.90e+02
3.70e+02
3.50e+02
3.30e+02
3.10e+02
2.90e+02
2.70e+02
2.49e+02
2.29e+02
2.09e+02
1.89e+02
1.69e+02
1.49e+02
1.29e+02
3.77e+02
3.73e+02
3.70e+02
3.66e+02
3.62e+02
3.59e+02
3.55e+02
3.51e+02
3.48e+02
3.44e+02
3.40e+02
3.37e+02
3.33e+02
3.29e+02
3.25e+02
3.22e+02
3.18e+02
3.14e+02
3.11e+02
3.07e+02
3.03e+02
c Fluent Inc. September 27, 2006 15
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
1.53e+00
1.45e+00
1.38e+00
1.31e+00
1.23e+00
1.16e+00
1.08e+00
1.01e+00
9.35e-01
8.61e-01
7.87e-01
7.13e-01
6.39e-01
5.64e-01
4.90e-01
4.16e-01
3.42e-01
2.68e-01
1.94e-01
1.20e-01
4.56e-02
5.13e-02
4.88e-02
4.62e-02
4.36e-02
4.11e-02
3.85e-02
3.59e-02
3.34e-02
3.08e-02
2.82e-02
2.57e-02
2.31e-02
2.05e-02
1.80e-02
1.54e-02
1.28e-02
1.03e-02
7.70e-03
5.13e-03
2.57e-03
0.00e+00
16
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
2.28e+01
2.17e+01
2.05e+01
1.94e+01
1.82e+01
1.71e+01
1.60e+01
1.48e+01
1.37e+01
1.25e+01
1.14e+01
1.03e+01
9.12e+00
7.98e+00
6.84e+00
5.70e+00
4.56e+00
3.42e+00
2.28e+00
1.14e+00
0.00e+00
1.53e-02
1.45e-02
1.38e-02
1.30e-02
1.22e-02
1.15e-02
1.07e-02
9.94e-03
9.17e-03
8.41e-03
7.64e-03
6.88e-03
6.12e-03
5.35e-03
4.59e-03
3.82e-03
3.06e-03
2.29e-03
1.53e-03
7.64e-04
0.00e+00
c Fluent Inc. September 27, 2006 17
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
7. Define a custom field function for the ratio of static pressure to total inlet pressure.
Define −→Custom Field Functions...
(a) Select Pressure... and Static Pressure from the Field Functions drop-down lists
and click the Select button to add pressure (p) to the Definition field.
(b) Click the “/” button to add the division symbol to the Definition field.
(c) Enter 70727.32 (the SI value of total inlet pressure = 530.5 torr i.e., 70727.32 Pa)
using the calculator buttons.
You can use the DEL button to delete characters one at a time, if required.
(d) Enter p-ratio2 for New Function Name.
(e) Click Define and close the Custom Field Function Calculator panel.
8. Compare the results obtained for the simulation with available experimental data.
Plot −→XY Plot...
(a) Select Custom Field Functions... and p-ratio2 from the Y Axis Function drop-down
lists.
Custom Field Functions... is at the top of the upper Y Axis Function drop-down
list. Scroll up the list to select Custom Field Functions....
(b) Select symmetry.2 from the Surfaces selection list.
18
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
(c) Click the Load File... button and select the experimental data file EXP-410-ppo.xy.
The plot title of the experimental data file (Exp-410 DATA P/Po) will be displayed
in the File Data list.
(d) Click Plot (Figure 10).
1.0
FLUENT
0.9
Exp-410
0.8
0.7
p-ratio2 0.6
0.5
0.4
0.3
0.2
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16
Position (m)
p-ratio2
FLUENT 6.3 (2d, dp, dbns imp)
c Fluent Inc. September 27, 2006 19
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
iii. Enter moses-stein-ppo.xy for XY File and click OK to close the Select File
dialog box.
(f) Close the Solution XY Plot panel.
20
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
1.78e+00
1.69e+00
1.61e+00
1.52e+00
1.43e+00
1.35e+00
1.26e+00
1.17e+00
1.09e+00
9.99e-01
9.13e-01
8.26e-01
7.39e-01
6.53e-01
5.66e-01
4.79e-01
3.92e-01
3.06e-01
2.19e-01
1.32e-01
4.56e-02
c Fluent Inc. September 27, 2006 21
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
2. Define a custom field function for the ratio of static pressure to total inlet pressure.
Define −→Custom Field Functions...
(a) Enter p/70727.32 in the Definition field.
Select Pressure and Static Pressure from the Field Functions drop-down lists and
click the Select button to add pressure (p) to the Definition field.
(b) Enter p-ratio2 for New Function Name.
(c) Click Define and close the Custom Field Function Calculator panel.
3. Compare the results obtained for the condensing and non-condensing cases.
Plot −→XY Plot...
(a) Select Custom Field Functions... and p-ratio2 from the Y Axis Function drop-down
lists.
Custom Field Functions... is at the top of the upper Y Axis Function drop-down
list. Scroll up the list to select Custom Field Functions....
(b) Select symmetry.2 from the Surfaces selection list.
(c) Click the Load File... button and select the data file moses-stein-ppo.xy.
Click the Free Data button to remove previously loaded data files from the File
Data selection list before loading the new file.
(d) Click Plot (Figure 12).
1.0
0.7
0.6
p-ratio2
0.5
0.4
0.3
0.2
0.1
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16
Position (m)
p-ratio2
FLUENT 6.3 (2d, dp, dbns imp)
The jump in pressure for the condensing case is due to the initiation of droplet
nucleation, which is absent in the non-condensing (isentropic expansion) case.
22
c Fluent Inc. September 27, 2006
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
Summary
This tutorial demonstrated the set up and solution for a non-equilibrium steam condensa-
tion problem using the wet steam model in FLUENT. The results for the steam condensation
case were compared with available experimental data. The problem was set up to exclude
calculation of the wet steam equations (i.e., for non-condensation) and results for the con-
densing and non-condensing cases were compared.
c Fluent Inc. September 27, 2006 23