Flow Erosion Research Articles

Research on Erosion Characteristics of the Sleeve-Type Blowdown Valve in the Shale Gas Gathering and Transportation Station

Summary In the shale gas separation and sewage system, the separator removes the sand and sewage in the produced gas and releases the sand-carrying sewage to the blowdown pipeline. As an important throttling component in this system, the sleeve-type blowdown valve is severely eroded during operation. To address the problem, this paper carries out numerical research on the erosion characteristics during the real opening and closing process of the sleeve-type blowdown valve based on the FLUENT fluid simulation software. The change of the erosion at different velocities, sand mass flow, sand size, and shape coefficient is analyzed, and the main factor affecting the erosion characteristics is evaluated. Based on the above research, a novel blowdown valve is proposed, and the erosion characteristics and flow performance before and after improvement are compared. The results show that when the velocity is 8 m/s, the maximum erosion rate after improvement is 4.79×10−5 kg/(m2·s). Compared with the maximum erosion rate before improvement of 1.22×10−2 kg/(m2·s), the erosion inhibition rate reaches 99.6%. Moreover, the pressure loss is obviously reduced after the improvement, and the flow capacity is enhanced. This provides useful guidance for the improvement of blowdown valves and ensures the safe operation of shale gas production.

CFD ANALYSIS OF EROSION RATE IN OIL AND GAS PIPELINE INDUSTRY

Generally, crude oil has been transported via pipelines. But sand particles normally will be deposited on the bottom of pipeline due to the unsteady flow which it can lead to flow concern and erosion deformation of the pipe wall if not well treated. As well, this issue can cause the pipe destruction and burden the maintenance cost if not prevented in the early stage. Thus, this study was simulated to observe and predict the erosion rate at various parameter conditions with presence of solid particles in pipeline surface via COMSOL© Multiphysics 5.4 software. The model of k-ω turbulent and particle tracing were applied where several different potential impacting factors on the formation of erosion were investigated including fluids flow velocity, sand particle size, sand flow rates, pipe orientation and pipe diameter. The result simulation showed the area around the bend pipe had high predictions of erosion deformation where the maximum erosion rate of 732 mg/m2s was observed for 2-inch pipe diameter and 100mm particle size. Besides, it shows the erosion rate increased as the higher mass flow rate and fluid velocity. By applying these erosion models, it could be possible to foresee the maximum point of erosion deformation along the pipelines which can reduce maintenance cost and prevent flow assurance issues.

Correlation analysis and clinical significance of Musculoskeletal Ultrasound semi-quantitative grading with bone salt metabolism, Rheumatoid factor and Erythrocyte Sedimentation rate in patients with Rheumatoid Arthritis.

To determine the correlation and clinical significance of musculoskeletal ultrasound semi-quantitative grading with bone salt metabolism, rheumatoid factor and erythrocyte sedimentation rate (ESR) in patients with rheumatoid arthritis. This is a clinical comparative study. A total of 240 patients with rheumatoid arthritis admitted to Baoding NO.1 Central Hospital were selected according to the DAS28 score of rheumatoid arthritis, and were divided into four groups, with 60 cases in each group from May 2020 to May 2022. The differences and correlation of musculoskeletal ultrasound semi-quantitative grading, bone metabolism indicators, erythrocyte sedimentation rate and rheumatoid factor among the four groups were statistically analyzed. The scores of bone erosion, synovial hyperplasia, joint effusion and intrasynovial blood flow in Group-H were significantly higher than those in Group-R, L and M, with statistically significant differences(p=0.00). The procollagen Type-1 N-terminal propeptide(P1NP), bone-specific alkaline phosphatase(BALP) and osteoprotegerin(OPG) in Group-H were significantly lower than those in Group-R, L and M, with statistically significant differences(p=0.00); The tartrate-resistant acid phosphatase(TRAC) in Group-H was significantly higher than that in Group-R, L and M, with a statistically significant difference(p=0.00). The levels of RF and ESR in Group-H were significantly higher than those in Group-R, L and M, with statistically significant differences(p=0.00). Musculoskeletal ultrasound semi-quantitative grading is correlated with the level of bone salt metabolism, rheumatoid factor and ESR in patients with rheumatoid arthritis. It can be combined with laboratory examination to objectively judge the severity of the course of rheumatoid arthritis.

Future perspective on the flow disturbances in Francis turbine caused by sediment particles and implementation of protection technique

