Next Article in Journal
Suitability Mapping for Managed Aquifer Recharge: Development of Web-Tools
Next Article in Special Issue
CFD Investigations of Transient Cavitation Flows in Pipeline Based on Weakly-Compressible Model
Previous Article in Journal
The Interplay between Flow Field, Suspended Sediment Concentration, and Net Deposition in a Channel with Flexible Bank Vegetation
Previous Article in Special Issue
A Landscape Connectivity Approach for Determining Minimum Ecological Lake Level: Implications for Lake Restoration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Prediction of the Short-Term Trajectory of Microplastic Particles in Laizhou Bay

1
School of Environment, Beijing Normal University, Beijing 100875, China
2
The Key Laboratory of Water and Sediment Sciences of Ministry of Education, School of Environment, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(11), 2251; https://doi.org/10.3390/w11112251
Submission received: 30 August 2019 / Revised: 21 October 2019 / Accepted: 24 October 2019 / Published: 27 October 2019

Abstract

:
Microplastic particles are easily captured by microorganisms and enter the food chain, which poses a threat to ecological health. These particles are abundant in coastal areas because of the influence of anthropic activities and the interaction between the sea and land. Although much research on microplastics has been done, predicting the transportation of microplastic particles in coastal zones is still a challenge. In this paper, the trajectories of microplastic particles released from four river mouths around Laizhou Bay are investigated using the lattice Boltzmann method coupled with the Lagrangian particle-tracking method, involving inter-particle and particle-wall collisions. The trajectories of particles released from four river mouths are recorded within 30 days.

1. Introduction

Plastic debris is frequently found in rivers, lakes and oceans. The transportation of plastic debris needs the entrainment of water flows. Although its degradation process is extremely sluggish, plastic can be broken down into small particles during its transportation. Plastic particles with a diameter less than 5 mm are recognized as microplastics [1]. These particles can remain in the soil and sediment for hundreds of years and even longer [2]. In the meantime, microplastics may be captured by marine life and then transferred into the food chain. The microplastics themselves are harmful substances the are poisonous to living beings [3]. Hence, as a persistent pollutant, the study of microplastics has provoked much attention.
Currently, many microplastics are accumulated in coastal zones, released by land sources (tourist beaches, factory sewage and rivers) and marine sources (ship transport, fishing and marine cultivation). Browne et al. [4] tested the impact of the wind and the deposition mode on the spatial distribution of microplastics in the Tamar Estuary, the U.K. Law et al. [5] studied the concentration of the floating microplastics in the surface water of the North Atlantic. Dubaish and Liebezeit [6] studied the distribution of microplastics and black carbon particles in the surface water of Jade System, Germany. Nor and Obbard [7] investigated and categorized the microplastics extracted from sediment in several mangrove habitats in Singapore. Còzarzar [2] provided the worldwide distribution of microplastics and predicted the accumulation sites. In the following year, Còzar [8] studied the suspended plastics in the Mediterranean Sea and evaluated the pollution level. The existing research is focused on the sampling process, categorizing, identifying the source and assessing the ecological effect of the microplastics. However, the trajectories of microplastic particles need more attention. The area with more particle trajectories enhances the chance captured by the marine creature and should not be selected as a marine agriculture zone. These trajectories also provide an optimized area for the recycle process of microplastic particles.
The movement of microplastic particles in flows can be described as fluid-solid two-phased flow. Usually, this problem is simulated by the Eulerian or Lagrangian method. For the former method, those particles are treated as a continuum, which is not suitable for the cased with less quantity of particles. However, the Lagrangian one can demonstrate the trajectories of individual particles and dealing with the inter-particle collision and the particle-wall collision according to particles’ properties. Many research describe the particle-tracking method [9,10,11] and inter-particle collision [12,13,14,15,16]. Commonly, these models were applied to the sediment dynamics simulation as well [17,18]. The model proposed by Tanaka and Tsuji [12] is adopted in this paper.
The numerical method in this paper is the lattice Boltzmann method (LBM), developed into a general calculation in hydrodynamic modelling for its efficiency, simplicity and accuracy. The theory of the LBM in shallow water simulation shows a good performance for flow with complex boundary conditions and source terms [19,20,21,22]. This method can also deal the wetting-drying process on a dry bed with the wet-dry interface treatment [23]. Its inborn nature of both Eulerian and Lagrangian attributes [24] is benefit for coupling with Lagrangian particle-tracking method. These advantages make the LBM performing better in modelling the trajectories of microplastic particles.
This paper took Laizhou Bay as the study area and investigated the movement trajectory of microplastic particles. Laizhou Bay lies in the south of Bohai Sea, Shandong Province, China. In this area, there are many rivers flowing into the sea, such as the Yellow River, the Wei River and the Yu River [25]. This area has been severely exposed to anthropic activities, which produces microplastics by means of garbage dumping, fishing cultivation and local tourism [26]. Considering a great number of factories and domestic sewage flows into Laizhou Bay through rivers, this study focuses on predicting the trajectories of microplastic particles, which are discharged from the river mouth in this region, to determine the trajectory and the moving range of these microplastics.

2. Governing Equations

2.1. Shallow Water Equations

Considering that the vertical scale is far smaller than the horizontal scale in Laizhou Bay, this two-dimensional hydrodynamic field can be simulated based on shallow water equations [27], stated as
h t + ( h u j ) x j = 0 ,
( h u i ) t + ( h u i u j ) x j = g h 2 2 x i + ν 2 ( h u i ) x j x j + F i ,
where h is the water depth; t is time; u is the flow velocity; x is the Cartesian coordinate; g is the gravitational acceleration; and F i is the force term in the i direction, determined as
F i = g h z b x i + τ w i ρ τ b i ρ + E i ,
where z b is the bed elevation; ρ is the water density; τ w i is the wind shear stress and is governed by τ w i = ρ a C w u w i u w j u w j , in which ρ a is the density of air; C w is the resistance coefficient; u w i is the component of the wind velocity in i direction; τ b i is the bed shear stress and is given by τ b i = ρ C b u i u j u j , where C b is the coefficient of bed friction and is expressed as C b = g n b 2 / h 1 / 3 and n b is the Manning’s coefficient at the bed; E i is the Coriolis term, defined as
E i = f c h v , i = x , f c h u , i = y ,
in which f c = 2 ω s i n ϕ ; ω is the rotation of earth, given by ω 7.3 × 10 5 r a d / s ; ϕ is the latitude of the study area.

2.2. Microplastic Dynamics

Considering that the microplastic particles are very small and hydrophobic, these particles are treated as Lagrangian particles. They are assumed floating on the surface constantly, and the position of particles can be predicted by
d x m j d t = u m j ,
d u m j d t = 3 ν C D R e m 4 d m 2 S ( u j u m j ) ,
in which x m j is the position of microplastic particle j; Subscript m donates the microplastic particles; u m j is the velocity of microplastic particle j; u j is the flow velocity; d m is the diameter of microplastic particles and S is the ratio of microplastic particle density to fluid density. R e m is the Reynolds number of microplastic particles and C D is the drag coefficient, defined as below.
R e m = d m | u j u m j | ν ,
C D = 24 R e m , R e m < 1 ; 24 R e m ( 1 + 0.15 R e m 0.687 ) , 1 < R e m < 1000 ; 0.44 , R e m > 1000 .

3. Methods

3.1. Lattice Boltzmann Method

The lattice Boltzmann method is selected as the numerical method, by virtue of its efficiency and flexibility. The lattice Boltzmann method is a completely discrete mesoscopic model that treats the fluid as fractional particles in regular lattices, providing a promising approach to complex partial differential equations such as the shallow water equations [28].
To simulate the hydrodynamic field, the lattice Boltzmann method was used, and the lattice Boltzmann equation is given by
f α ( x + e α Δ t , t + Δ t ) f α ( x , t ) = 1 τ ( f α f α e q ) + W α 3 Δ t e 2 e α i F α ( x , t ) , α = 0 8 ,
in which f α is the distribution function of fictional particles in the α direction, f α e q is the local equilibrium function and τ is the relaxation time. W α is the weight coefficient is given by Equation (10), employing the D2Q9 lattice pattern as Figure 1. F α is the force component in α direction.
W α = 4 9 , α = 0 ; 1 9 , α = 1 , 3 , 5 , 7 ; 1 36 , α = 2 , 4 , 6 , 8 .
The lattice speed e is determined by
e = Δ x / Δ t .
The particle velocity vector e α is shown as
e α = ( 0 , 0 ) , α = 0 ; e c o s ( α 1 ) π 4 , s i n ( α 1 ) π 4 , α = 1 , 3 , 5 , 7 ; 2 e c o s ( α 1 ) π 4 , s i n ( α 1 ) π 4 , α = 2 , 4 , 6 , 8 .
The local equilibrium distribution function is expressed as
f α e q = h 5 g h 2 6 e 2 2 h 3 e 2 u i u i , α = 0 ; g h 2 6 e 2 + h 3 e 2 e α i u i + h 2 e 2 e α i e α j u i u j h 6 e 2 u i u j , α = 1 , 3 , 5 , 7 ; g h 2 24 e 2 + h 12 e 2 e α i u i + h 8 e 2 e α i e α j u i u j h 24 e 2 u i u j , α = 2 , 4 , 6 , 8 .
The macroscopic water depth and velocity can be calculated by
h ( x , t ) = α f α ( x , t ) ,
u i ( x , t ) = 1 h ( x , t ) α e α i f α ( x , t ) .
According to Equation (5), the position of microplastic particle at time t + Δ t can be obtained as
x m j t + Δ t = u m j t Δ t + x m j t ,
where x m j t + Δ t is the position of microplastic particle j at time t + Δ t , x m j t is the position of microplastic particle j at time t and u m j t is velocity of microplastic particle j at time t, given as
u m j t + Δ t = 3 ν C D R e p Δ t 4 d m 2 S ( u i t u m j t ) + u m j t .

3.2. Inter-Particle Collision

To investigate the influence extent of particles collision to trajectories, the inter-particle collision model is employed. The microplastic particles are assumed floating on the water surface and their trajectories could be affected by inter-particle collisions. When the concentration of particles is low, the binary collision model is adopted [12]. In this study, the particles are rigid, and the collision occurs instantaneously. This critical step of this model is obeying the law of conservation of momentum.

3.2.1. Collision Identification

The total number of particles is N in the computational field. In each computational step, N ( N 1 ) × 2 pairs of particles should be checked for the collision identification. Take particles i and j as an example to show the identifying progress; assume particle i is advanced to particle j and the relative motion of these two particles is demonstrated in Figure 2. The distance of these two particles is R 0 in time t and R Δ t in time t + Δ t . The relative position is R c = R 0 + k ( R Δ t R 0 ) , in which k is a factor of the nondimensional time scale. The occurrence of the collision requires that the following equation has two real roots, k 1 and k 2 ( k 1 < k 2 , 0 k 1 1 ).
| R 0 + k ( R Δ t R 0 ) | 2 = d m 2 ,
in which d m is the diameter of microplastic particles.

3.2.2. Post-Collision Quantities

After collision, the velocities of the collision pair can be obtained as
u m i = u m i + J / m m
u m j = u m j J / m m
where u m i and u m j are precollision velocities of particles i and j. J is the impulsive force, consisting of the normal J n and the tangential J t components.
J = J n n + J t t ,
J n = ( 1 + e ) c · n 1 / m i + 1 / m j ,
J t t = m i n [ μ J n , 2 | c f c | 7 ( 1 / m i + 1 / m j ) ] ,
in which c is the precollision relative velocity, c = u j m u i m , n is the normal unit velocity of the relative precollision position, n = R c / | R c | , t is the tangential unit vector of the slip velocity, and c f c is the tangential component of c .
The post-collision position of particles i and j is updated as
x m i = x m i + u m i k 1 Δ t + u m i ( 1 k 1 ) Δ t ,
x m j = x m j + u m j k 1 Δ t + u m j ( 1 k 1 ) Δ t ,
in which x m is the post-collision position, and x m is the precollision position.

3.3. Particle-Wall Collision

The particle-wall collision frequently occurs near riverbanks and seashores. To address the complex irregular boundary shape, the lattice Boltzmann method requires that the computational domain is covered with lattices, and then, the solid boundary is converted to a denoted value stored in the center of the lattices. To begin, considering the particles are in a meshless system, it is necessary to determine the lattice where the particle located to decide whether the particle will collide with wall or not. Taking the particle collision with lower boundary as an example, except for the particle on the solid boundary, there are four lattices surrounding a particle. Select the nearest lattice position and check its solid data. A solid lattice indicates that this particle collided with wall during the time ( t Δ t , t ) . Second, determine the particle position ( x m i , t Δ t ) at time t Δ t and the location of this collision ( ( x c o l l i d e , t + t w Δ t ) ), as shown in Figure 3. t w is the time before the collision, defined by t w = | x m i x c o l l i d e | / u m i , where u m i is the velocity before the collision. Finally, post-collision velocity is calculated by Equations (19), (21)–(23), where m j is equal to positive infinity [29]. The post-collision position at time t is
x m i = x c o l l i d e + u m i ( Δ t t w ) .

4. Case Study

4.1. Model Setup

Laizhou Bay (Figure 4) is one of the three bays in Bohai Sea, belonging to Shandong Province [30]. More than ten rivers, such as the Yu River and the Wei River, flow into this bay and carry sediment particles and nutrient substance into Laizhou Bay, building the habitat of natural animals and providing fishery resources [25]. However, with the increasing influence of anthropic activities, the sea water is polluted by different pollutants [26], one of which is microplastics. The degraded microplastic particles are transported with the sewage discharging into rivers and then flows into the sea.
The position and the bottom elevation of Laizhou Bay are demonstrated in Figure 5. This computational field is divided into 140 × 230 square lattices ( Δ x = 460 m). The north border of this area is set to the tidal boundary of 12 h period, as shown in Figure 6. The east water border is the open boundary, next to the shoreline, where the wet-dry boundary is used [31]. The most frequent wind is considered: SSW wind with an average speed of 2.5 m/s (with ρ α = 1.205 and C w = 0.0026 ). The initial water depth is 20m and the velocity is 0 m/s. The time step is 4s. The Manning coefficient of the bottom roughness is 0.01. Four groups of microplastic particles are released from the following four river mouths: the Yu River, the Di River, the Wei River and the Jiaolai River and numbered as Groups 1–4, respectively. They are released at the same time. Each group has 81 particles, equally distributed in a lattice (Figure 7). All particles are assumed floating on the water surface and the density of particles is 850 kg/m 3 . The initial velocity of these particles is 0 m/s.
This model is validated by the measured value of two sites in Laizhou Bay, Weifang Harbor and Laizhou Port [31]. In Figure 8, the results show a good agreement with the measured data, with less than 15 % error, demonstrating the well accuracy of this numerical model.
The hydrodynamic model is run for 24 h to obtain a well-developed flow field and then reset the computational time to 0. Considering that the hydrodynamic condition in this bay is varying with the tidal, the trajectories of microplastic particles are varying correspondingly. Hence, the release time of 324 microplastic particles at 4 river mouths are (1) 0 h, (2) 3 h, (3) 6 h and (4) 9 h, which are used to investigate the influence of the hydrodynamic condition on particle transportation.

4.2. Results and Discussion

4.2.1. Particles Input at t = 0 h

There are four groups of particles released from the four river mouths at t = 0 h in this section, respectively. The particles numbered 1 of each group are selected as representative, and their trajectory are recorded within 24 h. Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 demonstrate the trajectories of these particles in the following four time periods: 0–6 h, 6–12 h, 12–18 h and 18–24 h, respectively. x is the distance in the lateral direction and y is the distance in the longitudinal direction. The units are in meters. The origin is located at ( 118.86 E, 37.09 N). At time t = 0 (Figure 9), these particles are in the initial position and then move with the water flow. During this period, the transport distance of the particle in Group 4 is the largest, while Group 1 is the smallest. Figure 10 shows the amplified diagram of trajectories of 4 particles and the velocity fields at t = 0 , 3 , 6 h. At t = 0 h, the flow moves toward the offshore direction, then it moves eastward at t = 3 h and southwestward at t = 6 h. The movement direction of particles is changing along with the flow direction. The particle of Group 1 starts from (36,156, 17,296). At first, it moves northeastward to (36,828, 18,312) and then moves back to (36,560, 17,460). During the first six hours, it travels 2111 m, and the absolute distance is 436 m. Similarly, the particles in Groups 2 and 3 start from (42,136, 14,996) and (54,556, 14,996) and transport northeast to (43,338, 16,433) and (55,703, 17,187), respectively. Eventually, they reach (43,140, 15,220) and (55,959, 15,062). They move approximately 3102 m and 4613 m, and the absolute distances are 1028 m and 1404 m. The particle in Group 4 moves from (62,836, 15,456) northeast to (62,731, 18,285) and stays at (63,118, 15,264). Although its movement distance is the largest (5876 m), the absolute distance is only 341 m by the end of this period.
Unlike the previous situation, the traveling distances in Figure 11 are shorter, and the general tendency is towards the shore during 6–12 h. The particles in Group 1, 2 and 3 are transported from the final position of the first 6 h period to (35,261, 17,140), (41,537, 14,779) and (54,088, 14,829). Their transport distance is basically equal to the absolute distance as 1338 m, 1663 m and 1886 m. During this period, apparently, the particle in Group 4 collides on the shore at (63,213, 15,181) and then bounces back to (63,183, 15,230) at the end, moving approximately 184 m. The cumulative collision number of 81 particles in Group 3 is 431 during this period (Figure 13), while there is no collision that occurs between particles of other groups. Although collide happens, it merely influences changing trajectories of particles.
The trajectories of particles during 12–18 h and 18–24 h are similar to those during 0–6 h and 6–12 h (Figure 14, Figure 15, Figure 16 and Figure 17), while the transport distances are shorter. During 12–18 h, the transport distances in Groups 1–4 are 1302 m, 1899 m, 3110 m and 3510 m, respectively. No collision happens in this period. In the last six hours, the particles in four groups are moving 969 m, 671 m, 219 m and 123 m. The collision number during 18–24 h is 801, of which 191 happen in Group 3 and 610 in Group 4.
In general, the transport distances of particles range from 5720 m to 9828 m and the absolute distances are within the range of 366–2031 m during 24 h. The absolute distances of particles in Group 3 and 1 are 2031 and 1630 m, followed by that in Group 4 (514 m) and in Group 2 (366 m). The transport distance of particles in Group 3 is the farthest (9828 m), being 1.7 times as long as that in Group 1 (5720 m). Although the absolute distances of particles in Group 4 and 2 are short, their transport distances reach 9693 m and 7335 m. The transport distances of particles in Groups 1–4 near the releasing points are 2397 m × 1224 m, 1354 m × 1689 m, 2296 m × 2461 m and 594 m × 3104 m, respectively. Among them, the transport distance of particles in Group 3 is the largest while in Group 4 is the smallest. If there were microplastic particles released from these four river mouths within 24 h, it is reasonable to collect pollution particles within the scope of 2397 m × 3104 m near the seashore.

4.2.2. Particles Input at t = 3 , 6 , 9 h

The trajectories of particles in four groups released at t = 0 , 3 , 6 and 9 h within 24 h are demonstrated in Figure 18. It can be seen that the moving ranges of the particles released at t = 9 h are the largest, being 4593 m × 1952 m in Group 1, 2353 m × 2618 m in Group 2, 2010 m × 5670 m in Group 3 and 1270 m × 8242 m in Group 4. These ranges are over twice than those particles released at t = 0 , 3 and 6 h, respectively, of which the moving ranges are similar, being 2398 m × 1224 m in Group 1, 1355 m × 1690 m in Group 2, 2297 m × 2462 m in Group 3 and 595 m × 3105 m in Group 4; 1591 m × 1252 m in Group 1, 1555 m × 1668 m in Group 2, 1345 m × 2467 m in Group 3 and 325 m × 3105 m in Group 4; 2307 m × 1247 m in Group 1, 1957 m × 1700 m in Group 2, 2236 m × 2463 m in Group 3 and 421 m × 3104 m in Group 4, respectively. Hence, it can be perceived that the microplastic particles released at t = 9 h can do little beneficial to the recycling process.
The cumulative collision number of particles released at t = 3 h is the most, which is 1.01 × 10 5 , followed by those released at t = 9 h ( 8.48 × 10 4 ). The numbers of collision of particles released at t = 0 and 6 h are only 4.05 × 10 4 and 3.32 × 10 4 .

4.2.3. Particles Input at t = 0 h within 30 Days

Compared to the trajectories of particles within 24 h, it can be seen in Figure 19 that the moving ranges of particles within 30 days expands obviously. The transport distances of particles in Groups 1–4 near the releasing points are 4732 m × 1224 m, 2375 m × 1709 m, 6393 m × 2466 m and 3587 m × 3105 m, respectively. In x direction, the transport distances of particles in Group 3 is the largest, while it is a relatively narrow range in Group 2 (2.7 times less than that in Group 3). Moreover, the particles in Group 4 have the widest transport distances in y direction, 2.5 times larger than that in Group 1. Generally, the transport distance of particles in Group 3 is the largest. The scope of 6393 m × 3105 m can be the minimal collect scope of microplastic particles releasing within 30 days.

5. Conclusions

In this work, the movement trajectories of microplastic particles in Laizhou Bay are simulated based on the lattice Boltzmann framework coupled with the Lagrangian particle-tracking method. In addition to the current, inter-particle collisions and particle-wall collisions are also involved. After releasing, particles move forth and back with the changing tides. The trajectories of particles released at t = 0 , 3 , 6 , 9 h are recorded with tide boundary of 12 h period. The results show that the moving range of the particles releasing at t = 9 h is the largest, which is unfavorable for the recycle process of microplastics. The simulated results indicate that these particles normally drift within the scope of 4593 m × 8242 m near the releasing point. Collisions indeed occur with as many as 3.32 × 10 4 within 24 h, which has little effect on changing trajectories of particles. Even within 30 days, the transportation of these particles is still near the seashore and their moving ranges expand to 6393 m × 3105 m from the releasing point.

Author Contributions

Conceptualization, H.L., W.Y. and Y.D.; methodology, H.L. and Y.D.; software, W.Y. and Y.D.; validation, Y.D.; resources, H.L.; data curation, H.L.; writing—original draft preparation, Y.D; writing—review and editing, H.L. and W.Y.; supervision, H.L.

Funding

This research was funded by the National Key R&D Program (2018YFC1406401) and the National Natural Science Foundation of China (51779011).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arthur, C.; Baker, J.; Bamford, H. Proceedings of the International Research Workshop on the Occurance, Effects, and Fate of Mircroplastic Marine Debris; NOAA Technical Memorandum NOS-OR&R-30; National Oceanic and Atmospheric Administration: Silver Spring, MD, USA, 2009. [Google Scholar]
  2. Còzar, A.; Echevarría, F.; González-Gordillo, J.I.; Irigoien, X.; Úbeda, B.; Hernández-León, S.; Palma, A.T.; Navarro, S.; García-de Lomas, J.; Ruiz, A.; et al. Plastic debris in the open ocean. Proc. Natl. Acad. Sci. USA 2014, 11, 10239–10244. [Google Scholar] [CrossRef] [PubMed]
  3. Ashton, K.; Holmes, L.; Turner, A. Association of metals with plastic production pellets in the marine environment. Mar. Pollut. Bull. 2010, 60, 2050–2055. [Google Scholar] [CrossRef] [PubMed]
  4. Browne, M.A.; Galloway, T.S.; Thompson, R.C. Spatial patterns of plastic debris along Estuarine shorelines. Environ. Sci. Technol. 2010, 44, 3404–3409. [Google Scholar] [CrossRef] [PubMed]
  5. Law, K.L.; Morét-Ferguson, S.; Maximenko, N.A.; Proskurowski, G.; Peacock, E.E.; Hafner, J.; Reddy, C.M. Plastic accumulation in the North Atlantic subtropical gyre. Science 2010, 329, 1185. [Google Scholar] [CrossRef]
  6. Dubaish, F.; Liebezeit, G. Suspended Microplastics and Black Carbon Particles in the Jade System, Southern North Sea. Water Air Soil Pollut. 2013, 224, 1352–1360. [Google Scholar] [CrossRef]
  7. Nor, N.H.M.; Obbard, J.P. Microplastics in Singapore’s coastal mangrove ecosystems. Mar. Pollut. Bull. 2014, 79, 278–283. [Google Scholar]
  8. Còzar, A.; Sanz-Martín, M.; Martí, E.; González-Gordillo, J.I.; Ubeda, B.; Gálvez, J.A.; Irigoien, X.; Duarte, C.M. Plastic accumulation in the Mediterranean sea. PLoS ONE 2015, 10, e0121762. [Google Scholar] [CrossRef]
  9. Wiberg, P.L.; Smith, J.D. A theoretical model for saltating grains in water. J. Geophys. Res. 1985, 90, 7341. [Google Scholar] [CrossRef]
  10. Niño, Y.; García, M.; Ayala, L. Gravel saltation: 1. Experiments. Water Resour. Res. 1994, 30, 1907–1914. [Google Scholar] [CrossRef]
  11. Lee, H.Y.; Lin, Y.T.; You, J.Y.; Wang, H.W. On three-dimensional continuous saltating process of sediment particles near the channel bed. J. Hydraul. Res. 2006, 44, 374–389. [Google Scholar] [CrossRef]
  12. Tanaka, T.; Tsuji, M. Numerical simulation of gas-solid two-phase flow in a vertical pipe: On the effect of inter-particle collision. Gas-Solid Flows. ASME 1991, 121, 123–128. [Google Scholar]
  13. Mei, R.; Hu, K.C. On the collision rate of small particles in turbulent flows. J. Fluid Mech. 1999, 391, 67–89. [Google Scholar] [CrossRef]
  14. Yeganeh-Bakhtiary, A.; Shabani, B.; Gotoh, H.; Wang, S.S.Y. A three-dimensional distinct element model for bed-load transport. J. Hydraul. Res. 2009, 47, 203–212. [Google Scholar] [CrossRef]
  15. Lin, J.H.; Chang, K.C. Particle Dispersion Simulation in Turbulent Flow Due to Particle-Particle and Particle-Wall Collisions. J. Mech. 2016, 32, 237–244. [Google Scholar] [CrossRef]
  16. Fede, P.; Simonin, O. Direct Simulation Monte-Carlo predictions of coarse elastic particle statistics in fully developed turbulent channel flows: Comparison with deterministic discrete particle simulation results and moment closure assumptions. Int. J. Multiph. Flow 2018, 108, 25–41. [Google Scholar] [CrossRef] [Green Version]
  17. Niño, Y.; García, M. Using Lagrangian particle saltation observations for bedload sediment transport modelling. Hydrol. Process. 1998, 12, 1197–1218. [Google Scholar] [CrossRef]
  18. Chang, Y.; Scotti, A. Entrainment and suspension of sediments into a turbulent flow over ripples. J. Turbul. 2003, 4, 1–22. [Google Scholar] [CrossRef]
  19. Bhatnagar, P.L.; Gross, E.P.; Krook, M. A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems. Phys. Rev. 1954, 94, 511–525. [Google Scholar] [CrossRef]
  20. Liu, H.; Zhou, G.J.; Burrows, R. Lattice Boltzmann model for shallow water flows in curved and meandering channels. Int. J. Comput. Fluid Dyn. 2009, 23, 209–220. [Google Scholar] [CrossRef]
  21. Zhou, J.G. Lattice Boltzmann Methods for Shallow Water Flows; Springer: Berlin, Germany, 2004; Volume 4. [Google Scholar]
  22. Sun, T.; Li, W.; Dong, B. Numerical and experimental study on the motion characteristics of single bubble in a complex channel. Int. J. Comput. Fluid Dyn. 2015, 29, 346–356. [Google Scholar] [CrossRef]
  23. Liu, H.; Zhou, J.G. Lattice Boltzmann approach to simulating a wetting-drying front in shallow flows. J. Fluid Mech. 2014, 743, 32–59. [Google Scholar] [CrossRef]
  24. Chen, S.; Doolen, G. Lattice Boltzmann method for fluid flows. Annu. Rev. Fluid Mech. 1998, 30, 329–364. [Google Scholar] [CrossRef]
  25. Zhang, Y.; Li, H.; Wang, X.; Zheng, C.; Wang, C.; Xiao, K.; Wan, L.; Wang, X.; Jiang, X.; Guo, H. Estimation of submarine groundwater discharge and associated nutrient fluxes in eastern Laizhou Bay, China using 222 Rn. J. Hydrol. 2016, 533, 103–113. [Google Scholar] [CrossRef]
  26. Chang, Y.; Hu, B.X.; Xu, Z.; Li, X.; Tong, J.; Chen, L.; Zhang, H.; Miao, J.; Liu, H.; Ma, Z. Numerical simulation of seawater intrusion to coastal aquifers and brine water/freshwater interaction in south coast of Laizhou Bay, China. J. Contam. Hydrol. 2018, 215, 1–10. [Google Scholar] [CrossRef] [Green Version]
  27. Zhou, J.G. Enhancement of the LABSWE for shallow water flows. J. Comput. Phys. 2011, 230, 394–401. [Google Scholar] [CrossRef]
  28. Ding, Y.; Liu, H.; Peng, Y.; Xing, L. Lattice Boltzmann method for rain-induced overland flow. J. Hydrol. 2018, 562, 789–795. [Google Scholar] [CrossRef]
  29. Yamamoto, Y.; Potthoff, M.; Tanaka, T.; Kajishima, T.; Tsuji, Y. Large-eddy simulation of turbulent gas-particle flow in a vertical channel: Effect of considering inter-particle collisions. J. Fluid Mech. 2001, 442, 303–334. [Google Scholar] [CrossRef]
  30. Rong, Q.; Liu, J.; Cai, Y.Z.L.; Zhao, Z.; Yue, W.; Xia, J. Leaf carbon, nitrogen and phosphorus stoichiometry of Tamarix chinensis, Lour. in the Laizhou Bay coastal wetland, China. Ecol. Eng. 2015, 76, 57–65. [Google Scholar] [CrossRef]
  31. Liu, H.; Zhang, J.; Wang, H.; Ding, Y.; Yi, Y. Numerical modeling of the tidal wave run-up and the eelgrass habitat at the Laizhou Bay. Ecol. Model. 2017, 360, 378–386. [Google Scholar] [CrossRef]
Figure 1. Sketch of D2Q9 lattice pattern.
Figure 1. Sketch of D2Q9 lattice pattern.
Water 11 02251 g001
Figure 2. Sketch of binary collision.
Figure 2. Sketch of binary collision.
Water 11 02251 g002
Figure 3. Sketch of particle-wall collision.
Figure 3. Sketch of particle-wall collision.
Water 11 02251 g003
Figure 4. Map of Laizhou Bay.
Figure 4. Map of Laizhou Bay.
Water 11 02251 g004
Figure 5. The position and the bottom elevation of Laizhou Bay (referencing to the datum level of the Yellow Sea).
Figure 5. The position and the bottom elevation of Laizhou Bay (referencing to the datum level of the Yellow Sea).
Water 11 02251 g005
Figure 6. The tidal curve of the north boundary.
Figure 6. The tidal curve of the north boundary.
Water 11 02251 g006
Figure 7. The layout of 81 microplastic particles in a lattice.
Figure 7. The layout of 81 microplastic particles in a lattice.
Water 11 02251 g007
Figure 8. The validation of tidal data in Weifang Harbor (a) and Laizhou Port (b).
Figure 8. The validation of tidal data in Weifang Harbor (a) and Laizhou Port (b).
Water 11 02251 g008
Figure 9. The trajectories of 4 particles during 0–6 h.
Figure 9. The trajectories of 4 particles during 0–6 h.
Water 11 02251 g009
Figure 10. The amplified diagram of trajectories of 4 particles during 0–6 h.
Figure 10. The amplified diagram of trajectories of 4 particles during 0–6 h.
Water 11 02251 g010
Figure 11. The trajectories of 4 particles during 6–12 h.
Figure 11. The trajectories of 4 particles during 6–12 h.
Water 11 02251 g011
Figure 12. The amplified diagram of trajectories of 4 particles during 6–12 h.
Figure 12. The amplified diagram of trajectories of 4 particles during 6–12 h.
Water 11 02251 g012
Figure 13. The trajectories of 81 particles in Group 3 during 6–12 h.
Figure 13. The trajectories of 81 particles in Group 3 during 6–12 h.
Water 11 02251 g013
Figure 14. The trajectories of 4 particles during 12–18 h.
Figure 14. The trajectories of 4 particles during 12–18 h.
Water 11 02251 g014
Figure 15. The amplified diagram of trajectories of 4 particles during 12–18 h.
Figure 15. The amplified diagram of trajectories of 4 particles during 12–18 h.
Water 11 02251 g015
Figure 16. The trajectories of 4 particles during 18–24 h.
Figure 16. The trajectories of 4 particles during 18–24 h.
Water 11 02251 g016
Figure 17. The amplified diagram of trajectories of 4 particles during 18–24 h.
Figure 17. The amplified diagram of trajectories of 4 particles during 18–24 h.
Water 11 02251 g017
Figure 18. The trajectories of particles released at different time.
Figure 18. The trajectories of particles released at different time.
Water 11 02251 g018
Figure 19. The trajectories of 4 particles in 30 days.
Figure 19. The trajectories of 4 particles in 30 days.
Water 11 02251 g019

Share and Cite

MDPI and ACS Style

Ding, Y.; Liu, H.; Yang, W. Numerical Prediction of the Short-Term Trajectory of Microplastic Particles in Laizhou Bay. Water 2019, 11, 2251. https://doi.org/10.3390/w11112251

AMA Style

Ding Y, Liu H, Yang W. Numerical Prediction of the Short-Term Trajectory of Microplastic Particles in Laizhou Bay. Water. 2019; 11(11):2251. https://doi.org/10.3390/w11112251

Chicago/Turabian Style

Ding, Yu, Haifei Liu, and Wei Yang. 2019. "Numerical Prediction of the Short-Term Trajectory of Microplastic Particles in Laizhou Bay" Water 11, no. 11: 2251. https://doi.org/10.3390/w11112251

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop