Numerical Assessment of Shear Boundary Layer Formation in Sewer Systems with Fluid-Sediment Phases

: Numerical and empirical studies of soil slurry transport and deposition in urban stormwater sewers are few, presumably due to the di ﬃ culty of direct observation of soil slurry ﬂow in stormwater pipes. Slurry in a sewer system includes both suspended load and bedload, but few studies have attempted to demarcate these two components. A boundary layer is a crucial determinant of sediment transport capacity. Stormwater runo ﬀ enters the sewer in turbulent ﬂow, mostly mixed with soil slurry generated by rainfall. In this paper, we attempt analysis using ANSYS Fluent commercial CFD software. We describe the development of a numerical analytical methodology capable of predicting the ﬂow of soil slurry in stormwater pipes, and propose a method for estimating the sediment–ﬂow boundary layer. Using this model, we simulated stormwater runo ﬀ with a large content of soil slurry during a rainfall event. We investigated soil slurry transport and predict the formation of shear boundary layer by varying the inlet conditions (volume of soil slurry entering the stormwater sewer system) and by analyzing the ﬂow velocity ﬁeld and soil slurry volume fraction in the pipes under various experimental ﬂow conditions. Based on the shear and settling velocity of sediment particles, we propose criteria for the formation of a shear boundary layer in stormwater pipes.


Introduction
South Korea's total annual precipitation is high, but highly concentrated in the summer months, with storms producing excessively concentrated local rainfall events [1][2][3]. Heavy rains in the summer transport large amounts of soil and floating debris, such as twigs, driftwood, and pebbles, which can block drainage pipes and thus cause landslides or floods. Mountainous regions affected by high rates of sediment transport often suffer damage caused by sedimentation of soils and debris in drainage pipes. Sedimentation exceeding the rate of subsidence results in overflow, increasing the risk of damage due to landslides or floods [4,5]. To prevent such damage, drainage pipe design must be informed by an understanding of flow characteristics, taking account of the sedimentation of soil and debris in the drainage network. Slurry transport through the pipeline causes various problems; for example, deposition of slurry solids lowers the performance and efficiency of the entire drainage system [6][7][8][9].
During a heavy rainfall event, soil slurry is detached from the base surface by rainwater and flows with it into the drainage system. Debris, including soil sludge, undermines the discharge capacity of a in different conditions of flow velocity and soil slurry volume, and, based on the interrelationship between the latter, we developed a numerical analysis model that simulates sedimentation of soil slurry in stormwater pipes.

