A comparison of Eulerian and Lagrangian methods for vertical particle transport in the water column
-
Published:2023-09-19
Issue:18
Volume:16
Page:5339-5363
-
ISSN:1991-9603
-
Container-title:Geoscientific Model Development
-
language:en
-
Short-container-title:Geosci. Model Dev.
Author:
Nordam TorORCID, Kristiansen Ruben, Nepstad Raymond, van Sebille ErikORCID, Booth Andy M.
Abstract
Abstract. A common task in oceanography is to model the vertical movement of particles
such as microplastics, nanoparticles, mineral particles, gas bubbles, oil
droplets, fish eggs, plankton, or algae. In some cases, the distribution of the
vertical rise or settling velocities of the particles in question can span a
wide range, covering several orders of magnitude, often due to a broad
particle size distribution or differences in density. This requires
numerical methods that are able to adequately resolve a wide and possibly
multi-modal velocity distribution. Lagrangian particle methods are commonly used for these applications. A
strength of such methods is that each particle can have its own rise or
settling speed, which makes it easy to achieve a good representation of a
continuous distribution of speeds. An alternative approach
is to use Eulerian methods, where the partial differential equations
describing the transport problem are solved directly with numerical methods.
In Eulerian methods, different rise or settling speeds must be represented
as discrete classes, and in practice, only a limited number of classes can be
included. Here, we consider three different examples of applications for a
water column model: positively buoyant fish eggs, a mixture of positively
and negatively buoyant microplastics, and positively buoyant oil droplets
being entrained by waves. For each of the three cases, we formulate a model
for the vertical transport based on the advection–diffusion equation with
suitable boundary conditions and, in one case, a reaction term. We give a
detailed description of an Eulerian and a Lagrangian implementation of these
models, and we demonstrate that they give equivalent results for selected
example cases. We also pay special attention to the convergence of the model
results with an increasing number of classes in the Eulerian scheme and with the
number of particles in the Lagrangian scheme. For the Lagrangian scheme, we
see the 1/Np convergence, as expected for a Monte Carlo method,
while for the Eulerian implementation, we see a second-order (1/Nk2)
convergence with the number of classes.
Funder
Norges Forskningsråd Horizon 2020 H2020 European Research Council Nederlandse Organisatie voor Wetenschappelijk Onderzoek
Publisher
Copernicus GmbH
Reference73 articles.
1. Abuhegazy, M., Talaat, K., Anderoglu, O., and Poroseva, S. V.: Numerical
investigation of aerosol transport in a classroom with relevance to COVID-19,
Phys. Fluids, 32, 103311, https://doi.org/10.1063/5.0029118, 2020. a 2. Benson, D. A., Aquino, T., Bolster, D., Engdahl, N., Henri, C. V., and
Fernandez-Garcia, D.: A comparison of Eulerian and Lagrangian transport and
non-linear reaction algorithms, Adv. Water Resour., 99, 15–37,
https://doi.org/10.1016/j.advwatres.2016.11.003, 2017. a 3. Billingsley, P.: Probability and measure, John Wiley & Sons, New York
Chichester Brisbane Toronto, ISBN 0471031739, 1979. a, b 4. Craig, P. D. and Banner, M. L.: Modeling wave-enhanced turbulence in the ocean
surface layer, J. Phys. Oceanogr., 24, 2546–2559,
https://doi.org/10.1175/1520-0485(1994)024<2546:MWETIT>2.0.CO;2, 1994. a 5. Cui, F., Boufadel, M. C., Geng, X., Gao, F., Zhao, L., King, T., and Lee, K.:
Oil Droplets Transport Under a Deep-Water Plunging Breaker: Impact of Droplet
Inertia, J. Geophys. Res.-Oceans, 123, 9082–9100,
https://doi.org/10.1029/2018JC014495, 2018. a
|
|