This study aims to propose an experimental and numerical framework to examine the sediment erosion in the two components of Francis turbine (FT) namely guide vanes and runner blades. This study presents a critical review on the problems in hydroturbines due to the flow of sediment particles. Initially, a comprehensive literature survey on various numerical studies on flow disturbances and erosion in FT has been presented. Based on the literature survey, a research methodology has been proposed which can be implemented to protect the FT. The experimental methodology can be used as a validation to the numerical model for sediment erosion testing. Experimental methodology includes the protection of guide vane and runner materials with the help of suitable thermal spray coating. Numerical methodology was proposed to test the sediment erosion efficiency of FT while guide vane and runner are protected with hard thermal spray coatings.

Multiscale CFD Simulation of Multiphase Erosion Process in a Connecting Pipe of Industrial Polycrystalline Silicon Unit

Severe erosion phenomena often occur in industrial polycrystalline silicon units, leading to hydrogen leakage accidents and affecting long-term operation. It is favorable to use a computational fluid dynamics (CFD) simulation with the dense discrete phase model (DDPM) and the sub-grid energy-minimization multi-scale (EMMS) drag model to improve the prediction accuracy of complex multiphase erosion phenomena in a connecting pipe of an industrial polycrystalline silicon unit. Furthermore, the effect of droplet the specularity coefficient on boundary conditions is thoroughly considered. The predicted erosion behaviors are consistent with industrial data. The effects of operations parameters were discussed with three-dimensional CFD simulation, including droplet size and hydrogen volume fraction on erosion behaviors. The results indicated that the non-uniform multiphase erosion flow behavior near the wall can be simulated accurately with the EMMS drag model in a coarse mesh. A suitable droplet specularity coefficient such as 0.5 can also improve the accuracy of erosion position. Small liquid droplets, such as those of 30 μm size, will follow the gas phase better and have a lower erosion rate. The inertia effect of large droplets, such as those of 150 μm size, plays a dominant role, resulting in obvious erosion on the elbow walls. The erosion range and thinning rate enlarge with the increase in hydrogen volume fraction. A few silicon solid particles, such as 0.01% volume fraction, change local flow behaviors and probably cause the variation of local erosion positions. The process of erosion deformation first circumferentially extended and then accelerated at the local center position deeper.

Impact characteristics and erosion mechanism of solid particles in a centrifugal pump

To study the impact and erosion mechanism of solid particles in centrifugal pumps, the standard k–ε turbulence model and SIMPLE algorithm are adopted in this paper. Based on the discrete phase model (DPM) of the Lagrange method and McLaury erosion model, the flow impact characteristics and erosion mechanism of solid particles impacting the surface of flow passage components in a single-stage centrifugal pump were numerically simulated, and the test data were compared with the numerical simulation results of the external characteristics of the pump in clean water. The results show that the erosion mechanism of the pressure surface of the blade is mostly the impact erosion caused by high-speed particles with large impact angles, and the impact angle and impact velocity are larger near the tail of the pressure surface. The impact angle of solid particles on the shroud and hub is relatively small, but the erosion mechanism is still impact erosion. The erosion mechanism of the volute wall is mostly the cutting friction erosion caused by the low-velocity particles with small impact angles, and it is only impacted by the particles with large angles near the volute tongue, which is impact erosion. Overall, the average impact angle and impact velocity of the particles on the pressure surface of the blade are higher than those on the volute, so the erosion of the pressure surface is more serious than that of the volute in theory. The research results have certain theoretical reference value for improving the wear resistance of a centrifugal pump.

Ablation resistance of Ti2AlC ceramics at 1600 °C in nitrogen plasma flame

In this study, the ablation behaviors of high-purity and dense Ti2AlC ceramics were investigated using a nitrogen plasma flame at 1600 °C. The macromorphology, microstructure, and phase compositions of the Ti2AlC ceramics were systematically investigated at ablation times of 60–180 s. At the longest ablation time of 180 s, the linear ablation rate and mass ablation rate are 3.09 μm/s and −0.473 mg/s, respectively, which indicate noteworthy ablation resistance. The ablated cross-sections of the samples can be classified into three layers: a surface oxide layer, transition decomposition layer, and matrix layer. The surface oxide layers formed at different ablation times are composed of Al2TiO5 and TiO2, and the transition decomposition layer comprises mainly TiC and Al2O3 with pores and cracks. The high ablation resistance of the Ti2AlC ceramics is primarily due to the formation of the surface oxide layer, which effectively prevents oxidation and plasma flow erosion.