Numerical Analysis
Two-phase flows involving solid particles can cause an unstable nature, promoting mixing and energy dissipation which can cause transport problems. The flow properties of unsteady flow are very important for high concentration slurry transportation. There are two main approaches to numerical modeling that account for these unstable flows [33][34][35], such as a fluid flow containing solid particles: the Eulerian approach, in which solid particles are assumed to be a continuous, fluid-like medium; and the Lagrangian approach, in which Newton's law of motion is applied to solid particles. In the Eulerian approach, the different phases treated as a continuum can be further classified into a single-fluid model that applies continuity equations, momentum equations, and particle component conservation equations by handling a mixture of soil and water as a fluid, and a two-phase separation model that applies the laws of conservation to the soil particles and the fluid. The concept of phase volume fraction was introduced, since the volume of phases cannot be occupied by other phases [31][32][33]. In this study, a two-phase flow in which fluid and soil are introduced was assumed. It was assumed that in the two-phase flow in which fluid and soil were introduced, flow conditions did not change and no deformed flow occurred. In this process, the full condition was considered and the sediment transport characteristics of the pressure pipeline were examined. This method can reduce calculation time and carry out an optimal solution through repeated calculations [36][37][38][39].
In this study, we carried out numerical analysis using ANSYS Fluent Workbench 13.0 [40]. ANSYS Fluent (FLUENT Inc.) is a program for modeling fluid flow using two-dimensional (2D) and three-dimensional (3D) finite-volume numerical procedures based on a graphical user interface (GUI) for hydrodynamics (fluid dynamics, heat transfer), and chemical reactions. Capable of performing comprehensive tasks including preprocessing, calculation, and postprocessing required for computational fluid dynamics (CFD) analysis, ANSYS Fluent comprises the Gambit (preprocessing), Tgrid (fluent meshing), and FLUENT (calculation and postprocessing) programs. The aforementioned 2D numerical models are unstructured grid-calculation programs that support tetrahedron, hexahedron, prism, pyramid, and polyhedral grids using the finite-volume method (FVM) to perform modeling of various complex shapes, and to provide numerical analysis for field-specific applications. In this study, boundary conditions of inlet velocity were taken as the inlet condition for soil and water while the interior and outflow boundary conditions were employed as mixture phase. During this process, the flow analysis considering the mixture was modeled in 2D. The numerical results are represented as contours and vectors for the soil-water discharges cases. As mentioned above, the 2D interior and outlet soil-water velocities were calculated from the 3D experimental area.
This study used the volume-of-fluid (VOF) method to analyze the behavior of multiphase fluids dealing with the pressurized pipe flow. For the calculation of turbulent flow, turbulence models such as the k-ε model and the large eddy simulation (LES) model were used to calculate water level and flow velocity for each fluid volume of a finite-volume network of water bodies such as ports, estuaries, rivers, or sewers, by simulating both steady and unsteady flow conditions [41].
The Eulerian-Eulerian method obtains a flow rate field and a soil particle distribution field by solving the continuity Equations (1) and (2), the momentum Equations (3) the standard k-ε turbulence equation. A conservative equation was derived at each step of this processing by combining all phases. Interphase bond was achieved through pressure and momentum exchange coefficients between the phases. Therefore, even in two-phase conditions, the sum of the volume concentrations must where x i (i = 1, 2, 3) represents the orthogonal coordinates, u i the orthogonal (directional) element of the velocity vector that represents average flow velocity, ρ is fluid density, p is pressure, and F i is the sum of influential forces, such as gravity, deflecting force, and centrifugal force. For each dispersed phase to be expressed as a continuum, a high-volume concentration is required. The method lends itself well to flow analysis of phases with high density or concentration [42][43][44][45].
In Equation (4), τ ij represents a viscous stress tensor representing Newtonian flow fluid, including the velocity of the deformed stress tensor; and i j represents structural correlation [46][47][48]: The k-ε turbulence model, which we selected for the purpose of this study, is based on the Boussinesq assumption that the Reynolds stress generated during time integration of Navier-Stokes equations is proportional to the gradient of the mean flow rate. This assumption takes into account the eddy viscosity coefficient (u T ) in Equations (5)- (7): where C u corresponds to the empirical constant, and k and ε are the turbulent kinetic energy and dissipation rate. Respectively. These can be calculated and using the transport Equations (5)- (7), where C u , σ k , σ ε , C 1ε , C 2ε , and C 3ε are empirical constants [41,49].

Analysis Model Validation
For 3D model analysis, configuration of the specifications of the sewer system to be analyzed should be preceded by the validation of the model to be developed, tailored to the analysis target. This is because validation with only the original model is insufficient, since it makes correcting of errors in case of upscaling labor-intensive and time-consuming. In this study, we validated a model constructed on the basis of earlier research in order to run a downscaled experimental model and reflect on the results when errors were detected.
We used the specifications determined in an experimental study with soil slurry particles measuring 1.4 to 2.0 mm by Matousek [50] to compare with the results of our numerical analysis. We validated our model by applying the same diameter, density, and composition ratio of water and soil particles of the soil slurry, and the specifications and gradients of the sewer system as that earlier study. Regarding particle-size distribution, our analysis was performed by injecting the same quantity of soil particles sized 1.4, 1.53, 1.68, 1.84, and 2.0 mm. For verification, we identified the concentration profiles conforming to those obtained experimentally from measurements at the deposition sites in each section Water 2020, 12, 1332 5 of 15 of the sewer system. We then checked the time of sediment deposition over the entire section of the sewer system, and identified the sediment distribution pattern.
Song et al. [32] also compared the results of numerical analysis of their study in a sewer system with these earlier experimental results [50], using the same flow conditions (volume concentration: 26%; flow velocity: 2.5-3.0 m/s), with soil particles sized 1.4-2.0 mm. The results of this comparison are presented in Figure 1. Extending this experiment, the present study added a flow velocity condition of 3.5 m/s ( Figure 2). The patterns of our modeling results are generally in good agreement with the experimental results. The turbulence model for our numerical analysis was validated by the observation that the volume concentration-dependent internal distribution profiles showed the same tendency in both flow velocity conditions.

