Zhuojin Pipe Fitting Co., LtdNumerical Simulation of Erosion Effects in Tee Fittings Under Asymmetric Inlet Conditions in Natural Gas Pipelines
Literature Overview
This paper by Fu Wencong, published in the Journal of Petroleum and Natural Gas Science and Engineering in 2025, presents a comprehensive numerical simulation study on erosion behavior in tee fittings under "two-in-one-out" operating conditions typical of natural gas gathering and transmission systems. The research employs a multiphase mixture model combined with a discrete phase model (DPM) to simulate gas-liquid-solid three-phase flow, coupled with the E/CRC erosion model to quantify erosion patterns. The study is particularly relevant to engineers working on natural gas pipeline integrity, as tee fittings are critical points of material loss in production systems where free water and solid particles coexist.
Core Technical Approach and Methodology
The simulation framework integrates three key components: a mixture multiphase flow model for the continuous gas-liquid phase, a discrete phase model for tracking solid particle trajectories, and the E/CRC (Erosion/Critical) erosion model for material removal prediction. This combination allows the researcher to capture the complex interactions between the continuous fluid phase, dispersed liquid droplets, and solid particles that collectively drive erosion in tee geometries.
The study systematically varies three independent parameters: inlet velocity differential between the main run and branch, inlet liquid volume fraction differential, and pipe diameter differential between the two inlet ports. Each parameter combination generates a distinct flow field topology that governs the erosion pattern and rate on the tee interior surfaces.
Key Findings and Technical Interpretation
| Parameter Variable | Effect on Flow Field | Erosion Pattern Change | Erosion Rate Trend |
|---|---|---|---|
| Main vs. branch inlet velocity differential | Alters dominant flow pattern and impingement zone | Changes erosion morphology and spatial distribution | Non-linear variation depending on which inlet dominates |
| Branch inlet liquid volume fraction increase | Increases high-velocity region on lower wall | Erosion shape transitions from elliptical to bell-shaped | Increases with higher liquid fraction |
| Branch pipe diameter increase | Modifies lower wall velocity distribution and particle dispersion | Erosion area expands | Decreases despite larger affected area (inverse relationship) |
The most striking finding is the inverse relationship between erosion area and erosion rate as branch pipe diameter increases. This occurs because a larger branch diameter distributes the flow energy over a wider impingement zone, reducing the local kinetic energy density while spreading the erosion footprint. This has direct implications for tee sizing decisions in production systems where the branch-to-run diameter ratio can be adjusted during design.
Connection to Engineering Practice
In natural gas production systems, tee fittings fabricated from carbon steel (typically X65 or X70 per API 5L) or low-temperature grades (API 5L X65 LT) are routinely subjected to erosive three-phase flow. The study's findings directly inform several practical decisions:
- Material selection: When branch inlet liquid volume fractions exceed certain thresholds, the bell-shaped erosion pattern concentrates wear on the lower downstream wall, potentially requiring erosion-resistant lining or hardfacing in that specific zone.
- Geometry optimization: The inverse relationship between branch diameter and erosion rate suggests that oversizing the branch connection may be a viable strategy to reduce erosion rate, even if it increases the total affected surface area.
- Inspection planning: The transition from elliptical to bell-shaped erosion patterns with increasing liquid fraction provides a basis for developing condition-based inspection intervals. Engineers can correlate observed liquid loading levels with predicted erosion morphology to prioritize ultrasonic thickness measurements (UT/PAUT) on critical zones.
- Erosion allowance in wall thickness: The E/CRC model outputs provide quantitative erosion rate predictions that can be incorporated into API 579/ASME FFS-1 fitness-for-service assessments, enabling more accurate remaining life calculations.
Reflections on Limitations and Future Directions
While the study provides valuable parametric insights, several limitations merit consideration in engineering application. The E/CRC model, while widely used in industry, has known limitations in predicting erosion rates at very high impact angles and in the presence of thick liquid films. The study does not address the effect of particle size distribution, which is critical in real production systems where particle size spectra can span from sub-micron to millimeter-scale. Additionally, the simulation assumes steady-state conditions, whereas actual production systems experience significant flow regime transitions between slug, bubble, and annular flow patterns.
For practical implementation, I recommend cross-referencing the simulation predictions with field data from pigging campaigns and inline inspection (ILI) results. The predicted erosion shapes can serve as validation targets for UT/TOFD scan patterns during in-service inspections. Engineers should also consider coupling these erosion predictions with corrosion rate data, as the combined effect of erosion-corrosion synergy often exceeds the sum of individual mechanisms by a factor of 2-5 in sour gas environments.
This study represents a meaningful contribution to the understanding of tee erosion under asymmetric operating conditions, and its parametric findings can be directly applied to improve design margins, inspection strategies, and remaining life assessments for tee fittings in natural gas production and gathering systems.