Water–soil interactions: Unravelling the processes and stages involved in the wetting of water repellent soils

The water repellent behaviour of soils is a widely studied phenomenon given its implications for infiltration, runoff, erosion and preferential flow. However, the principles underlying the eventual penetration of water into affected soils remain poorly understood. Theoretical considerations of the energetics and kinetics involved as a water drop makes contact with a water repellent soil surface and eventually penetrates into the soil suggest three distinct stages in the overall process. These stages are 1) adhesional wetting as soil and water first make contact, followed by 2) a kinetic barrier transitional stage in which molecular reorganisation of organics on soil reduces the water-soil contact angle to allow the water drop to sit deeper over soil particles of initial contact such that there is contact with particles in directly underlying soil layers, and finally 3) branching interstitial wetting as water penetrates into the bulk soil. Studies presented here of optical microscopy, mass of soil initially wetted, penetration time through layers of soil of different thicknesses, and time-dependent measurements of contact angle, volume of water penetrated, and mass of soil wetted, all give results consistent with this model. However, only for highly water repellent soils can distinct stages in wetting be clearly resolved experimentally, presumably because only these soils have a high enough kinetic barrier in the transitional stage for good separation between stages. For less water repellent soils, while the general time dependent behaviour remains consistent with the model, the distinction between the three stages is not so easy to resolve experimentally. The roles of contact angle, particle size distribution and drop size in determining the rates of these stages is considered, and the implications of the model for understanding soil water repellency are discussed.

Temporal and Spatial Evolution Laws of Freezing Temperature Field in the Inclined Shaft of Water-Rich Sand Layers

This study investigated the distribution and evolution characteristics of the temperature field during the freezing and excavation of inclined shafts, with the freezing open-excavation section of Shengfu Mine’s main inclined shaft (located in Shaanxi Province) as the project background. Utilizing field-measured data and the finite element software COMSOL Multiphysics, a 3D freezing temperature-field numerical calculation model was constructed to examine the temporal and spatial evolutions of the temperature field during the construction of the inclined shaft. The findings showed that after 88 days of freezing, the average temperature of the frozen wall in the open-excavation section was below −12 °C. The frozen wall thickness in the sidewalls of different layers exceeded 4 m, and the thickness at the bottom plate exceeded 5 m, meeting the excavation design requirements. For the same freezing time, the average temperature of the frozen wall in the fine sand layer was 0.28 to 2.39 °C lower than that of the frozen wall in the medium sand layer, and its effective thickness was 0.36 to 0.59 m greater than that in the medium sand layer. When the soil was excavated, and the well side was exposed, a phenomenon known as “heat flow erosion” occurred in the soil at the well-side position, causing the well-side temperature to rise. Nevertheless, this increase was generally limited, and when continuous cooling was applied, the well side could maintain a very low negative temperature level. Consequently, there was no spalling phenomenon. The effective thickness of the frozen wall during excavation did not decrease, with the average temperature remaining below −10 °C. Consequently, there was no large-scale “softening” of the frozen wall during excavation, thus ensuring construction safety. The numerical calculation model in this paper can be used to predict the development law of the freezing temperature field of the water–rich sandy layers in Shengfu Mine and adjust the on–site cooling plan in real time according to the construction progress. This research provides valuable theoretical insights for the optimal design and safe construction of freezing inclined-shaft sinking projects.

Experimental Investigation of Levee Erosion during Overflow and Infiltration with Varied Hydraulic Conductivities of Levee and Foundation Properties in Saturated Conditions

This study investigated erosion during infiltration and overflow events and considered different grain sizes and hydraulic conductivity properties; four experimental cases were conducted under saturated conditions. The importance of understanding flow regimes during overflow experiments including their distinct flow characteristics, shear stresses, and erosion mechanisms in assessing the potential for levee failure are discussed. The failure mechanism of levee slopes during infiltration experiments involves progressive collapse due to piping followed by increased liquefaction and loss of shear stress, with the failure progression dependent on the permeability of the foundation material and shear strength. The infiltration experiments illustrate that the rate of failure varied based on the permeability of the foundation material. In the case of IO-E7-F5, where the levee had No. 7 sand in the embankment and No. 5 sand in the foundation (lower permeability), the failure was slower and limited. It took around 90 min for 65% of the downstream slope to fail, allowing more time for response measures. On the other hand, in the case of IO-E8-F4, with No. 8 sand in the embankment and No. 4 sand in the foundation (higher hydraulic conductivity), the failure was rapid and extensive. The whole downstream slope failed within just 18 min, and the collapse extended to 75% of the levee crest. These findings emphasize the need for proactive measures to strengthen vulnerable sections of levees and reduce the risk of extensive failure.

