Valve Flow Analysis
Theory and Physics
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$.
$$ \sigma = \frac{p_2 - p_v}{p_1 - p_2} $$
$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.
$$ F_L = \sqrt{\frac{p_1 - p_2}{p_1 - p_{vc}}} $$
$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.
Coffee Break Trivia
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.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out in an unstable, spluttering manner, but after a while, it becomes a steady flow, right? This "period of change" is described by the temporal term. The pulsation of blood flow due to heartbeats, and the flow fluctuations each time an engine valve opens and closes are all unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning this term is set to zero. Since computational cost is significantly reduced, starting with a steady-state solution is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far corner of a room is also due to the "carrier" air transporting heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They are completely different! Convection is transport by flow, conduction is transmission by molecules. There is an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while they naturally mix, right? That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move in a "thick" manner. In low Reynolds number flows (slow, viscous), diffusion is dominant. Conversely, in high Re number flows, convection overwhelms and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push the plunger of a syringe, liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Flow arises where there is a pressure difference"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results become strange immediately after switching to compressible analysis, it might be due to confusion between absolute/gauge pressure.
- Source Term $S_\phi$: Warmed air rises—why? Because it becomes lighter (lower density) than its surroundings, so it is pushed upward by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat generated by a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside" and are expressed by source terms. What happens if you forget the source term? If you forget to include buoyancy in a natural convection analysis, the fluid won't move at all—a physically impossible result where warm air doesn't rise in a room with the heater on in winter.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (molecular mean free path ≪ characteristic length).
- Newtonian Fluid Assumption: Linear relationship between shear stress and strain rate (viscosity model required for non-Newtonian fluids).
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. Consider compressibility effects for Mach number ≥ 0.3.
- Boussinesq Approximation (Natural Convection): Consider density variation only in the buoyancy term, using constant density in other terms.
- Non-applicable Cases: Rarefied gas (Kn > 0.1), supersonic/hypersonic flow (shock wave capturing required), free surface flow (VOF/Level Set, etc. required).
Dimensional Analysis and Unit Systems
Variable SI Unit Notes / Conversion Memo Velocity $u$ m/s When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units. Pressure $p$ Pa Distinguish between gauge pressure and absolute pressure. Use absolute pressure for compressible analysis. Density $\rho$ kg/m³ Air: approx. 1.225 kg/m³@20°C, Water: approx. 998 kg/m³@20°C Viscosity Coefficient $\mu$ Pa·s Be careful not to confuse with kinematic viscosity coefficient $\nu = \mu/\rho$ [m²/s]. Reynolds Number $Re$ Dimensionless $Re = \rho u L / \mu$. Indicator for laminar/turbulent transition. CFL Number Dimensionless $CFL = u \Delta t / \Delta x$. Directly related to time step stability.
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$.
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.
- 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.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out in an unstable, spluttering manner, but after a while, it becomes a steady flow, right? This "period of change" is described by the temporal term. The pulsation of blood flow due to heartbeats, and the flow fluctuations each time an engine valve opens and closes are all unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning this term is set to zero. Since computational cost is significantly reduced, starting with a steady-state solution is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far corner of a room is also due to the "carrier" air transporting heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They are completely different! Convection is transport by flow, conduction is transmission by molecules. There is an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while they naturally mix, right? That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move in a "thick" manner. In low Reynolds number flows (slow, viscous), diffusion is dominant. Conversely, in high Re number flows, convection overwhelms and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push the plunger of a syringe, liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Flow arises where there is a pressure difference"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results become strange immediately after switching to compressible analysis, it might be due to confusion between absolute/gauge pressure.
- Source Term $S_\phi$: Warmed air rises—why? Because it becomes lighter (lower density) than its surroundings, so it is pushed upward by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat generated by a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside" and are expressed by source terms. What happens if you forget the source term? If you forget to include buoyancy in a natural convection analysis, the fluid won't move at all—a physically impossible result where warm air doesn't rise in a room with the heater on in winter.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (molecular mean free path ≪ characteristic length).
- Newtonian Fluid Assumption: Linear relationship between shear stress and strain rate (viscosity model required for non-Newtonian fluids).
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. Consider compressibility effects for Mach number ≥ 0.3.
- Boussinesq Approximation (Natural Convection): Consider density variation only in the buoyancy term, using constant density in other terms.
- Non-applicable Cases: Rarefied gas (Kn > 0.1), supersonic/hypersonic flow (shock wave capturing required), free surface flow (VOF/Level Set, etc. required).
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units. |
| Pressure $p$ | Pa | Distinguish between gauge pressure and absolute pressure. Use absolute pressure for compressible analysis. |
| Density $\rho$ | kg/m³ | Air: approx. 1.225 kg/m³@20°C, Water: approx. 998 kg/m³@20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity coefficient $\nu = \mu/\rho$ [m²/s]. |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Indicator for laminar/turbulent transition. |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time step stability. |
Numerical Methods and Implementation
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).
$$ \dot{m}^+ = C_{prod} \frac{3 \alpha_v \rho_v}{R_B} \sqrt{\frac{2}{3} \frac{\max(p_v - p, 0)}{\rho_l}} $$
$$ \dot{m}^- = C_{dest} \frac{3 \alpha_v \rho_v}{R_B} \sqrt{\frac{2}{3} \frac{\max(p - p_v, 0)}{\rho_l}} $$
$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
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 |
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.
| 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
Related Topics
なった
詳しく
報告