حل تحلیلی معادله ی انتقال جرم در شبکه ی رودخانه

در این پژوهش، با روش تبدیل لاپلاس حل تحلیلی معادله‌ی انتقال آلودگی در دامنه‌ی محدود برای شبکه‌ی رودخانه به ازای شرایط مرزی بالادست و پایین‌دست دیریکله و شرط اولیه‌ی صفر استخراج شد و شبیه‌سازی برای دو شبکه‌ی شاخه‌یی و حلقه‌یی با شرایط مرزی ثابت و متغیر انجام شد. بعد از نام‌گذاری گره‌ها، با تشکیل ماتریس‌های نحوه‌ی اتصال، مشخصات جریان و هندسی رودخانه برای هر شبکه به عنوان ورودی مسئله، ماتریس انتشار براساس تابعی از متغیر s لاپلاس ایجاد می‌شود که با حل دستگاه پیچیده‌ی ایجاد شده و استفاده از تبدیل لاپلاس معکوس، مقدار غلظت در هر گره محاسبه می‌شود. در نهایت برای اعتبارسنجی، حل تحلیلی با حل عددی مقایسه و سپس شاخص‌های خطای آماری محاسبه شد. در شبیه‌سازی‌های انجام شده، به علت تغییر دبی ورودی به هر شاخه و تغییر سطح مقطع جریان در هر شاخه، در محل‌های اتصال شاخه‌ها به یکدیگر تغییرات غلظت آلودگی به صورت افزایشی یا کاهشی رخ می‌دهد. نتایج حاکی از عملکرد مطلوب و توانایی بالای حل تحلیلی در مدل‌سازی دو شبکه‌ی مذکور و انطباق خوب آن با حل عددی است که به دلیل دقت و سرعت محاسبات بالا می‌تواند جایگزین حل عددی شود. همچنین از حل تحلیلی موجود می‌توان به عنوان ابزاری برای صحت‌سنجی سایر حل‌های تحلیلی در شبکه‌ی رودخانه استفاده کرد.

Analytical solution of mass transport equation in river network

In this study, by Laplace transform method, the analytical solution of the pollution transport equation in the limited domain for the river network to the upstream and downstream Dirichlet boundary conditions, and the initial condition of zero was extracted, and simulation was performed for two branch and loop networks with fixed and variable boundary conditions. After naming the nodes, by forming matrices of how to connect, flow characteristics and geometry of the river for each network as input to the problem, the diffusion matrix is created based on a function of the Laplace variable, which The value of the concentration in each node is calculated by solving the complex device created and using the inverse Laplace transform. Then, using the analytical solution extracted in a branch of the river for the pollution transfer equation, the analytical solution can calculate the value of pollution concentration at any desired location and time along with the river network. Finally, for validation, the analytical solution was compared with the numerical solution, and then the statistical error indices were calculated. The results indicate the optimal performance and high ability of analytical solution in modeling the two networks and its good adaptation to the numerical solution, which can replace numerical solution due to high accuracy and speed of calculations. Generally, due to common errors in numerical solutions such as numerical dispersion error, Roundoff error, Truncation error of Taylor expansion mathematical sentences, analytical solutions, if any, for the river network are recommended over numerical solutions. Also, in the performed simulations, due to the change of the inlet flow and the flow crosssection to each, changes in the pollution concentration occur in the areas where the branches connect to each other, increasing or decreasing. The proposed analytical solution for river networks can model more complex river networks and can be considered a criterion for the validation of numerical solutions. Also, the existing analytical solution can be used as a tool to validate other analytical solutions in the river network.