Numerical 2D wave simulation

A variety of numerical simulations of 2D wave equation using Finite Difference Method.
Simulation-based analysis of wave basic phenomena: diffraction and interference.

The Wave Equation

$u_{tt}=\nabla u \xrightarrow{} u_{tt}=u_{xx}+u_{yy}$
$\frac{\partial^2 u(x,y,t)}{\partial t^2}=\frac{\partial^2 u(x,y,t)}{\partial x^2}+\frac{\partial^2 u(x,y,t)}{\partial y^2}$

Initial conditions

$u_{t}(x,y,0)=g(x,y)$

To run all the simulations I have used the same IC

$u(x,y,0)=e^{-1000(x^2+y^2)}$
$u_{t}(x,y,0)=0$

This IC best simulates a puntual wave source caused by a droplet in a water surface for instance.

Boundary conditions

I have define a squared domain with boundaries sufficiently far away from the region of study that they will not have a particular influence. However I have run some tests with homogeneus Dirichlet conditions, which are the simplest ones.

Dirichlet conditions:

The Finite Difference Method

Based on numerical differentiation expressions, we can handle easily PDEs. Given that the basic wave equation consists of 2nd order partial derivatives, the same expresion can be used to simpify the whole equation:

$f''(x_0)\approx\frac{x_{-1}-2x_0+x_1}{h^2}$

Space-time discretization

In this problem, there are 3 variables involved: 2 space dimensions and time. So, we will perform a 3D discretization which consists of splitting every dimension in parts to perform numerical solution of the wave equation.

$x=x_0+i\Delta x$
$y=y_0+j\Delta y$
$t=k\Delta t$

For simplicity, mesh will be equal in both spatial dimensions:

$\Delta x=\Delta y=h$

Wave PDE algorithm

$u(x_i,y_j,t_k)=u_{i,j,k}$
$\frac{u_{i,j,k-1}-2u_{i,j,k}+u_{i,j,k+1}}{k^2}=\frac{u_{i-1,j,k}-2u_{i,j,k}+u_{i+1,j,k}}{h^2}+\frac{u_{i,j-1,k}-2u_{i,j,k}+u_{i,j+1,k}}{h^2}$
$\frac{u_{i,j,k-1}-2u_{i,j,k}+u_{i,j,k+1}}{k^2}=\frac{u_{i-1,j,k}+u_{i+1,j,k}+u_{i,j-1,k}+u_{i,j+1,k}-4u_{i,j,k}}{h^2}$
GENERAL ALGORITHM
$u_{i,j,k+1}=-u_{i,j,k-1}+2u_{i,j,k}+k^2\frac{u_{i-1,j,k}+u_{i+1,j,k}+u_{i,j-1,k}+u_{i,j+1,k}-4u_{i,j,k}}{h^2}$

Name	Name	Last commit message	Last commit date
Latest commit History 8 Commits 8 Commits
README.md	README.md
double_slit.m	double_slit.m
wave2d.m	wave2d.m

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Repository files navigation

Numerical 2D wave simulation

The Wave Equation

Initial conditions

Boundary conditions

The Finite Difference Method

Space-time discretization

Wave PDE algorithm

Single slit experiment. Difraction pattern

Double slit experiment. Difraction-interference pattern

About

Uh oh!

Releases

Packages

Uh oh!

Contributors

Uh oh!

Languages

Search code, repositories, users, issues, pull requests...

Folders and files

Latest commit

History

Repository files navigation

Numerical 2D wave simulation

The Wave Equation

Initial conditions

Boundary conditions

The Finite Difference Method

Space-time discretization

Wave PDE algorithm

Single slit experiment. Difraction pattern

Double slit experiment. Difraction-interference pattern

About

Resources

Uh oh!

Stars

Watchers

Forks

Releases

Packages 0

Uh oh!

Contributors

Uh oh!

Languages

Packages