Yıl 2018, Cilt 22, Sayı 2, Sayfalar 364 - 382 2018-04-01

Computer Code Development for the Numerical Solution of Two Dimensional Shallow Flow Equations on Unstructured Grid
İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi

Nuray Öktem [1]

170 261

In this paper a novel computer code is developed for numerical solution of two dimensional shallow flow equations on unstructured grid. The code is tested by running on some benchmark problems. The numerical method is based on a cell centered finite volume method applied on an unstructured triangular mesh. Weighted Averaged Flux (WAF) method is used for the computation of the fluxes in continuity and momentum equations. The most prominent property of the shallow flow equations is the discontinuities due to shock waves occuring on the free surface of the flow. Therefore, WAF is combined with HLLC Riemann solvers in order to make the method capable to describe the shock waves and also have the opportunity to use flux limiters in the computation of the intercell fluxes. Besides, the Total Variation Diminishing (TVD) theory is utilized for smoothing the oscillations which may occur due to the second order accuracy of the numerical solution. The ratio of the upwind to local variations involved in the interface flux limiter functions available in the literature for the structured grid systems is described here for the first time for the unstructured grids and it is made usable with the new gradient approach. Test results show that the numerical method and the code developed are successful and can be applied to real cases. 

 Bu çalışmada, iki boyutlu Sığ Akım Denklemleri’nin (SAD) sırasız ağ üzerinde sayısal çözümü için özgün bir yazılım geliştirilmiştir. Geliştirilen yazılım bazı ölçüt problemler için koşturularak test edilmiştir. Sayısal yöntem sırasız üçgensel bir çözüm ağı üzerinde hücre merkezli sonlu hacim yöntemine dayanmaktadır. Süreklilik ve momentum denklemlerindeki akı hesabı için ikinci derece doğruluklu Ağırlık Ortalamalı Akı (Weighted Average Flux, WAF) yöntemi kullanılmıştır. Sığ akım denklemlerinin en belirgin özelliği su yüzeyinde oluşabilecek şok dalgalarından kaynaklanan süreksizliklerdir. Sayısal yöntemin şok dalgalarını ayrıntılı bir şekilde tanımlayabilmesi ve hücreler arası akı hesabında akı sınırlayıcılarını kullanmaya olanak vermesi için WAF yöntemi HLLC Riemann çözücüleriyle birleştirilmiştir. Ayrıca, ikinci dereceden doğruluklu sayısal çözüm nedeniyle oluşabilecek sayısal salınımları söndürmek için de Toplam Salınım Azaltma (Total Variation Diminishing, TVD) teorisinden faydalanılmıştır. Literatürde sıralı çözüm ağları için mevcut olan ara-yüz akı sınırlayıcı fonksiyonlarının içerdiği rüzgar yönlü değişimlerin yerel değişimlere oranı, sırasız çözüm ağları için ilk kez tanımlanmış ve yeni gradyan yaklaşımlarıyla beraber kullanılabilir hale getirilmiştir. Test sonuçları, sayısal yöntemin ve geliştirilen yazılımın doğru çalıştığını ve gerçek problemlere uygulanabilirliğini göstermektedir.

  • E.F. Toro, "Riemann Problems and the WAF Method for Solving the Two-Dimensional Shallow Water Equations," Philosophical Transactions: Physical Sciences and Engineering, vol. 338, no. 1649, pp. 43-68, 1992.
  • S. J. Billett and E. F.Toro, "On WAF-Type Schemes for Multidimensional Hyperbolic Conservation Laws," Journal of Computational Pyhsics, vol. 130, pp. 1-24, 1997.
  • C. Zoppou and S. Roberts, "Numerical solution of the two-dimensional unsteady dam break," Applied Mathematical Modelling, vol. 24, pp. 457-475, 2000.
  • D.H. Kim, Y.S. Cho, A.M.ASCE, and W.G. Kim, "Weighted Averaged Flux-Type Scheme for Shallow-Water Equations with Fractional Step Method," Journal of Engineering Mechanics, 2004, 130(2): 152-160, vol. 130, no. 2, pp. 152-160, 2004.
  • W.D. Guo, J.S. Lai, and G.F. Lin, "Finite-volume multi-stage schemes for shallow-water flow," International Journal for Numerical Methods in Fluids 2008; 57:177–204, vol. 57, pp. 177-204, 2008.
  • K. Anastasiou and C.T. Chan, "Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes," International Journal for Numerical Methods in Fluids , vol. 24, pp. 1225-1245, 1997.
  • D. Pan and J. Cheng, "A second-order upwind finite-volume method for the Euler solution on unstructured triangular meshes," International Journal for Numerical Methods in Fluids, vol. 16, pp. 1079-1098, 1993.
  • J.Hou, Q. Liang, H.Zhang, and R. Hinkelmann, "An efficient unstructured MUSCL scheme for solving the 2D shallow water equations," Environmental Modelling & Software, vol. 66, pp. 131-152, 2015.
  • T.H. Yoon, F. ASCE, and S.K. Kang, "Finite volume model for two-dimensional shallow flows on unstructured grids," Journal of Hydraulic Engineering , vol. 130, no. 7, pp. 678-688, 2004.
  • C. Zoppou, M.ASCE, and S. Roberts, "Explicit schemes for dam-break simulations," Journal of Hydraulic Engineering, vol. 129, pp. 11-34, 2003.
  • P.O. Persson and G. Strang, "A simple mesh generator in matlab," SIAM, Society for Industrial and Applied Mathematics, vol. 46, no. 2, pp. 329-345, 2004.
  • E. F. Toro, Shock-capturing methods for free-surface shallow flows, Chichester: John Wiley & Sons, Ltd, 2001, pp. 15-24.
  • E.F. Toro, Riemann solvers and numerical methods for fluid dynamics, Berlin, Heidelberg: Springer-Verlag, 2009.
  • H.K. Versteeg and W. MalalaSekera, An introduction to computational fluid dynamics, the finite volume method, Pearson: Prentice Hall, 2007, pp. 321-329.
  • S. Tu and S. Aliabadi, "A slope limiting procedure in discontinuous galerkin finite element method for gasdynamics applications," International Journal for Numerical Analysis and Modeling, vol. 2, no. 2, pp. 163-178, 2005.
  • T. Hino, L. Martinelli and A. Jameson, "A finite volume method with unstructured grid for free surface flow simulations," in Sixth International Conference on Numerical Ship Hydrodynamics, Tokyo, Japan, 1993.
  • R.J. Fennema and M.H. Chaudhry, "Explicit methods for 2D transient free-surface flows," Journal of Hydraulic Engineering, ASCE, vol. 116, pp. 1013-1034, 1990.
  • F. Alcrudo and P. Garcia-Navarro, "A high-resolution godunov-type scheme in finite volumes for the 2D shallow-water equations," International Journal for Numerical Methods in FLuids, VOL. 16,489-505 (1993), vol. 16, pp. 489-505, 1993.
  • L.A. Monthe, F. Benkhaldoun, and I. Elmahi, "Positivity preserving finite volume Roe schemes for transport-diffusion equations," Computer Methods in Applied Mechanics and Engineering, vol. 178, pp. 215-232, 1999.
  • Youssef Loukili and Azzeddine Soulaimani, "Numerical tracking of shallow water waves by the unstructured finite volume WAF Approximaton," International Journal for Computational Methods in Engineering Science and Mechanics, vol. 8, pp. 1-14, 2007.
