/2D-Helmholtz-FEM-BEM-Coupling

Primary LanguagePythonMIT LicenseMIT

This is a Python implementation of FEM-BEM coupling for the Helmholtz equation

$$-\Delta u -k^2 u = \left\lbrace \begin{array}\end{array}\right.$$

Configuring

To configure, edit config/worldwide.py. Set the PATH variable to where you would like the graphs created to be saved.

Running

The files in the main directory will do the following when run:

helm.py             FEM-BEM coupling with 8 internal FEM nodes
helm-2.py           As above but with CIP FEM.
converging-helm.py  Draw graphs of the convergence to the solution
fix-helm.py         As above but with CIP FEM.
eigen.py            Find the eigenvalues of the system.