Cliffs. Elena Tolkova NorthWest Research Associates, Redmond, WA. NTHMP Benchmarking workshop: Tsunami currents Portland, OR Feb 9-10, 2015

Transcription

1 Cliffs Elena Tolkova NorthWest Research Associates, Redond, WA NTHMP Bencharking workshop: Tsunai currents Portland, OR Feb 9-, 25!

2 Shallow-Water Equations (atrix notation) Solution ethod: Titov and Synolakis, 995, 998 2D becoes a sequence of Ds Transition to Rieann Invariants 3 single convection probles

3 Difference Schee Taylor expansion: Maintain steady state: Possible difference schees: () j- j j+ centered difference or average between left and right cells (2)

4 Difference Schee Taylor expansion: Maintain steady state: VTCS (Titov & Synolakis 995): MOST, Cliffs Cliffs-2

5 MOST: sub-critical flow iplied wall in the edge wet node Cliffs: no iplications, irror ghost node wall in-between wet and dry nodes Wall in MOST, when coputing: in x-dir (rows) and in y-dir (coluns) Wall in Cliffs, any dir

6 E ( 5 J) Cliffs MOST 22 h decay ECliffs=3% EMOST Energy of 2 Tohoku tsunai Pacific-wide vs tie after EQ, log scale Surface elevation at Monterey Bay gage, local odel forced by the ocean-wide run: gage record MOST Cliffs hrs after quake MOST kept the correct signal level for 3 h, Cliffs - for twice as long c c hour after quake

7 Fro vertical wall to sloping beach flooding depth hin u Cliffs Birth Certificate: E. Tolkova. Land-Water Boundary Treatent for a Tsunai Model With Diensional Splitting. Pure and Applied Geophysics, Vol. 7, Issue 9 (24), pp Cliffs Code: GitHub, CSDMS,

8 Benchark Proble # Flow over suberged cone!

9 n=; dx=2.5 c (394 x 6); dt=.2 s aniation confr.ov

10 n=.; dx=2.5 c (394 x 6); dt=.3 s aniation confr.ov

11 n=.5; dx=2.5 c (394 x 6); dt=.3 s aniation con_fr5.ov

12 n= u (/s) v (/s) u2 (/s) v2 (/s) sec

13 n=. Cliffs, data v2 (/s) u2 (/s) v (/s) u (/s) sec

14 n=.5 Cliffs, data v2 (/s) u2 (/s) v (/s) u (/s) sec

15 Benchark Proble # 2 Paraeters lon x lat (cells) dx () dt (s) n 4 x x x

16 .5 Input Control Point c in Control Point Input hr

17 Control Point Hilo gage Tie after EQ (hrs)

18 Control Point Cliffs Prescribed Hilo gage 2 Cliffs data Tie after EQ (hrs)

19 .3 HA25, W E current.2. /s..2 Cliffs Cliffs 2 data HA25, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) 2 resolution

20 .3 HA25, W E current.2. /s..2 Cliffs Cliffs 2 data HA25, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) resolution

21 .3 HA25, W E current.2. /s..2 Cliffs Cliffs 2 data HA25, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) 5 resolution

22 HA26, W E current.5 /s.5 Cliffs Cliffs 2 data HA26, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) 2 resolution

23 HA26, W E current Cliffs Cliffs 2 data.5 /s HA26, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) resolution

24 HA26, W E current Cliffs Cliffs 2 data.5 /s HA26, S N current Cliffs Cliffs 2 data.5 /s Tie after EQ (hrs) 5 resolution

25 Maxial Current (/s) at resolution: 2 5 Cliffs Cliffs-2

26 Benchark Proble # 3 A Beacon Motoriki ADCP Tug B Sulphur P 837 x 747, dx=3, dt=.2 s, n=.25

27 .4 Input coputed by deconvolution.2.2 Response at Control Point Control P. Input hr Input Control Point.5 c in hr Input o Response =!.6.4 Control P. Input ControlPoint hr

28 .5 ABeacon: Tsunai only ABeacon: Tsunai+tide ABeacon: Tide only hour

29 .5 ABeacon: Tsunai only Cliffs data ABeacon: Tsunai+tide ABeacon: Tide only hour

30 .5 TugBerth: Tsunai only TugBerth: Tsunai+tide TugBerth: Tide only hour

31 .5 TugBerth: Tsunai only.5 Cliffs data TugBerth: Tsunai+tide TugBerth: Tide only hour

32 .5 SulphyrP: Tsunai only SulphyrP: Tsunai+tide SulphyrP: Tide only hour

33 .5 SulphyrP: Tsunai only Cliffs data SulphyrP: Tsunai+tide SulphyrP: Tide only hour

34 .5 Motoriki: Tsunai only Motoriki: Tsunai+tide Motoriki: Tide only hour

35 .5 Motoriki: Tsunai only Cliffs data Motoriki: Tsunai+tide Motoriki: Tide only hour

36 .5 ADCP: Tsunai only /s ADCP: Tsunai+tide /s ADCP: Tide only /s hour

37 /s.5.5 ADCP: Tsunai only Cliffs data ADCP: Tsunai+tide /s ADCP: Tide only /s hour

38 ADCPexact: Tsunai only /s ADCPexact: Tsunai+tide /s ADCPexact: Tide only /s hour

39 /s.5.5 ADCPexact: Tsunai only Cliffs data ADCPexact: Tsunai+tide /s ADCPexact: Tide only /s hour

40 Non-linearity? X velocity 2 /s 2 total su Y velocity total su /s hour total=tsunai&tide, su=(tsunai only) + (tide only)

41 Maxial Current

42 Benchark Proble # 4 Input: WG (5 ) thru left boundary WG3 B B4 B6 B9 Paraeters: x: ; dx: (depth following) y: ; dy:.3 dt =.5 s; n =.5

43 Benchark Proble # 4 aniation seaside.ov

44 Surface Elevation at WG 3.2 WG3, data Cliffs s

45 Water Surface Elevation: Cliffs, data.2 B.2 B B6.2 B s s

46 U-velocity: Cliffs, data B B4 /s /s B s B s

47 Moentu flux: Cliffs, data.8 B.8.6 B4 3 / c / c B s B s

48 No houses, Surface Elevation: Cliffs, data.2 B.2 B B6.2 B s s

49 No houses, U-velocity: Cliffs, data B B4 /s /s B s B s

50 No houses, Moentu Flux: Cliffs, data.8 B.8.6 B4 3 / c / c B B s s

51 Benchark Proble # 5 aniation conisland.ov

52 ().4.2 X=7.5, Y= X=7.5, Y=5. ().4.2 X=3., Y= X=3., Y=5. ().4.2 X=2., Y= X=2., Y=5. () ().4.2 X=25, Y=.4.2 X=25, Y= Tie (sec) X=25, Y= Tie (sec) WGs: Cliffs data

53 U&V velocities: Cliffs data /s 2 U.5 X=3, Y= V /s /s U U X=2, Y= X=2, Y= V V 2 Tie (sec).5 2 Tie (sec)

54 aniation Tohoku2AK.ov ;

55 Thank you!