Skip to content

Commit

Permalink
Fixed equations in nonreflectingBC. Refs #3.
Browse files Browse the repository at this point in the history
  • Loading branch information
cbolisetti committed May 15, 2020
1 parent 166ce56 commit e44d0f2
Showing 1 changed file with 12 additions and 4 deletions.
16 changes: 12 additions & 4 deletions doc/content/source/bcs/NonReflectingBC.md
Original file line number Diff line number Diff line change
Expand Up @@ -4,11 +4,19 @@

## Description

This boundary condition applies a Lysmer damper [!citep](lysmer1969finite) on a given boundary to absorb the waves hitting the boundary. To understand Lysmer dampers, let us consider an uniform linear elastic soil column and say a 1D vertically propagating P wave is traveling through this soil column. Then the normal stress at any location in the soil column is given by:
This boundary condition models a Lysmer damper [!citep](lysmer1969finite) on a given boundary to absorb the waves hitting the boundary. To understand Lysmer dampers, let us consider an uniform linear elastic soil column and say a 1D vertically propagating P wave is traveling through this soil column. Then the normal stress at any location in the soil column is given by:

$$
\sigma = E \epsilon = E \frac{du}{dx} = \frac{E}{V_p} \frac{du}{dt}= \rho V_p \frac{du}{dt}, $$
where $E$ is the Young's modulus, $\sigma$ is the normal stress, $\epsilon$ is the normal strain, $\rho$ is the density, $V_p = \sqrt{\frac{E}{\rho}}$ is the P-wave speed and $\frac{du}{dt}$ is the particle velocity. Note that for a 3D problem, the P-wave speed is $V_p = \sqrt{\frac{E(1-\ nu)}{(1+\nu)(1-2\nu)}$.
\begin{equation}
\label{eqn:wave}
\sigma = E \epsilon = E \frac{du}{dx} = \frac{E}{V_p} \frac{du}{dt}= \rho V_p \frac{du}{dt}
\end{equation}

where $E$ is the Young's modulus, $\sigma$ is the normal stress, $\epsilon$ is the normal strain, $\rho$ is the density, $V_p = \sqrt{\frac{E}{\rho}}$ is the P-wave speed and $\frac{du}{dt}$ is the particle velocity. Note that for a 2D or a 3D problem, the P-wave speed is

\begin{equation}
\label{eqn:vp}
V_p = \sqrt{\frac{E(1-\nu)}{(1+\nu)(1-2\nu)}}
\end{equation}

The stress in the above equation is directly proportional to the particle velocity which makes this boundary condition analogous to a viscous damper with damping coefficient of $\rho V_p$. So truncating the soil domain and placing this damper at the end of the domain is equivalent to simulating wave propagation in an infinite soil column provided the soil is made of linear elastic material and the wave is vertically incident on the boundary.

Expand Down

0 comments on commit e44d0f2

Please sign in to comment.