Numerical Analysis Conditions
The characteristics, and transport and deposition processes of sediments inside stormwater pipes, are distinct from those of sediments on riverbeds. Sediments on the riverbed contain sand and gravel particles of different sizes. Sediments flowing into a stormwater sewer through a manhole, by contrast, contain a limited range of particle sizes, and consequently show no significant variance in

Numerical Analysis Conditions
The characteristics, and transport and deposition processes of sediments inside stormwater pipes, are distinct from those of sediments on riverbeds. Sediments on the riverbed contain sand and gravel particles of different sizes. Sediments flowing into a stormwater sewer through a manhole, by contrast, contain a limited range of particle sizes, and consequently show no significant variance in Figure 2. Comparison of modeling verification results (this study) with experimental data [50], for flow velocity of 3.5 m/s.

Numerical Analysis Conditions
The characteristics, and transport and deposition processes of sediments inside stormwater pipes, are distinct from those of sediments on riverbeds. Sediments on the riverbed contain sand and gravel particles of different sizes. Sediments flowing into a stormwater sewer through a manhole, by contrast, contain a limited range of particle sizes, and consequently show no significant variance in density and unit weight values. Sediments in a sewer system thus have more uniform characteristics compared with those on the riverbed.
As illustrated in Figure 3, we simulated the soil slurry transport and sedimentation processes by configuring a stormwater sewer measuring 0.6 m by 10 m and consisting of 140,000 grid meshes; the Cartesian coordinate system under pipe flow was governed by the control equation. We conducted a 2D-model analysis, which had a high computation efficiency under the same inlet conditions. Using the inlet velocity of soil slurry and volume concentration and particle size of the soil as variables, we compared and analyzed the flow velocity and particle distribution profiles of the slurry mixtures in a total of 63 conditions. The calculation was performed up to a total of 500 s by applying the calculation time interval ∆t = 1 s, and the calculation was set to repeat 10 times the internal calculation (step length factor). Using very small values in repetitions of the computation could have led to a termination of the analysis during particle transport and deposition, the intervals used in the study reflected sufficiently leisurely to avoid this result. It was confirmed that the results did not change during the course of the analysis and that a constant value was derived. Table 1 presents the material properties of the input mixture of fluid and soil slurry and the boundary conditions of our numerical analysis. Table 1. Input and boundary conditions for the numerical analysis in this study. density and unit weight values. Sediments in a sewer system thus have more uniform characteristics compared with those on the riverbed. As illustrated in Figure 3, we simulated the soil slurry transport and sedimentation processes by configuring a stormwater sewer measuring 0.6 m by 10 m and consisting of 140,000 grid meshes; the Cartesian coordinate system under pipe flow was governed by the control equation. We conducted a 2D-model analysis, which had a high computation efficiency under the same inlet conditions. Using the inlet velocity of soil slurry and volume concentration and particle size of the soil as variables, we compared and analyzed the flow velocity and particle distribution profiles of the slurry mixtures in a total of 63 conditions. The calculation was performed up to a total of 500 s by applying the calculation time interval Δt = 1 s, and the calculation was set to repeat 10 times the internal calculation (step length factor). Using very small values in repetitions of the computation could have led to a termination of the analysis during particle transport and deposition, the intervals used in the study reflected sufficiently leisurely to avoid this result. It was confirmed that the results did not change during the course of the analysis and that a constant value was derived. Table 1 presents the material properties of the input mixture of fluid and soil slurry and the boundary conditions of our numerical analysis.

