Valve Flow Analysis
Valve Flow: Theoretical Foundations
Overview
Teacher! What is the purpose of performing flow analysis on valves?
CFD analysis of valves aims to predict flow coefficients (Cv/Kv values), evaluate cavitation characteristics, and identify pressure recovery downstream of the valve and noise sources. Each valve type, such as butterfly valves, ball valves, gate valves, and globe valves, has its own unique flow phenomena.
Governing Equations
Please tell me the basic flow characteristic equation for valves.
The flow characteristic of a valve is expressed by the flow coefficient $C_v$ (US system) or $K_v$ (European system).
$Q$ is the flow rate [US GPM], $SG$ is the specific gravity (water=1), and $\Delta p$ is the differential pressure [psi]. The relationship with $K_v$ in the SI unit system is:
Definition of $K_v$: At $\Delta p = 1$ bar and water temperature of 15°C, $K_v$ [m³/h] of water flows.
How do you determine Cv using CFD?
It is back-calculated from the pressure loss $\Delta p$ between the inlet and outlet and the flow rate $Q$ obtained from CFD using the above equation. There are evaluation methods conforming to the ISA/IEC 60534 standard.
Cavitation
How do you evaluate cavitation?
It is evaluated using the cavitation index $\sigma$.
$p_1$ is the upstream pressure, $p_2$ is the downstream pressure, and $p_v$ is the saturation vapor pressure. Cavitation begins when $\sigma$ falls below the critical value $\sigma_i$ (Incipient Cavitation Index).
| Valve Type | Typical $\sigma_i$ |
|---|---|
| Butterfly Valve (fully open) | 0.2~0.5 |
| Ball Valve (fully open) | 0.15~0.3 |
| Gate Valve (fully open) | 0.15~0.25 |
| Globe Valve | 0.5~1.5 |
The large σ for globe valves is because the pressure recovery at the valve body is small, right?
Exactly. Globe valves have a curved flow path, resulting in small pressure recovery, and the difference between the minimum pressure at the Vena Contracta (contracted flow section) and the downstream pressure is small. Consequently, cavitation is less likely to occur (larger $\sigma_i$) for the same $\Delta p$.
Pressure Recovery Coefficient
The pressure recovery coefficient $F_L$ is a valve-specific coefficient defined in IEC 60534.
$p_{vc}$ is the pressure at the Vena Contracta, right? Since CFD can directly read the pressure at the Vena Contracta, $F_L$ can be accurately determined.
Yes. In experiments, it is measured with downstream pressure taps, making it difficult to accurately identify the location of the Vena Contracta. With CFD, the point of minimum pressure within the flow path can be directly visualized.
Practical Considerations
- Upstream/Downstream Straight Pipe Sections: ISA/IEC standards recommend at least 10D upstream and 5D downstream.
- Wall Roughness: Cast valves have significant roughness (0.5~2 mm).
- Analysis for Each Opening: Evaluate not only fully open but also 25%, 50%, and 75% openings.
- Fluid Compressibility: Consider compressibility for gas valves when Mach number > 0.3.
History of Valve Flow Theory — Zhukovsky's Water Hammer Equation (1898)
The first mathematical description of water hammer theory was by the Russian hydraulic engineer Nikolai Zhukovsky. In his 1898 paper "On the Hydraulic Hammer in Water Supply Pipes," he derived the pressure rise ΔP = ρ×a×ΔV (a: pressure wave speed) due to rapid valve closure. This "Zhukovsky equation" is still used as a fundamental formula in piping design today. Interestingly, Zhukovsky also made significant contributions to aerodynamics (Zhukovsky transform for airfoil lift) during the same era—a rare researcher who established two fundamental theories in different fields of fluid mechanics. Modern CFD goes beyond the Zhukovsky equation, enabling complete analysis that includes the effects of valve body shape, pipe bends, and cavitation.
Computational Methods for Valve Flow
Details of Numerical Methods
Please tell me the specific implementation methods for valve CFD.
Mesh Strategy
The interior of a valve has a complex 3D shape with rapid changes in flow path cross-section. Mesh quality significantly affects the results.
| Region | Mesh Size | Remarks |
|---|---|---|
| Around Valve Seat | Diameter D/100~D/50 | Resolution of Vena Contracta |
| Valve Disc/Ball Surface | D/80~D/40 | Pressure distribution, Fluid Force |
| Seal Clearance (for small openings) | 1/5 of clearance or less | Minimum 5 cells |
| Upstream Straight Pipe | D/20 | Ensuring developed flow |
| Downstream Straight Pipe (separation region) | D/30~D/20 | Resolution of reattachment |
| Wall Prism Layer | y+ ≒ 1~30 | Match to turbulence model |
The clearance mesh for small valve openings seems particularly challenging.
Yes. For a butterfly valve at 10% opening, the clearance between the disc and the pipe wall is only a few millimeters. At least 5 cell layers must be ensured in this clearance. This is handled with an Inflation Layer (prism layer).
Boundary Conditions
| Boundary | Condition | Setting Value |
|---|---|---|
| Inlet | Pressure Inlet or Mass Flow | Upstream pressure or design flow rate |
| Outlet | Pressure Outlet | Downstream pressure |
| Valve Wall | No-Slip | Roughness setting (cast: 0.5~2 mm) |
| Pipe Wall | No-Slip | Roughness setting (steel pipe: 0.045 mm) |
How do you decide whether to use Pressure Inlet or Mass Flow Inlet?
For calculating Cv values, determining Δp at a constant flow rate yields higher accuracy. The combination of Mass Flow Inlet + Pressure Outlet is recommended. Conversely, to determine flow rate with fixed Δp, use Pressure Inlet + Pressure Outlet.
Turbulence Model
Valve flows involve separation, reattachment, and strong curvature effects, so SST k-omega is the most reliable.
| Valve Type | Recommended Model | Reason |
|---|---|---|
| Butterfly | SST k-omega | Separation behind disc |
| Ball | SST k-omega | Separation around sphere |
| Gate | Realizable k-epsilon | Relatively simple flow path |
| Globe | SST k-omega | Complex curved flow path |
Cavitation Model
Please tell me how to model cavitation in CFD.
The Schnerr-Sauer model or the Zwart-Gerber-Belamri model is widely used. It is combined with the VOF (Volume of Fluid) method to calculate bubble generation (evaporation) and collapse (condensation).
$C_{prod}$ and $C_{dest}$ are empirical constants, right? Are the default values okay?
Fluent's default values (Zwart: $C_{prod}=50$, $C_{dest}=0.01$, $R_B=10^{-6}$ m) produce reasonable results in many cases. However, calibration may be necessary for very high operating pressures or special fluids.
Solver Settings
| Parameter | Single-Phase Flow | Cavitation Analysis |
|---|---|---|
| Solver | Pressure-Based, Steady | Pressure-Based, Transient |
| Multiphase Flow Model | None | VOF (Mixture) |
| Pressure-Velocity Coupling | Coupled | Coupled |
| Time Step | - | CFL < 1 |
Numerical Methods for Valve Flow CFD — Incompressible, Compressible, Cavitation: Discriminating Between Three Physical Regimes
In valve flow analysis, accurately discriminating between three physical regimes based on the pressure ratio and flow velocity upstream and downstream of the valve and selecting the appropriate solver is fundamental to accuracy. ① Incompressible (Mach<0.3): Pressure-based solver is sufficient. Used for basic design of Cv values and pressure loss. ② Subsonic compressible (Mach 0.3~1.0): Density-based solver or low-Mach number correction is required. High-pressure steam valves fall into this regime. ③ Cavitation (liquid, local pressure < saturation vapor pressure): Schnerr