Changeset 8724


Ignore:
Timestamp:
Mar 5, 2013, 1:44:12 AM (12 years ago)
Author:
mungkasi
Message:

Updating report/results.tex for dam break dry.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/validation_tests/Tests/Analytical_exact/dam_break_dry/results.tex

    r8719 r8724  
    11
    2 \section{Dam Break Involving a Dry Area}
     2\section{Dam break involving a dry area}
    33
    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.
     4The 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.
    55
    66The 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}
     8u(x,0)=0, ~~v(x,y)=0, ~~\textrm{and}~~
    99h(x,0) = \left\{ \begin{array}{ll}
    10 h_1 & \textrm{if $x < x_0$}\\
    11 0 & \textrm{if $x > x_0$}\\
     10h_1 & \textrm{if $x < 0$}\\
     110 & \textrm{if $x > 0$}\\
    1212\end{array} \right.
    1313\end{equation}
    14 where $h_1>0$.
     14where $h_1>0$. The topography is a horizontal flat bed.
    1515
    16 The analytical solution at time $t>0$ is~\cite{Ritter1892, Stoker1948, Stoker1957}
    17 \begin{equation} \label{eq:h_sol_dry}
     16
     17The analytical solution~\cite{Ritter1892, Stoker1948, Stoker1957} at time $t>0$ is
     18\begin{equation}
    1819h(x) = \left\{ \begin{array}{ll}
    1920h_1 & \textrm{if $x \leq -t \sqrt{gh_1}$}\\
     
    2324\end{equation}
    2425which is the free surface and
    25 \begin{equation} \label{eq:u_sol_dry}
     26\begin{equation}
    2627u(x) = \left\{ \begin{array}{ll}
    27280 & \textrm{if $x \leq -t \sqrt{gh_1}$}\\
     
    3536
    3637\subsection{Results}
    37 For our test, we consider $x_0=0$ and $h_1=10$ in (\ref{eq:rdbp}).
     38For our test, we consider $h_1=10$ in (\ref{eq:db_dry_init}).
    3839The 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.
    3940
Note: See TracChangeset for help on using the changeset viewer.