Flow Characteristics of a Soil Slurry Mixture in the Stormwater Sewer System
For analysis of multiphase flows including fluid-soil phases, the total pressure drop in the pipe is basically related to four main components: frictional, hydrostatic, acceleration, and pressure drop. In many applications, the calculation of the pressure drop is the most complex component. The deposition of a sediment layer on the bottom of the pipe is a major cause of formation of a shear boundary layer in the soil slurry transport process, and thus serves as a major reference point for

Flow Characteristics of a Soil Slurry Mixture in the Stormwater Sewer System
For analysis of multiphase flows including fluid-soil phases, the total pressure drop in the pipe is basically related to four main components: frictional, hydrostatic, acceleration, and pressure drop. In many applications, the calculation of the pressure drop is the most complex component. The deposition of a sediment layer on the bottom of the pipe is a major cause of formation of a shear boundary layer in the soil slurry transport process, and thus serves as a major reference point for Water 2020, 12, 1332 7 of 15 facility maintenance. Friction velocity is a key parameter used in studies of erosion in waterbodies (streams, estuaries, etc.) and airborne particle transport. Friction velocity, defined as the square root of the ratio between wall shear stress over fluid density; it is decided by near wall flow conditions. Boundary layers are classified according to the relative value of shear velocity, and a sediment layer is formed by particles having a shear force exceeding that value. A major factor influencing the process of soil slurry transport and deposition is flow resistance. Sediments settle on the bottom of the pipe and thus increase the flow resistance and decrease the area of discharge, as shown in Figure 4.
Water 2020, 12, x FOR PEER REVIEW 7 of 16 facility maintenance. Friction velocity is a key parameter used in studies of erosion in waterbodies (streams, estuaries, etc.) and airborne particle transport. Friction velocity, defined as the square root of the ratio between wall shear stress over fluid density; it is decided by near wall flow conditions. Boundary layers are classified according to the relative value of shear velocity, and a sediment layer is formed by particles having a shear force exceeding that value. A major factor influencing the process of soil slurry transport and deposition is flow resistance. Sediments settle on the bottom of the pipe and thus increase the flow resistance and decrease the area of discharge, as shown in Figure  4. Earlier works have noted that the migration of sediments in a stormwater sewer occurs in three basic stages: entrainment, transport, and deposition [51,52]. Figure 5 presents the shear characteristics for each stage. The thickness of boundary layers is not uniform, and most of the flow segments are zones of turbulent flow. The demarcation of the boundary layers is determined by the relative magnitudes of shear velocity (u * ) and settling velocity (v ) of particles on the pipe bottom [53], as shown in Figure 4. Table 2. Flow conditions and sediment transport status [53].

Condition
Sediment transport status v /u * < 0.6 ∎ Particle travel in suspension 0.6 < v /u * < 2.0 ∎ Particle travel in saltation v /u * > 2.0 ∎ Particle travel as bedload Putting Figure 5 and Table 2 in context, different layers can be defined in the context of slurry transport and sedimentation: the term "bedload layer" refers to a layer on which interparticle adhesion force acts as a force to keep the particles together on the bottom (bed) of the pipe, making Earlier works have noted that the migration of sediments in a stormwater sewer occurs in three basic stages: entrainment, transport, and deposition [51,52]. Figure 5 presents the shear characteristics for each stage. The thickness of boundary layers is not uniform, and most of the flow segments are zones of turbulent flow.
Water 2020, 12, x FOR PEER REVIEW 7 of 16 facility maintenance. Friction velocity is a key parameter used in studies of erosion in waterbodies (streams, estuaries, etc.) and airborne particle transport. Friction velocity, defined as the square root of the ratio between wall shear stress over fluid density; it is decided by near wall flow conditions. Boundary layers are classified according to the relative value of shear velocity, and a sediment layer is formed by particles having a shear force exceeding that value. A major factor influencing the process of soil slurry transport and deposition is flow resistance. Sediments settle on the bottom of the pipe and thus increase the flow resistance and decrease the area of discharge, as shown in Figure  4. Earlier works have noted that the migration of sediments in a stormwater sewer occurs in three basic stages: entrainment, transport, and deposition [51,52]. Figure 5 presents the shear characteristics for each stage. The thickness of boundary layers is not uniform, and most of the flow segments are zones of turbulent flow. The demarcation of the boundary layers is determined by the relative magnitudes of shear velocity (u * ) and settling velocity (v ) of particles on the pipe bottom [53], as shown in Figure 4. Table 2. Flow conditions and sediment transport status [53].