Birincil Dil tr
Konular İnşaat Mühendisliği, Matematik ve İstatistik
Yayımlanma Tarihi Nisan 2018
Dergi Bölümü Araştırma Makalesi
Yazarlar

Yazar: Nuray Öktem
Kurum: Ankara Yıldırım Beyazıt Üniversitesi, Mühendislik ve Doğa Bilimleri Fakültesi, Matematik Bölümü
Ülke: Turkey


Bibtex @araştırma makalesi { saufenbilder328076, journal = {Sakarya University Journal of Science}, issn = {1301-4048}, eissn = {2147-835X}, address = {Sakarya Üniversitesi}, year = {2018}, volume = {22}, pages = {364 - 382}, doi = {10.16984/saufenbilder.328076}, title = {İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi}, key = {cite}, author = {Öktem, Nuray} }
APA Öktem, N . (2018). İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi. Sakarya University Journal of Science, 22 (2), 364-382. DOI: 10.16984/saufenbilder.328076
MLA Öktem, N . "İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi". Sakarya University Journal of Science 22 (2018): 364-382 <http://www.saujs.sakarya.edu.tr/issue/30829/328076>
Chicago Öktem, N . "İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi". Sakarya University Journal of Science 22 (2018): 364-382
RIS TY - JOUR T1 - İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi AU - Nuray Öktem Y1 - 2018 PY - 2018 N1 - doi: 10.16984/saufenbilder.328076 DO - 10.16984/saufenbilder.328076 T2 - Sakarya University Journal of Science JF - Journal JO - JOR SP - 364 EP - 382 VL - 22 IS - 2 SN - 1301-4048-2147-835X M3 - doi: 10.16984/saufenbilder.328076 UR - http://dx.doi.org/10.16984/saufenbilder.328076 Y2 - 2017 ER -
EndNote %0 Sakarya University Journal of Science İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi %A Nuray Öktem %T İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi %D 2018 %J Sakarya University Journal of Science %P 1301-4048-2147-835X %V 22 %N 2 %R doi: 10.16984/saufenbilder.328076 %U 10.16984/saufenbilder.328076
ISNAD Öktem, Nuray . "İki Boyutlu Sığ Akım Denklemlerinin Sırasız Ağda Sayısal Çözümü için Bilgisayar Yazılımı Geliştirilmesi". Sakarya University Journal of Science 22 / 2 (Nisan 2018): 364-382. http://dx.doi.org/10.16984/saufenbilder.328076