Changeset 8792

Timestamp:

Mar 30, 2013, 12:56:51 AM (12 years ago)

Author:

mungkasi

Message:

Automated report for rundown_mild_slope.

Location:

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/rundown_mild_slope

Files:

• plot_results.py (modified) (1 diff)
• produce_report.py (added)
• produce_results.py (modified) (1 diff)
• report.tex (added)
• results.tex (modified) (2 diffs)
• run_channel.py (modified) (2 diffs)

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/rundown_mild_slope/plot_results.py

@@ -46 +46 @@
 #--------------------
 pyplot.clf()
-pyplot.ion()
 
 line, = pyplot.plot( (p2.x[v].min(),p2.x[v].max()) ,( (p2.stage[:,v]-p2.elev[:,v]).max(),(p2.stage[:,v]-p2.elev[v]).min() ) )

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/rundown_mild_slope/produce_results.py

@@ -15 +15 @@
     run_validation_script('run_channel.py')
     run_validation_script('plot_results.py')
+    run('pdflatex', 'report.tex')
+    run('bibtex', 'report')
+    run('pdflatex', 'report.tex')
+    run('bibtex', 'report')
+    run('pdflatex', 'report.tex')
+    run('bibtex', 'report')     
 
 def clean():

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/rundown_mild_slope/results.tex

@@ -5 +5 @@
 
 \section{Shallow flow down a mild slope}
-This case simulates very shallow flow running down a mild slope topography. It represents an idealisation of the rainfall-runoff problem, which will often involve very shallow flows down such a topography. This case has an analytical solution (the steady-uniform solution, bed slope = friction slope).   
+This case simulates very shallow flow running down a mild slope topography. It represents an idealisation of the rainfall-runoff problem, which will often involve very shallow flows down such a topography. This case has an analytical solution, and in particular, we consider the steady-uniform solution with the values of bed slope and friction slope are the same.   
+
+Suppose that we are given a one dimensional domain. The steady state conditions with a contant water depth everywhere make the shallow water equations to the single identity
+\begin{equation}
+z_x = - S_f.
+\end{equation}
+Here $q=uh$ is the momentum or water discharge and $S_f$ is the symbol for the force of bottom friction involving Manning's coefficient $n$. We take 
+\begin{equation}
+S_f = n^2 \frac{q|q|}{h^{10/3}}.
+\end{equation}
+If $q$, $n$, and $z_x$ are given, then the analytical solution is
+\begin{equation}
+u(x)= \left[- n^{-2} q^{4/3} z_x\right]^{3/10},
+\end{equation}
+\begin{equation}
+h(x)= \frac{q}{u}\,.
+\end{equation}
 
 \subsection{Results}
+For our test, we consider a square dimensional domain with length and width 100. We take $q=0.2$, $n=0.03$, and $z_x=-0.1$.
+The topography is
+\begin{equation}
+z(x, y)= -0.1 x\,.
+\end{equation}
+The initial condition is $u=v=0$ and
+\begin{equation}
+w(x,y,0)= -0.1 x + 0.01\,.
+\end{equation}
+
+
+Some simualtion results are as follows.
 Figures~\ref{fig:depthdownchan} shows the steady state depth in the downstream direction. There should be a good agreement with the analytical solution, at least away from the boundaries.  
+Figures~\ref{fig:xvelscrosschan} and~\ref{fig:yvelscroschan} show the steady state $x$- and $y$-velocities, along a slice in the cross slope direction (near $x=50$). The $x$-velocities should agree well with the analytical solution, and the $y$-velocities should be zero. 
 
 \begin{figure}[h]
@@ -17 +46 @@
 \end{center}
 \end{figure}
-
-Figures~\ref{fig:xvelscrosschan} and~\ref{fig:yvelscroschan} show the steady state $x$- and $y$-velocities, along a slice in the cross slope direction (near $x=50$). The $x$-velocities should agree well with the analytical solution, and the $y$-velocities should be zero.  
+ 
 
 \begin{figure}[h]

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/rundown_mild_slope/run_channel.py

@@ -26 +26 @@
 from anuga.utilities.argparsing import parse_standard_args
 alg, cfl = parse_standard_args()
-alg = '2_0'
 domain.set_flow_algorithm(alg)
 domain.set_CFL(cfl)
@@ -60 +59 @@
 Br = anuga.Reflective_boundary(domain) # Solid reflective wall
 domain.set_boundary({'left': BdIN, 'right': Bt, 'top': Br, 'bottom': Br})
+
+
+#------------------------------------------------------------------------------
+# Produce a documentation of parameters
+#------------------------------------------------------------------------------
+parameter_file=open('parameters.tex', 'w')
+parameter_file.write('\\begin{verbatim}\n')
+from pprint import pprint
+pprint(domain.get_algorithm_parameters(),parameter_file,indent=4)
+parameter_file.write('\\end{verbatim}\n')
+parameter_file.close()
+
 #------------------------------------------------------------------------------
 # Evolve system through time