Changeset 9216

Timestamp:

Jun 24, 2014, 10:48:53 AM (10 years ago)

Author:

steve

Message:

Moved cross sectional calculations to plot_results

Location:

trunk/anuga_core/validation_tests/analytical_exact/trapezoidal_channel

Files:

• numerical_channel_floodplain.py (modified) (2 diffs)
• plot_results.py (modified) (5 diffs)
• produce_results.py (modified) (1 diff)
• project.py (added)
• results.tex (modified) (1 diff)

trunk/anuga_core/validation_tests/analytical_exact/trapezoidal_channel/numerical_channel_floodplain.py

@@ -20 +20 @@
 alg = args.alg
 
-
-floodplain_length = 800.0 # Model domain length
-floodplain_width = 14.0    # Model domain width
-floodplain_slope = 1./300.
-chan_initial_depth = 0.65  # Initial depth of water in the channel
-chan_bankfull_depth = 1.0  # Bankfull depth of the channel
-chan_width = 10.0          # Bankfull width of the channel
-bankwidth = 2.             # Width of the bank regions -- note that these protrude into the channel
-man_n=0.03                 # Manning's n
-l0 = 1.00 # Length scale associated with triangle side length in channel (min_triangle area = 0.5*l0^2)
+from project import *
+
+# floodplain_length = 800.0 # Model domain length
+# floodplain_width = 14.0    # Model domain width
+# floodplain_slope = 1./300.
+# chan_initial_depth = 0.65  # Initial depth of water in the channel
+# chan_bankfull_depth = 1.0  # Bankfull depth of the channel
+# chan_width = 10.0          # Bankfull width of the channel
+# bankwidth = 2.             # Width of the bank regions -- note that these protrude into the channel
+# man_n=0.03                 # Manning's n
+# l0 = 1.00 # Length scale associated with triangle side length in channel (min_triangle area = 0.5*l0^2)
 
 assert chan_width < floodplain_width, \
@@ -181 +182 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-for t in domain.evolve(yieldstep=10.0, finaltime=500.0):
+for t in domain.evolve(yieldstep=10.0, finaltime=1000.0):
     if myid == 0 and verbose: print domain.timestepping_statistics()
 

trunk/anuga_core/validation_tests/analytical_exact/trapezoidal_channel/plot_results.py

@@ -3 +3 @@
 import numpy
 
+from project import *
+
+filename = 'channel_floodplain.sww'
 
 # Time-index to plot outputs from
-p2 = util.get_output('channel_floodplain.sww')
+p2 = util.get_output(filename)
 p=util.get_centroids(p2, velocity_extrapolation=True)
 v = (p.x>6.0)*(p.x<8.0)
@@ -88 +91 @@
 pyplot.plot(p.x[v1],Analytic_Stage,'o', label='analytical stage')
 pyplot.plot(p.x[v1],p.elev[v1],'o', label='bed elevation')
-pyplot.legend(loc=10)
+pyplot.ylim([-4,2])
+pyplot.legend(loc=8)
 pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Downstream channel regions (95 to 105m)' +'\n')
 pyplot.xlabel('Cross-channel distance (m)')
@@ -117 +121 @@
 pyplot.plot(p.x[v1],Analytic_Stage,'o', label='analytical stage')
 pyplot.plot(p.x[v1],p.elev[v1],'o', label='bed elevation')
+pyplot.ylim([-4,2])
 pyplot.legend(loc=10)
 pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Central channel regions (495 to 505m)' +'\n')
@@ -147 +152 @@
 pyplot.plot(p.x[v1],Analytic_Stage,'o', label='analytical stage')
 pyplot.plot(p.x[v1],p.elev[v1],'o', label='bed elevation')
+pyplot.ylim([-4,2])
 pyplot.legend(loc=10)
 pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Downstream channel regions (695 to 705m)' +'\n')
@@ -152 +158 @@
 pyplot.ylabel('Generic scale (m or m/s)')
 pyplot.savefig('fig4downstream_channel.png')
+
+
+
+
+
+print '#======================================================================'
+print '# Extract some cross section info'
+print '#======================================================================'
+
+from anuga.shallow_water.sww_interrogate import get_flow_through_cross_section
+
+polyline0 = [ [floodplain_width, 10.0], [0., 10.0]]
+polyline1 = [[floodplain_width, floodplain_length-300.0], [0., floodplain_length-300.0]]
+polyline2 = [[floodplain_width, floodplain_length-1.0], [0., floodplain_length-1.0]]
+        
+time,Q0  = get_flow_through_cross_section(filename, polyline0, verbose=True)
+time,Q1  = get_flow_through_cross_section(filename, polyline1, verbose=True)
+time,Q2  = get_flow_through_cross_section(filename, polyline2, verbose=True)
+
+pyplot.figure(figsize=(12.,8.))
+pyplot.plot(time, Q0, label='10m')
+pyplot.plot(time, Q1, label='500m')
+pyplot.plot(time, Q2, label='799m')
+pyplot.plot([0,time[-1]], [Qin,Qin], label='Input Q')
+pyplot.ylim([0,7])
+pyplot.legend(loc=10)
+pyplot.title('Cross sectional flow across transect at 10m, 500m and 799m')
+pyplot.xlabel('Time (sec)')
+pyplot.ylabel('Discharge (m^3/sec)')
+pyplot.savefig('cross_section_10_500_790.png')
+
+
+
+

trunk/anuga_core/validation_tests/analytical_exact/trapezoidal_channel/produce_results.py

@@ -1 @@
-#--------------------------------
-# import modules
-#--------------------------------
-from anuga.validation_utilities.fabricate import *
-from anuga.validation_utilities import run_validation_script
-from anuga.validation_utilities import typeset_report
 
-# Setup the python scripts which produce the output for this
-# validation test
-def build():
-    run_validation_script('numerical_channel_floodplain.py')
-    run_validation_script('plot_results.py')
-    typeset_report()
 
-def clean():
-    autoclean()
+import anuga
+from anuga.validation_utilities import produce_report
 
-main()
+args = anuga.get_args()
+
+produce_report('numerical_channel_floodplain.py', args=args)
+
+

trunk/anuga_core/validation_tests/analytical_exact/trapezoidal_channel/results.tex

@@ -27 +27 @@
 \end{figure}
 
-Table~\ref{tab:trapztab} shows the discharge computed at a number of cross-sections in the channel, at a number of time-steps on the way to near steady-state. By the end of the simulation they should all be essentially the same. Large variations may suggest mass conservation errors (small variations are probably due to the interpolation that occurs in the routine:
-\begin{equation*}
-\textrm{compute\_flow\_through\_cross\_section}.
-\end{equation*}
+Figure~\ref{fig:hydrographs} show the hydrographs through various cross sections showing the flows limiting to the 
+expected inflow $Q$. 
+
+\begin{figure}
+\begin{center}
+\includegraphics[width=\textwidth]{cross_section_10_500_790.png}
+\caption{Hydrographs over a number of cross-sections.}
+\label{fig:hydrographs}
+\end{center}
+\end{figure}
+
+%Table~\ref{tab:trapztab} shows the discharge computed at a number of cross-sections in the channel, at a number of time-steps on the way to near steady-state. By the end of the simulation they should all be essentially the same. Large variations may suggest mass conservation errors (small variations are probably due to the interpolation that occurs in the routine:
+%\begin{equation*}
+%\textrm{compute\_flow\_through\_cross\_section}.
+%\end{equation*}
 
 
-\DTLloaddb{dischargeout}{../analytical_exact/trapezoidal_channel/discharge_outputs.txt}
+%\DTLloaddb{dischargeout}{../analytical_exact/trapezoidal_channel/discharge_outputs.txt}
 %\DTLloaddb{dischargeout}{discharge_outputs.txt}
-\begin{table}
-\caption{Discharge through cross-sections at a number of $x$-position, at different instants in time}
-\label{tab:trapztab}
-\DTLdisplaydb{dischargeout}
-\end{table}
+%\begin{table}
+%\caption{Discharge through cross-sections at a number of $x$-position, at different instants in time}
+%\label{tab:trapztab}
+%\DTLdisplaydb{dischargeout}
+%\end{table}
 
 \endinput