Next Article in Journal
Stochastic Equations of Hydrodynamic Theory of Plasma
Next Article in Special Issue
Proposed Approach for Modelling the Thermodynamic Behaviour of Entrapped Air Pockets in Water Pipeline Start-Up
Previous Article in Journal
Advanced RBF Methods for Mapping Aerodynamic Loads onto Structures in High-Fidelity FSI Simulations
Previous Article in Special Issue
Flow Modeling of a Non-Newtonian Viscous Fluid in Elastic-Wall Microchannels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correction Factors for the Use of 1D Solution Methods for Dynamic Laminar Liquid Flow through Curved Tubes

Department of Mechanical Engineering, University of Saskatchewan, Saskatoon, SK S7N 5A9, Canada
Fluids 2024, 9(6), 138; https://doi.org/10.3390/fluids9060138
Submission received: 15 May 2024 / Revised: 29 May 2024 / Accepted: 1 June 2024 / Published: 6 June 2024
(This article belongs to the Special Issue Pipe Flow: Research and Applications)

Abstract

:
The modeling of transient flows of liquids through tubes is required for studies in water hammer, switched inertance hydraulic converters, and noise reduction in hydraulic equipment. While 3D gridded computational fluid dynamics (CFD) methods exist for the prediction of dynamic flows and pressures in these applications, they are computationally costly, and it is more common to use 1D methods such as the method of characteristics (MOC), transmission line method (TLM), or frequency domain methods. These 1D methods give good approximations of results but require many orders of magnitude less computation time. While these tubes are typically curved or coiled in practical applications, existing 1D solution methods assume straight tubes, often with unknown deviation from the curved tube solution. This paper uses CFD simulations to determine the correction factors that can be used for existing 1D methods with curved tubes. The paper also presents information that can be used to help evaluate the expected errors resulting from this approximation.

1. Introduction

Modeling of fluid flows through curving pipes has a long history, both as examples of theoretical calculus solution methods as well as for practical engineering. The interested reader is referred to two review papers [1,2] for details, but a brief overview is presented here. Consider the curved tube shown in Figure 1, with the radius of curvature defined as r b , the inner radius of the circular cross-section of the rigid tube as r, and the curve angle is θ . Boundary conditions include the pressure, p ( t ) , and volumetric flow rate, q ( t ) , at the inlet and outlet (denoted with subscripts A and B, respectively). The pressure is assumed to be constant across the inlet or outlet. Given the initial pressure and velocity fields within the tube and two of these four transient boundary conditions, the problem is to determine the other two boundary conditions.
Although others had previously noted the effect of curved passageways on an open channel flow and air through pipes, steady-state solutions to this problem typically start with the 1927 works of Dean [3], who established much of the theoretical groundwork and nomenclature used in the subsequent research, including this paper. He defined what would become known as the Dean Number, as
D = R e r r b   1 2
where R e is the Reynolds number and r / r b is the curvature ratio. This number characterizes the nature of the flow via the ratio of the square root of the product of the inertial and centrifugal forces to the viscous forces [2].
The steady-state flow solution results in a secondary vortex flow perpendicular to the main streamwise flow, with the strength and number of these vortices related to the Dean number. It has also been noted that these vortices can have a significant effect on the flow for some distance downstream of the pipe bend [2]. For the sake of this paper, we will not concern ourselves with the internal details of the flow pattern but will concentrate on the inlet and outlet flows.
Although theoretical solutions exist for some simplified situations, the most reliable relations come from empirical correlations. The early part of the 20th century resulted in much of the careful experimental data still in use today. This includes references [4,5], which resulted in Hasson’s empirical correlation for steady-state flow, formulated as
λ λ s 0.0969 2 D 1 2 + 0.556
where λ / λ s is the ratio of the steady-state resistance to the flow of a curved pipe relative to the resistance of a straight pipe with the same parameters ([6] via [2]). Note that the solution for the laminar flow through a straight pipe is the classical Poiseuille flow solution:
λ s = p A p B q A = p A p B q B = 8 π μ L π r 2
where μ is the fluid dynamic viscosity and L is the tube length. Refer to [1] for a more comprehensive list of suitable correlations for various situations. The above correlation is for laminar flow, noting that transition to turbulent flow has been observed to occur at higher Reynold numbers than for straight pipes.
While the steady-state solution for the flow through curved pipes is relatively well understood, the dynamic response due to wave propagation under the influence of viscous friction is much less studied. Starting with straight tubes, analytical solutions taking into account the frequency-dependent nature of viscous friction exist in the frequency domain [7,8,9]. These solutions can be calculated very quickly at a very low computational cost but are only applicable for problems where the entire system is linear and time is invariant. In situations where there are other parts of the system model that require nonlinear dynamic models (common in hydraulic system design), it is common to use either the method of characteristics (MOC) or the transmission line method (TLM). The method of characteristics [10,11,12] divides the tube into an evenly spaced 1D grid along its length. The flow and pressure at any grid point can then be calculated, moving forward through time with a time step selected so that the propagating wave moves exactly one spatial grid distance in each time step. This method is somewhat more computationally expensive than the frequency-domain method, but allows for parameters that change in time, as well as along the length of the tube, and also may be connected to arbitrary dynamic models at its inlet and outlet. However, the required equations have not yet been developed for curved pipelines, and it appears that this may not be theoretically possible.
Another method of modeling dynamic responses is using the transmission line method. Developed by a number of researchers over recent years [13,14,15], this method uses a four-port model (flow and pressure at inlet and outlet) and applies a number of parameterized linear transfer functions and delays to connect the four ports, ignoring the internal pressures and flows. The parameters of these transfer functions are selected to give the desired overall input–output characteristics, typically to minimize the error with respect to an analytical solution. This is an approximate method, but the errors can be small, and it has the advantages of being relatively computationally efficient and integrating well with dynamic modeling software packages such as Simulink [16], Hopsan [17], etc. This method was originally designed for simple, straight, rigid-walled tubes, but the author of the present paper has collaborated with others to expand the concept to allow for straight tapered tubes, tubes with elastic pipe walls with varying thickness, as well as arbitrary cross-sections [18,19].
There has been relatively little work carried out, either experimentally or theoretically, on the dynamic response of general flow through curved tubes. One series of experimental results is available [20,21,22], but the apparatus had a very large bend radius and was designed to accentuate fluid–structure interactions, so it is difficult to separate this effect from the purely fluid response. Some work has been performed on specific applications, such as blood flow through the circulatory system [23], but these results are sufficiently specialized with respect to the pressure signal limitations and fluid properties to be inapplicable in most other general applications.
Some recent works have also used 3D computational fluid dynamics (CFD) to solve the transient flow of fluids through curved and other complex geometries, such as [24,25,26,27,28]. Unfortunately, these methods can be unstable and require many orders of magnitude more computation time and resources to complete than 1D methods.
The main objective of the present paper is to invest time in CFD results to determine the correlation for correction factors that will enable the use of high-speed 1D solution methods for curved tubes. This paper will also present information related to the expected errors when using these methods with or without these correction factors.

