Home

HAWASSI-VBM2 User Manual by © LabMath

1. sssss 5 78 Figure 5 41 Model setup of test case 6a 0 ccccceeeeeeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 5 78 Figure 5 42 Test case 6a Snapshot of the surface elevation esses eene 5 79 Figure 5 45 Domain settings of test case OD oer reoxu ses ii 5 80 Figure 5 44 Settings for mesh properties and model control for the test case 6b 5 81 Figure 2 45 Model Setup of test Case OD ud aana aka a a anaa paka ii 5 81 Figure 5 46 Test case 6b Snapshot of the surface elevation eese eene 5 82 1 5 Page ad HAWASSI List of Tables Table 4 1 Data format for user defined initial conditions eee 4 20 Table 4 2 User defined influx data format eese nennen nennen nne nn enne nnns 4 26 Table 4 5 Conste daba O iia UN ai ne a SD SLE saa 4 30 Table 4 4 Bathymetry data format esssssssssessesesssssseseseeeeeeeee nennen nnn nnne nennen sss nsns nenen anane 4 3 Table 4 5 Output HIes TOIdBTS rro AA ARA 4 40 l 6lPage dll 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 transport of people and goods half of all the world crude oil and increasingly more natural gas a
2. Figure 5 24 Settings for mesh properties and model control for the test case 2 2D signal Tl 1000 bathymetry amp numerical setup t 1200 bj E 500 gt 1000 F D 800 F E m 200 1 E 600 F s m 0 0 t 1000 0 1000 400 x m Amp Spec of 1D Signal 0 00 0 14 0 22 0 27 0 32 0 35 200 F OF 2 200 EBEN l T MEN E 1000 500 500 1000 3 o x m D 1 Dispersion Quality 0 1 02 03 0 1 2 3 4 5 f Hz kh Figure 5 25 Model setup of test case 2 5 66 Page ad HAWASSI Figure 5 26 The upper plot shows the snapshot of surface elevation and the lower plot shows contour line of the maximum wave height Observation o Look at the diffraction pattern especially behind the wall o Try to investigate the dependence of the diffraction pattern on wavelength o Use various influx waves also Jonswap influx type o Use various wave models 5 67 Page ad HAWASSI 5 4 Test Case 3 Refraction Test case 3 is a continuation of test case 1 now for a bathymetry that includes a shallow part that will show the refraction phenomenon i e the bending of waves towards the shallow area because of reduced speed over shallower areas The domain configuration settings for the mesh generator and model control are illustrated in Figure 5 27 Figure 5 28 and Figure 5 29 File Edit Model Setting Log Tool Help Physical properties Simulation Boundary Bathymetry Mesh data imported data Y p v 1 1 Show bat
3. Non experienced users are suggested to start using the software with the test cases that are provided in order to understand the input requirements and to explore the various possibilities of output from the simulation 2 2 Model features HAWASSI VBM2 accounts for the following physical phenomena of waves in 2 HD horizontal directions long and short crested waves e Wave propagation dispersion diffraction refraction and shoaling e Nonlinear wave wave interactions Features of the software include e The quality of dispersion is optimized for specific wave problem to be simulated which makes it possible to simulate deep ocean waves or very short waves kh 15 or more and infragravity waves e Various method for wave influx e Use of efficient damping zones and partially reflecting walls Facilities of the software include e GUI Graphical User Interface for making a depth dependent unstructured mesh in a given geometry with sea and or wall boundaries and for input of wave characteristics and model parameters e GUI for post processing of the output of the wave simulation e Time partitioned simulation is possible to reduce computer hardware requirements e Project examples with harmonic and irregular waves above bathymetry in simple as well as complicated geometries such as harbours Simulation results of the model have been compared with series of experiments in laboratories see the references in Section 6 2 lllPage dl
4. 3 843 12552 1892 43 ooo occ coc co cc coc Fe 579 1454 3076 518 14716 26213 37991 38124 22898 594 1144 26 cOooooocococcocoocooocsso 244 613 1295 3172 6198 11058 16021 16331 9672 3198 482 jk O Oooo occ cc c oc EI 51 128 271 662 1296 2317 3349 3410 2018 669 101 OOO OOO OOO S qc N 685lPage
5. Boundary types can be specified by using tools in Physical properties panel see Figure 5 8 The Damping Zone is used to dissipate absorb the outcoming waves The Partial Reflective Boundary PRB is used to have a controllable reflection on the boundary 5 3 lPage ad HAWASSI Physical properties delete influx Damping Zone damp 3 m Figure 5 8 Physical properties panel Activate and choose Damping Zone area in the following steps e Click Damping Zone button in Physical properties panel Figure 5 8 e Put four points that will become the area of the damping zone the first two points create a dashed line that indicates the beginning of the Damping Zone see Figure 5 9 e Put four points in the east position ee ee te E BAENA YO 500 0 500 1000 Figure 5 9 Creating Damping Zone in HAWASSI VBM2 552lPage ad HAWASSI e Editing the damping points can be done by activating physical properties module in the Edit menu then select M to edit the points by specifying the x y position To have a good damping boundary the width of the damping zone should cover at least two times the longest wave YA RO E 1400 Ih x y I 0 200 400 600 800 400 200 Figure 5 10 Creating 3 Damping Zones as transparent boundaries at the sea sides To select Partial Reflective Boundaries PRB
6. The Demo version of HAWASSI VBM2 has restricted functionality and facilities Only 1 parabolic vertical profile Only linear Non depth dependence mesh generator No internal flow calculations No comparison with external measurement data Full functionality and facilities under licence Licence for University Thesis Projects gt Research Licence for extending capabilities and or functionalities 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 1 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 9 Page dll HAW ASSI 2 Description of HAWASSI VBM2 2 1 Introduction HAWASSI VBM2 The purpose of this section is to provide the user relevant background information of HAWASSI VBM2 and to give some general advice in choosing the basic input for the computations HAWASSI VBM is a software package for the simulation of realistic waves in wave tanks 1HD 1 Horizontal Dimension oceanic and coastal areas harbor etc 2HD 2 Horizontal Dimensions The acronym of HAWASSI stands for Hamiltonian Wave Ship Structure Interaction HAWASSI VBM is a finite element imp
7. and investigate wave resonances there o Compare the wave disturbance caused by long crested monochromatic and irregular waves by short crested irregular waves linear and nonlinear simulations 5 79 Page ad HAWASSI 5 7 2 Breakwater configuration 2 The domain configuration as built in 2009 settings for the mesh generator and model control are illustrated in Figure 5 43 Figure 5 44 and Figure 5 45 File Edit Model Setting Log Tool Help v 1 1 Working Directory C Users DiditiDocuments M m model OVBM 1 profilegW Nonlinearity Initial cond Influx 1 Harmonic Time 0 1 300 Project Name TC6b RealisticHarbour v2t User s note Wave entering realistic har Simulation Boundary Bathymetry Mesh data y Depth imported data Show bathymetry A 1400 hJ E e q A 9 HK A oa f d a a 200 400 600 600 1000 1200 1400 Project Setting file is loaded Figure 5 43 Domain settings of test case 6b Physical properties Influx line 1 delete influx damp 3 0 2 Generate Mesh Preparation RUN STOP Post Processing 5 80 l Page y m 1400 F 1200 f 1000 800 600 400 sea E mesh properties Depth dependance Parameters Approx Tp s ppl max ds m shallowest depth m ad HAWASSI E model control Initial conditions Zero Time Signal Harmonic Influx 1 Direction deg 270
8. in the Desktop or in the Start menu The main GUI of HAWASSI VBM2 will appear as in Figure 5 1 Instruction of what to do in the next step is shown the comment box panel see Figure 5 2 After the software is opened select your Working Directory by clicking mm button It is advised to use CAUsersM UserName MDocumentAHAWASSI_VBM folder as the working directory see Figure 5 3 Test cases of HAWASSI VBM2 are stored in this folder named as TestCases_VBM2 548lPage ad HAWASSI File Edit Model Setting Log Tool Help semen C Show bathymetry WAKO bh Figure 5 1 Main GUI of VBM2 Nic please select your working directory Figure 5 2 Comment box panel showing short instruction for running the software wt Je Documentation Ji Output TestCases VBMI TestCases VBM2 gt j HAWASSI VBM Userlnput VEM Folder HAWASSI VBM Figure 5 3 Select working directory in My Document 5 49 Page ad HAWASSI 5 1 2 Selecting wave model HAW ASSI VBM provides two wave model the non dispersive Shallow Water Equation SWE and the dispersive Variational Boussinesq Model VBM By default the chosen model will be the linear Optimized VBM with one vertical profile see Figure 5 4 Help 1 Setting Log Tool Type SWE Model control Ctrl M VBM j Parabolic Profile Nonlin
9. the mesh generator will limit the smallest depth values with the Shallowest depth value The value should be a positive number If Depth dependence is chosen the user has to specify mesh size parameters in the Parameters panel Approx Tp for the approximate peak period of the waves to be simulated and ppl as abbreviation for point per wavelength denotes the number of points 1n the mesh to represent a peak wavelength Warning 1 Changing the project name after creating the mesh will cause the output of simulation be stored in the different output folder 2 Mesh size can become very small if the depth bathymetry is very shallow at some parts which may lead to longer computation times 3 Shallowest depth should be a positive value which means that the bathymetry data should only contain positive depth values and does not has topography land negative depth value If topography land is included in the bathymetry data the software will automatically only take the depth values and replaces the topography land data by Shallowest depth value 4 35 Page dll HAW ASSI 4 1 10 Advanced Settings HAWASSI VBM2 provides an additional GUI so called Advanced Settings to change given standard settings that are used in the simulation The Advanced Settings interface can be called by clicking the menu Setting gt Advanced Settings see Figure 4 34 and Figure 4 35 Physical parameters that
10. with respect to the depth the length should be more than 2 peak wavelengths Nonlinear adjustment Ares 2 Figure 4 22 Nonlinear adjustment for nonlinear model 4 28lPage ad HAWASSI 4 1 7 Simulation Time The time interval for the simulation can be specified in the box Time in the Model Control GUP Figure 4 23 tstart is the starting time of the simulation should be a non negative value dt is the output time discretization for the simulation should be a positive value and tend is the ending time of the simulation All values should be in seconds Time discretization dt is not the time discretization for calculation of the time integration HAW ASSI VBM uses automatic internal time stepping in the matlab odesolver tstart dt tend 0 05 200 Figure 4 23 Time control of simulation 4 1 8 Simulation Boundary Bathymetry and Boundary Conditions 4 1 8 1 Simulation Boundary The simulation boundary defines the boundary of the computational domain The user can create and or edit the boundary points by using available tools in the main GUI of HAWASSI VBM2 see Figure 4 25 When these tools are not active the user can activate them by clicking Edit gt Scatter module For loading an existing simulation boundary the user should click Simulation Boundary button in the main GUI the user should then locate the scatter data of coastline The data format of the coastline is illust
11. 18 18 22 28 20 29 32 EE RENE DEPRECOR aen 400 600 00 1000 1200 x m 0 3 e Pese z El E UE T ree Corre a Ga aa od po O 0 1 2 0 3 f Hz Figure 5 22 Model setup of test case 1 overview as result of Preparation 563 lPage ad HAWASSI Suggestions o Look at the log file to find the calculated wave length How many points will be in one wavelength for the given dx Estimate how many waves you will see in which interval with the expected amplitude o Runthe simulation and perform post processing Verify if your expectations are correct Perform other simulations for this wave case o See effect of changing grid size by changing ppl parameter in the mesh properties change period and choose sensible grid size o Change influx line position damping zone boundary values etc o Perform nonlinear simulations see effect of increasing the amplitude and explain o Change wave type to Jonswap Hs 1 Tp 12 Gamma 3 Estimate peak wave length Activate MTA when plotting animations and snapshots NOTE for comparison of different wave models the same initial signal has to be used it has to be loaded from a previous simulation since the parameters alone do NOT specify the Jonswap input signal because phases are added randomly o Check performance of linear simulation using SWE Observe that propagation is a pure translation o Now use VBM dispersion parabolic 1 and 2 profiles and observe differences
12. G File Sd Ka 2 19 9 D E3 2D signal bathymetry amp numerical setup E 200 A00 e 100 sim Oo D 500 1000 E 200 x m TA Amp Spec of 1D Signal 0 00 0 11 0 18 0 22 0 26 0 29 0 200 400 600 800 1000 1200 x m C amp V m s Dispersion Quality 0 0 1 2 3 4 5 f Hz kh e Numerical setup Domain of computation Damping zone Partial Reflective Boundary Embedded wave generation e Bathymetry Amplitude Spectrum Dispersion quality of the influx signal Velocity of Exact vs Model Figure 4 38 Overview of the numerical setup panel showing details of the simulation to be performed 0 00 0 11 0 18 0 22 026 0 29 Dispersion quality Velocity of Exact Airy theory versus Model Phase Velocity m s C amp V m s meee Group Velocity m s Exact Airy theory Dispersion Quality a 0 1 2 3 4 5 kh kh of the peak wave to be simulated Figure 4 39 Description of Dispersion Quality plot blue solid for Exact red dashed for Model set up 4 39 Page dll HAW ASSI 4 1 12 Output Files Output files will be stored in the folder WorkingDirectory Output ProjectName A description of all output files is given in Table 4 5 Table 4 5 Output files folders No Output Files Folders 1 INPUT ProjectName mat Input data for the simulation a result from Preparation 2 DATA ProjectName iter i mat
13. Output data from the i th step of the partitioned simulation a result from RUN If the Split Output file in the Advanced Settings GUI is activated then the data of simulation will be split into the number of Time Partition If the Split Output file is not active the iteration is only 1 3 LOG ProjectName log A log file of the case ProjectName If the ProjectName is not changed details of anew simulation will be added with to the existing log file INFLUX ProjectName txt The influx signal data txt or ASCII file that is used in the simulation only available when the Save influx signal in the Advanced Settings GUI was checked B HOTSTART ProjectName t time txt The last state of the surface elevation n and surface potential at the end of a time partition only available when the user had checked the Hotstart checkbox in the Advanced Settings GUI Is a sub folder that will contain animation gif files from the PostProcessing GUI Is a sub folder that will contain 3D animation gif files from the PostProcessing GUI Is a sub folder that will contain figure fig amp png files from the PostProcessing GUI 10 Is a sub folder that will contain frames png or jpg from Post Processing GUI 4 A40 l Page ad HAWASSI 4 2 Post Processing GUI After a simulation is finished the Post Processing GUI will automatically pop up and loads the simulation data The Post Processing
14. Tp 5 g Amplitude m 0 25 Figure 5 44 Settings for mesh properties and model control for the test case 6b bathymetry amp numerical setup A A 200 600 x m 200 400 900 2D signal 0 2 4 0 2 4 300 500 200 sim og 100 200 0 200 400 600 x m Amp Spec of 1D Signal 0 00 0 16 0 26 0 33 0 38 0 42 0 1 TETT AE se 2 6 z amp 0 05 gt 4 5 0 ikai 0 0 2 0 4 0 Figure 5 45 Model setup of test case 6b 5 8llPage ad HAWASSI t 00 04 18 hr min sec or t 259 87 sec idt n Li ta E wj JL LA E us we 1 X L cR Figure 5 46 Test case 6b Snapshot of the surface elevation Observation o Harbour lay out in Google Earth 7043 44 S 109001 32 E in gt 2009 o Different type of harbour waves depending on wave type and on main incoming direction Suggestion o Quantify the wave energy inside the fishing harbour and investigate wave resonances there o Compare wave disturbance simulations with long crested monochromatic and irregular waves with simulations of short crested irregular waves linear and nonlinear 5 82 Page dll HAW ASSI 6 References References to basic papers and application of VBM D Adytia Simulations of short crested harbour waves with variational Boussinesq modelling In Proceedings SOPE 2014 2014 912 918 D Adytia M Woran amp E van Groesen Effect of a possible A
15. amp E van Groesen Variational Boussinesq model for simulations of coastal waves and tsunamis Proceedings of the 5 International Conference on Asian Pacific Coasts APAC2009 13 16 October 2009 Singapore 9ed Soon Keat Tan Zhenhua Huang World Scientific 2010 ISBN 13 978 981 4287 94 4 Volume 1 ISBN 13 978 981 4287 96 8 pages 122 128 D Adytia A Sopaheluwakan amp E van Groesen Tsunami waveguiding phenomenon and simulation above synthetic bathymetry and Indonesian coastal area Proceedings of International Conference on Tsunami Warning Bali Indonesia November 12 14 2008 E van Groesen D Adytia amp Andonowati Near coast tsunami waveguiding phenomenon and simulations Natural Hazards and Earth System Sciences 8 2008 175 185 G Klopman M Dingemans amp E Van Groesen Propagation of wave groups over bathymetry using a variational Boussinesq model Proceedings nt Workshop on Water Waves and Floating Bodies eds Sime Malenica and Ivo Senjanovi Plitvice Croatia April 2007 pp125 128 G Klopman M Dingemans amp E van Groesen A variational model for fully non linear water waves of Boussinesq type Proceedings of 20 International Workshop on Water Waves and Floating Bodies Spitsbergen Norway 29 May 1 June 2005 6 83 Page ad HAWASSI Appendix I SWAN 1 owan standard spectral file version Data produced by SWAN version 41 01 Project projname Fun number runnun LOCATIONS locations in x y
16. can be changed in Advanced Settings GUI are the value for the gravitational acceleration g m s fluid density p kg m File Edit Model e Advanced Settings Figure 4 34 Calling Advanced Settings GUI In the Advanced Settings it is also possible to modify settings in the Time Stepping such as Time partition Hotstart Save signal influx and Split Output Data see Figure 4 35 Time partition is an option to divide the simulation calculation into several steps the default is 1 to reduce memory RAM requirements see Figure 4 36 This facility is useful especially when dealing with long time simulations or when solving a large problem When the Split Output Data option 1s checked the calculated data will be divided into the number of time partitions and stored as Output ProjectName DATA ProjectName iter mat with i denoting the i th partition The data of the i th iteration will be directly saved into the hard disk after the i th iteration is finished If it is unchecked by default the simulation data will be saved as one file named as Output ProjectName DATA iterl mat For saving the last state of the surface elevation n and surface potential o at the end of a specific time partition saving done as a txt ASCII file the option Hotstart should be checked is default The option Save signal influx is created to save the influx signal if it is a signa
17. compare outputs by combining them in a single plot Observe effects of changing Jonswap parameters o Perform non linear simulations increasing Hs S564lPage ad HAWASSI 5 3 Test Case 2 Diffraction Test case 2 1s an illustration of the diffraction phenomenon The configuration is chosen to mimic a semi analytical solution of the Helmholtz equation The contour lines that are a result of the VBM simulation shown in in Figure 5 26 illustrate the diffraction these results are comparable with the semi analytical solution see Shore Protection Manual 1984 page 124 The domain configuration setting for the mesh generator and model control are illustrated in Figure 5 23 Figure 5 24 and Figure 5 25 Influx line delete influx Nonlinearity No Initial cond Zero Influx 1 Harmonic Project Name TC2_diffraction User s note Diffraction phenomenon Figure 5 23 Domain settings of test case 2 e US Army Corps of Engineers 1984 Coastal Engineering Research Center Shore Protection Manual vol 1 5 65 Page y m ad HAWASSI E mesh properties Depth dependance Initial conditions Amplitude m Parameters A T E r Approx Tp s Width hump x v m zd 9s Location x y m Rotation angle deg Geometry max ds m Time 5ignal shallowest depth m Influx Direction deg Tp s Amplitude m
18. space 1 number of locations 22222 22 0 00 RFREQ relative frequencies in Hz 20 number of frequencies 0418 0477 0545 0622 0710 0810 0924 1055 1204 1375 1569 1791 2045 2334 2664 3040 3470 3961 4522 5161 b891 6724 f675 8 61 0000 CDIR spectral Cartesian directions in degr 12 number of directions 30 0000 60 0000 90 0000 120 0000 D 00000000000000000000000 6 84 Page 150 0000 180 0000 210 0000 240 0000 270 0000 300 0000 330 0000 360 0000 QUANT 1 VaDens m2 Hz degr 0 9900E 02 FACTOR 0 675611E 06 51 242 574 129 610 1443 273 1287 3054 665 3152 7463 1302 6159 14608 2328 10989 26020 3365 15922 37712 3426 16230 38440 2027 9612 22730 672 3178 7538 101 479 1135 2 11 26 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ad HAWASSI number of quantities in table variance densities in m2 Hz degr unit exception value 956 1288 1482 1481 1286 2402 3238 3 25 3724 3234 o084 6846 7872 7869 6837 12402 16712 24275 32688 43341 58358 62733 84492 63939 86109 37790 50909 12535 16882 1890 2542 2924 2923 2539 43 E a occ coc co cc co b of oo co 000 O 00 DO O I 19229 3 618 67109 97150 99010 908529 19440 66 coco oo cc Ka E Ka E Ka Ka kek 19221 16690 3 603 32644 67080 58281 97110 84380 98969 85995 98905 50841 19432 16870 66 of O DO o cco coco co ccc A 957 2406 5091 12419 24309 43401 62820 64027
19. the following steps should be done 1 Click PRB button in the Physical properties panel Figure 5 8 2 Select two points start and end point in the coastline data that will mark the boundary line of PRB The boundary line is in the counter clockwise direction from start to end point The reflection coefficient is in the range of O to 1 with 1for full reflection 5 533lPage dll HAWASSI WAWOR Figure 5 11 Creating PRB s at the land boundaries and harbour walls 5 1 6 Specifying the influx For influxing waves into the domain of computation the user needs to define an embedded influx line in the domain and provide the information for the time signal via the following steps 1 Click Influx line button in the Physical properties panel 2 Put two points that will create a line in the domain if the point is placed outside the domain it will be projected to the closest point on the domain boundary e For editing the influx line position activate the editing tools by in Edit menu gt Physical properties module then click p to edit the starting and ending points of the influx line by specifying the x y position 5 54 Page ad HAWASSI 400 200 0 200 400 600 800 Figure 5 12 Creating influx line 3 Define the influx signal by calling the Model Control GUI click Model menu gt Select Model Control or simply Ctrl M 4 Select Harmonic in the T
20. 0 f end Hz Figure 5 32 Settings for mesh properties and model control for the test case 4 bathymetry amp numerical setup 900 800 700 600 95 7 400 L 300 lb M ae eer MAD NEM nues E ie 200 400 600 600 x m 2D signal ng 800 600 E 400 ps t3 200 ex es Py lt 200 200 ah sim 00 t T 500 1000 Amp Spec of 1D Signal 0 000 110 180 220 260 290 32 C amp Y m s Dispersion Quality 0 01 02 03 012 34 5 6 f Hz kh Figure 5 33 Model setup of test case 4 5 72 Page ad HAWASSI t 00 02 39 hr min sec or t 159 80 sec mm d Figure 5 34 Test case 4 Snapshot of surface elevation in a simple harbour layout Observation o Compare wave disturbance from harmonic waves of different wavelength o Compare results of various wave models 1 and 2 profiles linear and nonlinear o Quantify the wave energy inside the harbour at one instant decide from that when starting from an initially flat sea the waves have completely occupied filled the harbour o Any standing waves Investigate buoy measurements at clever positions to support claims about standing waves Suggestion o Calculate on paper some simple standing wave patterns for the right side of the basin estimate the wavelengths that play a role and perform simulations to investigate the occurrence of such harbour resonances 5 73 Page ad HAWASSI 5 6 Test Case
21. 2 can be started as follows 1 Click the Preparation button this gives an overview to check the domain settings and the numerical setup based on the given input as illustrated in Figure 5 15 2 Click the activated RUN button After the calculation is finished there are several new data files created in the folder Output ProjectName containing detailed information of the performed simulation Details of names and contents of the files are described in Table 4 5 5 57 Page 1400 F bathymetry amp numerical setup 1200 p 1000 y m Bn 600 AUC 5 1 9 Post processing of the output RA 200 O 200 x m 400 Ee oar ROC EUER Te Tee mee DT A EE aN A Bn ad HAWASSI 0 signal 200 sim 0 10d 200 O 200 400 600 x m Amp pec of 1D Signal 0 00 0 16 0 26 0 33 0 38 0 42 e A E EMT z C amp V m s D 0 2 0 4 f Hz Figure 5 15 Overview of the numerical setup after Preparation is finished After the calculation is finished the Post processing GUI automatically pops up Figure 5 16 Post processing GUI The GUI can also be used to process existing output data by selecting the option Other in the Simulation Data panel and specifying the file through the dialog box Select a mat file To show the animation 2D or 3D for test case 6b the following steps have to be done mS STO Activate Animation panel in Post
22. 4 Figure 4 48 Setting land and wave color in 3D setting esses eee 4 45 Figure 4 49 Viewpoint with azimuth and elevation s anana naen 4 45 Figure 4 50 Bathymetry significant wave height Hs MTA Hamiltonian Energy and Wave disturbances M 4 46 Figure 4 51 Buoy setting option of Plotting sesanan aana aana aana eene 4 46 Figure 4 52 Option for filtering the extracted signal at specific Buoy location s 4 46 Figure 4 53 Snapshot setting option of Plotting sasapan anana nana Nak aaa Waa anaa aaa eS Na aaa kag KANA 4 47 Pour S L Mita CULO VBM a 5 49 Figure 5 2 Comment box panel showing short instruction for running the software 5 49 Figure 5 3 Select working directory in My Document ssssssseessssssseseeeeeeeeeeeee nennen nennen nnns 5 49 PISE 3A Select wave Mod ua 5 50 Figure 5 5 Load Simulation boundary data sess eene eene 5 50 Figure 5 6 Simulation boundary data is loaded sess eene 5 51 uti CUM Bro Tu babine AMA aa aa a ag aa BN Waaa a ANI a a ba anaa ba a aaa 5 51 Figure 5 5 Physical properties panel nana 5 52 Figure 5 9 Creating Damping Zone in HAWASSI VBM2 eeesssssssseseseeeeeseeeeseeeseeeeeeeeeeeeeeeeseeeeeeeessees 5 52 Figure 5 10 Creating 3 Damping Zones as transparent boundarie
23. 5 Simple harbour layout 2 Test case 5 is a harbour with a basin with an oblique basin all with constant depth The domain configuration setting for the mesh generator and model control are illustrated in Figures 5 39 5 40 5 41 A plot of a simulation 1s shown in Figure 5 38 Physical properties Influx line E 2 E Hh a delete influx Project Name TC5 SimpleHarbour2 User s note Wave entering a simple hai Figure 5 35 Domain settings of test case 5 5 74 Page y m EY mesh properties Depth dependance ad HAWASSI E model control Initial conditions Amplitude m Parameters Approx Tp s shallowest depth m Width hump x y m ppl Location x y m Rotation angle deg Geometry max ds m Time Signal Harmonic Direction deg Influx Tp 5 Amplitude m Spreading s cas 25 theta Filter f_0 f end Hz Figure 5 36 Settings for mesh properties and model control for the test case 5 2D signal pah yaa amp numerical setup l 1000 1200 0 E 500 1000 1 14 5 800 1000 200 600 sim 00 t 0 500 1000 1500 15 x m 400 Amp Spec of 1D Signal 0 000 110 180 220 26 0 290 32 200 H 0 3 15 5 y DH T 0 2 gt 200 3 16 0 1 o sa Ted Dis ersion Qualit x m 0 DISp y 0 0 1 02 03 0 1 2 3 4 5 6 f Hz kh Figure 5 37 Model setup of tes
24. Conditions and Time Signal information to be used for the simulation as shown in Figure 4 8 The GUI can be called manually by clicking Model Model control 4 18 Page ad HAWASSI EY model control File Edit Setting Log Type Initial conditions Zero Model control Ctrl M Time Signal Harmonic Influx 1 Direction deg Tp s Amplitude m Figure 4 8 Model Control 4 1 5 Initial Conditions In the Model Control GUI Initial conditions for the surface elevation n and the surface potential 0 have to be specified as shown in Figure 4 9 Choosing option Zero means that the surface elevation n and surface potential o have initial value zero flat water surface at rest Initial conditions can be specified with the option User defined the user will then be asked to click the file name of the prepared initial conditions through a dialog box The data format for the user defined initial conditions is illustrated in Table 4 1 The data should consist of four columns the first column specifies the discretization of the horizontal x coordinate the second specifies the discretization of the horizontal y coordinate the third and 4 9 Page ad HAWASSI fourth columns are the data of n and 6 respectively If only three columns are specified this will be interpreted as zero initial condition for no velocity EY model contral Initial conditions Fero User defined select init
25. GUI Figure 4 40 can also be called directly from the Main GUI by pressing the PostProcessing button To load data of previous simulations select other under the simulation data section in the GUI as illustrated in Figure 4 41 There are two main panels in the Post Processing GUI i e Animation and Plotting Each panel will be described in the next two subsections E Post Processing 2D v 1 1 Plotting Bathymetry Hamiltonian Enerqy Significant Wave Height Wave Disturbances Animation 4 LITA tinit tend Hat duis Time Xwest Xeast Ysouth Ynorth D Axis Lim Buoy s OCY save as GIF delay loop Signal Amp smoother save frames format ong Spectrum Energy 3D view 3D setting Xwest Xeast Y south north Snapshot t ID save plot s Figure 4 40 Post Processing GUI 4 4llPage dl HAWASSI Simulation Data uploaded file loaded loaded simulation data for Post Processing Figure 4 41 Load other simulation data in Post Processing GUI 4 2 1 Animation The Animation panel is to show animations of the wave elevation Information of domain and time interval for the animation is needed By default these fields are automatically filled based on the uploaded simulation data the whole area and whole time interval A restricted domain for the animation can be specified in Axis Lim see Figure 4 42 Animation Awest Xeast Y south Y
26. HAWASSI VBM2 User Manual HAWASSI 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 org VIA IF UE Copyright 02015 LabMath Indonesia dll HAW ASSI Contents A o a aaa Ta a E ag aga aaa da a aaa aan ag a ga aaa das aa baga aaa asana 1 7 RE JABATAN A DH srana a a e Eo o e E 1 9 2 Description of ITAWASSI VBMOA uti anan anana anna ene nana aaa nenen nenen naen 2 10 24 rene Ps Sich Ao iU Pr 2 10 2 2 ot Mi Dl PPP Po Un em aaa bles o o o t IP URN LES aa KA EUIS 2 11 2 3 Relation to other Boussinesq type wave models ooooooooooncnccccccnnnnnnnnnnnonononnnnnonoronnnnnnnonnnnnnnnos 2 12 2 4 X Units coordinate system and computational grid ooooononnnnnnccnonnnnnnnnnononononononanonnnnnnnnnnnnnnnnnos 2 12 3 Installing HAWASSI VBM software seeesseseeeeee aana nana nennen nennen nnn nnn nnne nnne nenne nne ene eee 3 13 3 1 System reduire MONIS 525 CO TET t A 3 13 22 Turststep Installing IVC Rusia air peace 3 13 3 3 Second step Installing HAWASSI VBM aaeeeeeenenanaaaa anana anan anana nana aaa nana anane nnn nn nennen nnns 3 13 A ATES Or AWAS SIN BM Zakaria aa P 4 14 4 1 Mat GUI sata eros olaaa EN NK EE 4 15 4 1 1 hora RAV Ue C VOI asasine ikan CO a ada ngga daka anang e ban
27. Wave Calculator that expects as x input the period frequency or wave length of a HAWASSI Calculator harmonic wave and the depth and will then calculate all wave relevant quantities by also specifying the amplitude the calculated steepness is added The calculator can be accessed by clicking Tool gt Wave Calculator INPUT Period s Figure 4 1 Wave Calculator Depth h m Amplitude a m 0 0503 125 m Relative wave length lambda h 2 5 k h 2252 Intermediate depth 0 0503 4 14 I Page ad HAWASSI 4 1 Main GUI Choosing the wave model characteristics wave parameters the domain and input signal and profile are managed in the main GUI Figure 4 2 Details of each input in the main GUI will be explained in detail in the next sub sections File Edit Model Setting Log Tool Help Simulation Boundary Bathymetry Constant v 1 1 Show bathymetry cuernos WAROM 1000 Preparation 400 500 600 RUN STOP HAWASSI VBM2 Post Processing Figure 4 2 The main GUI of HAWASSI VBM2 4 1 1 Working Directory To get started with HAWASSI VBM2 the user should first specify the working directory for the software This can be done by clicking m button and choose a location or folder Working Directory la Figure 4 3 Choosing working directory 4 15 Page dll HAW ASSI In the specified directory folder the software will create new fo
28. aa aa aa and ak ana a a aa baa ba apak 5 58 5 2 Test Case 1 Wave Propagation above a Flat bottom ooooncccccccncccccnononononnnonnnnnnnnnncnnnnnnnnnnnos 5 62 5 3 Test ase A on DE desee Va dni shed In 5 65 A MERI 6 Css ce Fa rre Rm 5 68 5 5 Test Case 4 Simple harbour layOuUt occcccccccccccnnnooonnnononnnnnnnnnnnnnnnnnnnnnononononannnnnnnnnnnnnnnnnnnos 5 71 2 0 Test Case Simple harbour lay Out 2 scscscacoudsnoccensaceccontsssssavacacaraed chee senencosteadsuseianacaiasandeaoanenacss 5 74 5 7 Test Case 6 Cilacap DarDOUE saananira ngana nana anapa na ga vase ko agan o AED pa aana ga ag aai 5 77 3l Breakwater CONSUMO aa angan trade sees a 5 77 laz Breakwate alero 018424040214 0 01390 ORAR PARO RE Oo UU o A 5 80 MEME Dioses 6 83 o E E 6 84 l2lPage dll HAW ASSI List of Figures Figure 2 1 Dispersion relation comparison of Boussinesq type models by Madsen amp S rensen 1992 Nwogu 1993 VBM with parabolic profile by Klopman et al 2010 with the Airy linear theory dispersion relation left as function of kh In the right plot a comparison of dispersion relations for VBM with 1 2 3 Airy profiles and parabolic profile with the Airy linear theory sssssssssseeeeseeeeeeeees 2 12 Beet Fe Waneh ISI EIE EE LT 4 4 Figure 4 2 The main GUI of HAWASSI VBMD sseseeeeeeeeeeeeee enne eene 4 15 Figure 4 3 Choosing working ireCtOry cccceeeeeeeeeseeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
29. able 4 4 below Table 4 4 Bathymetry data format loo pio Note Bathymetry information should cover the whole area inside the coastline 4 1 8 3 Boundary Conditions Two types of boundary that are used in the software are damping zone s and partial reflective boundaries PRB see Figure 4 27 and Figure 4 29 The purpose of the damping zone sponge layer is to absorb outgoing waves without reflecting or disturbing the waves in the rest of the computation domain The efficiency of the damping zone is determined by the length of the damping and its damping The length of damping zone should be at least 2 times the simulated peak wavelength a shorter length may create reflection Damping Zone Ri Create damping zone me wi Select 4 points to represent the area damping area pl pZ p3 p4 pl p2 p3 amp pd are coordinates pi xi WI pl amp p2 form a starting line of damping zone p3 amp p4 form a ending line of the damping z delete damp Figure 4 27 Damping Zone controller After the boundary of the domain has been made a damping zone can be created in the computation domain 1 Click Damping Zone button in the Physical properties panel Figure 4 27 2 Put four points that cover the area of the damping zone The first two points create a line that indicate the beginning of the damping zone area dash line Figure 4 28 4 3 Page dll HAWASSI sence i fiR 4 d
30. amp 1 delete damp pee td MG L o ek o gt gt 2 LE LES g Cross sectional view Damping area Width of damping Figure 4 28 Creating damping zone For editing the boundary point position select Edit Physical properties module Fig 4 24 Two tools for editing points will be activated as shown below move point by specifying the x y value IN y C M move point by dragging using mouse Advice To avoid reflection from the damping zone the width of damping zone should be at least 2 times the peak wavelength Warning Be sure that the influx line is not partly inside a damping zone A modified Sommerfeld boundary condition is used to model partial reflective boundary conditions in HAWASSI VBM2 For this boundary a reflection coefficient Cr e 0 1 has to be specified For Cr 1 the boundary acts as a fully reflecting wall for Cr 0 as a transparent boundary condition so that the out going wave will be fully transmitted Partially reflecting boundaries by setting for instance Cr 0 6 let the wave height of the outgoing wave be reflected for 60 and transmitted for 40 from the initial wave height this applies for all wave frequencies 4 32 Page dll HAWASSI Select an arc to be a Partial Reflective Boundary PRB Choose 2 points counter clockwise Assign a reflection coefficient ref 0 1 delete prb refl 1 is fully reflective boun
31. ct Animation gt 3D view To change the default setting of 3D output select Animation gt 3D view gt 3D setting see Figure 4 47 Animation D view 3D setting Figure 4 46 View the animation in 3D view There are three main panels in the 3D setting interface 1 e Color Scale and Frame Angle The color panel is used to set the color for topography and bathymetry along with the color for wave The data will be spanned in the given color data Set the color by specifying the value of red green and blue component The user can directly specify the value in the RGB box or through the color selection by clicking the button next to the RGB box To execute the color change click OK button Figure 4 48 4 A43 Page ad HAWASSI Eg GUI pp 3dsetting Color R G B Ema 1 00 0 95 0 07 1 00 0 69 0 39 Bathymetry Wave elevation 0 60 0 80 1 00 Wave 0 00 0 60 1 00 Scale Land wave limit Land Wave Frame Angle Azimuth Cancel Figure 4 47 3D setting of animation Using the Scale panel the data can be spanned in the given value range to execute the change click OK For changing the frame angle in 3D view select button O and then move the cursor to have a frame angle according to wish showing the azimuth and elevation value after pressi
32. dary hardwall refl 0 is fully transparent boundar Figure 4 29 Partial Reflective Boundary PRB controller T a 5 reflection coeff 0 1 os Counter Clockwise PRB line Figure 4 30 Creating Partial reflective boundary PRB The following steps are to specify a partial reflective boundary PRB condition at part of the domain boundary 1 Click PRB button in the Physical Properties panel Figure 4 29 2 Select two points at the coastline boundary line that will mark the area of PRB The area is ina counter clockwise direction from the first point until the second point The reflection coefficient is in the range of O to 1 1 for full reflection see Figure 4 30 If the user does not specify a boundary condition at part of the boundary line the software will treat that part as a fully reflecting boundary condition hard wall conditions Warning When PRB of HAWASSI VBM2 v 1 1 is used with Cr 0 fully transmitted not all wave frequencies may be transmitted perfectly short waves kh gt z may be partly reflected less than 10 4 33 Page ad HAWASSI 4 1 9 Mesh Generation HAWASSI VBM2 v 1 1 uses an unstructured triangular mesh The mesh generator is based on the algorithm of Distmesh A simple Mesh Generator in MATLAB that was introduced by Persson and Strang 2004 The mesh generator is designed to create a depth dependent mesh that is finer in shallower areas where the wave
33. dll HAW ASSI Figure 5 27 Domain settings Of test case 3 saanane eee 5 68 Figure 5 28 Settings for mesh properties and model control for the test case 3 sss 5 69 Figure 5 29 Model set p Of test CASE J sasae na aa A a a a a Pa TAA AE AA PA dun a AAS a aa a Ka KANA KN a NARA a a A a NE Pa aa 5 69 Figure 5 30 Test case 3 Snapshot of surface elevation shows the refraction phenomenon 5 70 Figure 5 31 Domain settings of test CAS A anaa era veo nuoc een pe ica 5 71 Figure 5 32 Settings for mesh properties and model control for the test case 4 sssssss 5 72 Figure 3 33 Model SCD Ot test Case RED 5 72 Figure 5 34 Test case 4 Snapshot of surface elevation in a simple harbour layout 5 73 Figure 5 35 Domain settings of test case 5 aane 5 74 Figure 5 36 Settings for mesh properties and model control for the test case 5 sss 5 75 dore 237 Modelspoor 6 SECUS I najana aa a aan PI cate Iac aga pa UU P cgo opa Mua o niodo es dua rap oU aaa pai 5 75 Figure 5 38 Test case 5 Snapshot of surface elevation in a rectangular basin with an internal oblique o cr MR cote 5 76 Figure 5 39 Domain settings of test case 6a Loue esee e esrec oi esset eic oa da sexe a cnt 5 77 Figure 5 40 Settings for mesh properties and model control for the test case 6a
34. e 1 Wave Propagation above a Flat bottom This introductory example of waves above a flat bottom is to illustrate wave influxing and effects of damping zones and partially reflecting boundaries Besides that various types of waves and the difference in propagation when different wave models are used can be investigated A simple domain configuration setting for the mesh generator and model control are illustrated in Figure 5 20 Figure 5 21 and Figure 5 22 File Edit Model Setting Physical properties Influx line Working Directory urs 4 48 I CAUsersDiditiDocumentsW z p delete influx Damping Zone damp 2 v Project Name TC1 IntroFlat User s note 600 1000 1200 1400 1600 Wave propagation above a Project Setting file is loaded Figure 5 20 Domain settings of test case 1 5 02 Page E mesh properties 200 dll HAWASSI E model control Initial conditions Zero Amplitude m Depth dependance Parameters Approx Tp s Width hump x y m ppl Geometry max de Im shallowest depth m Time Signal Direction deg Tp s Amplitude m Figure 5 21 Settings for mesh properties and model control for the test case 1 0 signal 400 1 300 200 1 bathymetry amp numerical setup ADO E 100 n 200 100 0 sim Oa 0 500 1000 x m Amp Spec of 1D Signal 100 10
35. earity Y 1 prof optimized 2 prof optimized 1 prof user defined 2 prof user defined Figure 5 4 Select wave model 5 1 3 Defining the Domain Start with defining the domain of computation or in this case the Simulation Boundary Simulation boundary files have been prepared for each test case To load the simulation boundary of Test case 6b click the Simulation Boundary button Figure 5 5 and locate the file location After the data is loaded the simulation boundary will be plotted in the GUI Figure 5 6 Simulation boundary data includes a list of two dimensional polylines that is always considered to be closed Simulation Boundary Load simulation boundary format data coastline xy in ascivbxt Figure 5 5 Load Simulation boundary data 5 50 Page dll HAWASSI 1500 1000 4 an 500 7 lt Ll a i bad m OF j 0 500 0 500 1000 Figure 5 6 Simulation boundary data is loaded 5 1 4 Reading bathymetry file After the simulation boundary is loaded the next step 1s to load the bathymetry data To that end click the Bathymetry pop up menu and select Scatter bathymetry data x y z see Figure 5 7 Locate and select the bathymetry data Simulation Boundary Bathymetry Consta iN Show bathymetry Select bathymetry type Constant bathymetry Scattered bathymetry data x y z Figure 5 7 Load bathymetry data 5 1 5 Assigning boundary types
36. ecified Buoy position inside the computation domain see Figure 4 51 For the spectrum the frequency band option is provided to filter the extracted signal information see Figure 4 52 The Snapshot option provides capability for the user to take a snapshot of the simulation at a time to be specified and within an interval to be specified see Figure 4 53 Platting OC locations where i the signal s will be extracted in the domain Te EY Din XL Y1 X2 Y2 XM YN Please provide the ey location EE t init t end Time frame of the signal frame of the signal Plot the signal amp buoy location s Ml Signal Please providemodity the time frame E Amp smother for the si Plot signalis at certain location s Iqnal Spectrum Energy Plot the spectrum tibuoy locationis select which spectrum to be plotted Ls BLA of for smoother function Plot the amplitude spectrum forthe spectrum Plot the energy spectrum default 3 Figure 4 51 Buoy setting option of Plotting f init end frequency band Buoy s CX Y D Filter the signalis Signal at the buoylsi location select a frequency band tin Hz Leave it blank to take all frec Spectrum Figure 4 52 Option for filtering the extracted signal at specific Buoy location s 4 46 Page ad HAWASSI Xvvest Xeast south north REM at N S zn Plot snapshot of simulation Time af the snapshot ix vj frame of the s
37. ed users a trial simulation with VBM parabolic profile will give an idea how the dispersive waves will propagate If the model s dispersion quality is too poor for the waves to be simulated the Optimized VBM with 1 profile and after that with 2 profiles can be tried Warning 1 The more Airy profiles are used the more computation time will be needed For rather simple waves simulations with 2 profiles will approximately cost 1 5 times more CPU time than needed for 1 Airy profile For complicated waves very steep broad spectrum the calculation time may be longer 4 l6lPage ad HAWASSI 2 The User defined option of VBM is recommended only for advanced users who know and understand the energy spectrum of the wave to be simulated and the effect of their choices 1 e how the choice of the frequencies in the Airy functions of VBM may affect the results Setting Help log Tool Type SWE Model control Ctrl M VBM Parabolic Profile Nonlinearity 1 prof optimized 2 prof optimized 1 prof user defined 2 prof user defined Figure 4 5 Choosing model type Airy user defined value Hz freg1 freg2 Figure 4 6 Pop up menu for Airy user defined value Hz 0 07 0 14 The option to include nonlinear terms in the mathematical model is given in the option under Model gt Type gt Nonlinearity gt weakly nonlinear see Figure 4 7 In the current version of HAWASSI VBM2 the weakly nonli
38. eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 15 Figure 44 Projectname amp USet S DOME tania 4 16 Disure 4 5 Choosing model typo di di iia 4 17 Figure 4 6 Pop up menu for Airy user defined value H2z eaaa nana aana 4 17 Figure 4 7 Specification of Linear or Nonlinear simulation sssseeeeeeeeee eaaa eea 4 18 Fisure des Model Control sola aaa 4 19 Figure 4 9 Initial condition for wave cleVatiOn cccccecceeeeseeseseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 20 Figure 4 10 Initial condition for single gaussian hump ooccccnnnnnnnncnnnnnononnnononnnnnnnnnnnnnnnnnnnnonnnnnnonnnnnnnnnos 4 21 Drsure 4 11 Single CASIO WW tr icon 4 21 Figure 4 12 AS 4 22 Ier rures Sd ls Bipolar MU Mie sana aa saraga aa a naa bana a Ag Gan aaa kaa a Nga aga agan aa aa a daa na aa AA A ja a aa naa aa aa aaa aa 4 22 Figure 4 14 Influx Ine controlleren BB AK AA EA NAN Ep AAN EE NG BB DAB Aa TONEN Aa aa aa ki 4 23 KD A AAA aana ng a D ag ats San ag aaa aaa agan E aaa a gan aa aaa ag aa dag aaa aa 4 23 Figure 4 16 Influx signal option cccccccccccccccccsssssessssesseeeeccececeeeeeeaaaaesssscseseeeeeeeeeeeeeseesaauaesesssseeeeeeeeeess 4 24 Figure 4 17 Signal parameters for harmonic signal oooonnnnnccnnnnnnncnnnnnnonononnnnnnnnnnnnnnnnnnnononononnnnnnnnnnnnnnss 4 24 Figure 4 18 Signal parameters for Jonswap option sessssseseeeeeeeeenenene nennen nenne nnn nnn nnns 4 25 Figure 4 19 Location of User defi
39. en nnns 4 35 Eretire 4 34 Calne Advanced Settings GU Lausana cas 4 36 Figure 4 35 Advanced Setiimes GU La eo 4 37 Figure 4 36 Time stepping parameters in the Advanced Settings sese 4 37 Figure 4 37 Buttons for running the simulation in the main GUI eese 4 38 Figure 4 38 Overview of the numerical setup panel showing details of the simulation to be performed4 39 l3lPage dll HAW ASSI Figure 4 39 Description of Dispersion Quality plot blue solid for Exact red dashed for Model set up 4 39 Figure 4 40 Post Processing GUI essssseeeessssesesssssseeeeeeenn nnn nnn nnnnnn nenne sens sensns ni i ren rn nnn nnne nnne 4 4 Figure 4 41 Load other simulation data in Post Processing GUI o oooooocoooonccccccnnnnnnnnnnnnnnononononononnnnnnnos 4 42 Figure 4 42 Domain for showing the animation ccccccccccccccccccceeaeesssseeeeeeccccceeeceseeesaaasssesseeseeeeeeeess 4 42 Figure 4 45 Animation MMe li ita 4 42 Figure 4 44 Saving animation as MEG deser e itin iii 4 43 Figure 4 45 Saving all frames in png pg format sesssssseeeeeeeneenenen nennen nennen nnn nnn nnne nennen 4 43 Figure 4 46 View the animation in 3D VICW cccceeeeeeeeeeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 43 Figure 4 47 3D setting of animation ssssseseeeesesseeseeeeeeeeeeeee nenne nnn hh nn nnn nennen sensns nenen ren rn nnn nnne enne 4 4
40. er by providing the Hamiltonian form of the dynamic equations HAWASSI 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 7 Page 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 agree that Yayasan AB is not liable for any loss or damage arising out of the use of the Softwa
41. hymetry Influx line Working Directory Xx 1 C Users DiditiDocuments M m mode OVBM 1 profil Nonlinearity No Initial cond Zero Influx 1 Harmonic Damping Zone Time 0 1 400 Project Name TC3 Reffraction 3000 User s note 500 1000 1500 2000 2500 Reffraction Phenomenon T Project Setting file is loaded Post Processing Figure 5 27 Domain settings of test case 3 5 68 Page y m ad HAWASSI EY mesh properties Yee EOS tel conditions Zero Parameters E H Approx Tp s ppl max ds m Time Signal shallowest depth m lec WA TEH Direction deg Tp 5 Amplitude m Figure 5 28 Settings for mesh properties and model control for the test case 3 2D signal 1000 0 5 bathymetry amp numerical setup o E 500 i 0 5 1 gt 1000 1000 lt ET 400 500 200 sim 00 1000 2000 3000 500 b x m Amp Spec of 1D Signal 0 00 0 13 0 18 0 22 0 26 0 2 e S P 0 15 E Son TET t 0 500 1000 1500 2000 2500 3000 E E x m 0 1 nk ki O 0 05 ed Salone D Dispersion Quality 0 of 02 0 2 4 6 8 f Hz V Figure 5 29 Model setup of test case 3 5 69 Page t 00 05 59 hr min sec or t 359 90 sec Figure 5 30 Test case 3 Snapshot of surface elevation shows the refraction phenomenon Observation o Observe the amplitude changes above slope o Reflections are more pr
42. i ajak naa A 4 15 4 1 2 Project INI de Users INODG akas sasa aba San aga ma aa a aa gana Wana b a ang estu Ei niens snc pobladas tiara 4 16 4 1 3 Wave Model Characters Susi ii oi 4 16 4 1 4 A Um m 4 18 4 1 5 Tarte dE panacea 4 19 4 1 6 TAG T e T a panasaran aaa aa oec 4 23 4 1 7 DUA TON MING APA OPA E o Aa AA a 09 890 9609 090 800 00 8 800 0 0 05099805 4 29 4 1 8 Simulation Boundary Bathymetry and Boundary Conditions 4 29 4 9 Nes CIE PN DT nie aan ane saa E aa ENG agek a ai pesada dn cum UE Dad UNE 4 34 O 234 0 SO SENIN aia kaag ag aha eed tbi Strobes saeeeeseabauentatagcasesntobes eenasdsdangacatacciessesovonsaeens 4 36 4LIL FAMINE SIAN cuna ico 4 37 AI a PEO A 4 40 da O AAA 4 4 4 2 1 AAAA O ET KR BEN 4 42 4 2 2 A a 4 46 A 5 48 5 Example of set up simulation and post processing of harbour waves eesssese 5 48 5 1 1 Gettine started the software sacos aana nag ak gana oir ca 5 48 l 1lPage dll HAWASSI 5 1 2 Selec ne wave model tail 5 50 5 1 3 Delfino tbe DOM osa A aaa EEN Ga AG TENGE oa AR ADO NE oirlo 5 50 5 1 4 Reading bath ymetty eucariotas praia 5 51 5 1 5 Assigning boundary CY PCS iria ET 5 51 5 1 6 Specs UE III EE 5 54 5 1 7 ee NS II E aaa aa 5 56 5 1 8 Running the simulation corales ka kaga E nayana ga kani gada aia aa 5 57 5 1 9 Post processes OL fhe OUI sseni aa aa aa ban kaa TE
43. ial conditions Single Gaussian hump Bipolar hump 1 Bipolar hump z Cancel Figure 4 9 Initial condition for wave elevation Table 4 1 Data format for user defined initial conditions Foo pro Other options of initial conditions are Single Gaussian hump bipolar hump 1 and bipolar hump 2 for these options additional parameters have to be specified such as the amplitude the width the center location of the initial profile and the rotation angle see Figure 4 10 Figure 4 13 4 20 Page ad HAWASSI E model contral Initial conditions Single Gaussian hump Amplitude m The amplitude m af initial profile the width m of initial profile in x amp y direction Width hump xyi m j HH Location x y m NE location m of initial profile in x amp y direction the rotation angle degree of initial profile In Mautical convention clackwise from Morth Rotation angle deg s Figure 4 10 Initial condition for single gaussian hump WIDTH Y Y LOCATION Y LOCATION X Figure 4 11 Single Gaussian hump 4 21 Page dll HAWASSI AMPLITUDE rr WIDTH LOCATION X Figure 4 12 Bipolar hump 1 AMPLITUDE LOCATION Y LOCATION X Figure 4 13 Bipolar hump 2 4 22 Page ad HAWASSI 4 1 6 Time Signal In HAWASSI VBM2 when dealing with a signaling problem an embedded wave influxi
44. ime Signal panel then set the following signal parameters Direction 270 nautical counterclockwise from north Tp 9 Amplitude 0 25 Time gt 0 1 300 5 55 Page ad HAWASSI EY model contral Time Signal Harmanic Influx 1 Direction deg 270 Tp s 9 Amplitude m Spreading s T n eta Filter f_0 f end Hz Figure 5 13 Set influx signal parameter in Model Control 5 Click the OK button to close Model Control GUI 5 1 7 Creating mesh Now the set up of the domain boundary conditions and wave influx is done the next step 1s to create the mesh by clicking Generate Mesh button After the button is pressed a pop up GUI Mesh properties will ask for parameters for the mesh generator The parameters should be set with the following information e Approx Tp 9 e ppl 10 e max ds 200 e shallowest depth 1 e check Depth dependence 5 56 Page dll HAW ASSI The mesh generation process will start after the Save button is clicked A progress bar indicating the progress of the unstructured grid generation is shown in the lower part of the main GUI Error Reference source not found with a result as shown In Figure 5 18 1400 1200 1000 B B00 400 400 200 0 200 LIE ZEN a Figure 5 14 Generated mesh 5 1 8 Running the simulation After the mesh has been created and model data have been specified the simulation with HAW ASSI VBM
45. iple for the boundary value problem in the fluid domain By restricting the set of competing functions in the minimization an approximation of K 7 is obtained The variational derivative 04K y On is the corresponding consistent approximation of the Dirichlet to Neumann operator e The approximate Hamilton 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 In VBM the interior flow is approximated by using a linear combination of vertical Airy profiles characterized by the values of wave numbers xm The choice of these values determines the dispersion 2 JOlPage dll HAW ASSI relation Given the spectrum of the influx signal or initial profile of the case under investigation the values Of Km are optimized for best performance for the relevant frequency interval Therefore VBM can have excellent tailor made dispersive properties deep water waves can be simulated just as well as infragravity waves Numerical Implementation A Finite Element method using piece wise linear splines can deal with the first order differentiations that appear in the approximate Kinetic Energy In addition to two scalar dynamic equations for and y a system of elliptic equations has to be solved for the amplitudes of the Airy functions Advice
46. l HAWASSI 2 3 Relation to other Boussinesq type wave models HAWASSI VBM is a phase resolving model where individual waves wave components in the energy spectrum are resolved with their phases and amplitudes This type of model is typically used for studying wave propagation in a small area such as in a harbor or near coasts Other Boussinesq type models that are adopted by commercial software are based on Madsen amp Sgrensen 1992 MIKE21 Boussinesq Wave by DHI and Nwogu 1993 BOUSS2D SMS by Aquaveo Both models have dispersion accuracy up to kh 3 14 which implies that it is possible to simulate waves with wavelength more than 2 times the depth The dispersive quality of HAWASSI VBM compared with these other two Boussinesq models is illustrated in Figure 2 1 N k y o g Airy YBM 1 profile 5 VBM 2 profiles 7 VBM 3 profiles VBM parabolic profile 6 2 4 6 8 10 12 i 5 10 15 20 25 30 kh kh Airy VBM parabolic profile t Madsen amp Sorensen 92 Nwogu 93 a 0 39 Figure 2 1 Dispersion relation comparison of Boussinesq type models by Madsen amp Sgrensen 1992 Nwogu 1993 VBM with parabolic profile by Klopman et al 2010 with the Airy linear theory dispersion relation left as function of kh In the right plot a comparison of dispersion relations for VBM with 1 2 3 Airy profiles and parabolic profile with the Airy linear theory 2 4 Units coordinate
47. lder Output when the directory does not contain the folder yet If the folder already exists it will keep and use the folder Al output files will be stored in the Output folder 4 1 2 Project Name amp User s Note After specifying the working directory the user should specify Project Name and if wanted User s note The User note will be printed in the log file which is stored as Output LOG ProjectName log Project Name User s note Figure 4 4 Project name amp User s note 4 1 3 Wave Model Characteristics HAWASSI VBM2 comes with two main versions the non dispersive model Shallow Water Equation SWE and various variants of the dispersive model the VBM There are 5 sub options of VBM as choices which vertical profile s to be used which will determine the dispersive quality of the model as illustrated in Figure 2 1 The user can select the model to be used under Model 2 Type 2 SWE or VBM see Figure 4 5 In the Optimized VBM the software will calculate automatically the optimized wave number x to be used in the Airy profile based on the problem to be solved signaling problem and or initial value problem In the User defined VBM the user can freely choose frequencies related to the wave numbers to be used as input in the Airy profile When this option 1s selected there will be a pop up menu asking for the custom values in Hz as shown in Figure 4 6 Suggestion For inexperienc
48. lementation of the Variational Boussinesq Model VBM Presently the code is for simulation of wave structure interactions coupled wave ship interaction is foreseen in future releases VBM is a Boussinesq type model first introduced by Klopman et al 2005 that is derived via the variational formulation for surface water waves The model has been further developed to have tailor made dispersion properties based on the problem to be solved with accuracy up to kh 15 or more see Adytia amp Groesen 2012 The interior fluid motion is modelled by a combination of a few Airy type depth profiles this makes it possible to optimize the dispersion properties depending on the specific case to be simulated Nonlinear effects are accounted for in a weakly nonlinear way that is sufficient for most applications The model is called the Optimized Variational Boussinesq Model OVBM which is the mathematical model behind the HAWASSI VBM Underlying Modeling Methods HAWASSI V BM is based on the following principles e The free surface dynamics of the irrotational flow of inviscid incompressible fluid is governed by a set of Hamilton equations for the surface elevation 7 and the potential at the surface e By approximating the kinetic energy functional K 7 explicitly as an expression in y and the simulation of the interior flow can be avoided the Boussinesq character of the codes e The way of approximating K 7 is based on Dirichlet s princ
49. length becomes shorter such that roughly the same number of points per wave length are taken in the whole domain An illustration of a depth dependent mesh is shown in Figure 4 31 2000 2000 1500 1500 1000 1000 500 1000 1500 2000 2500 1000 1500 2000 2500 Figure 4 31 Illustration of depth dependence mesh Colorbar indicates the water depth m After the user has designed loaded simulation boundary influx line and bathymetry the mesh can be created by clicking Generate Mesh button in the main GUI see Figure 4 32 A pop up dialog box so called mesh properties will appear for asking mesh parameters as shown in Figure 4 33 Generate Mesh Generate unstructured triangular grid Prepa based on domain configuration paints in the influx line will become grid ELI SIOP Figure 4 32 Generate mesh 4 34lPage al HAWASSI E mesh properties Depth dependance Parameters Approx Tp s ppl Geometry max ds m shallowest depth m Figure 4 33 Mesh properties dialog box The depth dependent mesh is activated when the check box Depth dependence is checked otherwise the mesh size will only be determined by max ds value maximum allowable grid size under the Geometry panel Shallowest depth is the value of the shallowest depth in the bathymetry As an illustration when the user loads a bathymetry data that contains topography land the depth is negative value
50. ling problem as a txt ASCIT file in the folder Output ProjectName the option is checked by default To save all information that have been configured in the Advanced Settings GUI click Save to cancel all changes click HE 4 36 Page ad HAWASSI EY Advanced Settings Parameters g m s 2 p kg m 3 Time Stepping Time partition 4 Hotstart E Split Output Data Save influx signal A Save surface potential phi Figure 4 35 Advanced Settings GUI E Advanced Settings Time Stepping Divide the simulation calculation into several steps Put a natural number Time partition Save the last calculation eta phi 7 Fl spit o Et end ee Y Hotstart split Output lids tie E aa Split the output simulation data AZ into the number of time partition Output data will be DATA ProjectName iter timepartition mat A Save signal influx F Save surface potential phi Save surface potential phi after the simulation Figure 4 36 Time stepping parameters in the Advanced Settings 4 1 11 Running Simulation After all parameters and input have been set the user can start the simulation Three main steps can be distinguished 1 Preparation After the Preparation button is pressed the software extracts all input data from the Main GUI the Model Control GUI and the Advanced Settings GUI The input data will be checked and processed An overview of the i
51. nak Krakatau explosion in the Jakarta Bay Proceedings Basic Science International Conference 2012 Malang Indonesia K1 5 ISBN 978 979 25 6033 6 D Adytia M Ramdhani amp E van Groesen Phase resolved and averaged Wave Simulations in Jakarta Harbour Proceedings 6th Asia Pacific Workshop on Marine Hydrodynamics APHydro2012 Johor Baru Malaysia 3 4 September 2012 pp 218 223 Ivan Lakhturov Optimization of Variational Boussinesq Models PhD Thesis UTwente 9 November 2012 Didit Adytia Coastal zone simulations with Variational Boussinesq Modelling PhD Thesis UTwente 24 May 2012 I Lakhturov D Adytia amp E van Groesen Optimized Variational 1D Boussineq modelling for broad band waves over flat bottom Wave Motion 49 2012 309 322 D Adytia and E van Groesen Optimized Variational 1D Boussinesq modelling of coastal waves propagating over a slope Journal Coastal Engineering 64 2012 pp 139 150 D Adytia and E van Groesen The variational 2D Boussinesq model for wave propagation over a shoal International Conference on Developments in Marine CFD 18 19 November 2011 Chennai India RINA ISBN No 978 1 905040 92 6 p 25 29 G Klopman Variational Boussinesq modelling of surface gravity waves over bathymetry PhD Thesis UTwente 27 May 2010 Gert Klopman Brenny van Groesen Maarten W Dingemans A variational approach to Boussinesq modelling of fully non linear water waves Journal Fluid Mechanics 657 2010 36 63 D Adytia
52. napshot tin m Provide time of the snapshot pe te oss n s Figure 4 53 Snapshot setting option of Plotting 4 47 Page dll HAW ASSI 5 Test cases HAWASSI VBM2 provides 6 test cases The first case is meant for practicing the software such that the user gets a first idea how the software works in handling influx signals boundary conditions and different wave models SWE 4 VBM The second and third test cases are to show the characteristic wave phenomena of diffraction refraction and shoaling Understanding such wave phenomena is important as we go further into test cases 4 6 using a harbor layout with real bathymetry where all wave phenomena are expected to happen in a single domain For getting started with the practical use of the software we take the last test case of a realistic harbour as example after that the test cases are described in successive sections 5 1 Example of set up simulation and post processing of harbour waves This section will show how to prepare a domain generate a mesh and run the simulation in HAWASSI VBM2 software for the test case of a realistic harbour in Indonesia test case 6b This section can be read without knowing the information of Section 4 to which the reader can refer for further details especially for knowing how to prepare input data that are provided for this and other test cases 5 1 1 Getting started the software To start the software click the shortcut so called VBM
53. nd 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 an approximative way this was formulated accurately in the 1960 1970 s by Zakharov and Bro
54. ne Mone Harmonic Jansweap User defined SWAN 20 Spectrum Figure 4 16 Influx signal option For the option Harmonic the user needs to specify the main wave direction in Nautical convention clockwise from North direction the wave period Tp and the amplitude Amp see Figure 4 17 E model contral Time 5ignal Harmonic Influx 1 Direction deg 90 m main Wave direction to degree s in Nautical convention clockwise from North Amplitude m 1 Filter fF O f end Hz filter the influx signal if yes then specify a frequency band f f low f high if nat leave it blank Figure 4 17 Signal parameters for harmonic signal 4 24 Page dll HAWASSI For the option Jonswap the user has to specify the main wave direction in Nautical convention clockwise from North direction the peak wave period Tp the significant wave height Hs Jonswap parameters such as gamma and directional spreading s see Figure 4 18 E model contral Time 5ignal Jonswap Influx 1 main Wave direction to degree Direction deg in Nautical convention clockwise from North Tp s 10 b x neve period 5 Hs m 1 La Significant wave heig Y 13 bs JONSWAP paramater gamma i default value 3 3 Spreading 5 inf Ls ns cos 2s theta Input s in the directional speading function cos 2s theta theta is wave di
55. near VBM is used see Adytia amp Van Groesen 2012 for further details For relevant applications that are given in Test Cases the weakly nonlinear version is sufficiently accurate to describe wave phenomena with a good match when compared with experimental data When using the nonlinear simulation the length of a nonlinear adjustment zone has to be prescribed in number of wave lengths in order to prevent spurious oscillations see Fig 4 7 Suggestion It is suggested to do first a linear simulation before any nonlinear simulation Besides the fact that a linear simulation is faster than a nonlinear one more important 1s that the linear results represents in many cases already 80 90 of the wave characteristics to be simulated except for cases with very strong nonlinearities Having studied the output of the linear simulation the differences caused by nonlinear effects can be better investigated 4 17 Page ad HAWASSI Warning The nonlinear calculation will cost at least 2 times the CPU time of the linear calculation Just as for the choice of the number of Airy profiles the CPU time will vary depending on the complexity of the waves to be simulated setting Log Tool Help Type SWE Model control Ctrl M VBM d Nonlinearity v Linear weakly nonlinear Figure 4 7 Specification of Linear or Nonlinear simulation 4 1 4 Model Control An independent GUI so called Model Control is provided to specify Initial
56. ned influx signal ccccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 26 Figure 4 20 Signal SWAN 2D Spe Cid sso gag aaa NENGA r ipro tes Ka AANE ita SANE DB EE DA Da Pa sio sd beso aded aaa 4 27 Petre dp Time Tater val For SOBE SE osse cessi baee pou ouod Sob oa FPE PR Boo ep 20d BE Nan KERA donostia 4 28 Figure 4 22 Nonlinear adjustment for nonlinear modetl sese 4 28 Fisure 4 23 Time control Of simulation rra soberana 4 29 Figure 4 24 Load simulation boundary data eaaa nn 4 29 Figure 4 25 Tools for editing creating the boundary points sese 4 30 Figure 4 26 Simulation boundary and Bathymetry option cccccccnccccnnnnononnnononnnnnnnnnnnnnnnnnnnnononononaninnnnnnnos 4 30 Fig re 4 27 Damping Zone controller aas asarana kesana aba doc iet ANGANGGE KEANGGE xd eor EE P dI PEDES Aou E06 Ba NE RC Ax 4 3 I 15016425 Creating damp ie ZONE eise mae aa aaa an aa a ak da a aa aa a aaa aa Bag E Tama IUE O IU M DD Rae aana baa aa ka A 4 32 Figure 4 29 Partial Reflective Boundary PRB controller sese 4 33 Figure 4 30 Creating Partial reflective boundary PRB esses eene 4 33 Figure 4 31 Illustration of depth dependence mesh Colorbar indicates the water depth m 4 34 Poua PA CILE SEI i E aaa aa 4 34 Figure 4 33 Mesh properties dialog box esses nennen nennen nennen nnn nnn nn nenn
57. ng method is used see Liam et al 2014 Instead of generating influx from a boundary an influx signal is generated in the computational domain from a line in the interior of the domain Provided that the boundary of the domain is already made the following steps will create an influx line l Click Influx line button under the Physical properties panel Figure 4 14 2 Put two points as the ends of the influx line Figure 4 14 J Specify the distance between points ds Figure 4 15 Physical properties Influx line M Create influx lines influx Z 1 oo ds delete influx x Figure 4 14 Influx line controller To have a point on the shoreline boundary line Put the starting point outside the domain the point will be projected automatically on to nearest the shoreline ree The influx line will automatically be created after the distance between points ds 1s specified Ar ello on sd o nme o recen eee A A SS E a a aana Figure 4 15 Creating influx line 4 23lPage dll HAWASSI In the Model Control GUI a time signal dialog will be active if there is an influx line in the domain after the mesh is generated There are four type of influx signal that can be chosen by the user a harmonic signal a signal with JONSWAP type of spectra a user defined signal and a signal with 2D Spectrum input 2D SWAN spectrum format see Figure 4 16 E model control Time Signal Mo
58. ng the Apply button to save all setting changes Fig 4 48 4 44 Page ad HAWASSI 1 00 0 69 0 39 0560 0 80 1 00 NODO TT el TILT LL ajajajaje jaja ajajajaj E ajajaja 11 ODIO ajajaja EN ajajaja EN Figure 4 48 Setting land and wave color in 3D setting Z Viewpoint Center of Plot Box y Figure 4 49 Viewpoint with azimuth and elevation 4 45lPage dll HAW ASSI 4 2 2 Plotting The Plotting panel can produce plots of bathymetry bottom profile the significant wave height Maximum Temporal Amplitude MTA maximum wave height Hamiltonian or total wave energy and wave disturbances normalized significant wave height Hs with respect to the Hs in the influx location see Figure 4 50 a FORO Plot the bathymetry A Bathymetry Hamiltonian Enercy Plot the Significant Wave Height Hs m as functian af time CE Plot the Hamiltonian Enerc Significant Wave Height Wave Disturbances Plot the Wave Disturbances DS Wave disturbances is normalized significant Si Plot the Maximum Temporal IV wave height Hs with respect to Hs Amplitude MTA m Y ie Figure 4 50 Bathymetry significant wave height Hs MTA Hamiltonian Energy and Wave disturbances Data at a specific time or at a specific location can be obtained by the option Buoys and or Snapshot Buoy will give the time signal and or its spectrum of the elevation at the sp
59. ng the installer and following the instruction in the installation wizard 3 3 Second step Installing HAWASSI VBM After the installation of the MCR is done the installation of HAWASSI VBMI and 2 is performed by double clicking the installer of HAWASSI VBM setup HAWASSI VBM vl l exe and follow 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 HAWASSI VBM VBM1 or VBM2 from the shortcut on the Desktop In the Main GUI that appears under Help go to Licence 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 After the installation is finished the software can be accessed from the shortcut in the Desktop and the Start Menu Documentation and Test cases of the HAWASSI VBM are provided and placed in My Document HAWASSI_VBM 3 l3lPage ad HAWASSI 4 GUTs of HAWASSI VBM2 For ease of operation HAWASSI VBM2 software includes two interfaces or GUI s Graphical User Interfaces namely the Main GUI for the input model and a Post Processing GUI for the wave simulation output The GUT s act as an input output manager There is a simple
60. north X Axiz Lim X Y I space domain m for showing the animation Figure 4 42 Domain for showing the animation A restricted time interval for the animation can also be specified and the animation interval steps a time step that is a multiple m of dt the value m to be given under tdt Figure 4 43 Animation tinit tend zat show the animation every Time what number af dt Ds Ls select a natural number Starting and Ending time 3 forthe animation Figure 4 43 Animation time The animation can be saved as a moving GIF file format by specifying a delay between each frame and the looping of the animation Figure 4 44 the delay information is in seconds and the loop information has to be a natural number 4 42 Page ad HAWASSI For making animations in other formats saving frames either in PNG or JPG format can be used afterwards with any existing animation software Figure 4 45 If no specific name is provided in the animation ID field the animation frames will be saved with anim as the default name Animation gk The saved gif file will be run save as GIF delay loop Da for how many loop choose 0 1 2 delay time at each frame during the making of gif file n s Figure 4 44 Saving animation as GIF iW Save frames format ana pna Figure 4 45 Saving all frames in png jpg format In HAWASSI VBM2 the animation result can also be viewed in 3D Sele
61. nput and numerical setup will be shown to the user via a panel plot called Model Setup The panel shows an overview of domain bathymetry boundary conditions 4 37 Page ad HAWASSI influx location the influx signal and its amplitude spectrum the initial condition of the surface elevation n and the dispersion quality of the model compared to the exact Airy linear theory see Figure 4 38 The dispersion quality plot is as in Figure 4 39 RUN The RUN button is enabled after the Preparation phase is finished RUN will execute the time integration from the numerical setup When the STOP button placed under the RUN button is pressed the simulation will terminate and the simulation results until the termination time will be stored If the calculation is finished the software will automatically call the Post Processing GUI and load the simulation data directly to it Post Processing When the Post Processing button is pressed the Post Processing GUI will appear Details of the GUI will be described in Section 0 Generate unstructured triangular grid based on domain configuration points in the influx line will become grid poi Process the model setup all input data will be read Run the timestepping Figure 4 37 Buttons for running the simulation in the main GUI 4 38lPage al HAWASSI Signal influx Initial condition 2 influx location for surface elevation 7 E Figure 1 Model Setup DTD A 2o
62. onounced for steep slopes o Investigate effect of different wave lengths 5 JOlPage ad HAWASSI 5 5 Test Case 4 Simple harbour layout Test case 4 is designed to illustrate standing resonant harbour waves in a rectangular basin with influx from part of the boundary The influxed waves are mainly reflected from the opposite wall and leave the domain not affecting or hindered by the influx line but diffraction causes the waves to partly enter the right part of the basin where for specific wave influxes a pattern may develop standing waves The domain configuration settings for the mesh generator and model control are illustrated in Figure 5 31 Figure 5 32 and Figure 5 33 File Edit Model Setting Log Tool Help Influx line Working Directory A NGA M y Im C Users DiditiDocuments M delete influx Nonlinearity No Initial cond Zero Influx 1 Harmonic Time 0 1 200 TC4_SimpleHarbour User s note 1000 1200 Wave entering a rectangle a Figure 5 31 Domain settings of test case 4 5 71 Page E mesh properties Depth dependance Parameters Approx Tp s Geometry shallowest depth m ad HAWASSI Initial conditions Zero Amplitude m Width hump x y m Location x v m Rotation angle deg Time 5ignal Direction deg Tp s Amplitude m Spreading s cos 25 theta Filter f_
63. processing GUI Keep the values for time and space at their default values For running the 2D animation click the RUN button directly For running the 3D animation check the 3D view and if desired adjust the 3D view setting by clicking 3D setting button To start the animation click the RUN button 5 58 Page ad HAWASSI Post Processing 2D E rocessing v 1 1 DATA TCBb Realish Bathymetry _ Hamiitonian Energy Significant Wave Height Wave Disturbances Animation MTA tinit tend Fat Tie Chi Len Xwest Xeast Ysouth Ynorth frequency band 0 0 99667 338 1603 755 371 1459 3312 Buoy s QXY tint t end save as GIF delay 0 1 loop inf Signal 0 300 Amp A r save frames format pn Spectrum a Jl jia Energy C 30 view 3D setting Xwest Xeast Ysouth Ynorth C Snapshot Gt 338 1603 755 371 1458 save plot s Done variable stored Figure 5 16 Post processing GUI 5 59 Page ad HAWASSI B Gut pp_3dsetting Wave elevation 0 60 0 80 1 00 0 00 0 60 1 00 Scale Land wave limit Land Min Max Wave Min Max Figure 5 18 2D Animation output 5 60 Page ad HAWASSI t 00 04 19 hr min sec or t 259 87 sec Y a ri uw pet Lo 24 LN E Lua 4 e 4 d E Figure 5 19 3D animation output 561lPage al HAWASSI 5 2 Test Cas
64. rated in Table 4 3 Simulation Boundary Load simulation boundary format data coastline x y in ascivtxt Figure 4 24 Load simulation boundary data 4 29 Page ad HAWASSI zoom to the extent of figure T X move point by specifying the x y value 2 x 5 move point by dragging using mouse Ni add point to the existing polygon simulation boundary e creating polygon right click to close the sequence of points Pa delete point s in the existing polygon simulation boundary a ae break point in the existing polygon ER connect two polygons Clear figure Figure 4 25 Tools for editing creating the boundary points Table 4 3 Coastline data format 4 1 8 2 Bathymetry The software provides two options for choosing bathymetry type flat bottom or custom bathymetry that is defined in a scatter data x y z For flat bottom the depth should be specified in the Depth box see Figure 4 26 Depth Simulation Boundary Bathymetry Consta kad Show bathymetry Select bathymetry type Constant bathymetry Scattered bathymetry data x y z Figure 4 26 Simulation boundary and Bathymetry option 4 30 Page dll HAWASSI For the bathymetry option Scatter data x y z the data input has to be located through a pop up dialog box The data should consist of three columns the first and second column for x and y coordinate the third column for the depth should be gt 0 at x y see T
65. re Although much care 1s 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 l SlPage ad HAWASSI 1 Introduction This document is the Manual of HAWASSI VBM2 software that serves as a guide for using and running the software HAWASSI VBM2 simulates phase resolved waves in 2 Horizontal Directions long and short crested waves as appear in wave tanks coastal and oceanic areas harbors with flat and varying bathymetry and with partially reflecting walls and damping zones Section 2 describes briefly the mathematical background and dispersive qualities of the code together with some 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 The full description of the software regarding GUIs and input output parameters is given in detail in Section 4 Section 5 starts with a short tutorial for getting ready to use the software without the need to know all the information from section 4 in particular to show how to deal with test cases that are provided in successive subsections www hawassi labmath indonesia org DEMO version with restricted functionality
66. rection for long crested wave s inf for short crested wave s positive integer Filter f_0 f end Hz if yes then specify a frequency band f f low f high If not leave it blank Figure 4 18 Signal parameters for Jonswap option For the option User defined the user should locate the data file of the signal influx in ASCII file through a pop up dialog box User defined signal data need to have a format as show in Table 4 2 The first box has to be filled by 0 The other boxes in the first row are the data position with respect to the first point The first column is the time information The column below data position is the signal data at the corresponding point The example of user defined signal can be found in test cases folders under the name INFLUX_ Project_Name txt see Figure 4 19 After the simulation 1s finished the input signal of the simulation has been saved by default in the output folder Output Project_name INFLUX_ Project_name txt 4 25 Page ad HAWASSI EY model control EY Choose signal data please refer to the manual for the format data Search TCI AA a ern Di TestCases VBM2 4 TC1 Organize Mew folder f a e E EEE Mame Date E COASTUNE_TC1 ba 7 6 Ej HOTSTART TC1 t200s txt 1 6 2 E INFLUX_TC1 bd 1 6 2 INPUT TCl mat MESH Type Text Documen
67. s at the sea sides 5 53 Figure 5 11 Creating PRB s at the land boundaries and harbour walls sess 5 54 Pietro 5 12 Create ATIN O aaa ones paga aa Nak EE dencia ipod enc aaa 5 55 Figure 5 13 Set influx signal parameter in Model Control esses 5 56 Fiure os Oenerated MES aerea obio 5 57 Figure 5 15 Overview of the numerical setup after Preparation is finished sss 5 58 Proure 5 16 Post processio Glorias ovato Queer tado Dr sta esu Ta rl 5 59 Figure 3 17 OUT TOR SD SUN di di ii 5 60 Figure 5218 2D Animation QUIN 5 60 Figure 3 19 SD animation UU ceci 5 61 I19ure 5 20 Domain settings of TSE CASS lista a Ep 5 62 Figure 5 21 Settings for mesh properties and model control for the test case 1 ooooonninnnnninininnnnnnnnnnnnm o 5 63 Figure 5 22 Model setup of test case 1 overview as result of Preparation sese 5 63 Figure 5 25 Domain settings of test case Zeta idad ie 5 65 Figure 5 24 Settings for mesh properties and model control for the test case 2 sss 5 66 Figure 5 25 Model setup of test case 2 cccccecceesssssessseeseseeseseeseeeeeeeeeeeseeeeeeeeeseeseeeeeeeeeeeeeeeeeeeeeeeeeeeeees 5 66 Figure 5 26 The upper plot shows the snapshot of surface elevation and the lower plot shows contour line OL AJAG gna qe MAA d meses EDU a A a A M De EE Es NUN MM Lip aaa 5 67 l4lPage
68. s loaded Physical properties Influx line 1 delete influx Damping Zone damp 3 m a delete damp PRB refi prb s 0 2 delete prb Generate Mesh Preparation RUN STOP Figure 5 39 Domain settings of test case 6a 5 Tl Page ad HAWASSI Bl mesh properties mr Ese E B model control Yol cuu Initial conditions Zero Parameters Approx Tp s ppl max ds m Time Signal shallowest depth m Harmonic Influx 1 Direction deg 270 Tp s 8 Amplitude m 0 25 Figure 5 40 Settings for mesh properties and model control for the test case 6a bathymetry amp numerical setup 2D signal 1 1400 021 2 0 1200 021 an SW 3n 0 100 200 1000 s mi 0 Y 200 0 200 400 600 E 4 x m ih Amp Spec of 1D Signal 0 00 0 16 0 26 0 33 0 38 0 42 800 5 A A 8 B 600 T E E O PA 4 o 400 2 iri i Dispersion Quality 200 200 400 600 800 0 0 2 0 4 05 1 3 3 4A 5 x m f Hz kh Figure 5 41 Model setup of test case 6a 5 78 Page ad HAWASSI t2 00 04 19 hr min sec ort 259 87 sec Figure 5 42 Test case 6a Snapshot of the surface elevation Observation o Harbour lay out in Google Earth 7043 44 S 109001 32 E in 2006 o Different type of harbour waves depending on wave type and on main incoming direction Suggestion o Quantify the wave energy inside the fishing harbour
69. system and computational grid The HAWASSI VBM2Q expects all quantities that are given by the user to be expressed in S I unit m kg s meter kilogram second As a consequence the wave height and water depth are in m wave period in s etc HAWASSI VBM2 only operates in a Cartesian coordinate system 2 2 Page dll HAW ASSI 3 Installing HAWASSI VBM software The HAWASSI VBM installer includes both HAWASSI VBM2 as well as HAW ASSI VBMI there is a separate manual for VBM1 the installation process will automatically install both codes HAWASSI VBM software is programmed under the MATLAB environment therefore the software needs the MATLAB Compiler Runtime MCR installer MCR will install MATLAB Runtime Libraries on the computer so that compiled MATLAB applications can run on PC machines that do not have MATLAB installed The installation of HAW ASSI VBM can be done in two main steps the installation of MCR and the installation of HAW ASSI VBM itself 3 1 System requirements HAW ASSI VBM v 1 1 can run on Windows operating system with 64bit architecture The minimum memory RAM required is 2GB 4GB RAM or more is advised 3 2 First step Installing MCR HAWASSI VBM package v 1 1 requires MCR installer for MATLAB version R2013b 8 2 for Windows operating system 64bit The MCR installer can be downloaded from the MATLAB website http www mathworks com products compiler mcr after downloading the MCR install it by double clicki
70. t MESH TULmat size 397 KB Date modified 7 6 2015 1 49 PN m Y m7 F J HAWASSI VBM Je Documentation J Output TestCases VBMI TestCases VBM2 i Ta i Tc2 File name Figure 4 19 Location of User defined influx signal Table 4 2 User defined influx data format 0 5 s m m t nay a2 nam t ney en nem ippo ee ta nan mn Mam 426lPage ad HAWASSI Note The User defined signal data will be interpolated along the influx line Therefore for long crest wave the data can also be defined only at the start and at the end point If the provided input data is shorter than the length of the influx signal or the influx line so either in space or time the missing information will be substituted with zero s For the option SWAN 2D Spectrum the user should choose the 2D spectrum file The data format of SWAN 2D Spectrum can be seen in Appendix I and an example file of SWAN 2D Spectrum can be found in UserInput VBM2 InfluxSignal SW AN SPEC2D 02 sp2 After the data file has been chosen the user should specify the type of SWAN simulation mode that is used when generating SWAN 2D spectrum N for Nonstationary or S for Stationary see Figure 4 20 E model control Sm Time Signal SWAN 2D Spect Influx 1 Mode M or N N M Monstationary 5 Stationary SWAN SPEC2D 0 Figure 4 20 Signal SWAN 2D spec
71. t case 5 5 75 Page ad HAWASSI Figure 5 38 Test case 5 Snapshot of surface elevation in a rectangular basin with an internal oblique basin Observation o Compare wave disturbance of simulation with long and short crested waves 1 profile linear and nonlinear o Any standing waves Be sure the harbour is fully filled Suggestion o Make simulation with a hot start o Perform time partitioned simulation with 3 intervals 5T6lPage 5 7 Test Case 6 Cilacap harbour ad HAWASSI Test case 6a and 6b show two realistic breakwater configurations in a fisher harbour at Cilacap South coast of Jawa Indonesia The aim is to show effects that different breakwater configurations will have on the wave disturbance in the inner harbour 5 7 1 Breakwater configuration 1 The domain configuration as built in 2006 settings for the mesh generator and model control are illustrated in Figure 5 39 Figure 5 40 and Figure 5 41 File Working Directory C Users DiditiDocuments M a Edit Model Setting Log Tool v 1 1 Nonlinearity Initial cond Influx 1 Time Project Name TC6a_RealisticHarbour_v2 No Zero Harmonic 0 1 300 EH NANA User s note Wave entering a realistic hi Help Simulation Boundary Bathymetry Mesh data Y Show bathymetry m Depth imported data TTT 7 200 400 600 00 1000 1200 1400 1500 Project Setting file i
72. trum Time Interval for Soft Start For the option User defined or SWAN 2D spectrum the user has to specify a time interval for soft start in an input dialog see Figure 4 21 This time interval is used as time duration to ramp the signal such that the signal is smoothly influxed into the domain of computation The default value is 10s For the other option Harmonic and Jonswap the time interval for soft start 1s automatically set to be 3 times the peak period 4 27 Page dll HAWASSI Time interval for soft start 5 10 Figure 4 21 Time interval for soft start Nonlinear Adjustment zone When a signal is influxed into a nonlinear wave model the nonlinearity may produces undesirable spurious modes on the generated waves see Liam et al 2014 a phenomenon that is well known in wave maker theory The appearance of spurious modes can also be expected when using the generation method of HAW ASSI VBM2 In order to avoid the generation of spurious modes the nonlinear terms can be smoothly introduced into the computation domain in the downstream direction from the influx location HAWASSI VBM2 provides this facility the length of a so called nonlinear adjustment zone can be prescribed The length has to be specified related to the peak wavelength see Figure 4 22 Advice For nonlinear simulations the length of Nonlinear adjustment should be at least 2 peak wavelengths For influxing rather high waves

HAWASSI-VBM2 User Manual by © LabMath

Contents

Download Pdf Manuals

Related Search

Related Contents

HAWASSI-VBM2 User Manual by &copy; LabMath

Contents

Download Pdf Manuals

Related Search

Related Contents

HAWASSI-VBM2 User Manual by © LabMath