Comparing turbulent flow and bank erosion with controlled experiments in a field-scale meandering channel

Abstract Bank erosion commonly occurs in alluvial rivers, shaping landscapes and riparian habitats and impacting water quality and infrastructure. Several models have been proposed that link shear stresses to bank erosion. However, data to test these hypotheses for characteristic geometries of meandering channels are sparse and technically challenging to acquire. Here we present results from a controlled experiment in a naturalistic channel to isolate the relationships between turbulent flow and nascent bank erosion. We ran the experiments at the Outdoor StreamLab (St Anthony Falls Laboratory, University of Minnesota) and gathered high-precision, contemporaneous measurements of the turbulent flow field and topography near a standardized, erodible bank at five locations along a single meander. The measurements show that the rate of bank erosion varied both along the channel and vertically and, local bank erosion was not correlated with any single hydrodynamic parameter. Upstream of the meander apex, erosion correlated with the near-bank time-averaged streamwise velocity magnitude while downstream of the apex, bank erosion correlated more strongly with near-bank turbulence parameters and depth. These results support field measurements that suggest that fluid shear contributions to outer bank erosion reflect multiple components of turbulent flow structure in river meanders.

Stability analysis of riverbanks with a dual structure under water-root-soil coupling.

The collapse mechanism of dual-structure vegetation riverbanks at different water levels is unclear. A method for calculating the critical collapse width of a dual-structure vegetation bank under different failure modes that consider the variations in river and groundwater levels and the influence of vegetation roots is proposed. Combined with the influence of flow lateral erosion and slope toe accumulation, a calculation model of riverbank stability was established. The results show that shear failure is the main failure mode when the cohesive soil layer on a dual-structure bank is thick, and the critical collapse width of the bank with root soil is higher than that of the soil bank. The critical collapse width of the bank varied with the water level during different water level periods. Compared with a soil riverbank, a rooted soil riverbank can significantly prolong the bank collapse time. The collapse width of a soil bank without vegetation roots is smaller than that of a rooted soil bank, and the cumulative collapse width is related to calculation time. The greater the thickness of rooted soil, the slower the decay rate of bank stability under water flow erosion.

Are high-severity burns in Alpine beech forests related to eruptive fire behavior?

Abstract In the complex Alpine topography, canyon-like terrain morphology combined with steep slopes may cause sudden changes in fire behavior, for instance, turning fast-spreading surface fires into eruptive fires. This phenomenon was observed during the vegetation dormancy in mountain beech (Fagus sylvatica L.) stands, a forest type usually spared from wildfires. The aim of this study is to understand the impact of eruptive fire behavior on the post-fire trajectories of beech trees. To this end, over the course of 3 years, we monitored the post-fire dynamics of 100 beech trees subjected to different fire intensities in February 2020, including a small-scale eruptive fire behavior associated with a steep slope (~41.4° on average) and canyon-like terrain features. Dendrometric and sylvicultural parameters and fire-related damages such as the presence of bark blackening, scars, foliage loss, crashes of branches or whole trees, and parasitic attacks were assessed. The observed medium to high fire severity was highly dependent on the position of the tree within the eruptive fire sector, whereas the relationship with dendrometric and sylvicultural tree characteristics, such as the crown position or the stem diameter, was less clear. Due to the lethal temperatures generated by the heat convection within the eruptive fire sector, the exposed bud and cambium tissues of the beech trees partially or totally failed to develop in the first post-fire season. Beech crown and stem tissue killed by exposure to convective heat fluxes rapidly dries out, loses its elasticity, and becomes brittle. Starting from the second post-fire season, this greatly increases the risk of the crashing of these crown parts and/or entire trees. A detailed understanding of the factors leading to an eruptive fire behavior would highly benefit the fire risk assessment and the detailed planning of the firefighting tactics in mountain areas, eventually increasing the safety of firefighters and preventing post-fire natural hazards, such as rockfall, erosion, surface runoff, and debris flow.