Condition
Sediment transport status v /u * < 0.6 ∎ Particle travel in suspension 0.6 < v /u * < 2.0 ∎ Particle travel in saltation v /u * > 2.0 ∎ Particle travel as bedload Putting Figure 5 and Table 2 in context, different layers can be defined in the context of slurry transport and sedimentation: the term "bedload layer" refers to a layer on which interparticle adhesion force acts as a force to keep the particles together on the bottom (bed) of the pipe, making The demarcation of the boundary layers is determined by the relative magnitudes of shear velocity (u * ) and settling velocity (v s ) of particles on the pipe bottom [53], as shown in Figure 4.
Putting Figure 5 and Table 2 in context, different layers can be defined in the context of slurry transport and sedimentation: the term "bedload layer" refers to a layer on which interparticle adhesion force acts as a force to keep the particles together on the bottom (bed) of the pipe, making them migrate together; the term "transition layer" refers to a boundary layer atop an existing sediment layer, on which deposited particles are separated from the sediment layer under high turbulence stress and move along with the flow; and the "wash load layer" and "suspended load layer" are subjected only to turbulence stress during transport. Table 2. Flow conditions and sediment transport status [53].

Condition
Sediment Transport Status v s /u * < 0.6 Particle travel in suspension 0.6 < v s /u * < 2.0 Particle travel in saltation v s /u * > 2.0 Particle travel as bedload Since the frictional resistance to the flow is minimal, the particles contained in soil slurry have a very low settling velocity. Therefore, turbulence stress exists as a complete turbulent flow. Based on the flow direction of the sewer pipe, shear stress decreases linearly towards the upper end, where flow velocity becomes 0 and intense mixing between fluid and solid particles occurs. Despite the presence of vortices or eddies, such sections are characterized by steady flow velocity. As discussed above, several flow boundary layers emerge in the process of soil slurry transport, but it needs to be pointed out that clear demarcation between these flow systems is impossible.
Shear stress on the bottom of the pipe are the forces that carry and transport sediments. The important ratio in the mechanics of the shear boundary layer is the flow on the bed per unit area. The threshold value (critical boundary shear stress or critical erosion rate) is the value of the shear stress (threshold of motion) below which no entrainment can take place. Bottom (or bed) roughness is determined by the size of sediment particles accumulated on the bottom; it affects the flow rate distribution and sediment transport capacity in the pipe.
In summary, settling velocity affects sediment transport by exerting hydrodynamic force on the sediment particles, and thus causing the formation of boundary layers. With increase in the flow velocity in turbulent flow conditions, the force acting on the sediment particles exceeds the threshold value, and particles deposited on the bottom of the pipe reach the level of re-entrainment at a constant speed. This suggests that, in the critical condition, shear stress in the pipe dominates the sediment flow, showing high dependence on flow velocity.

