Skip to content

A simple Fortran program of Discontinuous Galerkin method(no Limiter now) solving 2D Euler Equation with the Isentropic Vortex initial value.

Notifications You must be signed in to change notification settings

L1ssandra/Simple-Fortran-DG

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

8 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

这是我在研0学习DG的过程中练习用的程序,比较简单,不过我觉得我现在写得已经很直观了,有需要的话可以自取

几点说明:

  1. MATLAB中,contour命令是等高线,而mesh命令是画surface的。所以,如果要改成画等高线图,请把contour的注释取消,然后把mesh注释掉;另外把axis(s)这一句也要注释掉,它是一个三维的画布。如果要看数值解的其他分量,请把Q1 = str2num(fileread('Q2.txt'))这一句中后面的“Q2”改成其他的。
  2. Fortran中,DG_2D_New是主程序,RK3DG是子程序,现在默认的是用80×80的网格解Euler方程,初值为等熵涡流。
  3. 关于程序的详细说明可参见知乎笔记:https://zhuanlan.zhihu.com/p/533958396

Reference:

[1] Chi-Wang Shu, Discontinuous Galerkin Methods: General Approach and Stability.

[2] Bernardo Cockburn, Chi-Wang Shu, The Runge–Kutta Discontinuous Galerkin Method for Conservation Laws V: Multidimensional Systems.

About

A simple Fortran program of Discontinuous Galerkin method(no Limiter now) solving 2D Euler Equation with the Isentropic Vortex initial value.

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published