Scavenging processes in multicomponent medium with first-order reaction kinetics: Lagrangian and Eulerian modeling

A process of the removal of dissolved elements in the ocean by adsorption onto settling particulate matters (scavenging) is studied analytically and using the Lagrangian and Eulerian numerical methods. The generalized model of scavenging in a multicomponent reactive medium with first-order kinetics consisting of water and multifraction suspended particulate matter was developed. Two novel numerical schemes were used to solve the transport–diffusion–reaction equations for transport dominated flows. The particle tracking algorithm based on the method of moments was developed. The modified flux-corrected transport method for the Eulerian transport–diffusion–reaction equations is a flux-limiter method based on a convex combination of low-order and high-order schemes. The flux limiters in the developed approach are obtained as an approximate solution of a corresponding optimization problem with a linear objective function. This approach allows the construction of the flux limiters with desired properties. The similarity solutions of the model equations for an idealized case of instantaneous release of reactive radionuclide on the ocean surface were obtained. It was found analytically that the dispersion of reactive contamination caused by reversible phase transition with increase of settling velocity, concentration of particulate matter and distribution coefficient can be much greater than that caused by diffusion, whereas an increase in the desorption rate results in a decrease of the dispersion caused by the phase transfer. The solutions using both numerical schemes are consistent with the analytical similarity solution even at zero diffusivity. The scavenging of the 239,240\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{239,240}$$\end{document}Pu that was introduced to the ocean surface due to the fallout from past nuclear weapon testing was simulated. The simulation results were in agreement with the observations in the northern Pacific. It was shown that even if the concentration of the 239,240\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{239,240}$$\end{document}Pu on the particulate matter does not exceed 2% of the total concentration, settling of particulate matter plays a crucial role in the vertical transport and dispersion of the reactive radionuclide. The importance of the scavenging by both the large fast-settling particles and small particles slowly settling and dissolving with depth due to the biochemical processes was demonstrated. For large particles, the “pseudodiffusivity” caused by phase transfer was 60 times greater than the diffusivity.