write docs

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