Growth promotion of Sargassum fusiforme by epiphytic microbes is dependent on the extent of interspecific interactions of the microbial community

Profound growth differences such as seedling length and biomass are often observed during the cultivation of Sargassum fusiforme despite the absence of detectable variance in abiotic factors that could have affected this process. This highlights the importance of biotic factors such as epiphytic microbiota in controlling seedling growth. Yet, how, and to what extent microbial activities can affect host growth in the presence of seawater flow and continuous erosion remains debatable. Particularly, the contribution of microbial network interactions to the growth of macroalgae remains poorly understood. This study aimed to compare the physicochemical properties of S. fusiforme seedlings via 16S rRNA gene Illumina sequencing-based profiling of the epiphytic microbial communities of seedlings with different lengths. Significantly different epiphytic bacterial communities were observed among S. fusiforme seedlings of different lengths. The result showed that community from longer seedlings maintained higher bacterial diversity with the taxa Gammaproteobacteria, Burkholderiales, Alteromonadales, Vibrionaceae, Ralstonia, Colwelliaceae, and Thalassotalea being selectively enriched. More importantly, microbial interspecific interactions, which were predominantly positive, were enhanced consistently in communities of the longer seedlings, indicative of reinforced prevalent and mutually cooperative relationships among the microorganisms associated with S. fusiforme seedlings of greater length. Furthermore, longer seedlings also displayed up-regulation of microbial functional potentials involved in N fixation and mineralization, P mineralization and transportation, and ion transportation compared with shorter ones. Lastly, stochastic processes dominated the community assembly of the epiphytic microorganisms. These findings could provide new insights into the relationship between microbial communities and growth in S. fusiforme seedlings and enable us to predict the community diversity and assembly of macroalgae-associated microbial communities. This could have important implications for linking microbial community diversity and network interactions to their host productivity.

CFD-DPM modelling of solid particle erosion on weld reinforcement height in liquid-solid high shear flows

Liquid-solid two-phase flow erosion can lead to great economic losses and safety concerns in the oil & gas storage and transportation industry. The solid particle erosion on weld reinforcement height (WRH) in liquid-solid high shear flows were systematically investigated by coupling the computational fluid dynamics (CFD) and the discrete phase model (DPM) method. The effects of liquid velocity, WRH height, particle flow rate, and particle size on the particle erosion in the WRH region were analyzed. Increasing erosion rate was observed with increasing liquid velocity. The erosion rate first increased and then decreased with increasing WRH height. The increasing erosion rate was observed with both increasing particle flow rates and particle sizes. The existence of defects on the WRH surface significantly increased the maximum erosion rates. This work can provide scientific guidelines for the corrosion protection and safe operation of long-distance oil & gas pipelines.

Flow field and erosion characteristics of the valve head of a vertical well inclinometer

During the drilling, the flow field characteristics of the vertical well inclinometer are complex. In a harsh underground environment, the valve head is prone to erode and wear. To overcome these disadvantages, the three-dimensional (3D) model of the valve head of a vertical well inclinometer was developed, and the flow field and erosion characteristics were simulated by using a numerical program. Results show that the maximum velocity and the maximum turbulence effect occur at the throttling inlet through which the drilling fluid passes the guide sleeve. By analyzing the particle trajectory of the discrete phase model, it is found that the particles above the valve head surface have obvious accumulation and rebound phenomena and that the erosion effect appears to some extent on the valve head surface. Analysis of the shear stress on the surface of the valve head shows that the stress concentration occurs in the valve head directly impacted by the drilling fluid and in the groove of the flow drilling fluid, respectively. The feasibility and accuracy of the simulation are verified by comparing the theoretical results with that of the field test. This investigation can be used to explain the failure of the valve head of the vertical well inclinometer and provide the design and optimization of the valve head with scientific support.

Impacts of Spontaneous Waterfall Development on Bedrock River Longitudinal Profile Morphology

AbstractRiver profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self‐form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1‐D numerical model to investigate how self‐formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km‐scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady‐state. Consistency between our model predictions and field observations of waterfall‐rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

Roles of soil surface roughness in surface–subsurface flow regulation and sediment sorting

Tillage-induced soil surface roughness (SSR) and rainfall interactions determine the process of water erosion. In this study, long-duration rainfall simulations were conducted to investigate the effect of SSR on the processes of surface–subsurface flow generation and soil loss, flow hydrodynamics, sediment sorting and transport on the steep red soil sloping farmland, southern China. Rainfalls with intensities of 30, 60, and 90 mm h−1 with duration of 180 min were applied to two soil surfaces (smooth and rough) with the slope gradient of 15°. Lesser surface flow and more subsurface flow were generated on the rough surface, which also accelerated the initial flow response and erosion stage advancement. The reduction effect of SSR on the sediment yield was significant under light rain or in the sheet erosion stage, while was limited when the rainfall intensity became stronger or entering the rill erosion stage. Differences in the spatial variability of soil crusts and the spatial distribution of rill networks were the main agents. Stream power took roughness factor into consideration well, thus was identified as the best hydraulic parameter for the prediction of soil loss rate. The investigation of particle size distribution of the eroded sediment showed the contents of sand, silt and clay changed from 15.9%, 64.9%, and 19.2% (average value of the original soil) to 22.9%, 53.1%, and 24% for the smooth surface, and to 37.8%, 42.1%, and 20.2% for the rough surface (average value of the eroded sediment). The selective and degrading effects of SSR on sediment particles were represented well. The rough surface was more inclined to transport clay particles than the smooth surface in the former period of erosion. Sand particles instead become more enriched with the extension of rainfall time. The results are beneficial for controlling the water erosion on agricultural sloping farmland.

Novel Strategy for Thermal Evaluation of Film-Cooled Blades Using Thermographic Phosphors at Simulated Engine Conditions

Abstract Practical strategy for the thermal evaluation of film-cooled blade is of great importance to the gas turbine community. Due to the physical or methodology limitations, it is difficult to evaluate the blade’s thermal performance at simulated engine conditions. The present study proposed novel focal-sweep-based phosphor thermometry for blade cooling inspection. While Mg4FGeO6:Mn (MFG) served as the temperature sensor to quantify the blade temperatures as well as simulated the thermal barrier coating (TBC) effect, the focal sweep method was adopted to overcome the optical constraints in cascade testing. The obtained MFG results of microstructures, jet impingement, and anti-erosion test demonstrated that the MFG phosphor is robust enough to simulate the thermal insulation effect of TBC and can withstand high-speed flow erosion. Furthermore, the proposed strategy clearly captured the blade temperature distributions (mainstream at T0,∞=∼850K) with high spatial resolution, which was then successfully remapped onto the three-dimensional twisted blade. Additional comparisons with the thermocouples demonstrated that the simulated TBC has a thermal insulation effect of about 68 K. This study addressed the common problems of phosphor thermometry in blade cooling evaluation, offering a practical strategy for future thermal diagnostics of the gas turbine.

Numerical modelling of multiphase flow hydrodynamics by CFD analysis methods in a double-elbow pipeline

Introduction. Prediction of multiphase flow patterns in pipelines of gas and oil industry is a complicated hydrodynamic process. Hydraulics of gas-liquid flows and erosion processes in pipelines with many bends in which gas, water, oil and air move are insufficiently studied. There are known studies of movement in single-phase flow conditions with one elbow, which is insufficient for design of modern pipeline systems. Therefore, there is a need for in-depth analysis of interaction of media in multiphase flow and transported sand particles. In this paper the effects of elbows on multiphase flow hydraulics and erosion of pipe sections are investigated using computational fluid dynamics (CFD) modelling tools for 100 seconds’ of process.
 
 Materials and methods. A liquid volume model (VOF) was used to simulate three-phase flow: air-water with solids in a pipe with two elbows. Turbulence effects have been accounted for with the RNG k–ε model. The model verification methods from previous studies have been used. For the numerical solution, the Ansys Fluent version 20.1 software package has been used. The results of CFD modelling of the volumetric gas content have shown a good agreement with the available experimental data.
 
 Results. The results showed that the change in the flow regime remains unchanged before the first elbow and after the se-cond elbow in the cork flow. However, for the initial churn flow regime, the flow regime varies at different segments of the flow area before and after each elbow.
 
 Conclusions. The preliminary churn flow in the upper vertical section was transformed to wavy stratified flow in the horizontal section between the two elbows and wavy annular flow in the vertical pipe after the second elbow. The flow pattern in slug flow maintained same after the first elbow, but the Taylor bubbles are lengthier in the horizontal section between the first and second elbows.