HAWASSI-AB User Manual by © LabMath
Contents
1. ccccccccceccccceeeeeeeesseseeeseeeeeeeeeeeeeeeeeaaas 5 25 5 3 1 2BOOIIrSlope MARIN_103001 Irregular wave above a slope 5 13 5 25 23 9 2 2BO02HarmbBar Harimornic Over Dat estre rris a Neko ESPERE Ee E oa PE voee YR ea ER YES CER V RES 5 25 5 3 3 2BOO3IrBar Irregular Wave over bar eeeeeeesssssssssesssseeeeeeene nennen nennen eene 5 25 34 Breaking waves above DatDymelb asasaran ai a nen aga pa e na KG NINEN GENE Pao ETE ON qva seo KAE o HE PLE Nin 5 26 5 4 1 BEOR mE DEO aaa sasar aa cate ated modius ca 10029 Ka sca a 50100529055 0 Soa Ka i ga aa an aa a aaa aa a 5 26 5 4 2 2BBrOO2IrBar Irregular wave Spilling breaker 5 cccccccssssssessessseseeeeeeeees 5 26 5 5 Run up of waves breaking and non breaking cccccceeccecceceeeeceaaaeseseeeeseseeeceeeeeeeeeeeeaaaas 5 27 5 5 1 3R001 Harm Harmonic Run up non breaking ccccsseeeeeecceeeeeeeeeeaaaseessseseeeeeeeeees 5 27 2 92 3RBrOOlIHarm Spilling Breaker Run up with interior flow 5 5 27 ET GN Ret 6 28 6 1 References to basic papers and applications sss 6 28 O KN urat o TOT 6 29 I 21Page dl HAW ASSI Preamble Waves are fascinating important and challenging The importance can be substantiated from some well known observations e Half of the world population lives less than 150 km from the coast e The sea is a relatively easy medium for transpo2. gt each of the components of the total pressure Be aware that the amount of data can become very large depending on the chosen discretization settings HAWASSI AB1 Internal Flow rec CD anan Working directory er b CA G2 Labmaih Indonesia Internal Flow Calculation Input data Other bi Load data Time interval s t coarse C Save data x interval m x coarse Calculate z interval m dz m Post Processing Input data Other Load data Pressure Velocity Acceleration Dynamic Mon Lin dynamic Total Horizontal Vertical Total C Horizontal Vertical Total O Plotting X Z axes t z axes x t axes x a xes z xes takes 0 Animation O Validation meas step Load measurement data Density plots LL atz Error plot Time mean plot setting xlim zlim clim save data t coarse GIF setting 1 100 inf save figures fig 4 IS Page dl HAW ASSI 4 4 Required lay out of user defined input files User input of data files for various purposes need to be prepared with an extension mat dat txt etc with a specified format as described below 4 4 1 Influx time signal A 2 column matrix time elevation first column the equidistant time s second column the corresponding elevation m 4 4 2 Initial wave profile A 3 column matrix space point x elevation first column the equidistant x value m covering the whole interval if only partially the data wi
3. All combinations are possible of choices for dispersion also user specified nonlinearity and breaking with various choices for each item as described below 2 9 Page dl HAW ASSI 2 3 5 1 Dispersion The main property of HAWASSI AB is that it can handle exact dispersion of small amplitude so called linear waves of any wave length owe to the Fourier character and implementation of AB As a consequence the continuity equation is exactly satisfied above flat bottom and in a very good approximation above varying bottom Besides that mainly for educational academic purposes to enable simulations for models with approximate dispersion other predefined or user defined dispersive models can be chosen gt Shallow Water SW dispersion gt KdV Korteweg de Vries dispersion KdV 1895 gt BBM Benjamin Bona Mahony dispersion BBM 1972 gt User specified dispersion to be specified in input panels For the given or user specified dispersion all nonlinear terms are calculated consistently above flat as well as above varying bottom A brief explanation is given referring to the table for the explicit expressions The exact dispersion relation o Q k is the correct expression for small amplitude waves For kd 0 i e rather long waves or rather shallow water this relation can be Taylor expanded leading to gt SW dispersion first order Taylor Q sw k which gives a translation of waves of any wave leng
4. Flat Depth 510 m Significant wave Height Hs 10 23 m Peak period Tp Po 12 656 8s Peak frequency nu 0 496 rad s Peak wave number kp 0 025 Peak wave length 200087 Lm Peak phase speed t 19 758 m s Steepness kp Hs 2 e 0 129 Relative wave length lambda h 0 49031 kp h bI Category Deep water 5 2 5 1FBr005Bor Undular breaking bore 5 Dynamic Model t HOZDEE Dispersion Model OmExact Bathymetry Flat Depth 1 m INITIAL VALUE PROBLEM HS4brF MTA amp profile time 420 5 1500 x m Irregular wave breaking Simulation of MARIN measurement HSS MTA A pees tines 15 40 1 F f Breaking undular bore at two different times Settings and results as Wei e a for higher amplitudes for which no validation data are available See also Testcase 1 F006Bor section 5 1 6 for the non breaking bore 5 24 Page da HAWASSI 5 3 Non breaking waves above non flat bathymetry 5 3 1 2B001IrSlope MARIN 103001 Irregular wave above a slope 5 13 Dynamic Model HS2B Dispersion Model OmExact Bathymetry Slope Max Depth 0 6 m Min Depth 0 3 m Slope 0 05 m Foot slope position 30 m Significant wave Height Hs 0 062 m Peak period Tp 1 701 s Peak frequency nu 3 694 rad s Peak wave number kp 1 769 Peak wave length 3 551 m Peak phase speed 2 088 m s Steepness kp Hs 2 0 055 Relative wave length lambda h 5 9186 kp h 1 062 Category Intermedia
5. Output ALL output of a simulation will be stored in this subfolder together with selected output of post processing A User Note gives the possibility to provide details or a short description of the specific simulation this text will be copied to the log file Input panels are separated to provide various details of the simulation to be executed Model Initial condition Wave Input Xspace Bathymetry Options there are two major options VVVVV ON o Internal Flow details can be given of the times at which in a post processing step interior flow properties have to be calculated see the Interior Flow GUI o ODE partition to reduce hardware requirements it is possible to split one simulation time interval in various consecutive time intervals the data of each subinterval will be stored at the end of that time interval so that memory requirements and size of data are restricted The option to collect the data in one file after finishing the simulation can be checked Clicking the Pre Proc pre processing button will prepare the input before the actual evolution simulation starts A pop up figure will summarise graphically the input including the geometric lay out and the quality of dispersion used in the computation compared to exact dispersion Warnings suggestions may be given to optimise the results of the simulation the input can be changed after which a new Pre Proc step is required The log file is available after Pre
6. and or to compare simulation results with experimental data or other simulations two formats are available gt Graphical output for time signals at specified positions calculated amplitude or energy spectra spatial values of the Energy of MTA and MWL height of Mean Water Level Hs Skewness Asymmetry and Kurtosis and a position vs time plot of breaking events Quantitative information for time signals at a specified position o thecorrelation of the signals for validation o the value or quotient of the Variance Skewness Asymmetry and Kurtosis of simulated and measured signal The definitions of these quantities 1s as follows for signals with zero mean and with denoting time averaging The correlation of the simulated signal s t and the measured signal m t is corr s m sam s m and for a signal s with H the Hilbert transform H Advar var 9 Sk s s ASA Ku s s 4 17 Page ad HAWASSI 4 3 Internal Flow GUI In order to calculate a posteriori internal flow properties it is needed that before the simulation is started this has been indicated in the Main GUI since some data during the calculation will be stored to be used for the interior flow calculations This storage will slow down the simulation and therefore the times of interest can be indicated Quantities that can be computed are gt the horizontal and vertical velocities and accelerations of the interior fluid motion
7. 41633 7 pages R Kurnia amp E van Groesen Spatial spectral Hamiltonian Boussinesq wave simulations Advances in Computational and Experimental Marine Hydrodynamics VOL 2 Conference Proceedings 2015 pp 19 24 ISBN 978 93 80689 22 7R Kurnia amp E van Groesen Localization for spatial spectral implementations of 1D Analytic Boussinesq equations Wave Motion to be published R Kurnia amp E van Groesen Design of Wave Breaking Experiments and A Posteriori Simulations Memorandum 2042 January 2015 http www math utwente nl publications ISSN 1874 4850 R Kurnia amp E van Groesen Localization in Spatial Spectral Methods for Water Wave Applications Proceedings ICOSAHOM 2014 Submitted W Kristina O Bokhove amp E van Groesen Effective coastal boundary conditions for tsunami wave run up over sloping bathymetry Nonlin Processes GeoPhys 21 2014 987 1005 Lie S Liam D Adyia amp E van Groesen Embedded wave generation for dispersive surface wave models Ocean Engineering 80 2014 73 83 R Kurnia amp E van Groesen High Order Hamiltonian Water Wave Models with Wave Breaking Mechanism Coastal Engineering 93 2014 55 70 Lie She Liam Mathematical modelling of generation and forward propagation of dispersive waves PhD Thesis UTwente 15 May 2013 A L Latifah amp E van Groesen Coherence and Predictability of Extreme Events in Irregular Waves Nonlin Processes Geophys 19 2012 199 213 E van Groesen amp I v
8. 5 437 1995 18 5 crest height on 75m depth kp h 14156 here for simulation in MARIN wave tank cic Intermediate depth simulation of MARIN measurement 5 1 3 IF003Irreg MARIN 224002F Irregular wave in deep water HS4F MTA amp proie time 643 Dynamic Model HS4F Dispersion Model OmExact Bathymetry Flat Depth 510 m Significant wave Height Hs 9o29 m Peak period Tp 13 990715 Peak frequency nu 0 452 rad s Peak wave number kp 0 021 Peak wave length 501 056 mJ Peak phase speed 21 679 m s Steepness kp Hs 2 0 103 Relative wave length lambda h 0 59031 Observe nonbreaking freak wave on deep kp h 10 644 Category Deep water water Simulation of MARIN measurement 5 21IPage ad HAWASSI 5 1 4 IF004Wall AB1 Wall Full wave reflection from wall 5 Dynamic Model gt HSIFWall sec iia pa accu Dispersion Model OmExact Bina aee ine tans MM Bathymetry Flat Depth 5 m LEE EE B GL B B gG g Signal type Harmonic nl Significant wave Height Hs 0 1 m Peak period Tp 2 s Peak frequency nu gt 3 141 rad s Peak wave number kp 1 006 Peak wave length 6 22 1m Peak phase speed 341239 m7 6 Steepness kp Hs 2 0 05 kangh Relative wave length lambda h 1 2494 Pn kp h 5 029 Full wave reflection of a linear harmonic Category Deep water wave from a wall Dynamic Model HS2FWall Dispersion Model OmExact Bathymetry Flat D
9. Proc and contains info about the waves to be simulated note that calculated data may slightly differ from input values because of the statics used to calculate wave length period etc The log file will be updated after finishing the calculation with info about the computation time Successive simulations under the same Project Name will be added successively in the log file but files of computations will be overwritten as warned on the GUI At the bottom of the screen there will be warnings suggestions when specifying input After RUN during the simulation a time indicator estimates the remaining time 4 15 Page 4 2 Post Processing GUI ad HAWASSI After a simulation is finished the Post Processing PP GUI will automatically pop up loaded with the simulation data The GUI can also be called directly from the Main GUI and selected data can be loaded 4 HAWASSI ABI Post Processing File Help Working directory CA Simulation data Other C Plotting Profile s MTA Bathy Energy MTA and MWL Breaking events 2 Animation O Validation Load validation data Amplitude Signals Spectra Energy Setting xlim zlim x COarse t coarse Project name Buoy s Signal Spectrum Quantitative info Best correlation lime shifting Quantitative info MTA time interval Save figures fig Lxhrnath Indparzia Amplitude Energy GIF setting 11100 inf MTA and MWL Hs Sk As K Save
10. an approximative way this was formulated accurately in the 1960 1970 s by Zakharov and Broer by providing the Hamiltonian form of the dynamic equations HAW ASSI software is based on these last findings with methods for making the principal description into a practical numerical modelling and implementation tool The first release of the software deals with wave propagation but the developers are in the process to extend the capabilities to include coupled wave ship interactions amongst others in later releases We sincerely hope that the use of the software just as the design of it has been will be fascinating and challenging for students and academicians as well as for practitioners from both groups we hope to receive comments and suggestions for further improvements and extensions in a way that can be profitable for both sides Let nature tell its secrets Listen to the physics in its mathematical language Restrain from idealization Only then models will serve us in abundance 1 31Page ad HAWASSI O Copyright of HAWASSI software is with LabMath Indonesia an independent research institute under the Foundation Yayasan AB in Bandung Indonesia The software has been developed over the past years in collaboration with the University of Twente Netherlands with additional financial support of Netherlands Technology Foundation STW and Royal Netherlands Academy of Arts and Sciences KNAW By downloading and using the software you agr
11. 1 80 o Weiea gt G Wei J T Kirby S T Grilli R Subramanya A fully nonlinear Boussinesq model for surface waves Part 1 Highly nonlinear unsteady waves Journal of Fluid Mechanics 294 1995 71 92 Only by testing with realistic data the software can be validated and improved 5 20lPage da HAWASSI 5 1 Non breaking waves above flat bottom 5 1 1 1F001Foc MARIN 202002 Strong focussing wave Dynamic Model HSZFE Dispersion Model OmExact Bathymetry Flat Depth 1 m Significant wave Height Hs 0 013 m Peak period Tp 1 304 s Peak frequency nu 4 BlT7 rad s Peak wave number kp 2 404 Peak wave length 24613 mJ Peak phase speed 24 003 m s Jaman makan SENA i i Steepness kp Hs 2 0 015 aid EU REN Relative wave length lambda h 2 6133 Focussing Wave using dispersion to kp h 2 404 generate high waves in wave tanks Category Intermediate depth Simulation of MARIN measurement 5 1 2 IF002Draup MARIN204001 Draupner Wave HEF ee ie ime Dynamic Model HS2F NG pan agr o en ml Dispersion Model OmExact LEE VC RR e s Ie sy vow Bathymetry Flat Depth 1 m Significant wave Height Hs 0 076 m Peak period Tp 2 061 Ss Peak frequency nu 3 048 rad s Peak wave number kp 1 156 a en 70 Peak wave length 2 437 m Environmental Freak Wave measured at Peak phase speed 2 0638 m s e Gas Na IN en JEN gt 544 the Draupner plattorm in the NorthSea Relative wave length lambda h
12. 3Foc7 Focussing Breaking Wave 3 HS3brF OmExact Dynamic Model Dispersion Model Bathymetry Significant wave Height Hs Flat Depth 2 132 m D 078 m Peak Peak Peak Peak Peak period Tp frequency nu wave number wave length phase speed kp 99s 157 rad s 041 038 m Steepness kp Hs 2 Relative wave length lambda h kp h Category HS3brF OmExact Dynamic Model Dispersion Model Bathymetry Significant wave Height Hs Peak period Tp Peak frequency nu Peak wave number kp Peak wave length Peak phase speed Steepness kp Hs 2 i Relative wave length lambda h kp h Category 0 wW ON AF WEF 1 3 1 6 3 034 m s 0 2 2 04 m y 22 1 9 Intermediate depth Plat Depth 2 132 m 356 m 634 s 845 rad s SUN 157 m 544 m s 269 9498 222 Deep water HS FOE pre re BL tl Strong focussing breaking wave Simulation of MARIN measurement HS3brF MTA amp profile time 83 6 T T T x mi Strong focussing breaking wave Simulation of TUD measurement HEXhrE KATA amp profile ima 200 T T A E Fu P 4 FL L7 L s Pu px i i Ti y T ra lia a amy UP SIT y if a PN N H ws i b uy Bichromatic breaking wave 5 23 Page ad HAWASSI 5 2 4 1FBr004Ir MARIN 223002F Irregular wave breaking 10 Dynamic Model gt HS4brF Dispersion Model OmExact Bathymetry
13. 9 1 1 IPage dl HAW ASSI De o c A 5 20 5 Non breaking waves above flat bottom essere eene nnne nns 5 2 5 1 1 1FOO1 Foc MARIN 202002 Strong focussing wave eeeeeeeeeen 5 21 D L2 IF002Draup MARIN204001 Draupner Wave ccccccccccccccceeeessessseesseeececeeeeeeeeeeaaas 5 21 5 1 3 IFOO3Irreg MARIN 224002F Irregular wave in deep water eeeeeesssss 5 21 5 1 4 IF004Wall ABI Wall Full wave reflection from wall 5 ese 5 22 5 1 5 I F005WallFreq Frequency dependent wall reflection 5 sss 5 22 5 1 6 LEOOGBOE Undulit DOG saanak aaa aaa ran Ke paja ea ga daa ei vao pa E AA aa oi ag Ta a gape aa gana d ga 5 22 52 Breaking waves above flat DOMON sccisctccniascincsvaccidtetavatwsaverbanwlansieaavaacsbbntadatravaeebaawlansimaawensoass 5 23 52 1 IFBrOO1Foc MARIN203001 Strong focussing wave eeeeeeee 5 23 9 2 2 IFBr002Foc TUD1403Foc7 Focussing Breaking Wave 3 ss 5 23 3 2 9 IFBr003BiB TUD1403Bi8 Bichromatic wave breaking eeeeeseeees 5 23 5 2 4 1FBr004Ir MARIN 223002F Irregular wave breaking 10 5 24 32 IEBrO05Bor Undular_breaking bore 5 iter trot obo Go uo Rey R Roe io uad suede tod eode 5 24 5 3 Non breaking waves above non flat bathymetry
14. HAWASSI AB User Manual dan AWASS by LabMath Indonesia ver 1 150829 mail address LabMath Indonesia Lawangwangi LMI Jl Dago Giri No 99 Warung Caringin Mekarwangi Bandung 40391 Indonesia e mail hawassi labmath indonesia org home page www hawassi labmath indonesia or VATE i WA Copyright 2015 LabMath Indonesia ad HAWASSI Contents Igisciil eet c E E A M 1 3 SEntrotl Ule DIOE 50 025 E 135020539800 01050 0205395029539 0988000800 0205020802900 009 088 0980 00 8 52088 9980056 1 5 2 Description of HAWASSI AB asarana garane suna aka a gag saves se ab anaa a wka ag a ai aaa naa aa ga a E EIER RE MUR baa 2 0 2 1 Introduction HA W ASSI DB uassedescoisesenztalVeseic nip ddemeta pus esc ns paseuqua nana anana nana nana nana nean 2 6 22 NS and SONI pat ahi onal lee lls Peemenemnemeenewereteertetr tetenmrreett nr DO GEK A aa a NE EKA Cetera et tne Ai KA 2 7 2 9 Model features and tapabilit es sasada wen sesanan saa ian ew dun x Fev tn Pha dU UE aaa naban HU naa nali aa aeei 2 8 2 941 TIS i RN EN TT a a kaa aaa aan ata aing naa a aia a aa aa aa ga a a angina naa a Lang aa ga angka ap sana 2 8 2 3 2 D3tb suci aid FUIS ce saa aana aaa al e da a baa KN ga baa aa aga aa agak aa ag a aba ia 2 8 2 3 3 Embedded wave influxing of various wave LYPES cccccccccccccceeecenaeeseeseseeeeeeeeeeeeeeeeeeaaaas 2 8 2 3 4 TAAT value ot OD STING asah an aaa aa a anaa aaa nga ga anan aa aaa aa aana aan aka SUD P ei Na
15. T T mtu Simulation of Beji amp Battjes Experiment 5 25 Page da HAWASSI 5 4 Breaking waves above bathymetry 5 4 1 2BBr001HarmBar 10 HB bell Bathe amp prolio dit lima 72 1 E xx xr Dynamic Model t HSA2DbrU Dispersion Model OmExact Bathymetry Under water bar Significant wave Height Hs 0 058 m Peak period Tp 25994 8 Peak frequency nu 2 476 rad s Peak wave number kp 1 305 Peak wave length 4 816 m Peak phase speed 1 898 m s Steepnest kp Hs 2 0 038 Relative wave length lambda h 12 0391 kp h 0 522 Breaking harmonic waves over bar Category Intermediate depth Simulation of Beji amp Battjes experiment 5 4 2 2BBr002IrBar Irregular wave Spilling breaker 5 Dynamic Model HSZbrU Hint Beky MITA k profile di ties ERU Dispersion Model OmExact Bathymetry Under water bar i at sea 0 2 Significant wave Height Hs 0 035 m hr ae re Peak period Tp 1 898 s Peak frequency nu 3 3 fads Peak wave number kp 1 806 Peak wave length 3 479 m Peak phase speed 1 832 m s Steepness kp Hs 2 0 032 Relative wave length lambda h 8 6966 RDR Ta Spilling breaker over bar Category Intermediate depth Simulation of Beji amp Battjes experiment 5 26 Page ad HAWASSI 5 5 Run up of waves breaking and non breaking 5 5 1 3R001Harm Harmonic Run up non breaking Dynamic Model HSZBR Dispersion Model OmExact Bathymetry Slope Shore Max Dep
16. aab aga apa aa 2 9 2 52 Model yersions evolution eQUALEODS a aana aka Aan aa da AN aka eges bees edu cac naja a Popes Ie ag daja daana 2 9 2 3 6 Internal TOW caleulattons 5 2er pete cime A KRING NE AG PKK oe BA KENEN WR KE ag ea cue ERES cob a aa a Bagan 2 11 p E oin eel Em 2 11 3 Jnstalhme HAW ASSIEA DB SOIOWObO cscs aces coxa sscduteecennscazesecceesasetsavessusscecamnienacenocteaesacavascadsiGecemmieraceie 3 12 3 1 SUIS TIMES EI S Meni raaa aaa ec ee a ties stereo a anga Baka teers EE 3 12 e TS SIS PE Tk ale NIS aaa aaa aa nga a baan E A Bah AKE Sa Eka 3 12 3 3 Second step Installing HAWASSI AB aaaneeeeeeeeeenan nana anane anana aana nana nnn nnn nnne eaaa 3 12 4 GUP s of HAWASSI AB ccsctdecesastescsoonstacdosnenscdasettecasssieecaeonsteebostacsedasedluegeosteudemonstie bocracsedaseadeeaeoweode 4 13 4 1 EE D EEA anak aaa Ba 4 14 Ad euis ura 6 D 4 16 4 3 wa isi On 6 UI Tr ey 4 18 4 4 Required lay out of user defined input files eeeeeeeesessssssssssseseeeeeeee nnne 4 19 4 4 1 re aa EO EUM 4 19 4 4 2 Initral AVS DEODHG E EE aaa a Aa E aa EE EN 4 19 4 4 3 B TG TV RENNES TN 4 19 4 4 4 External ASU pie BE dataserie baa ag aa anaa aa an ag a E aa a a aaa enak naa EEE 4 19 4 4 5 Interior flow External measurement data cesses eene 4 19 4 4 6 MUS eS ON Me AON daja asa tegen baana aaa aia an ad a a TEE 4 1
17. alls with various reflection properties The acronym HAWASSI stands for Hamiltonian Wave Ship Structure Interaction HAW ASSI AB is a spatial spectral implementation of the Analytic Boussinesq Model AB Presently the code is for simulation of wave structure interactions coupled wave ship interaction is foreseen in future releases Underlying Modelling Methods HAW ASSI AB is based on the following principles e The free surface dynamics for inviscid incompressible fluid in irrotational motion is governed by a set of Hamilton equations for the surface elevation 77 and the potential at the surface e With H 9 77 the Hamiltonian the sum of potential and kinetic energy the Hamilton equations are given by Zakharov 1968 Broer 1974 6 H 7 0 0 6 H 9 7 Here 0 denotes the time derivative and the variational derivative with respect to and similarly for 7 e By approximating the kinetic energy functional K 9 n explicitly as an expression in 7 and the simulation of the interior flow can be avoided the Boussinesq character of the code e The way of approximating K 9 77 is based on Dirichlet s principle for the boundary value problem in the fluid domain By restricting the set of competing functions in the minimization an approximation of K 9 n is obtained The variational derivative K 9 7 0 is the corresponding consistent approximation of the Dirichlet to Neumann operator e The approximate H
18. amilton system conserves the approximate positive definite total energy exactly avoiding sources of instability e The time dynamics is explicit no CFL conditions are required Time stepping is done with matlab odesolver code with automatic variable time step 2 6 Page ad HAWASSI In AB with exact dispersion the interior flow 1s approximated by using a nonlinear extension of the potential as given by the Airy theory of small amplitude waves and Taylor expansion of the kinetic energy leading to Hamiltonian consistent approximations Numerical Implementation Fourier Integral Operators FIO s multiply the Fourier Transform of a function by the symbol of the operator These are generalizations of Pseudo Differential Operators since for FIO the symbol will depend both on the wave number and on the spatial variable the spatial dependence is for nonlinear extensions and varying bottom FIO s are used in a spatial spectral implementation these are approximated by interpolation techniques to enable efficient Fast FT methods 13 Localization methods a difficult point in Fourier type implementations have been successfully implemented to deal with walls run up breaking waves etc 5 2 2 Units and computational grid HAW ASSI AB expects all quantities to be expressed in S I units m kg s meter kilogram second As a consequence the wave height and water depth are in m wave period in 5 etc 2 7 Page HAW ASSI 2 3 Model fea
19. an der Kroon Fully dispersive dynamic models for surface water waves above varying bottom Part 2 Hybrid spatial spectral implementations Wave Motion 49 2012 198 211 E van Groesen amp Andonowati Fully dispersive dynamic models for surface water waves above varying bottom Part 1 Model equations Wave Motion 48 2011 657 666 E van Groesen amp Andonowati Time accurate AB simulations of irregular coastal waves above bathymetry Proceedings of the Sixth International Conference on Asian and Pacific Coasts APAC 2011 December 14 16 2011 Hong Kong China World Scientific ISBN 978 98 1 4366 47 2 pp 1854 1864 E van Groesen T Bunnik amp Andonowati Surface wave modelling and simulation for wave tanks and coastal areas International Conference on Developments in Marine CFD 18 19 November 2011 Chennai India RINA ISBN 978 1 905040 92 6 p 59 63 L She Liam amp E van Groesen Variational derivation of KP type equations Physics Letters A 374 2010 411 415 E van Groesen Andonowati L She Liam amp I Lakhturov Accurate modelling of uni directional surface waves Journal of Computational and Applied Mathematics 234 2010 1747 1756 E van Groesen L She Liam I Lakhturov amp Andonowati Deep water Periodic waves as Hamiltonian Relative Equilibria Proceedings of Waves 2007 N Biggs e a eds Reading UK 23 27 July 2007 pp 482 484 E van Groesen amp Andonowati Variational derivation of KdV type of models fo
20. animation Save data Run 4 16 Page al HAW ASSI The working directory will be as selected in Main GUI after automatic pop up when a simulation is finished Else the directory can be selected There are several panels in the PP GUI Simulation data is automatically loaded with results after finishing a simulation case data of previous projects including Test Cases can be selected to be loaded using Other Plotting Profile and Buoy to plot multiple wave profile s time signal s at specified time s position s with various options to include spectra bottom MTA maximal temporal crests and troughs and quantitative information see below Animation with options to make a gif movie on specified x and f interval Validation to compare simulation results with experimental data or other simulations For comparing at a certain position a time signal as simulated with a signal from a measurement or a previous simulation there is the option Time shifting to shift the simulated signal automatically optimized for best correlation or with a user specified number of time steps Quantitative information is provided also see below Setting options for any of the above graphical animation output methods Be aware that coarsening may change the quality such as time traces or profiles spectra etc and will influence the validation results Quantitative information To analyse properties of simulated results
21. ditional input fields may appear that have to be filled out User input is accepted for various purposes to replace pre programmed choices this 1s the case for e influx signal ows INPUT e initial wave profile D ER e Bathymetry e dispersion function Perd v Depth h m 50 e measurement data for comparison with simulations v Amplitude a m 6 The required lay out of such files or formula s is described in Section 4 4 Calculate OUTPUT l There is a simple Wave Calculator that expects as input the period frequency TIUS or wave length of a harmonic wave and the depth and will then calculate all festive wave number en 252 Wave length lambda 125 m other wave relevant quantities by also specifying the amplitude the tese spess 7 solm lambda h 22 5 calculated steepness is added Category Intermediate depth Steepness k a 0 302 Dispersion Exact 4 I3 Page ad HAWASSI 4 1 Main GUI Choosing the wave model characteristics wave parameters the domain input signal and initial profile are all managed in the Main GUI An overview of the GUI with its main functionalities and input requirements is shown in the Figure below some of the ingredients are described thereafter HAWASSI ABI Main ba E tw File Modules Help Working directory ca ikal Labmalh Indpnesia Project name Note Model Aspace Dynamic HS1 i Interval m Influx position
22. e at most a fraction 2m with m the so called cutfrac in the input for the specified code Note This cut only applies to the terms in the nonlinear equation the spatial grid remains as determined by p i e 2 gridpoints 2 3 5 3 Breaking 10 Breaking of waves is modelled with an eddy viscosity method to initiate the breaking process the value of a kinematic initiation breaking criterion has to be specified the quotient of fluid velocity at the crest and the velocity of the crest U C usually in the interval 0 6 1 2 3 6 Internal flow calculations 5 As an option it is possible to calculate in a post processing step but indicated in the preparation step gt the horizontal and vertical velocities and accelerations of the interior fluid motion all components of the total pressure at a user defined grid in horizontal and vertical direction in a specified time interval 2 4 Software facilities Facilities of the software include to be described in Section 4 e GUI for input of wave characteristics and model parameters with efficient project management e GUI for post processing of the output of the wave simulation and comparison with data e GUI for internal flow calculations e Wave Calculator e Time partitioned simulation to reduce computer hardware requirements e Pre processing step with warnings suggestions for improved settings e 18 TestCases see Chapter 5 with examples of various kind several of which include measure
23. ee that Yayasan AB is not liable for any loss or damage arising out of the use of the Software Although much care is taken to arrive at trustful results of simulations with HAWASSI Yayasan AB cannot be held responsible for any result of simulations obtained with the software or consequential actions or calculations that are based on the results e g because of possible bugs wrong use of the software or other causes 1 4 Page dl HAW ASSI 1 Introduction This document is the Manual of HAWASSI AB software that serves as a guide for using and running the software HAWASSI AB simulates phase resolved waves in Horizontal Direction 1HD long crested waves as are generated in wave tanks to simulate on scale coastal and oceanic waves above flat and varying bathymetry and with partially reflecting walls and damping zones Section 2 describes briefly the mathematical background and the capabilities of the code such as the various dispersive and nonlinear properties together with the features of the software it is advised to read this Section before continuing to the rest of the manual Section 3 provides a step by step installation procedure of the software A condensed description to handle the software regarding GUIs and input output parameters is given in Section 4 Section 5 describes briefly the 18 TestCases that show capabilities of the code and its use More test cases will become available on www hawassi labmath indone
24. epth 25 m Signal type User defined Significant wave Height Hs 0 1 m Peak period Tp 12 362 8 Peak frequency nu 0 508 rad s Peak wave number kp e Oat E TEGA fece dci em Reflection of an irregular wave with Peak phase speed t 13 9355 m s Steepness kp Hs 2 0 002 frequency dependent reflection at a wall Relative wave length lambda h 6 8906 from full reflection of long waves to half kp h 0 312 for short waves shown in the spectrum Category Intermediate depth plot red line 5 1 6 1F006Bor Undular bore HSA MTA amp prufile g time 50 Dynamic Model Lb HSGE Dispersion Model OmExact Bathymetry Flat Depth 1 m INITIAL VALUE PROBLEM K rmi Undular bore non breaking Settings and results as in Wei e a See also Testcase IFBrOO5Bor section 5 2 5 for the breaking bore 5 22IPage 5 2 Breaking waves above flat bottom ad HAWASSI 5 2 1 1FBr001Foc MARIN203001 Strong focussing wave Dynamic Model HS4brF Dispersion Model OmExact Bathymetry Flat Dep Signal type Significant wave Height Hs the 1 m User defined 0 038 m Peak period Peak Peak Peak Peak Tp frequency nu wave number wave length phase speed kp 446 s 346 rad s 998 145 m 176 m s Steepnest kp Hs 2 Relative wave length lambda h kp h Category PwWOoON WFR BE 038 1454 998 Intermediate 5 2 2 1FBr002Foc TUD140
25. ll be taken to have value 0 second column the prescribed elevation m third column the prescribed tangential velocity space derivative of the potential m s 4 4 3 Bathymetry A 2 column matrix space point x bathymetry first column the equidistant x value m covering the whole interval second column the corresponding bathymetry m 4 4 4 External measurement data Time signals at m measurement positions Matrix with m columns First row columns 2 to m 1 specify measurement position 0 position x 1 position x m Next rows time and elevation at the measurement positions time elevation 1 elevation m 4 4 5 Interior flow External measurement data Make a separate file for each horizontal measurement point with a name depending on the quantity that has been measured for instance data U XI mat for a measurement of horizontal U or vertical V velocity at position X7 Each data file has the following format Time signals at m measurement positions in the vertical direction matrix 24 T length m4 1 1 1 horizontal position x 7 2 water depth at position x 2 2 m 1 vertical positions of measurement 3 T end m 1 time u 1 u m 4 4 6 Dispersion relation The software can handle different dispersion relations Default is the exact dispersion but other dispersion relations for Shallow Water KdV and BBM dispersion are predefined and can be dealt with for simulations of n
26. m Breaking Grid size p Cutfrac k Dispersion Fart P elah Flat influx method Area Propagation Unis Depth Im Siope Initial condition Bottom friction Zero i E Wall A m xc m st dev Wave Input Wavetype Mone zo Fourier Bdy m Nonlin Adj ic A Im Tp s Y Options Filter E Internal flow Ramp Time interval Time step dt Time interval s ODE partition 7 Default sipartition Time step s Combine output files 4 I4 Page al HAW ASSI Opening File will show gt Open Project to go to a project that has been created before including the provided test cases from which data can be loaded to be inserted in the GUI gt Save Project Saves all data entered in the GUI after pre processing this info will be stored Clear clears the GUI from inserted input gt Quit HAWASSI Opening Modules will give possibility to activate the Post Processing GUI the Internal Flow GUI and the Wave Calculator Help contains info about the loaded version in About this manual in Documentation and Activation is used for loading the licence first use or renewing the licence The Working Directory can be chosen and specified the software will create a new folder named Output if the working directory does not contain this folder yet if the folder already exists it will keep and use it By specifying a Project Name the software will create a subfolder with that name under
27. ment data to compare with simulations e Comparison with experimental data that have been reported in various publications see the references in Section 6 and the examples of TestCases in Section 5 2 D IPage dl HAW ASSI 3 Installing HAWASSI AB software The HAWASSI AB installer will install the HAWASSI AB code including documentation HAWASSI AB software is programmed under the MATLAB environment The compiled MATLAB applications can be run on PC s that do not have MATLAB installed using the MATLAB Compiler Runtime MCR The MCR will install MATLAB Runtime Libraries on the computer The installation consists therefore of two main steps the installation of MCR and the installation of HAW ASSI AB 3 1 System requirements HAWASSI AB v 1 1 can run on Windows operating system with 64bit architecture The minimum memory RAM needed is 2GB for some test cases and extensive applications 4GB RAM or more 3 2 First step Installing MCR The HAW ASSI AB package v 1 1 requires MCR installer for version MATLAB R2014b for Windows operating system 64bit The MCR installer can be downloaded directly from the MATLAB website http www mathworks com products compiler mcr after downloading install the MCR by double clicking the installer and following the instruction in the installation wizard 3 3 Second step Installing HAWASSI AB After installing MCR the installation of HAW ASSI AB can be started by double clicking the AB installer setu
28. nical Physics 2 1968 190 194 6 29 Page
29. onlinear waves over bathymetry A user defined dispersion relation can be given through input panels in matlab formula style Needed are the gt dispersion relation c Q k d which should be defined as an odd function and gt the corresponding group velocity V dQ p dk user user 4 19 Page ad HAWASSI 5 Test cases HAWASSI AB provides 18 Test cases which are identified with a code of which the first letter has the following meaning F for various cases of non breaking waves above Flat bottom B for various cases of non breaking waves above non flat bathymetry Br for various cases of breaking waves above flat or varying bottom R for run up on a coast The basic properties of the test cases are listed with references to relevant publications in the next sections Acknowledgements We are very grateful to be allowed to use measurement data of gt MARIN Maritime Research Institute Netherlands Dr T Bunnik gt TUD Technical University of Delft Prof dr R H M Huijsmans Authors of publications o Beji amp Battjes gt S Beji J Battjes Numerical simulation of nonlinear wave propagation over a bar Coastal Engineering 23 1994 16 gt S Beji J Battjes Experimental investigation of wave propagation over a bar Coastal Engineering 19 1993 151 162 o Ting amp Kirby gt F C Ting J T Kirby Observation of undertow and turbulence in a laboratory surf zone Coastal Engineering 24 1994 5
30. p HAWASSI AB vl 1 exe and following the instructions in the installation wizard During the installation process a copyright and non liability agreement should be accepted to be able to proceed After the installation is finished start HAW ASSI AB from the shortcut on the Desktop In the Main GUI that appears under Help go to Activation and load licence lic Closing the software and starting again the licence will have been activated and the software can run for the licence period If a new version is downloaded and installed the same licence lic file will be valid for the new version until expiration time In the Help Menu of Main GUI the Documentation will show this manual Test Cases can be found in User My Documents HAWASSI AB1 Testcases 3 2 Page ad HAWASSI 4 GUTs of HAWASSI AB For ease of operation HAWASSI AB software includes three GUI s Graphical User Interfaces as input output managers Main GUI for providing Post Processing GUI for Internal Flow GUI for input for the simulation specifying output of the calculating interior flow simulation properties The GUI s will be described briefly in the next 3 sections The meaning of most required input fields needs no or little explanation the choices that can be made will be described The function of and required input format for input panels is indicated when the cursor is moved over it when an optional panel is checked ad
31. r surface water waves Physics Letters A 366 2007 195 201 6 28 Page dl HAW ASSI 6 2 Other references Beji amp Battjes BBM Broer KdV Ting amp Kirby Wei e a YAB LabMath Zakharov S Beji amp J Battjes Numerical simulation of nonlinear wave propagation over a bar Coastal Engineering 23 1994 16 S Beji amp J Battjes Experimental investigation of wave propagation over a bar Coastal Engineering 19 1993 151 162 T B Benjamin J L Bona amp J J Mahony On model equations for long waves in nonlinear dispersive systems Phil Trans Roy Soc London A272 1972 47 L J F Broer Approximate equationsn for long water waves Appl Sc Res 31 1975 377 395 D J Korteweg amp G de Vries On the change of form of long waves advancing in a rectangular canal and anewtype of long stationary waves Phil Mag 39 1895 422 F C Ting amp J T Kirby Observation of undertow and turbulence in a laboratory surf zone Coastal Engineering 24 1994 5 80 G Wei J T Kirby S T Grilli amp R Subramanya A fully nonlinear Boussinesq model for surface waves Part 1 Highly nonlinear unsteady waves Journal of Fluid Mechanics 294 1995 71 92 Y AB LabMath Water Wave Modelling amp Simulation with Introduction to HAW ASSI software V E Zakharov Stability of periodic waves of finite amplitude on the surface of a deep fluid J of Mechanics and Tech
32. rt of people and goods half of all the world crude oil and increasingly more natural gas and for intercontinental telecommunication through cables e Ocean resources of food and minerals are only at the start of discovery profits from wind parks and harvesting of wave energy in coastal areas 1s expanding Therefore a sustainable and safe development of the oceanic and coastal areas is of paramount importance Nowadays that means that for the design of harbours breakwaters and ships calculations are performed with increasingly more accurate and fast simulation tools Tools that are packaged in software based on the basic physical laws that describe the properties of waves the wave ship interaction the forces on structures etc HAWASSI software is aimed to contribute to extend the accuracy capability and speed of existing numerical methods and software using applied mathematical modelling methods that are at the basis A basis with a rich history that is fascinating and challenging Starting in the 18 century with Euler who generalized Newton s law for fluids in the 19 century Airy solved the problem to describe small amplitude surface water waves In that same century many renowned scientists like Scott Russel Stokes Boussinesq Rayleigh and Korteweg amp De Vries investigated the nonlinear aspects of finite amplitude waves As much as possible without the need to fully calculate the internal fluid motion started with Boussinesq in
33. s in the influxed wave when using a nonlinear wave model a nonlinear adjustment zone of length to be specified has to be applied in the adjustment zone a coefficient in front of the nonlinear terms in the Hamiltonian grows from O at the influx point to 1 at the end of the zone Typically especially for harmonic waves the required length will be at least 2 times the peak wavelength and substantially more for steeper waves on shallower water Wave types Any type of waves can be influxed from a user specified time signal The software makes it possible to specify parameters for harmonic waves and for irregular waves with Jonswap JS spectrum for irregular waves the phases are chosen randomly The parameterized influx signal will be stored for possible re use for comparison of different evolution models Any influx signal will start and end by default with a smooth ramp function the length of which can be specified as a number of periods 2 3 4 Initial value problems Instead of wave influxing data for an initial value problem initial elevation profile and initial potential can be user specified or chosen from predefined parameterized cases a Gaussian as a single hump and Nwave for an N wave shaped wave all with zero initial velocity The Gaussian 1s given by specifying the three parameters in n x 0 Aexp x x o and the N wave is the derivative with adjusted amplitude 2 3 5 Model versions evolution equations
34. sia org DEMO version with restricted functionality The Demo version of HAWASSI AB has restricted functionality and facilities Only exact dispersion Only linear and 2 order nonlinearity Only non breaking waves Partially reflective wall with reflection coefficient the same for all frequencies No internal flow calculations Comparison of demo simulations with AB simulations using full functionality instead of comparing with measurement data Full functionality and facilities under licence e Licence for University Thesis Projects e Research Licence for extending capabilities and or functionalities e Licence for companies commercial use tailor made on demand all proceeds will be used at Foundation Yayasan AB for improving extending the software Visit www hawassi labmath indonesia org for further information or send email to licence hawassi labmath indonesia org Users with limited experience in mathematical physical wave modelling may consult the service booklet 1 Water Wave Modelling amp Simulation with Introduction to HAWASSI software YAB LabMath 1 5 Page ad HAWASSI 2 Description of HAWASSI AB 2 1 Introduction HAWASSI AB This section provides background information of HAW ASSI AB about the basic scientific 1deas HAWASSI AB is a software package for the simulation of realistic waves in wave tanks 1HD i e long crested waves as can appear in oceanic and coastal areas with the option of reflections from w
35. te depth 5 3 2 2B002HarmBar Harmonic over bar HS2U OmExact Dynamic Model Dispersion Model Bathymetry Under water bar Significant wave Height Hs 0 027 m Peak period Tp Z026 8 Peak frequency nu 3 l0l rad s Peak wave number kp 1 675 Peak wave length 3 75 m Peak phase speed 14595 ay S Steepness kp Hs 2 0 023 Relative wave length lambda h 9 3752 kp h 0 67 Category Intermediate depth 5 3 3 2B003lIrBar Irregular wave over bar BB94SL Jonswap Low Freq Non breaking Dynamic Model HSZU Dispersion Model OmExact Significant wave Height Hs 0 021 m Peak period Tp 1 666 s Peak frequency nu 394742 rad s Peak wave number kp 2 109 Peak wave length 2 99 Int Peak phase speed l1 799 m s Steepness kp Hs 2 0 022 Influx position 5 7157 m Relative wave length lambda h 7 447 kp h 0 844 Category Intermediate depth W570 Daty MTA amp profile amp tima 7076 L A Sk Ms va MIN Ta Aine ji A ye yy ANI wise y Nala vem Non breaking irregular wave over slope Observe Freak Wave above deep part at x 71 t 282 8 Simulation of MARIN measurement WATU Uety MIA amp prolilo deo ems 190 Harmonic waves over under water bar Strong mode generation compare spectra at x 8 red and x 6 black below PEPP heres eevee 8 bees et tod u LI kasu Simulation of Beji amp Battjes Experiment HEAL Baty MTA amp profile iima 43 5
36. th 0 5 m Signal type Harmonic Significant wave Height Hs Q 006m Peak period Tp LO ts Peak frequency nu 0 62 rad s Peak wave number kp 0 281 Peak wave length 22 384 m Peak phase speed 2 1 T m s Steepness kp Hs 2 0 001 Relative wave length lambda h 44 7688 kp h 0 14 Category Shallow water Dynamic Model HS5Z2brUR Dispersion Model OmExact Bathymetry Slope Shore Max Depth 0 4 m Signal type Harmonic Significant wave Height Hs 0 125 m Peak period Tp 2 s Peak frequency nu 3 141 rad s Peak wave number kp 1 7 Peak wave length 3 695 m Peak phase speed 1 847 m s Steepness kp Hs 2 0 106 Relative wave length lambda h 9 2387 kp h 0 68 Category Intermediate depth WS BH Hathy MIA X profile i time AN T T Harmonic nonbreaking wave run up on 1 25 coast Harmonic wave run up on 1 35 coast spilling breaker Simulation including interior flow properties of Ting amp Kirby experiment 5 27 Page 6 dl HAW ASSI References 6 1 References to basic papers and applications 10 11 12 13 14 I5 16 17 18 19 E van Groesen amp Andonowati Hamiltonian Boussinesq formulation of Wave Ship interactions Part 1 Evolution equations submitted R Kurnia T van den Munckhof C P Poot P Naaijen R H M Huijsmans amp E van Groesen Simulations for design and reconstruction of breaking waves in a wavetank OMAE 2015 V012
37. th with the limiting speed c 4 gd gt KdV dispersion 3 order Taylor 9 QO pyy k note that short waves with kd v1 6 will travel in the opposite direction Remarks Since influxing uses properties of the group velocity Uni directional influxing in the KdV model will show the short waves running in the wrong direction corresponding to the dispersion relation Bi directional influxing will include these wrongly directed waves which is not corresponding to the original KdV dispersion relation for uni directional waves gt BBM dispersion Q T k same as KdV in 3 order but uni directional Note To avoid problems with too poor dispersion KdV and BBM uses exact dispersion for influxing In the table the explicit formulas and plots for the various cases are given Plots Disp relation amp Phase velocity Exact _ Q k _ nan m nh kd Dispersion relations dispersion Water KdV 0 OQ k Cok I L kd UPC 07Og k cok 1 Rd User Provide dispersion relation and group velocity specified in input panels using only wave number k and depth d as variables in matlab formula style i 99 2 JOlPage dl HAW ASSI 2 3 5 2 Nonlinearity The order of nonlinearity of the Hamiltonian System HS is specified by the number in the present version of HAWSSI AB HSm for orders m 1 linear 2 3 and 4 To avoid aliasing in the Fourier implementation the wave numbers have to be restricted to b
38. tures and capabilities HAW ASSI AB accounts for the following physical circumstances and phenomena of waves in 1 HD i e long crested waves 2 3 1 Geometry Simulation interval and Grid A uniform grid is defined by specifying an x interval ET oa and a grid with grid size dx Cr Keh N 1 where N 2 with p an integer Wave numbers k in Fourier space are defined in accordance with the spatial grid k N 2 1 N 2 xdk with dk 2z x Xun Fourier Boundary For the use of Fast Fourier Transformations all quantities except bathymetry are tapered to vanish near the end points this takes place is the so called Fourier Boundary the length of which can be specified The Fourier Boundary should be such that reflection of outgoing waves has to be prevented hence the Fourier Boundary also acts as a damping zone Walls 5 The position of a wall inside the simulation interval can be specified Depending on the reflection properties the following choices can be made gt Auniform partially reflecting wall by specifying the reflection coefficient in 0 1 for all wave lengths frequencies Afrequency depending non uniform partially reflecting wall by providing a reflection function with reflection coefficients depending on frequency in the input panel using only the frequency f as variable in matlab formula style 2 3 2 Bathymetry and run up The bathymetry can be user specified The software provides parameteri
39. zed bathymetries for flat bottom depth and linear sloping parts of the bottom including run up Bottom friction Bottom friction can be applied at a specified part of the bottom in the bottom friction formula C pur u Ju D n typical values for the friction coefficient are C 107 107 depending on Reynolds number and bottom 2 3 3 Embedded wave influxing of various wave types Wave influxing 9 In AB the wave influxing is done through a source in the continuity equation 2 8 Page al HAW ASSI The influx position and a time signal of the desired elevation at the influx position have to be specified A choice can be made between uni directional influxing for waves propagating in one direction in the direction of the positive unit or negative uni x axis or bi directional influxing for waves running symmetrically in both directions The spatial extent over which the influxing takes place can be chosen With a point influx the generation area is restricted to the immediate neighbourhood of the influx point using Dirac delta functions Then the desired time signal is modified into a much higher modified time signal to guarantee the correct waves being influxed A smoother influx better suited for high steep waves is area or area short influxing then the waves are generated over a broader interval better suited for high steep waves Nonlinear adjustment zone 9 In order to prevent spurious modulation
Download Pdf Manuals
Related Search