Measurements of the flow of a liquid–solid mixture/suspension through a segmented orifice | Scientific Reports

Scientific Reports volume 14, Article number: 269 (2024) Cite this article

725 Accesses

Metrics details

The paper attempts to solve the metrological problem that occurs when measuring the intensity of a flowing fluid with suspended solids with densities greater and less than the density of the fluid. The issue of the possibility of self-cleaning of a prototype variant of a segmented orifice from floating solid particles forming mixture/suspensions is discussed. For spherical particles of solids calculations have been made to allow for determining a borderline between their floating and entrainment by the flow, based on dimensionless numbers: Archimedes number and Reynolds number. Experimental tests and CFD simulations were conducted with a variable flow determined by Reynolds number for comparable segmental orifices with orifice module m = 0.102. Flow characteristics were plotted. Based on the results obtained from numerical simulations, positive influence of the inclination of skew segmental orifice downflow plane was presented. The results obtained from the study are a guideline for planning further studies to expand the knowledge of segmented orifices with inclined inflow plane.

The flow of fluid–fluid or fluid–solid mixtures often occurs in industrial processes. Therefore, knowledge of continuous phase of fluid (gas or fluid) in a mixture is desired when designing, manufacturing and using machines or equipment in the such processes. Flow of mixtures that occurs in technological processes can be categorized into: sedimentation, fluidization, pneumatic transport, and hydraulic transport. Source literature presents fluid mixtures in the form of water–oil or oil–water for different ratios in the studied mixture1.

However, considering only the fluid–solid mixture itself, depending on particle density of fluid and solids, we are dealing with unbounded movement of settling or entrained particles, including their forced movement on curvilinear paths. It is problematic to describe such complex problems with mathematical equations easily2. The occurring unbounded velocity of falling or entrainment depends on many factors. However, source literature mainly presents the problem of solids settling in a fluid. With contemporary state of knowledge and measurement instruments used, there is no problem in determining the velocity and trajectory of movement of individual particles. Yet, theoretical bases existing at this point are necessary for describing movement in the case of grouping of solids3. Such approach justifies further search for research methods and theoretical solutions in the case of movement of fluid–solid mixtures in horizontal pipes. If we limit the problem to hydraulic transport only, we will deal with fluid–solid type mixtures. This means of transport mainly occurs in construction and extractive industry4. It is used for transporting loose materials that do not dissolve in water, such as coal, ore or sand. In terms of transporting materials in horizontal pipes, literature presents a simplified model of behavior of particles in a solid–fluid mixture in transitory and turbulent flow (equation of motion, lift and drag coefficients)5.

Flow structure in a mixture depends on the ratio of fluid density to solid phase density, including its concentration, but also on the velocity of the flowing fluid. When solid particles are of similar density to the fluid, they move along with it with similar velocity to the flow, creating a uniform structure in its entire volume. Such flow type occurs for very small particles, usually smaller than 0.15 mm, and collisions among them are very rare. When particles of the solid phase are larger or their density is higher than the density of the fluid, the flow becomes asymmetric relative to the axis of the horizontal flow channel, the particles move in its bottom part6. When the velocity of the flowing fluid is reduced, the asymmetry increases, causing a flow of mixture with a moving bottom deposit. What happens is that solid particles roll on the bottom of the horizontal pipeline with much lower velocity than the flowing fluid7,8.

In industrial technological processes, where constant measurement of the volumetric flow rate is required, such mixtures create metrological problems even at low concentration of solid particles in the fluid (e.g. contamination of furnace oil with sand). For this reason, mechanical flow meters are useless for measuring, since rotating measurement elements become damaged due to concentration of solid particles in the metering mechanism. Among non-invasive flow meters, for measuring the volumetric flow rate of a fluid–solid suspension, we use ultrasonic and electromagnetic flow meters9,10. However, the most popularized, cheap and reliable instruments for measuring flow are constriction flow meters. They utilize sudden constriction of the pipeline that the flowing fluid encounters. In the constriction, rapid increase in flow velocity occurs, resulting in creating a difference in static pressure before and after the constriction. The following venturi are used as reducers for measuring the flow: venturi tubes, nozzles and orifices11,12. In spite of their flaws, among which are: constant pressure loss with fluid flow and low rangeability (4:1), orifices are popular in measuring flow, due to their simple structure and reliability in use. Among standardized, known structures of measuring orifices, the standard (centrical) orifice has been adopted as basic. It has a typical accuracy of 1.5–2%13. Conducted experimental tests has revealed that during the flow of contaminated fluids such as fluid–solid suspension, the geometry of the pipeline and inlet edge of the orifice changes through deposition of solids, causing a change in flow kinematics14. The example proves that when flow requirements for fluids, specified in the ISO 5167-1 standard15 cannot be met, there is a possibility to use non-standard orifices. They include such orifices, as: quadrant orifice (used particularly for flows with low Reynolds numbers), eccentric orifice, or segmental orifice. We use the latter two types of orifices mainly for the flow of liquids contaminated with solids in which the density of the contaminant particle is greater or less than the density of the fluid. In this case, we adjust the orifice opening so that the impurities can flow freely through the measuring orifice16,17,18. Particles of inclusions with density \({\uprho }_{{\text{cz}}}>{\uprho }\) (particle heavier than the fluid) are transported in the lower part of the pipeline, so the through hole of an eccentric or segmental orifice is placed so that it is tangent to the lower part of the pipeline. The inclusion particles with a density lower than that of the fluid \({\uprho }_{{\text{cz}}}<{\uprho }\), according to the buoyancy force, move in the upper part of the pipeline along with the fluid. In this situation, the through hole of the measuring orifice is placed in the upper part of the pipeline, which ensures the free flow of impurities in the measured fluid. This avoids the formation of deposits that restrict the flow in the pipeline in front of the venturi. Placing on the opposite side from the through hole of the orifice, impulse holes for taking pressure, eliminates the risk of plugging (clogging) of the holes, which can cause erroneous measurement of the difference of accumulated pressure14,15.