2. Materials and Methods

The objective of this paper is to provide correction factors so that existing 1D solution methods for transient pipe flows may be used to approximate the flows in curved pipes. The solution for a straight, rigid pipe with a laminar flow of Newtonian fluid and with a constant wave speed can be entirely defined by three parameters: the characteristic impedance, dissipation number, and time scale. The shape of the overall response is controlled by the dissipation number, while the scaling of the magnitude of the flow–pressure relationship is set by the characteristic impedance, and the scaling in time is controlled by the time scale. In this section, CFD methods are used to establish the factors required to correct these values for use with flow through curved tubes.
For the purpose of this paper, assume a laminar flow through a smooth, rigid tube. This ensures no structure–fluid interactions or compliance in the tube wall and ensures a linear, time-invariant dynamic system. Assume that the sonic speed of the fluid is large relative to the fluid velocities and assume small changes in pressure, which will allow for the wave propagation solutions to be superimposed over a variety of mean flows and to be calculated for arbitrary boundary conditions. The pressure distribution is constrained to be constant across the inlet and outlet, but the pressure distribution within the tube is not constrained, nor is the velocity profile constrained at any point.

2.1. Characteristic Impedance

The characteristic impedance of a tube is the ratio of the magnitude of a pressure wave traveling through a tube to the magnitude of the resultant flow wave. For a lossless straight tube, this is related to the ratio of the inertia per unit length to the compliance per unit length and is defined as
Z c s = ρ c A
where ρ is the fluid density, c is the wave speed, and A = π r 2 is the cross-sectional area of the tube (a subscript “s” is used to refer to quantities related to a reference straight tube in order to differentiate them from an arbitrarily curved tube under consideration). Note that this value is only related to the properties of the fluid and of the cross-section; there is no dependence on the length of the tube. This impedance would be measurable at the inlet of a tube of infinite length or at the inlet of a tube of finite length for the initial part of the response prior to the influence of returning reflected waves. This is shown in Figure 2 for a Method of Characteristics solution for the inlet flow, q A , in a straight tube of finite length subjected to a step in inlet pressure of Δ p A . The ratio q A Z c / Δ p A is equal to 1 in the middle portion (around t = L / c ) of this plot for a lossless tube (small dissipation number β ) . A similar method is used to measure the characteristic impedance of curved tubes, using response calculated using a CFD simulation.
The CFD responses are calculated using the OpenFOAM CFD package [29]. The “sonicLiquidFoam” solver was used, which is a compressible flow solver intended for studies of wave propagation in liquids with laminar flows. It uses an equation of state with a density described by
ρ = ρ 0 + ψ p p 0
where subscripts 0 denote reference conditions and ψ = 1 / c 2 is the (constant) derivative of density with respect to pressure. The default finite volume discretization of Gaussian integration is used.
The tube was meshed with the grid, shown in Figure 3, for half of the pipe (symmetry along the central plane is exploited to reduce the solution size). Following a similar study by Fries [25,26,27], the pipe’s cross-section is divided into a central core with equal cell spacing, while the outer 2/3 of the radius is graded, with smaller cells toward the pipe walls to allow for better resolution of near-wall effects. In each cross-section, there are a total of 290 cells. This cross-section is then extruded in either a straight line of total length L s or revolved about an axis over an angle of θ (with a centerline length of L s = θ r b ). The cells in the axial direction are spaced with a 1 mm spacing along the center line of the tube (approximately equal to the spacing in the radial and circumferential directions at the edge of the central core). OpenFOAM’s data structures are only concerned with the connections between cells, which permits the situation of θ > 2 π , which is physically impossible as the rotations would overlap in space, but allows for a longer length of tube and therefore more accurate solutions that minimize inlet effects. This can also be thought of as a helical coil in the limit as the helical pitch approaches zero. The baseline parameters for simulations are shown in Table 1.
The above tube was initialized with zero velocity and an initial pressure of 100 kPa. Boundary conditions were applied with a constant pressure of 100 kPa at the outlet and 200 kPa at the inlet, corresponding to an inlet pressure step of Δ P A = 100 kPa. A no-slip boundary condition was applied to the pipe walls and a symmetry plane boundary condition along the central plane. It should be noted that one could instead apply a flow boundary condition to the inlet, which can improve solution stability but would require some other method to determine the non-parabolic velocity profile of the inlet velocity.
The solution was then solved with a time step of 5 × 10 7 s for a period of 5 ms. This was selected as a time that would allow for the sonic wave to traverse the length of the straight tube approximately five times. (The change in the effective length of the curved tube will change the number of times that the wave actually transits the tube).
The volumetric flow was calculated by integrating the normal velocity over the plane of the inlet and outlet (and doubled to take into account the other symmetric half of the tube). The characteristic impedance was then calculated as the inlet pressure step divided by the flow, measured at the first time that the outlet flow equals the inlet flow ( t 1 in Figure 4). The characteristic impedance correction factor is defined as
K Z c = Z c Z c s
where Z c s is the geometrical characteristic impedance (Equation (4)) and Z c is the impedance measured from the CFD results. The above was then repeated for a range of curve radii.
These calculations were performed on a desktop computer with two Xeon 4114 processors (a total 20 physical cores) with 128 GB of RAM. ESI OpenFOAM version v2312 was used, running on Ubuntu Linux 22.04.4.