Analysis of the Flow Boundary Layer in Stormwater Pipes
In the case of rivers, analysis of discharge capacity can be based on various observations. This enables estimation of the influence of hydraulic parameters. Most of the previous studies of sedimentation have focused on riverbed data, and settling velocity has been used only as a concept associated with sediment transport. In the case of stormwater drainage, however, it is unlikely that the characteristics of the inlet soil slurry, that determine settling velocity, remain constant. It is reasonable to assume the value of each of these characteristic varies. Likewise, the application of a uniform settling velocity is unlikely to properly reflect the phenomenon of sediment transport in a stormwater sewer system.
The main characteristics that determine the boundary layer during the flow of turbid water in sewer pipes are the vertical distribution and the boundary layer of shear velocity on the bottom of the pipe. In this context, by investigating the shear velocity due to the sediment layer, we can measure the change in the shear velocity on the bottom of the pipe. Shear flow velocity and shear stress at the bottom represent the force accompanying and transferring the deposits. The shear stress along the critical boundary or the critical erosion speed is the limiting condition of deposit transfer and shear stress (limiting point of motion) at the start of the yield. The value of roughness at the bottom depends on the size of sediment particles there, and influences changes in the flow velocity distribution and the transfer capability of the soil slurry [17,27,32]. The shear flow velocity (also known as the friction velocity or shear-stress velocity), u * , of a flow is defined by the relation Equation (8), where τ is the shear stress in an arbitrary layer of fluid and ρ is the density of the fluid. This quantity has the dimensions of velocity and is frequently encountered in the study of boundary layers. For turbulent flow it is approximately constant in the region near to a wall. The flow velocity in the boundary layer is commonly expressed in terms of the friction velocity and the height, z When the flow velocity distribution on a given cross-section is known, shear velocity can be estimated using Equation (9), based on the Law of the Wall, widely used for a relatively thin layer (z/H < 0.2) close to the bottom: Water 2020, 12, 1332 where u is the mean flow velocity at the upper distance z relative to the bottom, z 0 is the distance from the bottom where flow velocity tends to 0, u * is the shear velocity, and κ is von Karman's constant (~0.41) [54,55]. Accordingly, the shear velocity and shear stress during transport and deposition can be calculated by considering the change caused by the continuous increase in velocity occurring on the boundary layer on the bottom of the pipe. The shear force is the frictional resistance generated along the pipe wall during flow. Therefore, by applying the change owing to an increase in continuous velocity occurring along the boundary of the bottom of the conduit, the shear flow velocity and shear stress in transfer deposition can be calculated. If y increases discontinuously at a certain calculation point, ∂u/∂z in Equation (10) increases infinitely [56][57][58][59][60][61]. This frictional resistance is proportional to the flow velocity and inversely proportional to the cross-sectional area of the pipe. The shearing stress operating on the pipe plane at height z in a flow depth h and θ in the channel slope is given. Expressed as the shear rate or shear stress, τ, it can be calculated from Equation (10) in Newtonian fluid [49,58]: where µ is the dynamic viscosity of the fluid. If y increases discontinuously at a given calculation point, however, ∂u/∂z in Equation (10) increases ad infinitum. The shear stress at the point where the change in flow velocity tapers off can be defined as the limiting shear stress. Figure 6 shows the results of numerical analysis of soil slurry with 0.5 mm particle size. Figure 7 below illustrates the shear velocity distribution profiles during the transport and deposition of soil slurry with particle diameter of 0.5 mm. It shows that the change in shear velocity has converged values on the bottom of the pipe, and that shear velocity does not change beyond a certain distance from the wall, demonstrating that shear friction exerts a greater force on the bottom, but does not extend to the pipe center.
Water 2020, 12, x FOR PEER REVIEW 9 of 16 u u * = 1 κ ln z z (9) where u is the mean flow velocity at the upper distance z relative to the bottom, z is the distance from the bottom where flow velocity tends to 0, u * is the shear velocity, and κ is von Karman's constant (~0.41) [54,55]. Accordingly, the shear velocity and shear stress during transport and deposition can be calculated by considering the change caused by the continuous increase in velocity occurring on the boundary layer on the bottom of the pipe. The shear force is the frictional resistance generated along the pipe wall during flow. Therefore, by applying the change owing to an increase in continuous velocity occurring along the boundary of the bottom of the conduit, the shear flow velocity and shear stress in transfer deposition can be calculated. If y increases discontinuously at a certain calculation point, ∂u/ ∂z in Equation (10) increases infinitely [56][57][58][59][60][61]. This frictional resistance is proportional to the flow velocity and inversely proportional to the cross-sectional area of the pipe. The shearing stress operating on the pipe plane at height z in a flow depth h and θ in the channel slope is given. Expressed as the shear rate or shear stress, τ, it can be calculated from Equation (10) in Newtonian fluid [49,58]: where µ is the dynamic viscosity of the fluid. If y increases discontinuously at a given calculation point, however, ∂u/ ∂z in Equation (10) increases ad infinitum. The shear stress at the point where the change in flow velocity tapers off can be defined as the limiting shear stress. Figure 6 shows the results of numerical analysis of soil slurry with 0.5 mm particle size. Figure  7 below illustrates the shear velocity distribution profiles during the transport and deposition of soil slurry with particle diameter of 0.5 mm. It shows that the change in shear velocity has converged values on the bottom of the pipe, and that shear velocity does not change beyond a certain distance from the wall, demonstrating that shear friction exerts a greater force on the bottom, but does not extend to the pipe center.  The settling velocity of sediment particles, also called the terminal or fall velocity, is one of the key variables in the study of sediment transport and is important in understanding suspension, deposition, mixing and exchange processes. Based on the results of numerical analysis, we used Equation (11) to calculate the settling velocity (v s ) in the thin layer on the bottom of the pipe [62]: where (ρ s − ρ)/ρ is the specific gravity (submerged) of the sediment in the fluid, ρ s is the density of the sediment (kg/m 3 ), ρ is the density of water (kg/m 3 ), g is the gravitational acceleration (m/s 2 ), d p is the particle diameter (mm), and C D is the coefficient of drag. The coefficient of drag varies depending on the degree of freedom of motion of the particles; however, given that its value usually converges to 0.44 in a turbulent flow, the latter value was used for this calculation. Table 3 presents the input settling velocity by particle size at (g d p /C D ) as calculated from Equation (10). Table 4 shows the shear velocity values calculated from two-phase analysis. Table 5 presents the results from the analysis of the shear boundary layer under all 63 conditions of shear velocity, as listed in Table 4. The settling velocity of sediment particles, also called the terminal or fall velocity, is one of the key variables in the study of sediment transport and is important in understanding suspension, deposition, mixing and exchange processes. Based on the results of numerical analysis, we used Equation (11) to calculate the settling velocity (v ) in the thin layer on the bottom of the pipe [62]: where (ρ − ρ)/ρ is the specific gravity (submerged) of the sediment in the fluid, ρ is the density of the sediment (kg/m 3 ), ρ is the density of water (kg/m 3 ), g is the gravitational acceleration (m/s 2 ), d is the particle diameter (mm), and C is the coefficient of drag. The coefficient of drag varies depending on the degree of freedom of motion of the particles; however, given that its value usually converges to 0.44 in a turbulent flow, the latter value was used for this calculation. Table 3 presents the input settling velocity by particle size at (g d /C ) as calculated from Equation (10). Table 4 shows the shear velocity values calculated from two-phase analysis. Table 5 presents the results from the analysis of the shear boundary layer under all 63 conditions of shear velocity, as listed in Table 4.  10%

Discussion
Looking at these results ( Table 5), most of the sediments exceeding the value of 2.0 were not transported, but were deposited on the bottom, as has also been reported in earlier works [53]. In our study, however, at the mean value less than 1.0 also tended to be deposited in the stormwater sewer system. From the abrupt increase in shear velocity (0.2-0.9 m/s), it can be inferred that the corresponding boundary layer is one of alternately recurring transport and deposition due to particle saltation, with a proportion of sediments traveling in suspended load. Sediment particles with a settling velocity value less than 0.2 m/s were subjected to frequent changes in shear velocity, resulting in turbid flow in suspended load. Table 6 provides mathematical formulas for the shear velocity conditions of the turbulent boundary layer with respective sediment transport/deposition status. Our work thus verified that shear velocity (u * ) and settling velocity (v s ) can be used to easily distinguish between the transport or deposition status of particles with respect to the flow characteristics of soil slurry in stormwater pipes. Adhesion and consolidation may be defined as the process by which the sediment density increased. Sediments change their physical properties after deposition due to the stress from the overburden (gravitational compaction) and as a result of settling reactions involving coagulation. It should be noted that if the initial sediment concentration exceeds a certain value, there will be no initial free settlement stage, and the sediment settlement directly enters the hindered settlement or self-weight consolidation settlement stage [59,63]. This allows the prediction of the boundary layer for each particle transport status. The results presented in this paper are estimations derived from the numerical analysis of various inlet conditions of soil slurry entering the stormwater sewer system. Analysis of the settling velocity, with the particle size used as a major parameter, permitted the quantification of the sediment transport status in sediment layers within sewer pipes. The key advantage of the boundary layer investigation approach described in this study is that it can predict the thickness of the boundary layer for transport of soil slurry in both low-volume concentration conditions, and, also, in high-volume concentration conditions, involving transport of large quantities of soil slurry. Furthermore, unlike earlier numerical models of riverbeds, results from this study are applicable to stormwater sewers.
The parameters selected in this study to investigate the shear boundary layer fall within the range of the parameters proposed by earlier workers [53]. This confirms that the method proposed in this study produces reliable results. However, this study is limited by the fact that it is based on the single scenario that a large quantity of soil slurry enters the sewer system in the initial hours of a rainfall event.
This study is expected to serve as the basis for designing drainage pipelines and for mitigating clogging or destruction of drainage pipelines due to disasters (e.g., landslides) caused by localized torrential rainfall. However, prior to examining sewer flow patterns for various types of soil particles and floating debris that enter the drainage channel, it is necessary to analyze various conditions by integrating numerical analyses into hydraulic model experiments. Validation of the proposed methodology through its application to many different rainfall event scenarios and through comparison with direct observations in in actual sewer systems are, therefore, required.

Conclusions
In this study, we investigated turbulent flow in stormwater pipes, and performed numerical analysis to determine the effects of soil slurry transport and deposition. A turbulence model was developed, and shear velocity was calculated using flow velocity and sediment volume concentration as parameters. Based on the analysis results, we developed a new approach to identify the turbulent shear boundary layer related to soil slurry transport in stormwater pipes. We conducted numerical analysis to understand the turbulent flow in sewer pipes where transport and deposition take place, and analyzed the effects of soil slurry transport and deposition. Applying a turbulence model, we investigated the mean flow velocity, volume distribution, shear velocity, and shear stress. The following conclusions were drawn from this analysis: (1) In conditions of slow inlet velocity of slurry particles and large particle diameter, the soil slurry traveled in suspended load closer to the bottom; in conditions of high inlet velocity and smaller particle diameter, the soil slurry flow showed a higher suspended load transport.
(2) When slurry solids are deposited near the inlet, the flow on the sediment boundary layer is accelerated, resulting in entrainment and transport of slurry particles. When slurry enters the inlet slowly and particle diameter is large, deposition of a larger quantity of solids impedes slurry transport.
(3) The overall flow rate distribution was weakened due to the drag caused by the similar volume fraction in the sewer channel with sedimentation characteristics, and the flow velocity gradient in the pipe increased due to sedimentation.
Based on the analysis of the flow velocity field in the pipe, the boundary layers were expressed as the relationship between shear velocity and settling velocity by particle diameter. Since the mean flow velocity and volume concentration in the pipe under different conditions exhibited similar patterns, we averaged the correlation between settling velocity and shear velocity. It is clear that there are no quantitative results in the literature on the standard high Reynolds number k-ε model of turbulent pipe flow. The information given by Launder and Spalding [33] about this model is only qualitative. In order to compensate for these limitations, the turbulence intensity of the pipe core for the fully developed pipe flow was estimated through a sufficient review of various prior studies and reflected in the analysis. In order to fill this gap, in this paper we present accurate spectral solutions of this model. The results obtained are encouraging and tend to prove that the standard k-ε model can be relevant for pipe flows.
It is expected that this research methodology can be improved and refined by checking the boundary layer model proposed here against the direct observations in actual stormwater pipes. This will be the subject of a follow-up study.