|
Article information
2017 , Volume 22, Special issue, p.113-124
Tsydenov B.O.
A numerical study of impurity propagation in a freshwater lake on the basis of water turbidity distribution
A mathematical 2.5D model for simulating the hydrodynamic processes of pollutant propagation on an example of Kamloops Lake (British Columbia, Canada) is presented. Water turbidity is considered as an indicator of pollution. The non-hydrostatic model includes the continuity, momentum, energy, salinity, and turbidity equations and takes into account the diurnal variability of the fluxes of shortwave and longwave radiation and latent and sensible heat. Turbulence closure of the system is performed with a two-parameter 𝑘 - 𝜔 Wilcox-type model and algebraic relations for the coefficients of turbulent diffusion. The convection-diffusion equations are solved by a finite volume method to satisfy the integral conservation laws. The numerical algorithm for finding the flow and temperature fields is based on a Crank-Nicholson difference scheme. The convective terms in the equations are approximated with a second-order upstream scheme, QUICK. The systems of grid equations at each time step are solved by an under-relaxation method. Numerical modelling has demonstrated that the river water with high values of turbidity first moves along the slope to the deeper parts, then spreads over the lake surface. The results of simulation have shown qualitative agreement with in-situ measurement described by John et al. (1976). It was also found that the variable heat flux affects the space-time distribution of turbidity.
[full text]
Keywords: lake modelling, numerical simulation, water bodies contamination, water turbidity, Kamloops Lake
Author(s):
Tsydenov Bair Olegovich
PhD.
Position: Senior Research Scientist
Office: Tomsk State University
Address: 634050, Russia, Tomsk, 36, Lenin Avenu
Phone Office: (3822)783593
E-mail: tsydenov@math.tsu.ru
SPIN-code: 1832-2785 References:
[1] Danilov-Danilyan, V.I. Water a strategic factor for the long-term development of the Russian economy. Vestnik RAN. 2009; 79(9):789798. (In Russ.)
[2] Holland, P.R., Kay, A., Botte, V. Numerical modelling of the thermal bar and its ecological consequences in a river-dominated lake. Journal of Marine Systems. 2003; 43(1-2):6181.
[3] Lawrence, S.P., Hogeboom, K., Le Core, H.L. A three-dimensional general circulation model of the surface layers of Lake Baikal. Hydrobiologia. 2002; (487):95110.
[4] Vasiliev, O.F., Bocharov, O.B., Kvon, V.I., Ovchinnikova, T.E., Kvon, D.V. Numerical Modeling of Thermal Bars in Deep Lakes. Proc. of the 3-rd Int. Conf. on Hydro-Science and Engineering. Berlin, Cottbus; 1998: 20.
[5] Tsvetova, E.A. Features of convection under ice in a deep lake. Interexpo GEO-Siberia. 2016; 4(1):110114.
[6] Blokhina, N.S. The influence of wind on the development of a thermal bar and currents in reservoirs of different depths during ice cover melting. Moscow University Physics Bulletin. 2015; 70(4):319325.
[7] Farrow, D.E. A model for the evolution of the thermal bar system. EJAM. 2013; 24(2):161177.
[8] Tsydenov, B.O. Chislennoe modelirovanie effekta vesennego termobara v glubokom ozere [Numerical modeling of the effect of the spring thermal bar in a deep lake]. Tomsk: TGU; 2013: 145. (In Russ.)
[9] Hurley, P. The air pollution model (TAPM) Version 2. Part 1 : technical description. CSIRO Atmospheric Research technical paper ¹ 55.
Australia: CSIRO; 2002: 51.
[10] Aleksandrova, M.P., Gulev, S.K., Sinitsyn, A.V. An improvement of parametrization of short-wave radiation at the sea surface on the basis of direct measurements in the Atlantic. Russian Meteorology and Hydrology. 2007; 32(4):245251.
[11] Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988; 26(11):12991310.
[12] Tsydenov, B.O., Starchenko, A.V. Application of the two-parametric 𝑘-𝜔 turbulence model for studying the thermal bar phenomenon. Tomsk State University Journal of Mathematics and Mechanics. 2014; 5(31):104113. (In Russ.)
[13] Chen, C.T., Millero, F.G. Precise thermodynamic properties for natural waters covering only limnologies range. Limnology and Oceanography. 1986; 31(3):657662.
[14] Tsydenov, B.O., Starchenko, A.V. To the selection of heat flux parameterization models at the water-air interface for the study of the spring thermal bar in a deep lake. Proceedings of SPIE. 2015; (9680):18.
[15] Orlanski, I. A simple boundary condition for unbounded hyperbolic flows. Journal of Computational Physics. 1976; 21(3):251269.
[16] Patankar, S. Numerical heat transfer and fluid flow. New York: CRC Press; 1980: 214.
[17] Leonard, B. A stable and accurate convective modeling procedure based on quadratic upstream interpolation. Computer Methods in Applied Mechanics and Engineering. 1979; 19(1):5998.
[18] Tsydenov, B.O., Kay, A., Starchenko, A.V. Numerical modeling of the spring thermal bar and pollutant transport in a large lake. Ocean Modelling. 2016; (104):7383.
[19] Tsydenov, B.O. Numerical modelling convective flows in a cavity. Proc. of the International Conference of Students and Young Scientists "Prospects of FundamentalSciences Development". Tomsk: TPU; 2009; (2):673 676. (In Russ.)
[20] John, B.E.St., Carmack, E.C., Daley, R.J., Gray, C.B.J., Pharo, C.H. The limnology of Kamloops Lake, B.C., B.E. Vancouver; 1976: 167.
Bibliography link:
Tsydenov B.O. A numerical study of impurity propagation in a freshwater lake on the basis of water turbidity distribution // Computational technologies. 2017. V. 22. XVII All-Russian Conference of Young Scientists on Mathematical Modeling and Information Technology. P. 113-124
|
