حل عددي موج سينماتيك با استفاده از روش ديفرانسيل كوادريچر

Numerical solution of kinematic wave equation using differential quadrature method

يكي از معادلات ساده و كاربردي كه در رونديابي هيدروليكي و هيدرولوژيكي استفاده مي شود، معادله موج سينماتيك است. با توجه به اينكه اين معادله فقط در شرايط خاصي حل تحليلي دارد، استفاده از روش هاي عددي براي حل آن رواج يافته است و يافتن يك روش عددي مناسب براي حل اين معادله، تلاش بسياري از محققين را به خود اختصاص داده است. روش ديفرانسيل كوادريچر (DQM)، يكي از روش هاي عددي است كه به دليل پايداري و كارايي، در زمينه هاي مختلف، مورد توجه محققين قرار گرفته است. در تحقيق حاضر از قابليت هاي اين روش براي حل عددي موج سينماتيك استفاده شده است. به اين منظور، روش DQ با ساختارهاي مختلف شبكه و با استفاده از توابع پايه اي مختلف، بكارگرفته شده است. براي ارزيابي روشDQ ، نتايج اين روش با نتايج حل تحليلي و داده هاي مشاهده اي مقايسه شده است. بررسي ها نشان داد، روش DQ نسبت به انتخاب تابع پايه اي حساسيت چنداني ندارد ولي در مورد توزيع نقاط شبكه، توزيع كسينوسي نسبت به توزيع يكنواخت بهتر عمل مي كند. به طور كلي روش DQ در حل معادله موج سينماتيك حتي با تعداد نقاط كم جواب هاي دقيقي ارايه مي-دهد و مدل عددي تهيه شده پايداري و كارايي خوبي دارد.

Kinematic wave equation is a simplification of the fully dynamic wave equation of the fluid flow in open channels. In this equation, the inertia and pressure forces are neglected compared to the gravity force and the friction. Due to the simplicity of kinematic wave, many researchers are still interested in this model in its range of applicability.Until now, various numerical methods such as Finite Difference (FD) and Finite Element (FE) have been used for solution of kinematic wave equation. Differential Quadrature Method (DQM) as an alternative numerical method has attracted the attention of some researchers because of its stability, accuracy and efficiency. In this study, the applicability of DQM for solution of the kinematic wave equation is investigated. For evaluation of this method, the results were compared with analytical solution and observed data.Results of this study show that the type of grid distribution affects the results significantly. The cosine grid distribution is better than uniform grid distribution as it was suggested by previous applications in the literature. DQM for simulating kinematic wave model is not sensitive to choice of the test functions. In general, DQM can be considered as another numerical method for kinematic wave equation using less number of grid points to achieve accurate results. It is also unconditionally stable. The main disadvantage of the kinematic wave equation is its failure to predict wave attenuation. This can be resolved by using diffusion wave model and the same numerical technique.