2.2. Time Scale

The second parameter required for the use of 1D solution methods is the time scale L / c . In this case, it is assumed that the wave propagation speed, c, is constant and assumed that changes in the response are due to changes in the effective length L. This change can be explained by the fact that some wave energy may preferentially travel from inlet to outlet following the shorter path length along the inside of the curve.
The same CFD results from Section 2.1 were used to determine this effect. In this case, we identify the time it takes for a wave to travel the length of the tube, reflect off the open boundary, and return. As shown in Figure 4, for a sample waveform, t 1 is defined as the first time that the outlet flow, q B , exceeds the inlet flow, q A . t 2 is defined as the first time that the outlet flow exceeds the inlet flow after t = 2   L s / c . The effective length can then be calculated using
L = t 2 t 1 c 2
and the correction factor is defined as
K L = L L s
where L s = θ r b is the centerline length.

2.3. Dissipation Number

The dissipation number for a tube reflects the damping induced by the fluid viscosity. This is analogous to the damping ratio of a dynamic system: a dissipation number near unity will reach the steady-state value with little overshoot, while a system with a low dissipation number will take many oscillations before reaching steady-state, as shown in Figure 2, for MOC solutions for varying dissipation numbers.
For a straight tube, the dissipation number is defined as
β s = ν L c r 2
which is equal to
β = λ 8 Z c
where λ is the steady-state laminar resistance of the tube; this second equation applies to curved tubes as well as straight tubes.
The estimate for the effective dissipation number of a curved tube is determined by using the same CFD model as described in Section 2.1, but with the end time extended to 10   L s / c . As before, the dynamic flow response to a step in pressure is calculated. The dissipation for a method of characteristics (MOC) solution (available for download from [12]) was determined for a straight tube that best fits this response. The characteristic impedance and effective length relations developed in the previous sections were used to compensate for these effects. The iterative, non-derivative global optimization scheme fminsearch, provided with the Matlab 2023a software package, was used to minimize the error between the MOC straight tube response and the CFD curved tube response by varying the MOC dissipation number. The normalized root mean squared error between the straight tube MOC solution, q B s , and the curved tube CFD solution, q B , was minimized, and defined as
E = Z c s Δ P A 1 t e n d   0 t = t e n d q B q B s 2   d t 1 / 2
where t e n d = 10   L s / c . The effective dissipation number correction factor can then be defined as
K β = β β s .
where β is the effective dissipation number and β s is the dissipation number for a straight tube of the same centerline length (Equation (9)).

