Construction of solutions and study of their closeness in $L_2$ for two boundary value problems for a model of multicomponent suspension transport in coastal systems

Three-dimensional models of suspension transport in coastal marine systems are considered. The associated processes have a number of characteristic features, such as high concentrations of suspensions (e.g., when soil is dumped on the bottom), much larger areas of suspension spread than the reservoir depth, complex granulometric (multifractional) content of suspensions, and mutual transitions between fractions. Suspension transport can be described using initial-boundary value diffusion–convection–reaction problems. According to the authors' idea, on a time grid constructed for the original continuous initial-boundary value problem, the right-hand sides are transformed with a “delay” so that the right-hand side concentrations of the components other than the underlying one (for which the initial-boundary value problem of diffusion–convection is formulated) are determined at the preceding time level. This approach simplifies the subsequent numerical implementation of each of the diffusion–convection equations. Additionally, if the number of fractions is three or more, the computation of each of the concentrations at every time step can be organized independently (in parallel). Previously, sufficient conditions for the existence and uniqueness of a solution to the initial-boundary value problem of suspension transport were determined, and a conservative stable difference scheme was constructed, studied, and numerically implemented for test and real-world problems. In this paper, the convergence of the solution of the delay-transformed problem to the solution of the original suspension transport problem is analyzed. It is proved that the differences between these solutions tends to zero at an $O(\tau)$ rate in the norm of the Hilbert space $L_2$ as the time step $\tau$ approaches zero.