Changeset 6735
- Timestamp:
- Apr 7, 2009, 8:53:39 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/publications/boxing_day_validation_2008/patong_validation.tex
r6734 r6735 95 95 96 96 \subsection{Validation data} 97 Eyewitness accounts detailed in~\cite{papadopoulos06} report that most people at Patong Beach observed an initial retreat of the shoreline of more than 100m followed a few minutes later by a strong wave (crest). Another less powerful wave arrived another five or ten minutes later. Eyewitness statments place the arrival time of the strong wave between 2 hours and 55 inutes to 3 hours and 5 minutes after the source rupture (09:55am to 10:05am local time). After the event (HOw long?) a survey mapped the maximum observed inundation at Patong beach. The inundation map is shown in Figure~\ref{fig:patongescapemap} and was obatined fromthe Thai Department of Mineral Resources \protect \cite{XXX}.97 Eyewitness accounts detailed in~\cite{papadopoulos06} report that most people at Patong Beach observed an initial retreat of the shoreline of more than 100m followed a few minutes later by a strong wave (crest). Another less powerful wave arrived another five or ten minutes later. Eyewitness statments place the arrival time of the strong wave between 2 hours and 55 inutes to 3 hours and 5 minutes after the source rupture (09:55am to 10:05am local time). After the event (HOw long?) a survey mapped the maximum observed inundation at Patong beach. The inundation map is shown in Figure~\ref{fig:patongescapemap} and was kindly provided by the Thai Department of Mineral Resources \protect \cite{XXX}. 98 98 99 99 \begin{figure}[ht] … … 143 143 144 144 The domain was discretised into approximately ...,000 triangles. The resolution of the grid was increased in certain regions to efficiently increase the accuracy of the simulation. The grid resolution ranged between a maximum triangle area of $...\times 10^5$ m$^2$ near the Western ocean boundary to $...$ m$^2$ in the small regions surrounding the inundation region in Patong Bay. Due to a lack of available data, friction was set to a constant througout the computational domain. A Manning's coefficient of 0.01 was chosen based upon previous numerical experiments conducted by the authors (FIXME: Citation Tom Baldock?? Or Duncan??). 145 In \cite{ Schoettle2007} values of Manning's coefficient in the range 0.007 to 0.030 is suggested for tsunami propagation over a sandy sea floor.145 In \cite{schoettle2007} values of Manning's coefficient in the range 0.007 to 0.030 is suggested for tsunami propagation over a sandy sea floor. 146 146 147 147 The boundary condition at each side of the domain towards the south and the north where no data was available was chosen as a transmissive boundary condition effectively replicating the time dependent wave height present just inside the computational domain. Momentum was set to zero. Other choices include applying the mean tide value as a Dirichlet type boundary condition but experiments as well as the result of the verification reported here showed that this approach tends to under estimate the tsunami impact due to the tempering of the wave near the side boundaries. FIXME(OLE): Should we include Nick's test example? … … 349 349 Geoscience Australia, in an open collaboration with the Mathematical Sciences Institute, The Australian National University, is developing a software application, \textsc{anuga}, to model the hydrodynamics of tsunamis, floods and storm surges. The open source software implements a finite volume central-upwind Godunov method to solve the non-linear depth-averaged shallow water wave equations. This paper investigates the veracity of \textsc{anuga} when used to model tsunami inundation. A particular aim was to make use of the comparatively large amount of observed data corresponding to the Indian ocean tsunmai event of December 2004, to provide a conditional assessment of the computational model's performance. Specifically a comparison is made between an inundation map, constructed from observed data, against modelled maximum inundation. This comparison shows that there is very good agreement between the simulated and observed values. The sensitivity of model results to the resolution of bathymetry data used in the model was also investigated. It was found that the performance of the model could be drastically improved by using finer bathymetric data which better captures local topographic features. The effects of two different source models was also explored. 350 350 351 different even types submarine mass failure generate larger events because of proximity more directional wave generation 352 353 even if data is available it is hard to access 351 different even types submarine mass failure generate larger events because of proximity more directional wave generation even if data is available it is hard to access 354 352 355 353 \begin{thebibliography}{7} … … 392 390 \bibitem{satake95} 393 391 Satake, K. (1995). Linear and nonlinear computations of the 1992 Nicaragua earthquake tsunami. Pure and Applied Geophysics, {\bf 144(3)}, 455-470. 394 395 396 \bibitem{Schoettle2007} 392 \bibitem{schoettle2007} 397 393 Schoettle E., and Sakimoto S. (2007), Modeling the Effects of Coral Reef Health on Tsunami Run-up 398 394 with the Finite-element Model ADCIRC, \url{http://istim.ce.nd.edu/2007/Posters/Schoettle_poster.pdf}. 399 400 395 \bibitem{synolakis05} 401 396 Synolakis, C., E. Okal, and E. Bernard (2005), The megatsunami of December 26 2004, {\em The Bridge, National Academy of Engineering Publications}, {\bf 35(2)}, 36-35.
Note: See TracChangeset
for help on using the changeset viewer.