3. Results

3.1. Validation

In order to have some confidence in the CFD solution, the CFD solution for a straight pipe was calculated with parameters in Table 1, which can be compared to the validated method of characteristics and frequency domain solutions. These results are shown in Figure 5, along with varied viscosity solutions to show the effect of dissipation number. The results have good general agreement. One area with some differences is during the rapid change in flow around t / ( L / c ) = 1, 2, 3, and 4 as the wave fronts are reaching the reflective boundaries. This may be due to the assumption in the 1D solutions that the tube is narrow with no waves moving in the radial direction. Since the 3D CFD solution allows some wave energy to move diagonally through the pipe, some energy will arrive later, leading to the wave front’s arrival being less sharp.
A mesh invariance test was also performed to ensure that the mesh had been selected appropriately. The mesh size was halved in each dimension, along with halving the time step to ensure a constant Courant number. The results are shown in Figure 6. The Courant number was also checked by maintaining the original grid spacing while reducing the time step, as shown in Figure 6. Again, there is good agreement, and one can have confidence in the time and spatial mesh.
Finally, the linearity of the transient solution was verified by doubling the driving pressure step magnitude, Δ p A . As shown in Figure 7 for curved tubes, the resulting flow doubles, but the shape is nearly exactly the same, with deviations only noticeable for extreme bend ratios and low dissipation numbers. This validates the linearity assumption and demonstrates that the solution does not rely on the Dean number (as it would for static solutions). The nonlinear static solution would be important for very large pressure waves, but the results presented here can be viewed as accurate for small variations about an operating point.

3.2. Characteristic Impedance

Examples of inlet flow wave forms for varying bend radii are shown in Figure 8. Note that as r / r b 1 (more extreme bends), the magnitude of the flow between the steps increases, meaning the characteristic impedance decreases. There is also a change in the timing of the steps, which will be addressed in the next section. The associated calculated characteristic impedance correction factors are shown in Figure 9. Data for double the tube radius is also included, which is nearly identical, demonstrating that the curve for the characteristic impedance correction factor is a function of the bend radius ratio, not the radius itself. This curve can be approximated by a fit of
K Z c = Z c Z c s 1 0.2239 r r b 3.276 .
For better accuracy, a higher order polynomial could be fit to the data in Table A1.

3.3. Effective Length

As shown in Figure 8, the curve ratio clearly affects the arrival time of the steps in flow (after they have traveled the length of the tube). The calculated effective length correction factor is shown in Figure 10, which is very similar to the correction factor for characteristic impedance. The correlation was fit to a similar form:
K L = L L s 1 0.2095 r r b 2.965
It is unclear at present whether the difference in values in this equation should be assumed to be due to simulation error and taken as identical to Equation (13) or whether this difference exists in reality.
It should be noted that the K L correction factor is not equal to the change attributable to following the shortest path length through the inside of the curve, indicating that the wave traveling along this shortest path is “pulled back” by the neighboring waves traveling longer paths.

3.4. Dissipation Number

Examples of the inlet flow responses are shown in Figure 11, Figure 12 and Figure 13 for a range of bend ratios and dissipation numbers, along with the associated 1D MOC solutions with the best-fit dissipation number. These demonstrate a relatively good fit, especially for the parts of the response between the steps with moderate rates of change and also for moderate bend ratios and dissipation numbers. The more extreme bend ratios exhibit an oscillation in the volumetric flow as the wave fronts reflect at the boundary. There is no fluid experimental data available to validate this effect for large bend ratios, but these oscillations are visible in some data from curved microwave antennas [30], which are governed by similar wave propagation equations. This may be due to the nonplanar wavefronts that are generated in the strongly curved tubes. An example of this nonplanar wave front is shown in Figure 14, with the white contour line showing that the wave has progressed faster in the middle of the tube than along the outer bend radius. Figure 15 shows a plot of the pressure distribution along a radial line. This shows that the assumption of constant pressure across the cross-section that is commonly made for straight tubes is not applicable here. It also shows the fact that constructive interference has occurred with the maximum pressure considerably above the inlet’s driving pressure.
The calculated dissipation number correction factor is plotted in Figure 16 and Figure 17, along with the RMS error in Figure 18. As shown in Figure 17, the correction factor can be significant for strongly curved tubes, especially if the dissipation number is low. In most cases, the corrected response provides a usefully accurate estimation of the flow, although it does not capture the oscillations immediately after each step for strongly curved tubes.
To provide an indication of the relative computation time required for the CFD and MOC solutions, the typical time required to calculate a single response in this section was 1600 s for CFD and 190 ms for the MOC, a speedup of over 8000 times. In addition, the CFD fully utilized the 20 cores of the computer, while the MOC code was not parallelized and ran on a single core, with associated speed improvements expected, it was to run on multiple cores.
Figure 11. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 1.20.
Figure 11. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 1.20.
Fluids 09 00138 g011
Figure 12. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 1.50.
Figure 12. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 1.50.
Fluids 09 00138 g012
Figure 13. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 2.00.
Figure 13. CFD (solid lines) and best-fit MOC (dashed lines) solutions for r b / r = 2.00.
Fluids 09 00138 g013
Figure 14. Example of the pressure distribution on the central plane after the wave has traveled clockwise from the inlet for 15 μ s , for an extreme bend radius of r b / r = 1.2 and β s = 0.1 . Note the nonplanar wave front, denoted by the white pressure contour at 150 kPa. Also, note the overpressure above the inlet pressure of 200 kPa, plotted along the green line in Figure 18.
Figure 14. Example of the pressure distribution on the central plane after the wave has traveled clockwise from the inlet for 15 μ s , for an extreme bend radius of r b / r = 1.2 and β s = 0.1 . Note the nonplanar wave front, denoted by the white pressure contour at 150 kPa. Also, note the overpressure above the inlet pressure of 200 kPa, plotted along the green line in Figure 18.
Fluids 09 00138 g014
Figure 15. The plot of the pressure distribution along the green line in Figure 17, showing the non-constant pressure distribution across the cross-section, which exceeds the inlet pressure of 200 kPa.
Figure 15. The plot of the pressure distribution along the green line in Figure 17, showing the non-constant pressure distribution across the cross-section, which exceeds the inlet pressure of 200 kPa.
Fluids 09 00138 g015
Figure 16. Measured dissipation number relative to the straight pipe dissipation number. Each “X” indicates a CFD simulation location.
Figure 16. Measured dissipation number relative to the straight pipe dissipation number. Each “X” indicates a CFD simulation location.
Fluids 09 00138 g016
Figure 17. Calculated effective dissipation number correction factor (relative to the straight pipe dissipation number).
Figure 17. Calculated effective dissipation number correction factor (relative to the straight pipe dissipation number).
Fluids 09 00138 g017
Figure 18. Residual error between corrected MOC solution and CFD solution after optimizing the dissipation number.
Figure 18. Residual error between corrected MOC solution and CFD solution after optimizing the dissipation number.
Fluids 09 00138 g018

3.5. Overall Effect

One additional piece of useful information is how much of an overall effect a curved tube has, and under what conditions a straight tube approximation can be used without correction factors. In this case, one can calculate the error between the CFD-calculated responses for curved and straight tubes with the same cross-section and centerline length. The resultant normalized root mean squared (RMS) error is shown in Figure 19. As one can see from this plot, it is possible to make the straight tube simplifying assumption for more extreme bends if the dissipation number is large. For example, for zero viscosity, the relative error reaches 1% at r b / r   = 16.2, while for β s = 0.1 the same error can be achieved with the more extreme bend ratio of r b / r = 3.5. Another feature that can noted from this graph is that the dissipation number makes little difference for β s less than 1 × 10 4 , which is similar to the effect in straight tubes for most applications of interest.

4. Conclusions

The results in the previous section can be used to provide guidance with regard to a number of important questions that may be relevant to those employing 1D models to speed up transient flow calculations relative to CFD solutions. The first is whether the curvature of the tube in a particular application can be ignored and safely modeled as a straight tube. This can be evaluated by referencing the error calculated in Figure 19, quantifying this effect for a range of curve ratios and dissipation numbers.
If it is determined that this level of error is unacceptable, one can then use the presented information to obtain better approximations using correction factors. In particular, the characteristic impedance can be corrected using Equation (13), the effective length corrected using Equation (14), and the dissipation number corrected using Figure 14 or Figure 15.
While these results are applicable for small pressure changes and laminar flows, it is left to future work to establish the exact bounds of applicability with regard to these factors and establish additional methods of operating outside of these bounds. It is also left to future work to experimentally verify the presented results.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

cwave propagation velocity
DDean number
Enormalized root mean squared error
K L length correction factor
K Z c characteristic impedance correction factor
K β dissipation number correction factor
Llength of pipe (or effective length)—subscript “s” for a straight pipe
ppressure. Subscript A for inlet, B for outlet
qvolumetric flow—subscript A for inlet, B for outlet
rtube inner radius
r b bend radius
R e Reynolds number
ttime—subscripts 1 and 2 denote times in Figure 4
Z c characteristic impedance—subscript “s” for a straight tube with the same dimensions
β dissipation number—subscript “s” for a straight tube with the same dimensions
θ bend angle
λ resistance to flow—subscript “s” for a curved pipe with the same centerline length
µdynamic viscosity
vkinematic viscosity
ρ density
ψ derivative of density with respect to pressure

Appendix A. Results Tables

Table A1. Results for the characteristic impedance correction factor, as plotted in Figure 9.
Table A1. Results for the characteristic impedance correction factor, as plotted in Figure 9.
rb/rKZc
r = 10 mmr = 20 mm
1.010.77230.7689
1.10.84770.8473
1.20.88650.8860
1.30.91020.9098
1.40.92670.9263
1.50.93880.9386
1.60.94800.9478
1.80.96130.9612
20.97020.9702
2.40.98110.9813
2.70.98600.9864
30.98950.9899
Table A2. Results for the effective length correction factor, as plotted in Figure 10.
Table A2. Results for the effective length correction factor, as plotted in Figure 10.
rb/rKL
r = 10 mmr = 20 mm
1.01N/A0.7765
1.10.84700.8505
1.20.88450.8870
1.30.90800.9098
1.40.92400.9255
1.50.93600.9377
1.60.94520.9472
1.80.95820.9600
20.96700.9683
2.40.97770.9790
2.70.98280.9835
30.98630.9870

References

  1. Sigalotti, L.D.G.; Alvarado-Rodríguez, C.E.; Rendón, O. Fluid Flow in Helically Coiled Pipes. Fluids 2023, 8, 308. [Google Scholar] [CrossRef]
  2. Berger, S.A.; Talbot, L.; Yao, L.S. Flow in Curved Pipes. Annu. Rev. Fluid Mech. 1983, 15, 461–512. [Google Scholar] [CrossRef]
  3. Dean, W.R. XVI. Note on the Motion of Fluid in a Curved Pipe. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1927, 4, 208–223. [Google Scholar] [CrossRef]
  4. White, C.M.; Appleton, E.V. Streamline Flow through Curved Pipes. Proc. R. Soc. Lond. Ser. A Contain. Pap. A Math. Phys. Charact. 1929, 123, 645–663. [Google Scholar]
  5. Adler, M. Strömung in Gekrümmten Rohren. ZAMM-J. Appl. Math. Mech. Z. Angew. Math. Mech. 1934, 14, 257–275. [Google Scholar] [CrossRef]
  6. Hasson, D. Streamline Flow Resistance in Coils. Res. Corresp. 1955, 1, 81. [Google Scholar]
  7. Muto, T.; Kinoshita, Y.; Yoneda, R. Dynamic Response of Tapered Fluid Lines: 1st Report, Transfer Matrix and Frequency Response. Bull. JSME 1981, 24, 809–815. [Google Scholar] [CrossRef]
  8. Tahmeen, M.; Muto, T.; Yamada, H. Simulation of Dynamic Responses of Tapered Fluid Lines. JSME Int. J. Ser. B Fluids Therm. Eng. 2001, 44, 247–254. [Google Scholar] [CrossRef]
  9. Zielke, W. Frequency-Dependent Friction in Transient Pipe Flow. J. Basic Eng. 1968, 90, 109–115. [Google Scholar] [CrossRef]
  10. Vítkovskỳ, J.; Lambert, M.; Simpson, A.; Bergant, A. Advances in Unsteady Friction Modelling in Transient Pipe Flow. In Proceedings of the 8th International Conference on Pressure Surges, The Hague, The Netherlands, 12–14 April 2000; BHR Group: Bedford, UK, 2000; pp. 471–482. [Google Scholar]
  11. Trikha, A. An Efficient Method for Simulating Frequency-Dependent Friction in Transient Liquid Flow. J. Fluids Eng. 1975, 97, 97–105. [Google Scholar] [CrossRef]
  12. Wiens, T. Method of Characteristics Solver for Transmission Lines. 2020. Available online: https://github.com/tkw954/MOC_solver (accessed on 10 May 2024).
  13. Krus, P.; Weddfelt, K.; Palmberg, J.-O. Fast Pipeline Models for Simulation of Hydraulic Systems. J. Dyn. Syst. Meas. Control 1994, 116, 132–136. [Google Scholar] [CrossRef]
  14. Johnston, N. The Transmission Line Method for Modelling Laminar Flow of Liquid in Pipelines. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2012, 226, 586–597. [Google Scholar] [CrossRef]
  15. Johnston, N.; Pan, M.; Kudzma, S. An Enhanced Transmission Line Method for Modelling Laminar Flow of Liquid in Pipelines. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2014, 228, 193–206. [Google Scholar] [CrossRef]
  16. Wiens, T.; van der Buhs, J. Transmission Line Models. 2020. Available online: https://github.com/tkw954/usask_tlm (accessed on 10 May 2024).
  17. Braun, R.; Krus, P. Multi-Threaded Distributed System Simulations Using the Transmission Line Element Method. Simulation 2016, 92, 921–930. [Google Scholar] [CrossRef]
  18. Wiens, T. Transmission Line Modeling of Laminar Liquid Wave Propagation in Tapered Tubes. J. Fluids Eng. 2019, 141, 101103. [Google Scholar] [CrossRef]
  19. Wiens, T. Modeling Arbitrarily Shaped Liquid Pipelines Using a Segmented Transmission Line Model. J. Dyn. Syst. Meas. Control 2020, 142, 064502. [Google Scholar] [CrossRef]
  20. Ferras, D.; Manso, P.A.; Schleiss, A.J.; Covas, D.I.C. Hydraulic Transients in Straight and Coil Pipe Rigs. In Proceedings of the 12th International Conference on Pressure Surges, Dublin, Ireland, 18–20 November 2015. [Google Scholar]
  21. Ferras, D. Fluid-Structure-Interaction in Pipe Coils During Hydraulic Transients: Numerical And Experimental Analysis. In Proceedings of the 36th IAHR World Congress, The Hague, The Netherlands, 8 June–3 July 2015. [Google Scholar]
  22. Ferras, D.; Manso, P.A.; Covas, D.I.C.; Schleiss, A.J. Fluid–Structure Interaction in Pipe Coils during Hydraulic Transients. J. Hydraul. Res. 2017, 55, 491–505. [Google Scholar] [CrossRef]
  23. Gijsen, F.J.H.; Allanic, E. The Influence of the Non-Newtonian Properties of Blood on the Flow in Large Arteries: Unsteady Flow in a 903 Curved Tube. J. Biomech. 1999, 32, 705–713. [Google Scholar] [CrossRef] [PubMed]
  24. Azhdari, M.; Riasi, A.; Tazraei, P. Numerical Analysis of Fluid Hammer in Helical Pipes Considering Non-Newtonian Fluids. Int. J. Press. Vessel. Pip. 2020, 181, 104068. [Google Scholar] [CrossRef]
  25. Fries, C.; Manhartsgruber, B. Dynamics of Transmission Line Junctions: Comparison of CFD Results Against Measurements. In Proceedings of the 8th FPNI Ph.D Symposium on Fluid Power, Lappeenranta, Finland, 11–13 June 2014; American Society of Mechanical Engineers: Lappeenranta, Finland, 2014; p. V001T02A003. [Google Scholar]
  26. Fries, C.; Manhartsgruber, B. CFD Based Modeling of Wave Propagation in Liquid Transmission Line Junctions. In Volume 1A, Symposia: Advances in Fluids Engineering Education; Advances in Numerical Modeling for Turbomachinery Flow Optimization; Applications in CFD; Bio-Inspired Fluid Mechanics; CFD Verification and Validation; Development and Applications of Immersed Boundary Methods; DNS, LES, and Hybrid RANS/LES Methods; American Society of Mechanical Engineers: Incline Village, NV, USA, 2013; p. V01AT03A025. [Google Scholar]
  27. Fries, C. Coupling of Liquid Transmission Line Models with 3-Dimensional CFD in OpenFOAM. Ph.D. Thesis, Johannes Kepler Univsitat Linz, Linz, Austria, 2016. [Google Scholar]
  28. Riedelmeier, S.; Becker, S.; Schlücker, E. 3D CFD Simulation of Water Hammer Through a 90° Bend: Applicability of URANS, 3D Effects and Unsteady Friction; Volume 4: Fluid-Structure Interaction, American Society of Mechanical Engineers: Anaheim, CA, USA, 2014; p. V004T04A094. [Google Scholar]
  29. Chen, G.; Xiong, Q.; Morris, P.J.; Paterson, E.G.; Sergeev, A.; Wang, Y. OpenFOAM for Computational Fluid Dynamics. Not. AMS 2014, 61, 354–363. [Google Scholar] [CrossRef]
  30. Thomann, W.; Isele, B.; Russer, P. Characterization of a 90 Degrees Curved Microstrip Transmission Line in the Time-Domain and Frequency-Domain with 3D-TLM-Method and Measurements. In Proceedings of the IEEE Antennas and Propagation Society International Symposium 1992 Digest, Chicago, IL, USA, 18–25 June 1992; Volume 2, pp. 658–661. [Google Scholar]
