پديد آورندگان :
عطائيان، آتوسا دانشگاه تربيت مدرس تهران - گروه مهندسي سازه هاي آبي , ايوب زاده، علي دانشگاه تربيت مدرس تهران - گروه مهندسي سازه هاي آبي , نبوي، عبدالرضا دانشگاه تربيت مدرس تهران - دانشكده مهندسي برق- الكترونيك , انخل گومز-لوپرا، سالوادور دانشگاه پليتكنيك كارتاخنا اسپانيا - گروه فيزيك كاربردي
كليدواژه :
انتقال جرم , نگهداشت موقت , مدارهاي الكتريكي معادل , مدل متشابه الكتريكي , روش احجام محدود , رودخانه
چكيده فارسي :
با توجه به آثار زيان بار آلاينده هاي ورودي به منابع آبي كه ناشي از فعاليت هاي انساني مي باشد، انجام مطالعات در زمينه پيش بيني غلظت اين آلاينده ها به منظور انجام اقدامات مقتضي براي كنترل، ضروري به نظر مي رسد. بدين منظور مدل هاي انتقال جرم متعددي براي پيش بيني غلظت محلول در آبراهه هاي طبيعي ارايه شده اند. اين مدل ها به ندرت داراي حل تحليلي بوده و عمدتا با استفاده از روش هاي عددي حل مي شوند. در اين تحقيق روشي تحت عنوان روش شبيه سازي مداري (NSM) كه بر پايه آنالوژي موجود بين معادلات ديفرانسيلي حاكم بر پديده هاي هيدروديناميكي و الكتريكي است، معرفي شده و كاربرد آن در حل مدل هاي مذكور مورد بررسي قرار گرفته است. مراحل اجراي اين روش شامل استخراج مدل الكتروآنالوژيكال و طراحي مدار الكتريكي معادل و در نهايت شبيه سازي مدل مداري با استفاده از كد مناسب در يك نرم افزار تحليل مدارهاي الكتريكي است. در اين مقاله ابتدا NSM با مدلسازي حالتي از معادله انتقال نگهداشت موقت كه حل تحليلي دارد، صحت سنجي شده و سپس دقت و كارايي آن در مقايسه با روش عددي احجام محدود (FVM) در حل مدل نگهداشت موقت تودرتو، برآورد شده است. نتايج مدلسازي ها حاكي از مطابقت بسيار عالي بين دو روش NSM و FVM با شاخص هاي خطاي نزديك به صفر است. حال آن كه پياده سازي شرايط مرزي در NSM ساده تر بوده و انعطاف پذيري بالاتري دارد. به علاوه، زمان محاسباتي مورد نياز NSM براي مثال هاي مورد مطالعه كمتر از زمان محاسباتي مورد نياز FVM مي باشد. بنابراين NSM به عنوان جايگزيني دقيق و كارا براي روش هاي عددي در حل معادلات انتقال توام يك بعدي پيشنهاد مي شود.
چكيده لاتين :
Human activities everyday release a huge amount of domestic, industrial and agricultural waste into water bodies and continuously change the ecosystem conditions in the world. Considering the harmful effects of these pollutants entering water resources, study about pollution transfer in streams and predicting the pollutant concentration at downstream points seem to be important. For this purpose, the well-known classical advection-dispersion equation (ADE) was presented as the first attempt for describing mass transfer and energy transfer in physical systems. This equation is useful for channels with relatively prismatic and uniform cross-sections. Experimental studies carried out in rivers show that ADE is no longer applicable for natural streams, especially mountain pool and riffle streams because of their irregular cross-sections. Afterward, some more accurate models, referred dead zone models or transient storage models, were suggested by several researchers for predicting solute concentration in natural rivers and calibrated using tracer approach. Such models cause more realistic concentration-time distributions which have lower picks and longer residence time. Solving such models, for which in most cases the analytical solution doesn’t exist, needs numerical methods –methods which usually deal with complexity and is time-consuming. In this study, we have applied Network Simulation Method (NSM) –a powerful and efficient computational method for simulating systems governed by differential equations based on the electric circuit concepts and the analogy between the governing differential equations of hydrodynamic and electrical phenomena– which according to the previous studies simulates desirably the transport of mass in natural streams, to solve two transient storage models. The method consists of two phases of designing and simulating. In designing phase, the system of differential equations corresponding to the prototype must be discretized spatially over the studied domain and then, for each term of the discretized equations the equivalent electrical devices are chosen. These electrical elements are connected based on the algebraic sign of the terms to satisfy Kirchhoff’s current low. Regarding the mathematical models, in most studies, electric potential and electric current are equivalent to the value of the unknown variable and its flux, respectively. The last step of designing the electro-analogical model is the implementation of initial conditions and boundary conditions of unknown variables using appropriate dependent and/or independent, voltage and/or current sources. Simulating this equivalent electrical network is performed through an appropriate electrical-computational circuit code, such as PSpice code. PSpice, which is a powerful circuit analysis software, uses the Newton-Raphson iterative algorithm to solve this set of nonlinear equations and performing the transient analysis. In this paper, NSM is firstly verified by simulating a transient storage transport model developed by Bencala and Walters (1983) for unsteady conservative solute transfer in pool and riffle streams. This model includes two equations for solute concentration in the main channel and in the storage zone and involves one storage zone. The analytical solution for this model has been presented in Laplace domain by Kazezyılmaz-Alhan (2008) considering a hypothetical channel with a constant cross-sectional area, flow velocity, and dispersion coefficient and for two types of upstream boundary conditions including a continuous injection and a pulsating solute injection. The results of this verification were desirable. Then, the accuracy and efficiency of NSM were compared with Finite Volume Method (FVM) –a widely used numerical method in computational fluid dynamics- through simulating an unsteady reactive solute transport using a nested two-storage zone transport model developed by Kerr et al. (2013). This model consists of three equations and involves two storage zones including the surface and hyporheic storage zones interacting together. The results of simulating a hypothetical solute transport problem with this nested model indicate a good match between these two methods with near-zero error indices. The computational time needed for NSM and FVM were 117 seconds and 505 seconds, respectively. So, NSM is much faster. Furthermore, the implementation of boundary conditions in NSM is direct, easier and more flexible. Therefore, NSM is proposed as a precise and efficient alternative for numerical methods in solving one-dimensional coupled differential equations of unsteady transport, simultaneously and providing benchmarks without complex mathematical calculations. Because of its analogical based concept, it can be used as a predicting and monitoring tool for transport phenomena instead of using troublesome physical hydraulic models to perform the water quality studies with less time, low expense and higher accuracy. Hence, in critical conditions, including a sudden spill of a high-hazardous contaminant in a specified point of the river or increasing the concentration of a chemical element to its maximum level, the monitoring and controlling measures at different parts of the river can be carried out with an acceptable accuracy and speed to improve the water quality.
