Sufficient conditions for convergence of positive solutions to linearized two-dimensional sediment transport problem

Introduction. The sediment transportation is one of the major processes that define the magnitude and back surface changing rate for water bodies. The most used prognostic studies in this field are based on the mathematical models that allow reducing, and in some cases, eliminating the expensive and often environmentally burdensome experiments. Spatially one-dimensional models are mainly used to predict changes in bottom topography. For actual coastal systems with irregular coastal configuration, the flow vector is generally not orthogonal to the tangent line for the coastline at each of its points. It also may not coincide with the wind stress vector. Therefore, to solve lots of practically important problems associated with the prediction of the dynamics of the back surface of water basins, it is necessary to use spatially two-dimensional models of sediment transportation and effective numerical methods of their implementation. Materials and Methods. The authors (A.I. Sukhinov, A.E. Chistyakov, E.A. Protsenko, and V.V. Sidoryakina) have earlier proposed a spatially two-dimensional model of sediment transport that satisfies the basic conservation laws (of material balance and momentum) which is a quasilinear equation of parabolic type. The linear difference schemes are constructed and studied; model and some practically important tasks are solved. However, the theoretical study on the proximity of solutions for the original nonlinear initial-boundary value problem and the linearized continuous task - on which basis a discrete model (difference scheme) was built - remained in the shadow. The study of the linearized problem correctness and the determination of sufficient conditions for positivity of solutions are of special interest because only positive solutions to this sediment transport problem have physical value within the framework of the considered models. Research Results. The investigated nonlinear two-dimensional model of sediment transport in the coastal zone of shallow reservoirs takes account of the following physically significant factors and parameters: soil porosity; critical value of the tangent stress at which load transport starts; turbulent interaction; dynamically varying of the bottom geometry; wind currents; and bottom friction. Linearization is carried out on the time grid - nonlinear coefficients of the parabolic equation are taken for the previous time grid step. Next, a chain of tasks connected by the initial data - final solutions of the linearized mixed Cauchy problems chain on a uniform time grid is constructed, and thus, the linearization for the initial 2D nonlinear model is carried out in large. Earlier, the authors have proved the existence and uniqueness of the solution to a linear tasks chain. Prior estimate of the proximity of the linearized problem chain solution to the initial non-linear task solution has been also obtained. The conditions of its solution positivity and their convergence to the nonlinear sediment transport problem are defined in the norm of the Hilbert space L₁ with the rate O(τ) where τ is a time step. Discussion and Conclusions. The obtained research results of the spatially two-dimensional nonlinear sediment transport model can be used for predicting the nonlinear hydrodynamic processes, improving their accuracy and reliability due to the availability of new accounting functionality of physically important factors, including the specification of the boundary conditions.

пространственно-двумерная модель транспорта наносов, прибрежная зона, нелинейная задача, линеаризованная задача, положительность решения, spatially 2D model of sediment transportation, coastal zone, nonlinear problem, the linearized problem, solution positivity

1. Marchuk, G.I., Dymnikov, V.P., Zalesny, V.B. Matematicheskie modeli v geofizicheskoy gidrodinamike i chislennye metody ikh realizatsii. [Mathematical models in geophysical hydrodynamics and numerical methods for their implementation.] Leningrad: Gidrometeoizdat, 1987, 296 p. (in Russian).

2. Sukhinov, A.I., Chistyakov, A.E., Alekseenko, E.V. Numerical realization of the three-dimensional model of hydrodynamics for shallow water basins on a high- performance system. Mathematical Models and Computer Simulations, 2011, vol. 3, iss. 5, pp. 562–574.

3. Leontyev, I.O. Pribrezhnaya dinamika: volny, techeniya potoki nanosov. [Coastal dynamics: waves, moving streams, deposits drifts.] Moscow: GEOS, 2001, 272 p. (in Russian).

4. Xiaoying Liu, et al. Predictive modeling in sediment transportation across multiple spatial scales in the Jialing River Basin of China. International Journal of Sediment Research, 2015, vol. 30, iss. 3, pp. 250–255.

5. Sukhinov, A.I., Chistyakov, A.E., Protsenko, E.A. Mathematical modeling of sediment transport in the coastal zone of shallow reservoirs. Mathematical Models and Computer Simulations, 2014, vol. 6, iss. 4, pp. 351–363.

6. Sukhinov, A.I., Sidoryakina, V.V. Sushchestvovanie i edinstvennost' resheniya linearizovannoy dvumernoy zadachi transporta nanosov. [Availability and uniqueness of solution to linearized two-dimensional sediment transportation problem.] Teoriya operatorov, kompleksnyy analiz i matematicheskoe modelirovanie: tezisy dokladov XIII mezhdunar. nauch. konf. [Theory of operators, complex analysis and mathematical modeling: Abstracts of XIIIth Int. Sci.Conf.] Vladikavkaz: YuMI VNTs RAN, 2016, pp. 184–185 (in Russian).

7. Sukhinov, A.I., Sidoryakina, V.V. O edinstvennosti resheniya linearizovannoy dvumernoy nachal'no-kraevoy zadachi transporta nanosov. [On the uniqueness of solutions of the linearized two-dimensional bottom deposit transportation problem.] Taganrog Institute Journal, 2016, no. 2, pp. 270–274 (in Russian).

8. Sukhinov, A.I., Sidoryakina, V.V., Sukhinov, A.A. Apriornaya otsenka resheniya dvumernoy zadachi transporta nanosov. [A priori estimate of solution to two- dimensional sediment transportation problem.] Aktual'nye problemy prikladnoy matematiki i avtomatizatsii: mat-ly mezhdunar. konf.; Nelokal'nye kraevye zadachi i sovremennye problemy analiza i informatiki: mat-ly XIV shkoly molodykh uchenykh. [Challenging issues of applied mathematics and automation: Proc. Int. Conf.; Nonlocal boundary-value problems and current problems of analysis and informatics: Proc. XIV School of young scholars.] Terskol, 2016, 363 p. (in Russian).

9. Sukhinov, A.I., Chistyakov, A.E. Adaptive modified alternating triangular iterative method for solving grid equations with a non-self-adjoint operator. Mathematical Models and Computer Simulations, 2012, vol. 4, iss. 4, pp. 398–409.

10. Sukhinov, A.I., Chistyakov, A.E., Protsenko, E.A. Matematicheskoe modelirovanie transporta nanosov v pribrezhnykh vodnykh sistemakh na mnogoprotsessornoy vychislitel'noy sisteme. [Sediment transport mathematical modeling in a coastal zone using multiprocessor computing systems.] Numerical Methods and Programming, 2014, vol. 15, iss. 4, pp. 610–620 (in Russian).

11. Sukhinov, A.I., et al. Sravnenie vychislitel'nykh effektivnostey yavnoy i neyavnoy skhem dlya zadachi transporta nanosov v pribrezhnykh vodnykh sistemakh. [Comparison of computational efficiency of explicit and implicit schemes.] Numerical Methods and Programming, 2015, vol. 16, iss. 3, pp. 328–338 (in Russian).

12. Godunov, S.K. Uravneniya matematicheskoy fiziki. [Equations of mathematical physics.] 2nd revised and enlarged ed. Moscow: Nauka, 1979, 392 p. (in Russian).