Awesome
DGFEM for Acoustic Wave Propagation
This repository implements a discontinuous Galerkin finite element method (DGFEM) applied to the linearized Euler equations and the acoustic perturbation equations. The solver is based on GMSH library and supports a wide range of features:
- 1D, 2D, 3D problems
- 4-th order Runge-Kutta
- High order elements
- Absorbing and reflecting boundaries
- Complex geometry and unstructured grid
For more information, a detailled report is available here.
| Auditorium | Isosurfaces | Bulk |
|---|---|---|
| <img src="https://raw.githubusercontent.com/pvanberg/DGFEM-Acoustic/master/assets/auditorium_source2_2.png" width="400" height="200" /> | <img src="https://raw.githubusercontent.com/pvanberg/DGFEM-Acoustic/master/assets/auditorium_source_iso1.png" width="400" height="200" /> | <img src="https://raw.githubusercontent.com/pvanberg/DGFEM-Acoustic/master/assets/auditorium_source_bulk1.png" width="400" height="200" /> |
Getting Started
Prerequisites
First, make sure the following libraries are installed. If you are running a linux distribution (ubuntu, debian, ...), an installation script is provided.
Gmsh
Eigen
Lapack
Blas
OpenMP
Installing
git clone https://github.com/pvanberg/MATH0471-DG.git
cd MATH0471-DG
mkdir build && cd build
cmake .. && make -j4
Running the tests
Once the sources sucessfully build, you can start using with the solver. It required two arguments: a mesh file created with Gmsh and a config file containing the solver options. Examples of mesh files and config files are given here.
cd bin
./dgalerkin mymesh.msh myconfig.conf
Minimal working example
2D propagation of an Gaussian initial condition over a square.
./dgalerkin ../../doc/2d/square.msh ../../doc/config/config.conf
Authors
- Pierre-Olivier Vanberg
- Martin Lacroix
- Tom Servais