Figure 1. Geometry and nomenclature for a curved tube, showing the curved tube in plan view and also the cross-section. The radius of the cross-section is defined as r , the bend radius as r b , and the bend angle as θ .
Figure 1. Geometry and nomenclature for a curved tube, showing the curved tube in plan view and also the cross-section. The radius of the cross-section is defined as r , the bend radius as r b , and the bend angle as θ .
Fluids 09 00138 g001
Figure 2. Inlet flow responses to a step in inlet pressure for a straight tube with an open (constant pressure) outlet for varying dissipation numbers. The (left) and (right) plots show the same data at different time-scale zoom levels. For straight tubes, this response is totally defined by the dissipation number, β , the characteristic impedance, Z c , and the time scale L / c , where L is the length and c is the wave speed.
Figure 2. Inlet flow responses to a step in inlet pressure for a straight tube with an open (constant pressure) outlet for varying dissipation numbers. The (left) and (right) plots show the same data at different time-scale zoom levels. For straight tubes, this response is totally defined by the dissipation number, β , the characteristic impedance, Z c , and the time scale L / c , where L is the length and c is the wave speed.
Fluids 09 00138 g002
Figure 3. Example meshing for a portion of a tube with an extreme bend radius. Note the inner core with equal grid spacing and the outer section with cells graded toward the outer wall. A symmetry plane bisects the tube.
Figure 3. Example meshing for a portion of a tube with an extreme bend radius. Note the inner core with equal grid spacing and the outer section with cells graded toward the outer wall. A symmetry plane bisects the tube.
Fluids 09 00138 g003
Figure 4. Sample CFD flow responses for inlet ( q A ) and outlet ( q B ), showing points used to measure K Z c and K L .
Figure 4. Sample CFD flow responses for inlet ( q A ) and outlet ( q B ), showing points used to measure K Z c and K L .
Fluids 09 00138 g004
Figure 5. Comparison of flow responses (inlet flow at left and outlet at right) for straight tubes for CFD (solid lines) and MOC (dashed lines).
Figure 5. Comparison of flow responses (inlet flow at left and outlet at right) for straight tubes for CFD (solid lines) and MOC (dashed lines).
Fluids 09 00138 g005
Figure 6. Mesh invariance test (left), showing little effect when the spatial mesh spacing is halved in each dimension while also halving the time step to maintain the same Courant number. (Right) the minimal effect of the Courant number is shown by reducing the time step (while maintaining a consistent spatial mesh).
Figure 6. Mesh invariance test (left), showing little effect when the spatial mesh spacing is halved in each dimension while also halving the time step to maintain the same Courant number. (Right) the minimal effect of the Courant number is shown by reducing the time step (while maintaining a consistent spatial mesh).
Fluids 09 00138 g006
Figure 7. Linearity test, showing the CFD flow response for a variety of bend ratios and dissipation numbers for the base case of Δ p A = 100 kPa (black), and for the case of Δ p A = 200 kPa (colored). The results are very similar and overlap for almost all cases. The only visible difference is for extreme bend ratios and small dissipation numbers (top right of the plot). The error between the two solutions, E, is calculated using Equation (11).
Figure 7. Linearity test, showing the CFD flow response for a variety of bend ratios and dissipation numbers for the base case of Δ p A = 100 kPa (black), and for the case of Δ p A = 200 kPa (colored). The results are very similar and overlap for almost all cases. The only visible difference is for extreme bend ratios and small dissipation numbers (top right of the plot). The error between the two solutions, E, is calculated using Equation (11).
Fluids 09 00138 g007
Figure 8. Flow responses to a step in inlet pressure for β = 0 , while varying bend ratio.
Figure 8. Flow responses to a step in inlet pressure for β = 0 , while varying bend ratio.
Fluids 09 00138 g008
Figure 9. Characteristic impedance correction factor for a curved tube relative to a straight tube of the same dimensions and wave speed. The same data is plotted in each graph, with the exponent applied to the left plot’s x-axis to show linearity.
Figure 9. Characteristic impedance correction factor for a curved tube relative to a straight tube of the same dimensions and wave speed. The same data is plotted in each graph, with the exponent applied to the left plot’s x-axis to show linearity.
Fluids 09 00138 g009
Figure 10. Effective length correction factor for a curved tube relative to a straight tube of the same dimensions. This also shows the same K Z c data as Figure 9 for comparison.
Figure 10. Effective length correction factor for a curved tube relative to a straight tube of the same dimensions. This also shows the same K Z c data as Figure 9 for comparison.
Fluids 09 00138 g010
Figure 19. Error between the straight and curved CFD flow solutions for a range of dissipation numbers and bend ratios.
Figure 19. Error between the straight and curved CFD flow solutions for a range of dissipation numbers and bend ratios.
Fluids 09 00138 g019
Table 1. Baseline simulation parameters.
Table 1. Baseline simulation parameters.
rTube inner radius10 mm
rbBend radius20 mm
LsCenterline length1 m
µDynamic viscosity1 Pa s
ρ Density1000 kg/m3
cWave propagation velocity1000 m/s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wiens, T. Correction Factors for the Use of 1D Solution Methods for Dynamic Laminar Liquid Flow through Curved Tubes. Fluids 2024, 9, 138. https://doi.org/10.3390/fluids9060138

AMA Style

Wiens T. Correction Factors for the Use of 1D Solution Methods for Dynamic Laminar Liquid Flow through Curved Tubes. Fluids. 2024; 9(6):138. https://doi.org/10.3390/fluids9060138

Chicago/Turabian Style

Wiens, Travis. 2024. "Correction Factors for the Use of 1D Solution Methods for Dynamic Laminar Liquid Flow through Curved Tubes" Fluids 9, no. 6: 138. https://doi.org/10.3390/fluids9060138

APA Style

Wiens, T. (2024). Correction Factors for the Use of 1D Solution Methods for Dynamic Laminar Liquid Flow through Curved Tubes. Fluids, 9(6), 138. https://doi.org/10.3390/fluids9060138

Article Metrics

Back to TopTop