Changeset 8724
- Timestamp:
- Mar 5, 2013, 1:44:12 AM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/validation_tests/Tests/Analytical_exact/dam_break_dry/results.tex
r8719 r8724 1 1 2 \section{Dam Break Involving a Dry Area}2 \section{Dam break involving a dry area} 3 3 4 The dam break problem involving a dry area was solved analytically by Ritter~\cite{Ritter1892} as well as Stoker~\cite{Stoker1948, Stoker1957}. The analytical solution exhibits a rarefaction fan as a parabolic curve. As water moves, it involves wetting process over the dry area. This is a one-dimensional-type problem.4 The dam break problem involving a dry area was solved analytically by Ritter~\cite{Ritter1892} as well as Stoker~\cite{Stoker1948, Stoker1957}. The analytical solution exhibits a rarefaction fan as a parabolic curve. As water moves, it involves wetting process over the dry area. 5 5 6 6 The initial condition is 7 \begin{equation} \label{eq: rdbp}8 u(x,0)=0 ~~,v(x,y)=0, ~~\textrm{and}~~7 \begin{equation} \label{eq:db_dry_init} 8 u(x,0)=0, ~~v(x,y)=0, ~~\textrm{and}~~ 9 9 h(x,0) = \left\{ \begin{array}{ll} 10 h_1 & \textrm{if $x < x_0$}\\11 0 & \textrm{if $x > x_0$}\\10 h_1 & \textrm{if $x < 0$}\\ 11 0 & \textrm{if $x > 0$}\\ 12 12 \end{array} \right. 13 13 \end{equation} 14 where $h_1>0$. 14 where $h_1>0$. The topography is a horizontal flat bed. 15 15 16 The analytical solution at time $t>0$ is~\cite{Ritter1892, Stoker1948, Stoker1957} 17 \begin{equation} \label{eq:h_sol_dry} 16 17 The analytical solution~\cite{Ritter1892, Stoker1948, Stoker1957} at time $t>0$ is 18 \begin{equation} 18 19 h(x) = \left\{ \begin{array}{ll} 19 20 h_1 & \textrm{if $x \leq -t \sqrt{gh_1}$}\\ … … 23 24 \end{equation} 24 25 which is the free surface and 25 \begin{equation} \label{eq:u_sol_dry}26 \begin{equation} 26 27 u(x) = \left\{ \begin{array}{ll} 27 28 0 & \textrm{if $x \leq -t \sqrt{gh_1}$}\\ … … 35 36 36 37 \subsection{Results} 37 For our test, we consider $ x_0=0$ and $h_1=10$ in (\ref{eq:rdbp}).38 For our test, we consider $h_1=10$ in (\ref{eq:db_dry_init}). 38 39 The following figures show the stage, $x$-momentum, and $x$-velocity at several instants of time. We should see excellent agreement between the analytical and numerical solutions. The wet/dry interface is difficult to resolve and it usually produces large errors. 39 40
Note: See TracChangeset
for help on using the changeset viewer.