From 26806930940273bec8c28a9757e22f2ce5fbdac3 Mon Sep 17 00:00:00 2001
From: SEOGI KANG <sgkang09@gmail.com>
Date: Wed, 17 Sep 2014 09:02:23 -0700
Subject: [PATCH] write docs

---
 docs/source/api_Acous.rst | 41 +++++++++++++++++++++++++++++++--------
 1 file changed, 33 insertions(+), 8 deletions(-)

diff --git a/docs/source/api_Acous.rst b/docs/source/api_Acous.rst
index a3b3791..68a3bb5 100644
--- a/docs/source/api_Acous.rst
+++ b/docs/source/api_Acous.rst
@@ -18,13 +18,16 @@
 Acoustic wave equation
 **********************
 
+Backgrounds
+-----------
+
 Time domain acoustic wave equation in 1st order form:
 
 .. math::
 
-        \rho\frac{\partial \phi}{\partial t} = \nabla \cdot \vec{u}+\frac{\partial s}{\partial t}\delta(\vec{r}-\vec{r}_s)
+    \rho\frac{\partial \phi}{\partial t} = \nabla \cdot \vec{u}+\frac{\partial s}{\partial t}\delta(\vec{r}-\vec{r}_s)
 
-        \mu^{-1}\frac{\partial \vec{u}}{\partial t} = \nabla \phi
+    \mu^{-1}\frac{\partial \vec{u}}{\partial t} = \nabla \phi
 
 where \\(\\rho\\) is the density, \\(\\mu\\) is the adiabatic compression modulus, \\(\\delta(\\vec{r}-\\vec{r}_s)\\) represents point source, and \\(\\phi\\) is the pressure.
 
@@ -39,13 +42,13 @@ By applying Fourier transform to frequency domain we have:
 
 .. math::
 
-        -\rho\omega^2\phi - \div\mu\grad\phi = -\omega ^2 s \delta(\vec{r}-\vec{r}_s)
+    -\rho\omega^2\phi - \div\mu\grad\phi = -\omega ^2 s \delta(\vec{r}-\vec{r}_s)
 
 where \\(\\omega = 2\\pi f\\) is angular frequency. Assuming contant \\(\\mu\\) and \\(\\rho\\) we have Helmholtz equation:
 
 .. math::
 
-        (\grad^2 + k^2)\phi = \mu^{-1}\omega^2 s \delta(\vec{r}-\vec{r}_s)
+    (\grad^2 + k^2)\phi = \mu^{-1}\omega^2 s \delta(\vec{r}-\vec{r}_s)
 
 where \\(\ k = \\omega\\sqrt{\\frac{\\rho}{\\mu}}=\\frac{\\omega}{v} \\) is the wave propagation constant.
 
@@ -53,12 +56,35 @@ where \\(\ k = \\omega\\sqrt{\\frac{\\rho}{\\mu}}=\\frac{\\omega}{v} \\) is the
 Discretization of problem
 -------------------------
 
-Artificial Boundary condition
------------------------------
+To compute the solution of time or frequency dependent partial differential equation (PDE)  as discussed above, we need to discretize those equations in both time and space. Our domain in real situation (i.e., earth) can be considered as infinite in space. However, in discrete space where we compute the solution of those PDE should be finite. Therefore, we need to implement artificial boundary conditions to get rid of out going wave, which propagates beyond our domain of interest. We first implement sponge boundary condition then perfectly matched boundary condition.
+
 
 Sponge boundary
 ===============
 
+A fundamental idea of sponge boundary condition is adding damping term to acoustic wave equation:
+
+.. math::
+
+    \rho\frac{\partial^2 \phi}{\partial t^2} + c\frac{\partial \phi}{\partial t}- \div\mu\grad\phi = \frac{\partial^2 s}{\partial t^2}\delta(\vec{r}-\vec{r}_s)
+
+Similarly in frequency domain we have:
+
+.. math ::
+
+    -\rho\omega^2\phi +\imath\omega c\phi- \div\mu\grad\phi = -\omega ^2 s \delta(\vec{r}-\vec{r}_s)
+
+
+The first order form of time domain equation can be written as:
+
+.. math ::
+
+    \rho\frac{\partial \phi}{\partial t} +c\phi- \nabla \cdot \vec{u}= +\frac{\partial s}{\partial t}\delta(\vec{r}-\vec{r}_s)
+
+    \mu^{-1}\frac{\partial \vec{u}}{\partial t} = \nabla \phi
+
+Using staggered grid, we discretize acoustic wave equation in time  (CITE):
+
 PML boundary
 ============
 
@@ -67,8 +93,7 @@ PML boundary
     :file: examples\refraction.html
 
 
-Backgrounds
-===========
+
 
 
 Notebooks