Changeset 8773

Timestamp:

Mar 20, 2013, 5:55:11 PM (12 years ago)

Author:

mungkasi

Message:

Adding and modifying automated report for Carrier--Greenspan transient.

Location:

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/carrier_greenspan_transient

Files:

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

trunk/anuga_core/source/anuga_validation_tests/Analytical_exact/carrier_greenspan_transient/numerical_cg_transient.py

@@ -107 +107 @@
 
 #------------------------------------------------------------------------------
+# 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
 #------------------------------------------------------------------------------

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

@@ -44 +44 @@
 #Plot the stages##############################################################
 pyplot.clf()
-pyplot.ion()
 pyplot.plot(p2_st.x[v2], p2_st.stage[0,v2], 'b.', label='numerical')
 pyplot.plot(p2_st.x[v2], p2_st.stage[1,v2], 'b.') 
@@ -61 +60 @@
 #Plot the momentums##########################################################
 pyplot.clf()
-pyplot.ion()
 pyplot.plot(p2_st.x[v2], p2_st.xmom[0,v2], 'b.', label='numerical')
 pyplot.plot(p2_st.x[v2], p2_st.xmom[1,v2], 'b.') 
@@ -77 +75 @@
 #Plot the velocities#########################################################
 pyplot.clf()
-pyplot.ion()
 pyplot.plot(p2_st.x[v2], p2_st.xvel[0,v2], 'b.', label='numerical')
 pyplot.plot(p2_st.x[v2], p2_st.xvel[1,v2], 'b.')

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

@@ -10 +10 @@
 def build():
     run_validation_script('numerical_cg_transient.py')
-    run_validation_script('plot_results.py')   
+    run_validation_script('plot_results.py')  
+    run('python', 'produce_report.py')   
 
 def clean():

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

@@ -1 +1 @@
 
-\section{Carrier-Greenspan transient solution}
+\section{Carrier--Greenspan transient solution}
 
-The Carrier-Greenspan transient solution for flows on a sloping beach. References: Carrier and Greenspan, Journal of Fluid Mechanics, 1958.
+A transient solution for flows on a sloping beach was proposed by Carrier and Greenspan~\cite{CG1958}. The water moves to the shore at an early time, then it becomes still when time is large.
+
+Consider the dimensionless shallow water equations, as presented in the Carrier--Greenspan periodic solution.
+
+The analytical solution is:
+\begin{equation}
+w = - \frac{u^2}{2} + \epsilon {Re} 
+\left[1- 2 \frac{5/4 - i\lambda}{\left\{(1-i\lambda)^2 + \sigma^2 \right\}^{3/2}}
++ \frac32 \frac{(1-i\lambda)^2}{\left\{ (1-i\lambda)^2 + \sigma^2 \right\}^{5/2}} \right],
+\end{equation}
+\begin{equation}
+u = \frac{8\epsilon}{a} {Im} \left[ \frac{1}{\left\{(1-i\lambda)^2 + \sigma^2 \right\}^{3/2}}
+- \frac34 \frac{1-i\lambda}{\left\{ (1-i\lambda)^2 + \sigma^2 \right\}^{5/2}}    \right],
+\end{equation}
+where
+\begin{equation}
+t = \frac12 a\lambda -u\,, \quad c = \frac14 a\sigma\,,
+\end{equation}
+in which $c=\sqrt{gh}$ is the wave propagation speed.
+Here $\sigma \geq 0$ and we take $a=1.5(1+0.9\epsilon)^{1/2}$. Carrier and Greenspan~\cite{CG1958} observed that the waves do not break if $\epsilon$ is very small, namely $\epsilon \leq 0.23$. Setting $\sigma=0$ into this solution, we get the motion of the shoreline.
+
+The initial condition is given by setting time $t=0$ in this analytical solution. Note that this analytical solution is defined in the dimensionless space. To implement this in the numerical test, we just need to scale it back to the dimensional space.
 
 \subsection{Results}
 
-
-We should see excellent agreement between the analytical and numerical solutions.
+We consider $\epsilon=0.2$. The following three figures show the stage, $x$-momentum, and $y$-momentum at several instants in time. We should see excellent agreement between the analytical and numerical solutions.
 
 \begin{figure}[h]