Effect of sea-dykes on tsunami run-up Kao-Shu Hwang Chen-Chi Liu Hwung-Hweng Hwung Tainan Hydraulics Laboratory National Cheng-Kung University Tainan TAIWAN.
Outline Introduction tsunami research in THL after the 24 Indian Ocean tsunami The motivation and the objectives of present study Preliminary results Ongoing experiments Future works 2
run up sensor 3.5 m 1.5 m wave generator 1:6 197 m 6 1 MEASURE 2 6 PROCESS 3 wave gages No. 51~78 spacing 2. m 88 m plane view of the flume wave gages No. 79~96 spacing.4~.8 m GENERATE 4 5 AD/DA 5. m wave gages No. 4~47 spacing 4.~6. m 15m 6 1:2 side view of the flume.4 m 1:4 5.2 m 1 2 3 4 5 6 : : : : : : : configuration file instruments database measured signals processing results AD/DA instruction file wave board control signals command file IPC1 MNDAS Experimental set-up 3
4
Surface wave profiles after breaking.6.4.2-2 -1 5-1 -5 5 1 15 -.2 -.4 -.6 5
R/h vs. H/h.5 c o tβ = 1 9 r u n u p la c o tβ = 2 E B M c o tβ = 1 9 c o tβ = 2.4.3. 8 5, S y n o la k is ( 1 9 8 7 ) w, 4, 6, L i( 2 3 ). 8 5, S y n o la k is ( 1 9 8 7 ), h = 3. 5 m, r a m p tr a je c to r y R /h c o tβ = 4, h = 2. 5 m, r a m p tr a je c to r y c o tβ = 6, h = 2. 5 m, r a m p tr a je c to r y.2 c s c s c s.1.1 H /h.2 o o o o o o tβ = 6, lita r y -w tβ = 6, lita r y -w tβ = 6, lita r y -w h av h av h av = e = e = e 2.5 tr a 3.2 tr a 1.5 tr a m, je c to r y m, je c to r y m, je c to r y.3 6
.2.4 x*=-.7.1.2 -.1 -.2 x*=-1..1.2 -.1 -.2 x*=-1.4.1.2 breaking at x*=-1. -.1 -.2 x*=-1.8.1.2 η* η* and u* (cotβ=2, H/h=.12) -.1 -.2 x*=-2.8.1.2 u* -.1 -.2 2 3 t* 4 5 7
.6 x*=3.8.3..4.9.6..3 -.2. -.4.9.2 x*=1.8.6..3 -.2. -.4.9.2 x*=.9.6 η*.2 x*=2.8..3 -.2. -.4 x*=..9.2.6..3 -.2. -.4 2 3 4 5 6 η and u* above shoreline (cotβ=2, H/h=.12) u* 7 t* 8
.2.4 x*=-1.8.1 -.1 -.2 x*=-2.7.1.2 -.1 -.2 x*=-4.5.1 η*.2 -.1 -.2 x*=-6.3.2 -.1 -.2 4 5 6 t* 7 breaking at x*=-6.3.2.1 η and u* (cotβ=4, H/h=.12) u* 8 9
.3.6 x*=-1.8.2.1.4.2 -.1 -.2 x*=-6.3.2.4.1.2 -.2 -.1 x*=-8.3.2.4.1.2 -.1.4.2.1 η* -.1 breaking at x*=-11.8 -.2 x*=-1..2 η and u* (cotβ=4, H/h=.23) -.2 x*=-11.8.2.4.1 -.1 -.2 u*.2 2 4 6 t* 8 1 1
.1.4.8 x*=7.3.6.2..4.2 -.2.1.4.8 x*=5.4.6.2..4.2 -.2.1.4 x*=3.6.8.6.2..4.2 -.2.1.4.8 x*=1.8.6.2..4 η* η and u* above shoreline (cotβ=4, H/h=.23).2 -.2.1.4.8 x*=..6.2..4.2 -.2. -.4 4 6 8 1 u* 12 t* 11
About 1/3 coastline is protected coast length 1,141 km with 369.8 km protected by sea-dykes or revetments design wave conditions:??????????????????????????? West coast??? H1/3 :4-6 m, T1/3:8-1 s East coast H1/3 :1-14 m, T1/3:1-16 s crown levels of sea dykes or????????????????????????????????? revetments West coast???????????? E.L. +4~ +6 m East coast E.L. +8~+1 m 12
A typical coastal landscape of Taiwan Seafloor bathymetry source: NCOR, TAIWAN
A report about the role of sea-dykes (revetments) in the 24 tsunami attack.. However, on Male Island, which is the national capital of the Republic of Maldives, revetments and other structures have been constructed around the island. The island was relatively unaffected by the Sumatra tsunami because these revetments helped protect the island from inundation. National Institute for Land and Infrastructure Management, Ministry of Land, Infrastructure and Transport, Japan 14
State of tsunami wave runup with and without revetment National Institute for Land and Infrastructure Management, Ministry of Land, Infrastructure and Transport, Japan 15
The objectives Tsunami run-up and overflow when there are sea-dykes (revetments) protecting the coast Hydrodynamics around the sea dykes (revetments) and the overflow field in a tsunami attack 16
Preliminary results by taking advantage of consulting projects that originally aimed to investigate the stability of sea-dykes under storm wave conditions
s o lita r y w a v e -fo r m x / h = -1 6 7 η /η m ax 1..5. -6. -4. -2.. θ 2. 4. 6. The flume 2 m 2 m 2 m The wave maker Servo-electrical piston type Stroke: 1. m Max..4 m solitary wave height under 1. m water depth 18
A typical sea-dyke profile of the east coast 19
The mean overflow rate 25 m ean overtopping rate (m ^3/s/m ) 2 15 1 Y am am oto (25) 5 w ith detached breakw ater w ithout detached breakw ater 1 2 3 4 5 6 over flow height (m ) 2
The mean overflow rate h H q =.55 gh h (m 3 /m/s) Where h is the overflow height by Takanashi and Yamamoto (25) 21
Computation domain 3 dx =.5 m B oundary particle = 648 Fluid particles = 887 2 1 On fixed boundary, boundary force is used to prevent penetration. -1 1 X (m) 2 3 22
23
24
25
3 2 1-1 1 15 2 25 3 26
Ongoing experiments
Considerations in test program tsunami heights the tsunami heights (H) are the same order to the sea-dyke height above the sea level (Hc) The sea-dyke front and rear slopes The slope behind the sea-dyke Hc H 28
Disaster cases of coastal structures by tsunamis (Yamamoto et. al.,26) Date of disaster Coastal name Structure type Damage level Type of destructive force Crown height (m) Tsunami height (m) Mean overflow rate (m3/m/s) Hokkaido Prefectur e (1993 July 12) Hirahama, Taisei-cho Revetment (crown covered with concrete) Complete destructio n 6. 8. 4.87 Umikurimae, Okushiricho Self-supported dike Complete destructio n 5.3 8.7 1.79 Monai, Okushiricho Self-supported revetment Complete destructio n 6.5 21. 95.7 Aonae, Okushiricho Self-supported revetment Complete destructio n By incident tsunami pressure By incident tsunami pressure By incident tsunami pressure By tsunami pressure from back coast 4.5 12.4 38.23 29
Date of disaster Coastal name Structure type Damage level Type of destructive force Crown height (m) Tsunami height (m) Mean overflow rate (m3/m/s) Akita Prefectur e facing the Nihon sea (1983 May 26) Kodomari fishery harbor Block-type revetment (crown not covered) Complete destruction By return flow 4. 5.3 2.55 Todoroki fishery harbor Self-supported revetment Partial destruction By return flow 3.7 4..28 Iwadate fishery harbor Self-supported revetment Partial destruction By incident tsunami pressure 4.65 7. 6.12 Iwadate fishery harbor Self-supported revetment Partial destruction By incident tsunami pressure 5.65 8.8 9.81 Minehama Self-supported revetment Partial destruction By incident tsunami pressure 5.3 6.9 3.52 Kotohama Reinforced concrete revetment (the crown was covered) Complete destruction Destruction of front face and outflow of back filling 5.25 8.5 1.9 3
2m NDV.5m.76m.5m 1 :2 1 m.1 8m Piston-type W avem aker Synchronizer A D C ard Am plifier A D system W avem aker system A solitary wave propagate on an inclined beach & a sea-dyke The run-up heights, mean overflow rate The pressure distribution on the structure The flow field 31
.8 d2 d18 d16 d2 Dike d18 Dike d16 Dike Synolakis (1987) R/h.6.4.2.. The run-up heights.1.2 H/h.3.4.5 33
3 2.5 2 1.5 1.5 3 2.5 2 1.5 1.5 d2h71t25 (cm) H (cm) P (cm/s) 12 1 8 6 4 2 6 V 7 8 9 1 Time(1/5 s) 11 12 Time series of the overflow depth, pressure and velocity on top of the sea-dyke 34
M e a n o v e r f l o w r a te ( m ^ 3 /m /s ). 8 T a k a n a s h i a n d Y a m a m o to ( 2 5 ) d2 d18 d16. 6. 4. 2... 1. 2 T s u n a m i o v e r f lo w h e ig h t ( m ). 3. 4 The overflow 35
Future works To use the super tank for larger scale tests on 36
Thank you! http://www.thl.ncku.edu.tw
38
Seafloor bathymetry around TAIWAN (source: NCOR, TAIWAN) 39
The worries Can it happen here? How bad can it be if it can happen here? How about the defense capability of the existing coastal defense works? 4
h 堤重 weight of the passive ground sea wall pressure 被動土壓 H static water pressure incident tsunami pressure q =.55 gh h (m 3 /m/s) Where h is the overflow height by Takanashi and Yamamoto (25) u = 1.1 gh Where H is the tsunami height by Iizuka and Matsutomi (2) 41
h H 42
.8 d2 d18 d16 d2 Dike d18 Dike d16 Dike R/h.6.4.2...1.2 H/h.3.4.5
South China.6 Sea Tsunami workshop 27 d2 d18 The run-up height d16.55 d 2 D ik e d 1 8 D ik e d 1 6 D ik e R /h.5.45.4.35.3.2.25.3.35 H /h.4.45 44
青草崙海堤改善斷面示意圖 青草崙海岸地形示意圖 A typical sea-dyke profile of the west coast 45
Design Wave condition VS. solitary waves sea level: E.L. +2.2 m crown level: E.L. +5. m Design wave H1/3=6. m, T1/3=12. s overflow:.45 m3/s/m Solitary wave H=6. m overflow:.987 (2.41) m3/s/m H=4.8 m.673 (1.237) m3/s/m H=3.6 m.23 (.827) m3/s/m H=2.4 m.173 (.67) m3/s/m tolerance:.5 m3/s/m H=1.2 m.162 (.297) m3/s/m 46
η u (z u (z b = 5 c m b =.5 c m ) ) 3. 2. 1.. -1. -2. -3. 1.. η (m ) x = 4m -1. 1. 3. 2. 1.. -1. -2. -3.. η (m ) x = 3 m -1. 3. 2. 1.. -1. -2. -3. 1. x = m. η(m ) u ( m /s ) u ( m /s ) u ( m /s ) 1. 1 s q r t( g H ) -1. -1 t1 ( s ) 1 2 47