Source literature contains few articles on flow through segmental orifices. One of them presents simulations of large vortexes (the SBES model) during flow through segmental orifices with different constriction. The obtained results of flow were verified by measurement using the ultrasonic flow meter and the Laser Doppler Anemometric flow meter (LDA)19. The influence of distortion caused by sudden change in flow direction due to inserting a 90° elbow in the flow system before the segmental orifice was presented in the form of experimental tests. The influence of distance between the elbow and the segmental orifice on the accuracy of flow rate measurement20. Whereas, article21 presents flow study of mixtures with various concentration, in the form of emulsion that combines water and special oil dissolved in. Study was conducted for two constriction coefficients \(\beta\) equal to 0.3 and 0.5. It was found that the value for flow coefficient decreases whenever the constriction coefficient increases. In another article the authors presented experimental and numerical tests for a prototype flow meter with an segmental orifice inclined by an angle of 60° against the axis of the pipeline for constriction \(\beta =0.5\). The authors assessed the value of flow coefficient in a developing turbulent flow22.

Metrological problems are encountered when measuring the stream of the liquid phase contaminated with a mixture of solid particles of different densities. Particles of solids with a higher density (\({\uprho }_{{\text{cz}}})\) than the density of the liquid (\(\uprho\)) flow freely through the orifice constriction set at the bottom of the pipeline. On the other hand, particles with a density lower than that of the liquid (\({\uprho }_{{\text{cz}}}>{\uprho }\)) accumulate in the upper part of the pipeline in front of the orifice disrupting the flow and the measurement of the differential pressure value being an indirect measurement in determining the flux of the flowing liquid–solid mixture. By inclining the inflow plane of the segmental orifice by angle \(\alpha^\circ\) it was assumed that the new velocity distribution will allow the flowing fluid to entrain more solid particles with density of \({\rho }_{cz}<{\rho }\) from area located before the measurement orifice. The presented assumption was defined as ‘self-purging process’ of the measurement orifice from deposits that settle before it.

In order to solve the problem in the top part of the pipeline, the article proposes a flow meter based on an area reducer in the form of a skew segmental orifice that enables self-purging. For such a prototype of a measurement orifice, numerical study of ‘self-purging’ have been conducted for a flowing mixture of fluid and solids with lower density than the fluid.

Experimental tests and CFD numerical simulations have been conducted for a segmental orifice with module \(m=0.102\) made according to standard17 and on a prototype segmental orifice wich inflow plane inclined by angle α in line with the flow direction (Fig. 1).

The studied object: (a) segmental orifice (m = 0.102); (b) segmental orifice (m = 0.102) with inclined plane [own work].

Results of conducted experimental measurements have produced flow characteristics through a segmental orifice and segmental orifices with inflow plane inclined by angle \(\alpha =70^\circ\). They have been presented as mass flow of the flowing fluid in a function of build-up value of static pressure on a measurement orifice during flow.

To theoretically determine the flux of a flowing incompressible, non-viscous fluid, we use the stream continuity equation in conjunction with Bernulli's equation. The stream continuity equation for a flowing incompressible fluid is given by the relation:

The point of constriction of the orifice \(({F}_{1}>{F}_{2})\), there is an increase in the velocity of the fluid (\({\vartheta }_{2}>{\vartheta }_{1}\)). According to the principle of conservation of energy (Bernoulli's law), when there is an increase in kinetic energy (increase in velocity), there is a decrease in potential energy (decrease in pressure), so:

After the simplifications adopted and the transformations carried out, the expression for the theoretical volume flux \({q}_{V}\) was obtained:

However, in practice, the value of the flux of the flowing fluid determined from the above relationship differs from the actual flux. This is a consequence of the simplifications that were adopted in the mathematical attempt to describe the physical phenomena occurring during the flow of fluid through the measuring venturi. For this reason, the flow coefficient C (proportionality factor) was introduced into the equation determining the actual value of the flowing fluid stream, the equation takes the form:

Results obtained from experimental tests have been subjected to further calculative analysis to determine type B uncertainty of the conducted measurements of flow and pressure build-up on the orifice. The relative standard uncertainty of the measurement of the pressure pile-up at the measuring orifice is equal to:

From relation (5), the relative expanded uncertainty of the measurement of the pressure buildup \(\Delta p\) at the measuring orifice was determined as \(U(\Delta p)/\Delta p\) (with the assumed confidence level \(p=\mathrm{0,95}\) and the expansion factor \({k}_{p}=2\)), which can be written in the form of the relation:

We determine the relative standard uncertainty of the measurement of the mass flux flowing through the measuring orifice from Eq. (7):

where:

the relative standard uncertainty of the flux measurement is:

the relative standard uncertainty of the determination of water density, it was determined using the table of water density as a function of temperature from the relationship:

In the measurement series carried out, the temperature of the flowing water changed in the range up 2 °C, which corresponds to the relative standard uncertainty of its determination

From relations (7–9), the assumed relative expanded uncertainty of the measurement of the mass flux \({q}_{m}\) flowing through the measuring orifice was determined as \(U({q}_{m})/{q}_{m}\) (with the assumed confidence level and expansion factor \({k}_{p}=2\)), which can be written in the form of the relation:

The final stage will use graphic velocity distribution, from numerical simulations, before the measurement orifice, for different angles of inclination of the downflow plane, to analyze the velocity of entrainment of a single particle. Based on the analysis, a reduction coefficient for adverse, blind spot in the area of the pulsive opening before the measurement orifice will be determined.

Numerical tests used in fluid mechanics (CFD) has become a very popular tool used in the development new industrial measurement systems. The CFD numerical simulation has been performed in ANSYS FLUENT 2020R1 software where turbulent models \(k-\omega BSL\) (2 equations) and \(TransitionSST\) (4 equations) have been subjected to selection analysis. Both presented models are an extension of model k–ω standard23,24.

Model k–ωBSL has been selected to conduct numerical flow tests in the work25 where validation from several turbulent models for similar flow conditions was performed. The adopted model k–ωBSL, allows for using the benefits of the known calculation model k–ω in internal area of the boundary layer, whereas in the external part modelling is performed in model k–ε which, excluding boundary layer, does a much better job. The model has been described with Eqs. (11, 12):

Equation (11) is identical to that in model k–ω standard, while in the energy transport equation \(\omega\) (12) the term \({D}_{\omega }\) has been added, which is a component of cross diffusion that allows to connect model \(k-\omega\) with \(k-\varepsilon\).

Source literature reveals that the model \(TransitionSST\) is the most accurate and reliable in calculating flow, compared to known models from \(k-\omega\) group. It allows for a smooth transition from the turbulent model \(k-\omega\) into the turbulent model \(k-\varepsilon\) which is predisposed for numerical calculations in the main flow. To determine the course of kinematic energy transport \(\omega\) in the model \(TransitionSST\), an equation from the model \(k-\omega standard\) has also been used (13):

The equation for turbulent kinematic energy k has been modified to26,27 (Eq. 14):

In this model, Eqs. (13, 14) have been complemented with optional equations for discontinuity transport \(\gamma\) (Eq. 15):

It includes transitioning in the boundary layer from laminar flow into turbulent flow, taking account of, apart from the former, transitory momentum thickness transport (Eq. 16):

The above models have a calculative mechanism of transition from model k–ω in the internal area of the boundary layer to the model k–ε in its external part. However, they differ in the degree of solving this phenomenon. Detailed information and methods for determining individual parameters in the presented Eqs. (11–16) can be found in the Ansys Fluent 2020R1 software manual26.

The necessary computational grid is an important component determining the results and time necessary for conducting numerical simulations28,29. Insufficient quality of the computational grid in the conducted numerical simulations will result in erroneous results in the solution of equations during iteration and/or insufficient convergence in results. Whereas, a grid built of individual cells which are too small (fine-grained structure) excessively prolongs the time necessary to generate it and the duration of calculations themselves30,31,32,33

Model domain is being divided into smaller, individual volumes, forming a MOV (finite volume method) type computational grid. To conduct numerical tests, a grid has been developed using the Mosaic Meshing Technology structure34. A significant refinement of the grid near the walls has been proposed, allowing for presenting flow phenomena in the boundary layer (transition of the laminar flow into turbulent flow). Vortexes begin to occur in this area, causing energy loss, and areas with accumulation of contaminating particles begin to appear. For the grid that represents boundary layer, elements in the shape of polyhedral prisms have been selected. They transit through buffer layers (including polyhedral elements) into cubical cells in the form of cuboids in the central part (core) of the pipeline.

The adopted boundary conditions determine movement inside the domain and the selection of correct parameters influences the process of calculations. For the conducted tests, a schema has been selected where the outlet plane is described with a condition Pressure Outlet (140 kPa of absolute pressure). The inlet plane to the measurement channel is described with a condition Velocity Intel, implementing planes with velocity profiles that correspond to the developed mass flows. Border walls of the models have been described by the parameter wall that does not consider wall roughness. The fluid adopted for numerical calculations is water with constant temperature of 20 °C, density at 998.2 kg/m3 and dynamic viscosity of 0.001003 kg/ms.

The values of flow rate of fluid flowing through the tested orifices has been selected in such a way that they coincide with the scope of obtained values of mass flow at the measuring station. Five velocity profiles have been created, corresponding to the following mass flow rates: 0.25 kg/s, 0.35 kg/s, 0.45 kg/s, 0.55 kg/s and 0.70 kg/s

The settings of the computational method were read in a standard configuration suitable for the Transition SST turbulnetic model using the Coupled scheme, with the double precision solver enabled. During the iteration, the value of the flux of the flowing fluid on the inlet and outlet planes was monitored. The condition for completing the iteration was for the residual monitor to reach a residual of 0.001 for the continuity and 5 model equations describing the Transition SST model.

For the purpose analysis of the research problem it has been assumed that solid particles that form inclusions in the passing fluid are spherical in shape, with density \({\rho }_{cz}\) and diameter \({d}_{cz}\) (Fig. 2). For the purpose analysis of the research problem it has been assumed that solid particles that form inclusions in the passing fluid are spherical in shape, with density \({\rho }_{cz}\) and diameter \({d}_{cz}\) (Fig. 2). Using the Archimedes criterion number, a set of particles suspended in a liquid with density \({\rho }\) and dynamic viscosity was determined \({\mu }\). This set is described by the following equation:

where \({L}_{V}\) is the ratio of spherical particle volume \({V}_{cz}\) to its total value \({P}_{{c}_{CZ}}\):

Forces acting on a spherical particle submerged in fluid [own work].

After simplifications, the relationship of Archimedes' number (17) can be written in the form of Eq. (19) in which the parameters defining the size and mass of the inclusion particle, as well as the density and viscosity of the flowing fluid, are "sewn in”:

Using this correlation of the physical properties of a particle (diameter \({d}_{cz}\) and density \({\rho }_{cz}\) suspended in a fluid characterized by dynamic viscosity \(\mu\) and density \(\rho\), the size of a particle suspended in a flowing stream of flowing fluid is defined in a dimensionless way.

A solid particle in the fluid can move on forced curvilinear paths resulting from velocity distribution as well as free falling/floating. Therefore, the manner of movement of particles in a certain area is dependent on the velocity of the fluid. To determine the phenomenon of self-cleaning of the measuring orifice from impurity particles (with a density \({\rho }_{cz}\) lower than the density of the fluid \(\rho\)) located in front of the measuring orifice, the test was limited only to velocities resulting from free falling/floating, depending on the density of the examined particle (\({\rho }_{cz}\))1,35. Sense of vector of particle velocity is determined by the difference between fluid density (\(\rho\)) and density of a solid (\({\rho }_{cz}\)) submerged in the fluid. In a case when \({\rho }_{cz}-\rho <0\) the fluid displaces the submerged solid, therefore the velocity vector is sensed in the opposite direction to gravitational force.

The values of the velocity of a freely moving particle are calculated from the balance of the forces acting on the particle of a solid body immersed in a fluid according to Eq. (20) The values of the forces of gravity \({F}_{G}\), buoyancy force \({F}_{W}\) and drag force \({F}_{D}\) which are illustrated in the figure (Fig. 2) are described by equations:

where

By substituting the formulas for the individual components of \({F}_{G}\), \({F}_{W}\) and \({F}_{op}\), under the force balance Eq. (20), the equation was obtained:

The drag coefficient \({C}_{D}\) existing in Eq. (24) is a function of number \({Re}_{cz}.\) Therefore, it is dependent on the diameter of particle \({d}_{cz}\) and its floating velocity \({v}_{u}\) as a result of free, unbounded movement36. The Reynolds number was calculated from the relationship (25):

Table 1 shows the equations for determining the CD value depending on the Reynolds number \({Re}_{cz}\) value for the flowing particle.

After rearranging the force balance Eq. (24), the value of particle velocity moving freely in the fluid has been determined in the following form:

When calculating the lift velocity \({v}_{u}\), one should perform a sequence of looped calculations, and for the first approximation, the value of the resistant coefficient \({C}_{D}=1\) should be taken. From the obtained value of the lift velocity, the value of the resistant coefficient \({C}_{D}\) based on Table 1 should be determined.

By generating a graph of the function of the criterion numbers \(Ar=f({Re}_{cz})\), one obtains a set of points forming the limiting curve of the particles that are carried away and float in the fluid. Knowing the equations determining the Archimedes number (19) and Reynolds number (25) of a particle, it is possible to determine the minimum fluid velocity that allows the floating particles to be carried away and transported through the segmental orifice and segmental orifice with an inclined inflow plane.

From the CFD numerical simulations performed, vector maps of the velocity of the flowing fluid through the measuring orifice under study were generated. They define the area upstream of the measuring orifice bounded by the pipeline plane, the inflow plane of the measuring orifice and a boundary line composed of points where the fluid velocity \(v\) is less than the calculated velocity that allows to carry away the particles \({v}_{u}\). By inclining the inflow plane of the segmented orifice, a new velocity distribution is created, with a decreasing area in which the fluid has a velocity \(0\le v{\le v}_{u}\). With an increase in the velocity of the flowing fluid (Reynolds number increases), the area of accumulation of particles in front of the measuring orifice is reduced by their entrainment by the flowing fluid stream.

To create a computational grid in a 3D model a volume closed inside the measurement pipeline has been adopted. The pipeline is composed of the tested measurement orifice and straight segments before and after the measurement orifice. It has been assumed the length of the inlet segment before the orifice will be 175 mm (3.5·D) long, and the outlet segment behind the orifice will be 425 mm (8.5·D) long.

The selection of numerical grid and turbulent model was made using a 3D model of the measurement pipeline with a segmental orifice with module \(m=0.102\). The tests were conducted for velocity profiles corresponding to mass flow 0.25 kg/s and 0.70 kg/s. 10 grids were generated in total (5 grid sizes for each defined mass flow) with refinement in the boundary layer. The grids have been marked with letters in alphabetical order—from the smallest A (ca. 1 million elements) to the largest E (ca. 4.5 million elements). The sizes are presented in Fig. 3.

The size of the tested grids [own work].

Simulation calculations have been made for all test grids (Fig. 4) using two turbulent models \(k-\omega BSL\) and \(TransitionSST\). On this basis, swelling values on the segmental orifice \(\Delta {p}_{CFD}\) have been determined for selected mass flows. In order to compare the value for differential pressure obtained from the CFD simulation with the theoretical value, engineering calculations of pressure build-up \({\Delta {p}}_{PN}\) have been made according to standard17 using the Eq. (27):

Results comparison of the tested numerical grids and turbulent models [own work].

In this equation, the value of mass flow \({q}_{m}\) is represented by a value in boundary conditions of numerical simulations as \({{q}_{m}}_{CFD}\) in the form of defined velocity profiles. Relative error of the obtained differential pressure on measurement orifice \({\delta }_{{\Delta {p}_{CFD}}}\) has been calculated. Relative error \({\delta }_{{\Delta {p}_{CFD}}}\) of comparison of swelling pressure has been calculated from the following dependence (28):

where: \(\Delta {p}_{PN}\) is the value of differential pressure on the measurement orifice, calculated according to standard17. \(\Delta {p}_{CFD}\) is the value of differential pressure determined from CFD numerical simulation.

The results obtained from the conducted validation were presented in the form of a graph (Fig. 4). When it comes to grid A with a specified mass flow \({q}_{m}=0.70\) kg/s and turbulent model k–ωBSL numerical calculations were not completed in the iteration process. The obtained relative error for the same computational grid with the same specified mass flow for model \(TransitionSST\) deviates from the results and trend of the other tested grids. For this reason grid A was rejected from further considerations.

With acceptable error of simulation and time of numerical calculations, computational grid C has been selected. It is composed of 2.14 million elements for flow \({q}_{m}=0.70\) kg/s and 1.84 million elements for flow \({q}_{m}=0.25\) kg/s Figure 5 presents a part of the selected grid in a cross-section and longitudinal section of the pipeline in the measurement orifice area.

The structure of the numerical grid for segmental orifice with module \(m=0.102\) for mass flow \({q}_{m}=0.70kg/s\) [own work].

Data collected from the conducted experimental tests and CFD simulations was used to create charts and graphs that represent the tested flow parameters of the segmental orifice and the segmental orifice with inclined inflow plane. Table 2 presents the values of differential pressure \({\Delta p}_{CFD}\) from CFD numerical simulation, obtained on a measurement orifice with specified mass flows \({q}_{{m}_{CFD}}\). Simulations were conducted according to the established boundary conditions, using computational grid C and turbulent model \(TransitionSST\).

Experimental tests were conducted according to assumptions presented earlier. Apart from measurement of swelling (pressure difference on the measurement element) for individual flow rates of the fluid passing through the measurement pipeline, fluid temperature was also recorded before and after the series of measurements. For further calculations, arithmetic means from measurement series were used. Tables 3 and 4 compile data obtained from experimental tests and extended type B relative uncertainties of flow \({q}_{m}\) and differential pressure \(\Delta {p}\) measurements.

In a hydraulic measurement station (Fig. 6), using an eccentric pump with constant rate of delivery (1), the fluid in the measurement unit is moved from the main tank with absolute pressure at ca. \(140\) kPa. Before pumping the fluid into the measurement system, the fluid enters a deaerating vessel (2). Forced by the pump, the fluid flows out of the bottom port (2a) whose end is located inside the tank, at \(3/4\) of its height from the bottom tank end. In the top tank end there is a port (2b) from which fluid flow with possible air bubbles exits, via adjustable side bleed, through replaceable glands (3). This allows for gradual adjustment of the measurement volume flow within the range of \(0.25{{\text{dm}}}^{3}/{\text{s}}<{q}_{V}<0.65{{\text{dm}}}^{3}/{\text{s}}\). The final outlet port (2c) is also located in the top tank end but it is submerged inside the deaerator at the depth of ca. \(3/4\) of the entire vessel height. This allows for collecting fluid with no air bubbles and forcing it into the hydraulic measurement station. The fluid, separated from possible gas fractions, flows through the measurement pipeline (4) and then returns to the main tank, creating a closed system of fluid flow.

Measurement station—block diagram [own work].

The measurement pipeline (Fig. 7) is composed of reducing tested flow meter, connected to a differential pressure transducer (5), electromagnetic flow meter acting as standard (6), and straight segments made of stainless steel, with internal diameter of D = 50 mm. The segment before the tested flow meter with segmental orifice is 2.05 m (\(41\cdot D\)) in length, and the one behind the measurement orifice is 1.55 m (\(31\cdot D\))—both lengths remain within the scope recommended by source literature20. In order to retain axial symmetry when replacing the tested flow meter with a skew segmental orifice and straight segments of the pipeline at the point of connection, fitting flanges with a centering lock have been used.

Measurement system of a segmental flow meter with a flow meter acting as standard [own work].

In order to determine the dependence of function \({q}_{v}=f(\Delta p)\) simultaneous measurements must be made for the value of swelling \(\Delta p\) obtained on the tested orifice and the flow \({q}_{v}\) of fluid flowing through it at the same time.

The static pressure of fluid on the tested segmental orifice was measured at pressure input points located at the corner, on the opposite side of the flow-through opening, in the top part of the flow meter. Differences between high pressure area (before the orifice) and low pressure area (behind the orifice) were measured with a programmable APR-2000/ALW differential pressure transducer with an output current signal from 4 to 20 mA. The transducer for tests was programmed for measurement range of \(\Delta p=12.8\) kPa with time constant of \(t=5\) s with a limiting error \(\Delta p=0.15\%.\) The value of volume flow of the fluid passing through the tested orifice was measured with a PROMAG 30AT15 electromagnetic flow meter with output current signal from 4 to 20 mA. It is located behind the tested flow meter in a parallel DN15 pipeline. Its measuring range \({q}_{v}\) was set to \(3.6\) m3/h with the same time constant \(t=5\) s. The flow meter has a limiting error \(\Delta {q}_{v}=\pm (0.2\%\cdot {q}_{{v}_{mier}})\pm 0.05\%\cdot {q}_{v}\), and was used as standard in the measurement system. Fluid temperature was measured with an electronic thermometer graduated in intervals of \(0.1\) °C after prior calibration with standard laboratory thermometer. Temperature measurement was registered (recorded) at the beginning and end of each individual measurement series—arithmetic mean was considered in calculations. Temperature value is essential for determining fluid density, while determining mass flow and kinematic viscosity for calculating Reynolds number.

Flow charts have been presented in the form of graphs. Figure 8 presents flow charts obtained from CFD experimental tests and simulations for a segmental orifice, while Fig. 9 presents the same for a segmental orifice with plane inclined by angle \(\alpha =70^\circ\). Power trend lines have been plotted on the graphs for \({q}_{m}=f(\Delta p)\), obtaining their equations and alignment coefficients \({R}^{2}\).

Flow chart for a segmental orifice with module \(m=0.102\) [own work].

Flow chart of a segmental orifice with inflow inclined by angle \(\alpha =70^\circ\) and module \(m=0.102\) [own work].

Figure 10 presents a graphic image of results from numerical simulations for a segmental orifice (\(\alpha =90^\circ )\) and a segmental orifice with inclined downflow plane \((\alpha =70^\circ )\) with passing mass flow \({q}_{m}=0.70\) kg/s. Static pressure distribution before and after the orifice is presented as a spectrum of colors and velocity distribution as a vector spectrum of colors.

Static pressure distribution with velocity vectors for mass flow \({q}_{m}=0.7\) kg/s: (a) segmental orifice, (b) segmental orifice with inclined plane [own work].

Having analyzed the influence of inflow plane inclination, an example set of solids has been determined which contains solids whose particle diameter remains within a range of \({d}_{cz}=0.2\dots 1\) mm with density \({\rho }_{cz}\) from a range \(700\dots 950\) kg/m3 submerged in a fluid with parameters consistent with the fluid used in numerical calculations (\(\rho =998.2\) kg/m3, \(\mu =10.3\times {10}^{-4}\) Pa s). Using Eq. (20) for this data, velocity value \({v}_{u}\) for a particle was calculated, and then the Reynolds number \({Re}_{cz}\). When determining the value of floating speed \({v}_{u}\) a series of 10 computational iterations was conducted, approximating the value of floating velocity \({v}_{u}\). When calculating the first approximation of floating velocity \({v}_{u}\) the influence of drag coefficient was omitted and value \({C}_{D}=1\) was assumed, resulting in obtaining input value of floating velocity. On the basis of the paper36, where the most popular methods of determining the resistance coefficient are discussed, the relationships presented in Table 137. Based on the Reynolds number of the particle \({Re}_{cz}\) the value of drag coefficient \({C}_{D}\) was determined from a function. When determining the value of floating speed \({v}_{u}\) a series of 10 computational iterations was conducted, approximating the value of floating velocity \({v}_{u}\). The calculation algorithm has been presented in the form of a diagram (Fig. 11).

Diagram of the process of approximating particle movement velocity [own work].

The graph in Fig. 12 presents the value of floating velocity \({v}_{u}\) obtained in each following iteration during the conducted calculations for a particle with \(Ar=12.74\) number (Eq. 19). The horizontal axis is marked with subsequent numbers of calculation loop (iteration) where zero denotes the value of particle velocity with drag coefficient CD = 1.

The course of the obtained values of floating velocity for a particle with Archimedes number Ar = 12.74 in the process of iteration [own work].

The graph in Fig. 13 presents a set of point described with characteristic numbers function \(Ar=f({Re}_{cz})\). The points form a limit curve of entrainment of particles floating in the fluid. Knowing the equations that determine the Archimedes number (19) and Reynolds number (25) of the particle, it is possible to determine the minimum velocity of fluid which allows for entrainment and transport of particles floating in the fluid through the segmental orifice and segmental orifice with inclined inflow plane.

Boundary value of the Reynolds number where particles described with Archimedes number are entrained [own work].

Considering the calculated values of solid particle floating velocity (26), CFD numerical simulations have been used to determine areas of deposits ‘stockpiling’. The velocity of solid particle entrainment that occurs in this area is lower than floating velocity. The areas were being determined on the basis of vectorial distributions of velocity in the passing flow presented as an example in Fig. 14 for a passing flow \({q}_{m}=0.70\) kg/s with Archimedes number \(Ar=12.74\). Before the orifice plane, in its upper part, additional local vortexes occur (increase in fluid velocity) that entrain particles, but keep them in the area where the vortexes occur. Therefore, the surface of the vortexes was also added to the area of deposits ‘stockpiling’. The illustrations show that when the downflow plane of the segmental orifice is inclined, both the size of the occurring vortex and the area of the deposit decrease.

Area of vortex occurrence before the plane of the segmental orifice for flow \({q}_{m}=0.70\) kg/s [own work].

The conducted CFD numerical simulations of velocity distributions were used to determine areas \({A}_{\alpha^\circ }\) where velocity values of the passing fluid \(v\) remain within range \(0\le v{\le v}_{u}\) for a segmental orifice (\({A}_{\alpha =90^\circ })\) and segmental orifice with inclined inflow plane α: 80°, 70° and 60° (\({A}_{\alpha =n^\circ })\). Based on the obtained areas \({A}_{\alpha =90^\circ }\) and \({A}_{\alpha =n^\circ }\) relative percentage reduction of ‘stockpiling’ area was determined for a segmental orifice with a plane, in relation to a segmental plane. The reduction has been described with the Eq. (29).

This relation was acknowledged as quality index (\({\psi }_{\alpha =n^\circ }\)) of ‘self-purging’ of the segmental orifice as a result of inclining the inflow plane. The obtained area values are presented in Table 5, depending on the flow size of the passing fluid defined by Reynolds number.

The values of purging index \({\psi }_{\alpha =n^\circ }\) obtained from the calculations in a function of the number Re for module \(m=0.102\) have been presented in the form of a graph (Fig. 15).

The value of purging index \({\psi }_{\alpha =n^\circ }\) of the tested segmental orifices with module \(m=0.102\) [own work].

The article has presented the results of simulation tests in Ansys Fluent 2020R1 software as well as experimental tests on two orifices (a segmental orifice and a segmental orifice with downflow plane inclined by angle α) with module \(m=0.102\) in a DN50 pipeline. Scientific research has been conducted within the range of Reynolds number 6300 < Re < 17600 (area of developing turbulent flow).

Before conducting target CFD numerical calculations, an analysis was conducted of matching two turbulence models: k–ωBSL and \(TransitionSST\), for theoretical calculations for a segmental orifice with module m = 0.102 (\(\alpha =90\)°)17. Out of the two models, the \(TransitionSST\) model was characterized by smaller matching errors, with extreme mass flows (\(0.25\) kg/s and \(0.70\) kg/s) adopted for calculations.

For this reason, the turbulent model \(TransitionSST\) with a grid of 2.14 million elements was selected for further CFD numerical calculations. With acceptable error and time of numerical calculations, it allowed for accurate reflection of physical phenomena occurring in a passing fluid flow.

The obtained simulation results are convergent with conducted experimental tests for the segmental orifice (\(\alpha =90^\circ\)) as well as orifice with inclined downflow plane (\(\alpha =70^\circ\)) with module \(m=0.102\). Matching has also been achieved between points and power trend lines characterized by coefficient R2 equal to one.

In the proposed solution, the occurring inclinations of the orifice plane causes a stepless, smooth transition from circular cross-section of the pipeline to the area of flow-through opening in the shape of circular sector. When inclination angle \(\alpha\) is reduced, the length orifice plane increases, reducing the flow section. Therefore, the measurement value of swelling pressure on the orifice decreases, and, consequently, so does constant pressure loss.

The analysis involved the possibility of segmental orifice self-purging from solid particles suspended in the passing fluid \(({\rho }_{cz}-\rho <0)\), based on numerical simulations. For this purpose, calculations for spherical solid particles were made, determining the boundary line between particle floating and entrainment by the flow, based on non-dimensional numbers: the Archimedes number and the Reynolds number (Fig. 13).

The article presented a method to reduce the area of ‘stockpiled’ deposits (concentration of floating solid particles in fluid) before the segmental orifice, by inclining its downfall plane. As the inclination angle of the segmental orifice inflow plane increases, the area of ‘stockpiled’ deposits decreases. For an inclination angle \(\alpha =60^\circ\) with number Re ˃ 8800 it decreases by more than 50% compared to a segmental orifice with angle \(\alpha =90^\circ\).

As can be seen, the inclination of the downflow plane causes entrainment of a much larger amount of floating particles, reducing the area of ‘stockpiled’ deposits, which improves kinematics of the passing fluid flow.

Further research is planned which will allow to optimize the inclination angle of the segmental orifice downflow plane due to the occurrence of a similar deposit ‘stockpiling’ area on the outflow side of the orifice.

The datasets used during the current study are available from the corresponding author on reasonable request.

Archimedes number (–)

Surface area (mm2)

Drag coefficient (°)

Pipeline internal diameter (m)

Particle diameter (sphere) (m)

Drag force (N)

Gravity force (N)

Displacement force (N)

Cross-sectional area of the pipeline (m2)

Area of the measurement orifice opening (m2)

Gravitational acceleration (m/s2)

Orifice constrictions (m)

Characteristic dimension (–)

Orifice module (–)

Particle (sphere) volume (m3)

Particle (sphere) surface area (m3)

Static pressure (Pa)

Mass flow (kg/s)

Volume flux (dm3/s)

Liczba Reynoldsa (–)

Orifice plane inclination angle (°)

Dynamic viscosity coefficient (Pa s)

Particle floating velocity (m/s)

Average velocity of the flowing fluid in the pipeline cross-section (m/s)

Average velocity of the flowing fluid in the cross-section of the orifice opening (m/s)

Fluid density (kg/m3)

Solid density (kg/m3)

Turbulent viscosity (m2/s)

Differential pressure (Pa)

‘Self-purging’ quality index (%)

Effective diffusivity coefficients

Coefficient of generating turbulence kinetic energy as a result of mean velocity gradients

Coefficient of generating dissipation velocity in the transport of turbulence kinematic energy

Coefficient of dissipation of turbulence kinetic energy

Coefficient of dissipation in the transport of turbulence kinematic energy

Coefficient of displacement effect in turbulence kinematic energy transport equation

Constant definable coefficients

Derivative coefficients \({G}_{k}\) and \({Y}_{k}\) caused by separation

Coefficients of transition source, discontinuity destruction/relaminarization

Constant coefficient in model (Reynolds number)

Charles, M. E., Govier, G. W. & Hodgson, G. W. The horizontal pipeline flow of equal density oil-water mixtures. Can. J. Chem. Eng. 39, 27–36. https://doi.org/10.1002/cjce.5450390106 (1961).

Article CAS Google Scholar

Chong, J. S., Christiansen, E. B. & Baer, A. D. Rheology of concentrated suspensions. J. Appl. Polym. Sci. 15, 2007–2021. https://doi.org/10.1002/app.1971.070150818 (1971).

Article CAS Google Scholar

Łukasiewicz, E. & Rząsa, M. Modelling of sedimentation for polydisperse mixture. in WIT Transactions on Engineering Sciences vol. 115 3–11 (Computational and Experimental Methods in Multiphase and Complex Flow, 2017). https://doi.org/10.2495/MPF170011.

Baxter, J. W. & Bumby, J. R. Proceedings of the institution of mechanical engineers, Part I. J. Syst. Control Eng. https://doi.org/10.1243/PIME (1995).

Article Google Scholar

Zandi, I. Advances in Solid–Liquid Flow in Pipes and Its Application (Elsevier, 1971). https://doi.org/10.1016/C2013-0-02346-8.

Book Google Scholar

Parzonka, W. & Vočadlo, J. Méthode de la caractéristique du comportement rhéologique des substances viscoplastiques d’après les mesures au viscosimètre de Couette (modèle nouveau à trois paramètres). Rheol. Acta 7, 260–265. https://doi.org/10.1007/BF01985787 (1968).

Article Google Scholar

Erdal, A. & Andersson, H. Numerical aspects of flow computation through orifices. Flow Meas. Instrum. 8, 27–37. https://doi.org/10.1016/S0955-5986(97)00017-4 (1997).

Article Google Scholar

Tebowei, R. Computational fluid dynamics (CFD) modelling of critical velocity for sand transport flow regimes in multiphase pipe bends. (2016).

Andruszkiewicz, A., Wędrychowicz, W., Synowiec, P. & Piechota, P. Analysis of the possibility of using an aperture in the pipeline as an orifice for continuous measurement of gas mass streams. Przeglad Elektrotechniczny 95, 41–43. https://doi.org/10.15199/48.2019.11.10 (2019).

Article Google Scholar

Liptak, B. G. Process Measurement and Analysis. Society.

Kremlewski, P. P. Raschodimiery i scetciki kolicestva, Izd. Politechnika (2002).

Kabza, Z. Pomiar strumieni płynów (przewodnik). (Politechnika Opolska, studia i monografie, 1996).

Golijanek-Jędrzejczyk, A., Mrowiec, A., Hanus, R., Zych, M. & Świsulski, D. Uncertainty of mass flow measurement using centric and eccentric orifice for Reynolds number in the range 10,000 ≤ Re ≤ 20,000. Meas. J. Int. Meas. Confederation. 160, 107851. https://doi.org/10.1016/j.measurement.2020.107851 (2020).

Article Google Scholar

Levin, B. M. & Lopatin, A. N. Measurement of the discharge of suspended matter-bearing flows in hydromechanization systems using segmental orifice plates. Meas. Tech. 36, 36–39. https://doi.org/10.1007/BF00980014 (1993).

Article Google Scholar

PN-EN ISO 5167-1:2005. Pomiary strumienia płynu za pomocą zwężek pomiarowych wbudowanych w całkowicie wypełnione rurociągi o przekroju kołowym—Część 1: Zasady i wymagania ogólne. (Polski Komitet Normalizacyjny, 2005).

Kasprzak, D. & Mrowiec, A. Analysis of the possibilities of measurement with eccentric orifice of flow medium for small Reynolds numbers. Pomiary Automatyka Robotyka 20, 25–28. https://doi.org/10.14313/PAR_220/25 (2016).

Article Google Scholar

PN-93M-5950. Pomiar strumienia masy i strumienia objętości płynów za pomocą zwężek pomiarowych. (Wydawnictwo Normalizacyjne ALFA - WERO Sp. z o.o., 1994).

Golijanek-Jędrzejczyk, A. et al. The assessment of metrological properties of segmental orifice based on simulations and experiments. Meas. J. Int. Meas. Confederation. 181, 109601. https://doi.org/10.1016/j.measurement.2021.109601 (2021).

Article Google Scholar

Straka, M., Fiebach, A., Eichler, T. & Koglin, C. Hybrid simulation of a segmental orifice plate. Flow Meas. Instrum. 60, 124–133. https://doi.org/10.1016/j.flowmeasinst.2018.02.006 (2018).

Article Google Scholar

Straka, M., Koglin, C. & Eichler, T. Segmental orifice plates and the emulation of the 90°-bend Segmentblenden und die Nachbildung des 90°-Krümmers. Tech. Mess. https://doi.org/10.1515/TEME-2019-0120 (2020).

Article Google Scholar

Buhidma, A. & Pal, R. Flow measurement of two-phase oil-in-water emulsions using wedge meters and segmental orifice meters. Chem. Eng. J. Biochem. Eng. J. 63, 59–64. https://doi.org/10.1016/0923-0467(95)03076-X (1996).

Article CAS Google Scholar

Mrowiec, A., Heronimczak, M. & Golijanek-Jędrzejczyk, A. Skośna kryza segmentowa do pomiaru przepływającego strumienia płynu. Metrologia naukowa, normatywna i przemysłowa: wybrane zagadnienia 295–308 (2020).

Lu, Z., Piro, M. H. A. & Christon, M. A. Mesh and turbulence model sensitivity analyses of computational fluid dynamic simulations of a 37M CANDU fuel bundle. Nucl. Eng. Technol. https://doi.org/10.1016/j.net.2022.06.004 (2022).

Article Google Scholar

Li, J., Zhong, C., Pan, D. & Zhuo, C. A gas-kinetic scheme coupled with SST model for turbulent flows. Comput. Math. Appl. 78, 1227–1242. https://doi.org/10.1016/j.camwa.2016.09.012 (2019).

Article MathSciNet Google Scholar

Golijanek-Jędrzejczyk, A. et al. A numerical and experimental analysis of multi-hole orifice in turbulent flow. Meas. J. Int. Meas. Confederation. https://doi.org/10.1016/j.measurement.2022.110910 (2022).

Article Google Scholar

Fluent Theory Guide Fluent2020R1. https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/main_page.html (2022).

Moshfeghi, M., Song, Y. J. & Xie, Y. H. Effects of near-wall grid spacing on SST-K-ω model using NREL Phase VI horizontal axis wind turbine. J. Wind Eng. Ind. Aerodyn. 107–108, 94–105. https://doi.org/10.1016/j.jweia.2012.03.032 (2012).

Article Google Scholar

Pawłucki, M. & Kryś, M. CFD dla inżynierów. Praktyczne ćwiczenia na przykładzie systemu ANSYS Fluent. 288 (2020).

Golijanek-Jędrzejczyk, A., Świsulski, D., Hanus, R., Zych, M. & Petryka, L. Uncertainty of the liquid mass flow measurement using the orifice plate. Flow Meas. Instrum. 62, 84–92. https://doi.org/10.1016/j.flowmeasinst.2018.05.012 (2018).

Article Google Scholar

Jurkowski, S. & Janisz, K. Analiza wpływu parametrów siatki obliczeniowej na wynik symulacji przepływomierza. Autobusy Technika Eksploatacja Systemy Transportowe 20, 129–134 (2019).

Google Scholar

Czapla, B. & Rząsa, M. Analiza wpływu gęstości siatki numerycznej na wyniki symulacji ścieżki wirowej von karmana. Inform. Control Meas. Econ. Environ. Protection 7, 66–69. https://doi.org/10.5604/01.3001.0010.4842 (2017).

Article Google Scholar

Karcz, J. & Kacperski, L. An effect of grid quality on the results of numerical simulations of the fluid flow field in an agitated vessel. in 14th European Conference on Mixing Warszawa 10–13 (2012).

Karkoulias, D. G., Tzoganis, E. D., Panagiotopoulos, A. G., Acheimastos, S. G. D. & Margaris, D. P. Computational fluid dynamics study of wing in air flow and air-solid flow using three different meshing techniques and comparison with experimental results in wind tunnel. Computation 10, 34. https://doi.org/10.3390/computation10030034 (2022).

Article CAS Google Scholar

Krishna, Z. et al. Ansys mosaic poly-hexcore mesh for high-lift aircraft configuration. in 21st Annual CFD Symposium 1–11 (2019).

Nguyen, C. T. et al. Viscosity data for Al2O3–water nanofluid—Hysteresis: Is heat transfer enhancement using nanofluids reliable?. Int. J. Thermal Sci. 47, 103–111. https://doi.org/10.1016/j.ijthermalsci.2007.01.033 (2008).

Article CAS Google Scholar

Kelbaliyev, G. I. Drag coefficients of variously shaped solid particles, drops, and bubbles. Theor. Foundations Chem. Eng. 45, 248–266. https://doi.org/10.1134/S0040579511020084 (2011).

Article CAS Google Scholar

Kelbaliyev, G. & Ceylan, K. Development of new empirical equations for estimation of drag coefficient, shape deformation, and rising velocity of gas bubbles or liquid drops. Chem. Eng. Commun. 194, 1623–1637. https://doi.org/10.1080/00986440701446128 (2007).

Article CAS Google Scholar

Download references

Calisia University, Polytechnic Faculty, 62-800, Kalisz, Poland

Marcin Heronimczak & Andrzej Mrowiec

Opole University of Technology, Faculty of Electrical Engineering, Automatic Control and Informatics, 45-758, Opole, Poland

Mariusz Rząsa

Department of Biosystems Engineering, Poznan University of Life Sciences, 60-637, Poznan, Poland

Krzysztof Koszela

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Conceptualization, M.H. and A.M.; Methodology, M.H. and A.M.; Software, H.M.; Validation, M.H., A.M. and M.Rz.; Formal analysis, H.M. and A.M.; Investigation, H.M. and A.M.; Data Curation, M.H. and A.M.; Writing - Original Draft, H.M. and A.M.; Writing—review and editing, H.M. and A.M.; Visualization, M.H.,A.M. and M.Rz.; Supervision, M.Rz and K.K.; Project Administration, A.M. All authors have read and agreed to the published version of the manuscript.

Correspondence to Marcin Heronimczak.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Heronimczak, M., Mrowiec, A., Rząsa, M. et al. Measurements of the flow of a liquid–solid mixture/suspension through a segmented orifice. Sci Rep 14, 269 (2024). https://doi.org/10.1038/s41598-023-50737-6

Download citation

Received: 11 July 2023

Accepted: 24 December 2023

Published: 02 January 2024

DOI: https://doi.org/10.1038/s41598-023-50737-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative