Changeset 7584

Timestamp:

Dec 9, 2009, 8:37:58 AM (15 years ago)

Author:

jgriffin

Message:

Added energy run-up calculations and breaking criteria to Madsen run-up

Location:

misc/tools/event_selection/wave_energy

Files:

• run_up.py (modified) (2 diffs)
• wave_energy.doc (modified) (previous)
• wave_energy.py (modified) (8 diffs)

misc/tools/event_selection/wave_energy/run_up.py

@@ -1 +1 @@
 """ This code implements run up solutions to the shallow water wave
 equations.
+
+References:
+Madsen and Fuhrman 2008. Run-up of tsunamis and long waves
+in terms of surf-similarity. Coastal Engineering, 55, p209.
 
 Creator: Jonathan Griffin, Geoscience Australia
@@ -8 +12 @@
 def Madsen(depth, slope, height, period, gravity = 9.81):
 
-    freq = (2*np.pi)/period
-    a = (depth*np.power(freq,2.0))/(gravity*np.power(slope,2.0))
-    print 'a',a
-    b = np.power(a,0.25)
-    print 'b',b
-    run_up =  2*np.power(np.pi,0.5)*b*height
+    freq = (2*pi)/period
+    a = (depth*power(freq,2.0))/(gravity*power(slope,2.0))
+    #print 'a',a
+    b = power(a,0.25)
+    #print 'b',b
+    run_up_madsen =  2*power(np.pi,0.5)*b*height
 
+    R_break = gravity*(power(slope,2.0))/(power(freq,2.0))
+    run_up = np.min(run_up_madsen, R_break)
     return run_up
 
-def surf_sim(depth, slope, height, wavelength, gravity = 9.81):
+#def surf_sim(depth, slope, height, wavelength, gravity = 9.81):
+# Add Madsen's surf similarity here
 
-    
+
+
+def energy_conserv(crest_integrated_energy,theta = 3.,alpha = 2.,gravity = 9.81):
+    E = crest_integrated_energy
+    print 'crest_integrated_energy', crest_integrated_energy
+    a = 1/tan(alpha*pi/180.)
+    b = 1/tan(theta*pi/180.)                                
+ 
+    F = a-b
+#    print 'F',F
+    R3 = 6.*E/(1000.*gravity*F)
+#    print 'R3', R3
+    R = power(R3,(1./3.))
+
+    return R
 
 import numpy as np
+from numpy import tan, sin, cos, power,pi
 if __name__ == '__main__':
     sys.exit(main())

misc/tools/event_selection/wave_energy/wave_energy.py

@@ -25 +25 @@
 index = 2872
 filebasename = 'sts_gauge_' + str(index) +'.csv'
-filename = join(filepath, filebasename)
+#filename = join(filepath, filebasename)
 
 filepathlist = []
@@ -44 +44 @@
 max_energy_inst_list = []
 run_up_list = []
+energy_run_up_list = []
+max_stage_list = []
 counter = 0
 
@@ -89 +91 @@
         x_vel = x_momentum[i]/(100+stage[i])
         y_vel = y_momentum[i]/(100+stage[i])
-        velocity = numpy.sqrt(x_vel*x_vel,y_vel*y_vel)
+        velocity = numpy.sqrt(x_vel*x_vel+y_vel*y_vel)
         # Kinetic energy KE = 1/2 roh h u^2
         KE = (1./2.)*(numpy.power(x_vel,2)+numpy.power(y_vel,2))*(1000)*(100+stage[i])
@@ -110 +112 @@
             if temp_max_energy > max_energy:
                 max_energy = temp_max_energy
+                max_energy_time = time[i] - start_time
             if temp_max_stage > max_stage:
                 max_stage = temp_max_stage
@@ -122 +125 @@
             if temp_max_energy > max_energy:
                 max_energy = temp_max_energy
+                max_energy_time = time[i] - start_time
             if temp_max_stage > max_stage:
                 max_stage = temp_max_stage
@@ -142 +146 @@
     # Add maximum values to list        
     max_energy_inst_list.append(max_energy_inst)
-    max_energy_list.append( max_energy)
+    max_energy_list.append(max_energy)
+
+    # get wave period from max time
+##    locator = time.index(max_time)
+##    time_count = 0
+##    for i in range(locator,len(time)):
+        
 
     # Account for both halves of the wave
@@ -151 +161 @@
     # call run-up module
     #####################################################################
-    Run_up = run_up.Madsen(100,0.01,max_stage,max_time)
-    print 'runup', Run_up
+    Run_up = run_up.Madsen(100,0.017,max_stage,max_time)
+    print 'Madsen runup', Run_up
     run_up_list.append(Run_up)
+    
+    energy_run_up  = run_up.energy_conserv(max_energy,1.0,0.8,9.81)
+    energy_run_up_list.append(energy_run_up)
+    print 'energy run up', energy_run_up
+    max_stage_list.append(max_stage)
 
     
@@ -170 +185 @@
 
 # Print results for each file and calculate maximum value across all events
+outFilename= filepath +'/wave_energy_summary.csv'
+outFile = file(outFilename,'w')
+outFile.write('filename, max crest-integrated energy (J.s/m^2), max instantatneous energy J/m^2, Run up (Madsen), Energy run up , max stage (m)\n')
 for filename in filepathlist:
+    outFile.write(filename+',')
+    outFile.write(str(max_energy_list[counter])+',')
+    outFile.write(str(max_energy_inst_list[counter])+',')
+    outFile.write(str(run_up_list[counter])+',')
+    outFile.write(str(energy_run_up_list[counter])+',')
+    outFile.write(str(max_stage_list[counter])+'\n')
     print filename
-    print 'max crest-integrated energy:', max_energy_list[counter],'J.s/m^2'
-    print 'max instantatneous energy:', max_energy_inst_list[counter],'J/m^2'
-    print 'Run_up:', run_up_list[counter], 'm'
+    print 'max crest-integrated energy:', max_energy_list[counter],'J/m^2'
+    print 'max instantatneous energy:', max_energy_inst_list[counter],'J/s*m^2'
+    print 'Run_up (Madsen 2009):', run_up_list[counter], 'm'
+    print 'Run_up (Energy 2009):', energy_run_up_list[counter], 'm'
     
     if max_energy_list[counter] > total_max_energy: