Home

HAWASSI-VBM1 User Manual by © LabMath

Contents

1. Time Signal Tp Amp i BC East Damping Harmonic finit f end Fitter freq band Hz Generation aa Project Name Method Area Bi direction al User s note Pari Eon RUN t_start dt t_end Post Processing Ee STOF HAWASSLVBM1 please select your working directory Figure 4 2 The main GUI of HAWASSI VBM1 4 14 Page HAWASSI 4 1 1 Working Directory To get started with HAWASSI VBM1 software the user should first specify the working directory for the software This can be done by clicking lawu button and choose a location or folder Working Directory Figure 4 3 Choosing working directory In the specified directory folder the software will create new folder 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 HAWASSI VBM1 comes with two main versions the non dispersive Shallow Water Equations SWE and various variants of the dispersive model the VBM There are 7 sub options of the VBM that deter
2. 1500 Figure 5 8 Plotting panel 5 45 Page HAWASSI A signal and or its spectrum at a specific location can be plotted by activating the Buoy sub panel The user has to specify the location for extracting the signal see Figure 5 9 and the length of the time interval Plots of wave profiles over a specified space interval are obtained through the Snapshot sub panel see Figure 5 10 mn Figure 3 Signal Buoys File f init fend N frequency band 0 2 0005 E EINARLEE Buoys X o n x 500m t init t end Signal 0 269 9377 Spectrum E Energy 100 150 200 t 5 Amplitude Spectrum x 500m 0 15 02 0 25 f Hz Figure 5 9 Plotting signal and or its spectrum Plotting By Figure 2 Snapshot _ Xwest Xeast File V Snapshot t 100 2000 2000 D MTA TTNG GETENG E save plotis Snapshot t 200s 1500 Figure 5 10 Plotting snapshot at specific time 5 46IPage da HAWASS 5 1 2 Comparison with experimental data test cases 3 6 For test cases 3 to 6 the output can be compared with experimental data from MARIN To do so the Post Processing GUI provides a Validation panel for this purpose To illustrate how to use the facility we take test case 5 the Focusing Wave Group as example Steps to be done are as follows 1 In the Post Processing GUI activate Validation panel by clicking
3. Please providetmoadifty the time frame Ms MA Amp smocther for the siqnal Spectrum Energy Plot the spectrum buoy locationis select which spectrum to be plotted b a of for smoother function Plot the amplitude spectrum FH forthe spectrum Plot the energy spectrum default 3 Figure 4 46 Option for plotting signal s at specific Buoy location s fint f_end frequency band 0 2 b E Buoys X Filter the signals at the buoy s location Signal Select a frequency band in Hz Leave it blank to take all frequencies Spectrum Figure 4 47 Option for filtering the extracted signal at specific Buoy location s x frame of the snapshot tin m Awrest Meast Y Snapshot t Plot snapshot of simulation Time ofthe snapshot Plot the Maximum Temporal Provide time of the snapshot t t1 t2 5 Amplitude MTA together iin s with snapshot plot Figure 4 48 Option to plot a Snapshot 4 2 3 Validation HAWASSI VBM1 provides a Validation panel to compare results of simulation with experimental data Data of the experiment are loaded by clicking Load File see Figure 4 49 The data should have the following format in columns Except from a first header row the first column contains the discretized time s and the other columns contain the wave elevation n at specified location s The locations are specified in the first row above the corresponding column t
4. ccccceessssssssseeseeeeeeeeeeeeeeseeeeeeeeeeeeeeeeeeeeeeeeeeeseeeeeeeeeeeeeeeeeeees 4 27 Figure 4 32 Time Stepping in Advance Setting ccccccccccccccccceeceeaeeeeesseseeeecceceeeeeeeseeaaaaaessseseeeeeeeeees 4 28 Figure 4 33 Buttons for running the simulation in the main GUL cece ceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 29 Figure 4 34 Overview of the numerical setup panel showing details of the simulation to be performed4 29 Figure 4 35 Description of Dispersion Quality plot ccccccccccsssssssssseeseeeeccceeeceeeeseeaaaesessseseeeeeeeees 4 30 Poured GADING S aana gag eure agan aaa aga anana a fr aga jaan Ehret rennet tcc her Shere Jaa a gag a naban bag ese eer ester 4 30 Pie 4 57 Post Processing GU asasaran aaa kana anapa aa aa kaag ec sees teen aaa papa AE aa aa Bana aa aa ank galan naka 4 32 Figure 4 38 How to load simulation data in Post Processing GUI cece eceeeeeeeeeeeeeeeeeseeeeeeeeeeeeeeeeees 4 33 1 3 Page HAWASSI Figure AAN on NE cc caryieg agan bana nag saat kani nen a aa aa gan gkah Ag aa Ban naen EE an ka kad TEETER 4 33 Figure 4 40 Domain for showing the animation cccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 33 Figure 4 41 Option to include Maximal Temporal Amplitude MTA in the animation 00 4 33 Figure 4 42 Saving animation as GIF os sasate aaa aa wau ng aa aa rer ner a ag da aa e a aaa naa a a a atest aana a ag aa a Tee 4 34 Figure 4 43 Sav
5. gt Open Project then select 5 41 Page 4 HAWASSI TC1_Intro_Flat mat is the input setting for Test case 1 After the project is loaded the Main GUI will be filled with data as in Figure 5 3 Load Project data GG HAWASSLVBML b TestCases VBM1 TO uh ts re New Project Open Project Ctrl Q_ Oa Arrange by Folder Name Ctri Q Tc1_Intro_Flatmat MATLAB yee MATLAB Data Size 1 30 KB di Documentation Date modified 2 9 2015 11 57 AM J Output Ji TestCases_VBM1 Wb TO W TE File name MAT files mat Figure 5 2 Open the project of Test Case 1 Log Tool Help Working Directory C Users Didit Documents MATLAB Wave model ______ Geo bath amp boundary Xwest dx Xeast Xspace 400 4 2000 Dispersion Optimized VBM 1 profile 4 Nonlinearity Linear z Bathymetry dvalue Hz si TYPE Flat bottom 4 Depth m 30 Initial conditions ______ Boundary Zero v BC West Damping Time Signal location Tp Amp Harmonic a o 2 f_init f_end Fiter freq band Hz Generation Method BC East Damping Area Uni Project Name TC1_intro_flat User s note Wave propagation above flat bottor rreman STOP Project file is loaded Figure 5 3 The Main GUI when the input project of test case 1 is loaded 4 Click Prepar
6. ilah Signal 0 300 sawe as GIF delay 0 1 inf Amp Spectrum Energy Awest Xeast Snapshot t 400 2000 MTA Save frames format png Animation ID save plotis 6 Validation Case File measurement position s Load File Spectrum smoother Energy Figure 4 37 Post Processing GUI 4 32 Page 4 HAWASSI Simulation Data uploaded file data loaded lt a DATA TC1_Intro_Fla data loaded Other Simulation data for Post Processin Figure 4 38 How to load 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 of the animation are needed By default these fields are automatically filled based on the uploaded simulation data The time interval for the animation can also be specified to make an animation with a time step that is a multiple m of dt the value m has to be given under dt Figure 4 39 a Animation tinit tend Hct show the animation every Time what number of dt K bs select a natural numbers starting and Ending time 5 for the animation Figure 4 39 Animation time The spatial interval can be specified in Xspace see Figure 4 40 To include the Maximal Temporal Amplitude MTA maximal wave elevation during the simulation the option MTA should be checked as shown in Figure 4 41 a Animation Awest gt seast Space A space domain m for showing the anim
7. 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 98 1 4287 94 4 Volume 1 ISBN 13 978 981 4287 96 8 pages 122 128 e 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 E van Groesen D Adytia amp Andonowati Near coast tsunami waveguiding phenomenon and simulations Natural Hazards and Earth System Sciences 8 2008 175 185 e G Klopman M Dingemans amp E Van Groesen Propagation of wave groups over bathymetry using a variational Boussinesq model Proceedings Int Workshop on Water Waves and Floating Bodies eds Sime Malenica and Ivo Senjanovic Plitvice Croatia April 2007 pp125 128 e 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 56 Page HAWASSI 6 2 Other references e P A Madsen and O R S rensen A new form of the Boussinesq equations with improved linear dispersion characteristics Part 2 A slowly varying bathymetry Coast Eng 18 183 204 1992 e O Nwogu Alternative form o
8. 0 02 0 3 measurement data from Generation area Uni Project Name TC4_ir12s MARIN Nonlinear adjustment 2 User s note TC4_Irr_103001_20_ 40min o Make a simulation with hot start at 8 min So RUN 1200 1 2400 STOP 5 53 Page 5 6 Test Case 5 Focusing Wave Group A focussing wave group is designed in a wave tank to get at a certain position a high amplitude wave by constructive interference of waves with different wave length and therefore different speed the slowest daHAWASS waves are influxed first and will be caught up by the longer waves that are influxed later as a consequence of the effect of dispersion When all waves are in phase a maximal amplitude will result Suggestions o Dispersion is essential to get the correct positioning of the various influxed waves compare for instance with a linear SWE simulation o Compare linear and nonlinear simulations see spectra at different positions in particular also at the focusing point o Compare result of simulation with the measurement data from MARIN Lakhturov Adytia amp Van Groesen 2012 File Setting Log Tool Help Working Directory CAUsers DiditiDesktop _ VBMi Wave model Dispersion Optimized VBM 2 profiles Nonlinearity Weakly nonlinear Initial conditions Time Signal location User defined 4 10 fint f end Fitter freq band Hz 0 15 1 5 Generation 7 4 1 Method Cea cca Nonlinear adjustment 1 t
9. The input data will be checked and processed An overview of the input and numerical setup will be shown to the user via a panel plot called Numerical Setup The panel shows an overview of domain bathymetry boundary conditions influx location the influx signal and its amplitude spectrum the initial condition of the surface elevation 7 and the dispersion quality of the model compared to the exact Airy linear theory see Figure 4 34 The dispersion quality plot is described in Figure 4 35 2 RUN The 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 pushed the simulation will terminate and the simulation results until the time when the calculation is terminated are stored When a calculation is finished the software will automatically call the Post Processing GUI and load the simulation data directly to it 3 Post Processing When the button is pushed the Post Processing GUI will appear Details of the GUI will be described in Section 4 2 4 28 Page Preparation Process the model setup HAWASSI Figure 4 33 Buttons for running the simulation in the main GUI mn Figure 1 Numerical Setup File Gus Bathymetry m KATONE DE JDamping zone Embedded wave generatigqa 500 0 500 100 200 t s Amp Spec of influx signal C amp V m s 1500 20
10. a boundary an influx signal is generated in the computation domain from a point in the interior of the domain The software provides three options for types of the influx signal to be used a harmonic signal a signal with JONSWAP type of spectrum and a user defined signal see Figure 4 12 With a choice from the three options the user should specify the location of the signal and parameters related to the chosen signal see Figures 4 13 and 4 14 Time Signal Mone a Mone select signal influx Harmonic Jone vap User detined Figure 4 13 Influx signal option For the option Harmonic the user needs to specify the wave period Tp and the amplitude Amp 4 19 Page HAWASSI For the option Jonswap the peak wave period Tp the significant wave height Hs and the steepness parameter Gamma y needs to be specified see Figure 4 14 Time Signal location Tki Amp FEFYETEETTTTETEETTITEETETTETTETTTTEETETTTT x location m of the signal influx Provide The peak period Tp 5 Amplitude m Figure 4 14 Harmonic influx signal Time Signal location Tp Hs gamma x location mm of the signal influx Provide The peak period Tp s Significant wave height Hs m Garmma Figure 4 15 Jonswap influx signal 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 The data should consist of two columns wit
11. and or its spectrum in the Validation panel 4 36 Page HAWASSI 4 3 Post Processing GUI for Internal Flow For the wave dynamics HAWASSI VBM1 only calculates surface variables the elevation and the potential at the surface of the fluid Using the expression for the internal flow that is the basis of the governing equations of HAWASSI VBM1 it is possible to calculate the internal flow properties 4 3 1 Interior Calculation To use Internal flow as a post processing step before the start of the simulation the user should have checked the Internal Flow panel in the GUI in the Advanced Settings GUI shown in Figure 4 31 see Figure 4 51 4 Advanced _Settings BJ vevi Tool Help Internal Flow Advanced Settings N Pressure E Vertical Velocity Horizontal Velocity T Horizontal Acceleration D Vertical Acceleration Time Interval Specify the time to save the input for internal flow calculation Figure 4 51 Activating internal flow calculation in the Advanced Settings 4 3 2 Post Processing of Internal Flow After the simulation is finished the Post Processing GUI for Internal flow can be accessed from the Post Processing GUI for the wave dynamics by clicking Other gt Interior Calculation see Figure 4 52 and Figure 4 53 GUI_PP Figure 4 52 Calling Post Processing GUI for Internal Flow 4 37 Page MI HAWASSI L lt Internal Flow oan l Interior Flow Calcu
12. be longer 2 The User defined option of VBM is only for advanced users who know and understand the energy spectrum of the wave to be simulated and the effect of their choices i e how the choice of the frequencies in the Airy functions of VBM may affect the results Wave model Dispersion Optimized WBM 2 profiles p54 SWE Non Dispersive YEM Parabolic profile Optimized WBM 1 profile Optimized WBM 2 profiles Nonlinearity Optimized VBM 3 profiles VBE 1 profiles user defined WBM 2 profiles user defined VBM 3 profiles user defined Figure 4 5 Choosing model type Wave model Dispersion VBM 2 profiles user d i Nonlinearity linear fregi freq Airy UserDefined value Hz Provide custom value for Airy profile in Hz Figure 4 6 Airy user defined value Hz The option to include nonlinear terms in the mathematical model is given in the option Nonlinearity In the current version of HAWASSI VBM1 the weakly nonlinear VBM is used see Adytia amp Van Groesen 2012 for further details For relevant applications that are given in Test Cases the weakly nonlinear 4 16 Page HAWASSI version is sufficiently accurate to describe wave phenomena with a good match when compared with experimental data 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
13. flOW eee eeeeeeeeeeeeeeeeeeeeeeeeeeees 4 39 Figure 4 56 Panel of Post Processing in the Post Processing GUI for Internal Flow 00 eee 4 39 Figure 4 57 Sub panel Plotting and Animation in the Post Processing GUI for Internal Flow 4 40 Figure 4 58 Animation of dynamic pressure of internal flOW 0cceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 40 Figure 5 1 Select Working Directory sasasaas aana adan anana sasadan abaan gana naak kah anae eaaa naak ah awane ea aaa dah aaa ANA waana 5 41 Piour 5 2 Open the project OF West Case sasanak ga aa NANANG KANG ENENGE DEGAN NE SAI NGGEN GA GAN nseri ONONE BA NGANAH 5 42 Figure 5 3 The Main GUI when the input project of test case 1 is loaded aaaaeeeeeeeeenanaanan anan ene 5 42 Figure 5 4 Preparation process and plot of the numerical setup for test case 1 eaaa aananane 5 43 Figure 5 5 Progress bar during the time integration of the simulation cccceeeeeeeeeeeeeeeeeeeeeeeeees 5 43 Figure 526 Post Process io pW aaa penegen a EEN aa ja Na aa GE E aaa a Ta aa E Ng NG Aa EN GAN ag KB agi a aaa ai 5 44 Figure 5 7 Calling animation in the post processing GUI 000aaaaaaaaananeenenanananaan anana nana nana aana eee 5 45 Fei Ae POC AS My aana sasaka anaa ANG a aa aaa TA AR E a aaa a Ak a a aa a Keeani a ga aaa aia aa a aaa ah bah 5 45 Figure 5 9 Plotting signal and or its spectrum cccceeeeee
14. getting started directly with the software www hawassi labmath indonesia org DEMO version with restricted functionality The Demo version of HAWASSI VBM1 has restricted functionality and facilities Only 1 parabolic vertical profile Only linear No internal flow calculations No comparison with external measurement data Full functionality and facilities under licence gt Licence for University Thesis Projects gt Research Licence for extending capabilities and or functionalities gt 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 8 Page HAWASSI 2 Description of HAWASSI VBM1 2 1 Introduction HAWASSI VBM1 The purpose of this chapter is to provide the user with relevant background information of HAWASSI VBM1 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 HAWAS
15. needed for calculating internal flow of VBM It is only calculated when the internal flow has been activated Internal Flow data It is calculated from the Post Processing GUI of internal flow Is a sub folder that will contain animation gif files from the PostProcessing GUI SE Is a sub folder that will contain figure fig amp png files from the PostProcessing GUI JA KA ahaaa ANNA Is a sub folder that will contain frames png or jpg from Post Processing GUI 3 4 5 7 4 31 Page MI HAWASSI 4 2 Post Processing GUI for Wave Dynamics The Post Processing GUI Figure 4 37 can 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 38 After a simulation is finished the Post Processing GUI will automatically pop up and then the data of the recently finished simulation will have been loaded automatically There are three main panels in the Post Processing GUI 1 e Animation Plotting and Validation Each panel will be described in the next subsection 4 2 1 4 2 2 and 4 2 3 respectively Module 1 0620 file DATA_TC1_intro_ flat Bathymetry Hamiltonian Energy Significant Wave Height Wave Disturbances o MTA tint t end Fat Time 0 300 4 t init fT end Xwest Xeast frequency band 0 2 ASpace 400 2000 Buoy s x t init t end
16. start dt t_end 0 0 1 150 Geo bath amp boundary Xwest dx Meast Aspace 10 0 075 120 Bathymetry TYP Flat bottom Depth m 1 Boundary BC West Damping BC Fast Damping kd Project Name TC5_FWG User s note TCS_FWG101013_24prof 5 54 1Page da HAWASS 5 7 Test Case 6 New Year Draupner Wave One of the highest waves ever measured in real seas at a fixed measurement point is done at the Draupner platform in the North Sea In a huge storm at January 1 1995 a video camera pointing downwards recorded a wave of approx 18m crest height the definite proof of the appearance of freak or rogue waves in nature The input signal used here has been designed and used at MARIN to reproduce this scaled wave with very broad spectrum in the laboratory gd vem File Setting Log Tool Help Suggestions Working Directory o Compare linear and GAliseralliiditWeskiopt VEM nonlinear simulations Wave model Geo bath amp boundary and conclude that a Xwest dx Xeast both shou r very Dispersion Optimized WEM 2 profiles Y Wace Ps high wave linear oe Weakly nonlinear gt Bathymetry dispersive focusing Type Flat bottom as in Test Case 5 z Depth m 1 seems to be the basic Initial conditions ingredient to which Boundary nonlinear effects add BC West Damping to increase the wave Time Signal height even more location ae d in o Compare result of the User de
17. the data of the bathymetry should be gt 0 see Table 4 3 Table 4 3 File format for User defined bathymetry Warning HAWASSI VBM1 v 1 01 does not have the capability for run up This means that the software can only take a positive value for the depth the bottom profile should be below the Mean Sea Level MSL Two types of boundary that are used in the software are damping zone and partial reflective boundary PRB see Figure 4 28 and Figure 4 29 The purpose of the damping zone sponge layer is to absorb the outgoing wave so that it will not reflect and disturb the waves in the rest of the computation domain The efficiency of the damping zone is determined by the length of the damping The length of damping zone should be at least 2 times the simulated peak wavelength a shorter length may create reflection Advice To avoid reflection from the damping zone the length of damping zone should be at least 2 times the peak wavelength Warning Do not place the influx location inside the damping zone 4 25 Page 4 HAWASSI Boundary width Bl West Damping wa Damping Width of the damping zone m at least 2 peak wavelengths Figure 4 28 How to set boundary condition in the West or left location Boundary width Bl East Damping a Damping Width of the damping zone m at least 2 peak wavelengths Figure 4 29 How to set boundary condition in the East or right location A modified Sommerfeld boun
18. the given dx Estimate how many waves you will see in which interval with the expected amplitude o Run the simulation and perform post processing Verify if your expectations are correct o Perform other simulations for this wave case o See effect of changing dx change period and choose sensible dx o Change influx position damping zone boundary values etc 5 49 Page HAWASSI 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 input signal because phases that 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 compare outputs by combining them in a single plot o Observe effects of changing Jonswap parameters o Perform non linear simulations increasing Hs 5 50 Page MI HAWASSI 5 3 Test Case 2 Wave Propagation above a Sloping Bottom This test case continues Test Case 1 now for waves above varying bottom Suggestions All suggestions from Test Case 1 apply here The main effects to be investigated observed are now related
19. 00 1500 2000 Initial wave profile 0 1000 2000 x m 0 000 130 180 220 260 290 32 0 Dispersion Quality 0 2 4 6 8 10 12 kh Numerical setup Domain of computation Damping zone Partial Reflective Boundary Condition with reflection coefficient Refl Embedded wave generation Bathymetry profile Signal influx influx location Initial condition for surface elevation n Dispersion quality Velocity of Exact vs Model Amplitude spectrum of the influx signal Figure 4 34 Overview of the numerical setup panel showing details of the simulation to be performed 4 29 Page MI HAWASSI Dispersion guality Velocity of Exact Airy theory versus Model amb D Phase Velocity m s 2 gt oS O Group Velocity m s Exact Airy theory Model kh of the peak wave to be simulated Figure 4 35 Description of Dispersion Quality plot 4 1 10 Output Files Output files will be stored in the folder WorkingDirectory Output ProjectName as illustrated in Figure 4 36 A description of all output files is given in Table 4 4 A J D TO intro flat InternalFlov lt Oei t VBMI y Output TCi intro_flatInternalFlow Organize Include in library Share with Burn New folder E3 Programs A TE Ta System and Security RR User Accounts and Family Safety E Recycle Bin J vem Anim_GIF m Figures di frames j DATA TCL intr
20. 1 0620 Bathymetry Hamiltonian Energy significant Wave Height Wave Disturbances MTA f_init f_end frequency 0 2 Buoys X t_ini t_end Signal 0 300 Amp Spectrum Energy Xwest Xeast 4400 2000 MTA Snapshot t save plotis Compare Signal Amp Spectrum Energy Figure 5 6 Post Processing GUI For showing an animation of the surface elevation the Animation panel should be activated In this panel the user can save the animation as a gif file or save them each frames by clicking save as GIF and save frames respectively see Figure 5 7 The animation will pop up after the RUN button 15 pressed 5 44 Page HAWASSI t 00 03 40 hr min sec or t 220 00 sec Figure 5 7 Calling animation in the post processing GUI 8 The Plotting panel can be activated for plotting some input and output parameters used in the simulation such as Bathymetry bottom profile Significant Wave height Maximum Temporal Amplitude MTA maximum wave height Hamiltonian Energy and Wave Disturbances normalized significant wave height Hs with respect to Hs at the influx location After the RUN button is pushed the plots will pop up see Figure 5 8 Plotting Balwmetre a Hamilton M By Figure 4 Maximal Temporal Amplitude MTA m nS AY Fite It ET EOCEEN Maximum Temporal Amplitude MTA Significant Wave Height C Wave Di MTA h
21. AT MAT VBM asia sss sarees acne ana Na aa E aaa KE ag na aga ana E a a Ba a aaa Ke ag aga aana na aaa akah 4 17 Figure 4 8 Initial condition for wave ecleVatiOn cccccceeseseesessseeseeeeseeeseeeeeeeseeeeeeseeeeeeeeeeeeeeeeeeeeeeeeees 4 17 Figure 4 9 Initial condition for single Gaussian hump ccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 18 Figure 4 10 Single Gaussian hump profile cccccseseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 18 Figure 4 11 Bipolar hump 1 positive hump at the right side 2 0 0 0 eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 19 Figure 4 12 Bipolar hump 2 positive hump at the left Side cc eeseesssseeeeeeceeeeeeeeeeeaaaeeeseseseeeeeeeess 4 19 Pidie tux Stet OPuOM css cscaiasiinssecneceascdecevanee ea n eE OREA eE anirno riria 4 19 Fienie 4 14 Harmonie mR SIS TI aaa ag a a aaa ga ENI AO ga raa aaa ga aa Na 4 20 Figure 4 15 Jonswap influx signal oc aesssuscassacecsuseiraventncneesievebenaadsenenccaneea yenwensossuhebeeasdsonwensaberseaeiuenenssabanntes 4 20 Figure 4 16 Time interval for soft start s of the user defined signal aoonananananan e nnnnne anana anana ene 4 21 Figures 4 1 7 Filter TING SIA en oe nm eR ROO a naa agan an a 4 21 Figure 4 18 Generation method for embedded influx cccssssssssesesseeeecccceeeeeeeeeeaaaesessseeeeeeeeees 4 21 Figure 4 19 Nonlinear adjustment sc2c 002susc5 0nsaxecontideneneniousientesce
22. HAWASSI VBM1 User Manual dH AWASS by LabMath Indonesia ver 1 150827 by The HAWASSI team 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 VALF Wa Copyright 2015 LabMath Indonesia HAWASSI Contents Pean O e E E E A E E KET nan A 1 6 e 5 as E a Najan aa Ga EE E E E E E 1 8 2 Desciphonol MAWAS Sl VIS MN kisenan asana abiane a a DB KAGEN EEE Dk aa a NE AG GEN Gg PG E AE D aaa S 2 9 2 1 Introduction HAWASSI VBM1 uu cccccccccccccsscceeeessssseeeeeeeccccccesesssaaaeeagsseeeeeeeeeeeceeeeeeseaaas 2 9 ma Mode bi e Sae jangan aaa an er en Rn ee 2 10 223 Relation to other Boussinesq type Wave models cccccccccececececceeeaeesssseseeeeeeeeeeeeeeeeeseqaas 2 11 2 4 Units and Computational grid 0 0 0 0 cccccccccccccccccccesssssssseseseeeeccceeeeceseseeaaaaaesssseseeeeeeeseeeeeeeaeaaas 2 11 3 Installing HAWASSI VBM1 software 00 0 0 ccccccccccsssssssseseeeeeececeeeeeeeeseeaaaaeesesseeseeeceeeeeeeeeeeeaaaas 3 12 3 1 DY 5 TAG TTS ests aa Ga SA ana a Aa baja a Aa anes Saja da a aa E a Ka aa Da aaa aa a Sanan 3 12 Be Firststep Mans ai MER saus anaa aa aa daa ga Kaga A e aaa aaa ENE aaa WE ag a a NR aa a WAN A e ERS 3 12 3 3 Second step Installing HAWASSI VBM1 i cccsssssessseseeeceeeeeeeeeeeesaaaaseeeessseseeeeeeeees 3 12 4 INS HAWA SIN BM aka
23. M1 software includes three interfaces or GUI s Graphical User Interfaces namely the Main GUI for the input model a Post Processing GUI for the wave simulation output and a Post Processing GUI for internal flow simulation The GUI s act as an input output manager HAWASSI Calculator There is a simple Wave Calculator that expects as input the period of a harmonic wave and the depth and will then calculate all wave relevant quantities by also specifying the amplitude the INPUT calculated steepness is added The calculator can Dispersion exact be accessed by clicking Tool gt Wave Calculator Period s Depth h m Amplitude a m Figure 4 1 Wave Calculator Frequency 0 696 rad s Wave number k 0 0503 Wave length lambda 125m Relatve wave length lambda h 2 5 k h 2 52 Intermediate depth 0 0503 13 9 ms 7 39 ms Exact 4 13 Page HAWASSI 4 1 Main GUI Settings for the wave model initial conditions time signal domain and boundary are managed in the Main GUI see Figure 4 2 Details of each input in the main GUI will be explained in detail in the next sub sections File Setting Log Tool Help Wave model i ia Geo bath amp boundary Xwest de Xeast Dispersion Optimized WBM 1 profile pe pr Mspace Nonlinearity Linear Bathymetry Type Flat bottom Depth m Initial conditions Bound ary BC West Damping
24. SI stands for Hamiltonian Wave Ship Structure Interaction HAWASSI VBM is a finite element implementation 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 VBM 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 7 and the simulation of the interior flow can be avoided the Bous
25. Validation then load the experimental data by clicking Load File as shown in Figure 5 11 Validation i Compare t init t_end Case File sps Signal l measurement position s Pak E eet Load File SA Amp Spectrum smoother Load the file of measurement data Energy to be compared with the simulation Figure 5 11 Load experimental data in the Validation panel 2 Select the file of experimental data in C Users UserName Documents HAWASSIVBM1 TestCases_VBMI TCS see Figure 5 12 When the validation is loaded the user is notified in the lower part of the GUI as illustrated in Figure 5 13 Bay GuLpP ra x Choose measurement data a G TestCases VBMI TCS Organize New folder di A Name J HAWASSI by LMI dJ HAWASSI VBMI di Documentation d Output TestCases VBMI ype Text Document b Ji Size 808 KB Ta Date modified 11 28 2014 6 48 PM E INFLUX_MARIN101013FWG bt 2 TC5_FWG101013 mat VALIDATION_MARIN1O1013FWG bet Wb Te Wb 13 bh ca di cs i TC6 m Figure 5 12 Select experimental data 5 47 Page HAWASSI Validation Tampane t init t end Load File Ss finit f end 2 IAI teapot Amp 0 200 0055 uploaded file C Users DiditDocumentsi HAV F aa T Energy Data validation loaded VALIDATION_MARINLOLO1SFW6 t ka Figure 5 13 Status bar when the experimental data is loaded The user can select Signal and or Spectrum
26. are specified these are interpreted as the x coordinate and the initial condition for 7 with the initial condition for considered to be zero no initial velocity Initial conditions Select initial conditions User defined Single Gaussian hump Bipolar hump 1 Bipolar hump z2 Figure 4 8 Initial condition for wave elevation 4 17 Page HAWASSI Table 4 1 Data format for user defined initial conditions 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 and the central location of the initial profile see Figures 4 8 till 4 11 Initial conditions Amp width cent x single n en hump i Select initial conditions Provide the 4rniplitude rr Width rr fcenter Location of the initial hump rm Figure 4 9 Initial condition for single Gaussian hump amplitude still water depth Figure 4 10 Single Gaussian hump profile 4 18 Page da HAWASSI z otis Sa width depth Figure 4 11 Bipolar hump 1 positive hump at the right side Z amplitude N cent_x still water width depth Figure 4 12 Bipolar hump 2 positive hump at the left side 4 1 5 Time Signal In HAWASSI VBM1 when dealing with a signalling problem an embedded wave influxing method is used Instead of generating influx from
27. ation Figure 4 40 Domain for showing the animation Animation A MTA b show the Maxirial Temporal Amplitude ChAT A during the animation Figure 4 41 Option to include Maximal Temporal Amplitude MTA in the animation 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 42 The delay information is in seconds and the loop information has to be a natural number For making animations in other formats saving frames either in PNG or JPG format can then afterwards be combined using any existing animation software Figure 4 43 If no 4 33 Page HAWASSI specific name is provided in the animation ID field the animation frames will be saved with anim as the default name Animation save as GIF delay loop The saved gif file will be run hs Tor how many loop choose 1 2 3 delay time at each frame during the making of git file ain a Figure 4 42 Saving animation as GIF Animation V save frames format png pny Figure 4 43 Saving all frames in png jpg format 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 F
28. ation button for getting an overview of the numerical setup to be simulated 5 42 Page HAWASSI nm Figure 1 Numerical Setup t_start dt t end File 0 0 5 600 SIG A S AAV Oe OF E Damping zone Embedded wave generatiqn 100 00 Optimizing 1 Vertical Profile of VBM completed 1000 1500 2000 mnnn 1000 1500 2000 Initial wave profile pE e 2000 E E w a Signal influx x m Rp ian a ie signal 0 000 120 180 240 260 290 32 7 Dispersion Quality aa a hak C amp V m s 3 Figure 5 4 Preparation process and plot of the numerical setup for test case 1 5 Run the simulation by clicking RUN button During the time integration process the progress of the calculation is shown by a progress bar in the lower level of the GUI as shown in Figure 5 5 Time Figure 5 5 Progress bar during the time integration of the simulation 6 The Post Processing GUI will automatically appear when the simulation process is done The GUI is illustrated in Figure 5 6 In the GUI the user can plot and animate the results of simulation 5 43 Page HAWASSI ation Data uploaded file DATA_TC1_intro_flat data loaded t init t end dt Time 0300 1 Kwest Xeast Xspace 400 2000 MTA sawe as GIF delay 0 1 loop inf save frames Animation ID Validation Case File measurement position s Load File
29. ces ea aana eter ee eco a aa a aaa aana aan 4 34 4 2 3 BY INO eee nee ae ee ee re ee 4 35 4 3 Post Processing GUI for Internal FlOW cc esessssssssseeeeceeceeeeeeeeeeaaaaeseessseeeeeeeeeeeeeeeeeeaqaas 4 37 4 3 1 Ian ag Cal CON AO ANAA TAN EA A A Er TAN RR TEN EA AE EA 4 37 4 3 2 Post Processing of Internal FLOW ccccccccccssssssssseeseeeeeeccceeeeceeeeaaaaaseeeseseeeeeeeeeeeeeeeaeaaas 4 37 5 HAWASSI Gettin started Test CAS CS a arai aii e e a e ian aaa gan aa dadaran nan 5 41 5 1 Mus at ons t0 use the GUT S scrcascadessesensancosouseenetarseadssandsnosccnscscatonsaemcaanivasanesdgboocenectiatenitoesaanas 5 41 5 1 1 ILA SAKA On Testcase Tessier sakara a ga aae EAEE Gana na aaa 5 41 5 1 2 Comparison with experimental data test cases 3 6 cccccccccccccccceceeeeeeseseessseseseeeeeeeees 5 47 5 2 Test Case 1 Wave Propagation above a Flat Bottom 000seeeeeeeeeeeeeeeeeenanananannnnn nane nane 5 49 5 3 Test Case 2 Wave Propagation above a Sloping Bottom seeeseeeeeeeeeenananaaannn nana nee 5 51 5 4 Test Case 3 Bichromatic Wave Propagation on a Sloping Bottom cccccceeeeeeeeeeeeees 5 52 5 5 Test Case 4 Irregular Wave Propagation above a Sloping Bottom cccccceeeeeeeeeeeeeees 5 53 30 TestCase 3 Focusing Wave TITO asasi daga aaa agan GN aa i isasi aiae iea iiaeia 5 54 5 7 Test Case 6 New Year Draupner Wave cccccsssssssssssseeeccccceeeeeeeeeesaaesesssee
30. dary condition is used to model partial reflective boundary conditions in HAWASSI VBM1 With this choice of boundary the user can specify a reflection coefficient Cp 0 1 With reflection coefficient Cpe 1 the boundary will act as a fully reflecting wall with Cr 0 as a transparent boundary condition so that the out going wave will be fully transmitted Partially reflected boundary can be illustrated by setting for instance Cp 0 6 to let the wave height of the outgoing wave be reflected for 60 and transmitted for 40 from the initial wave height this applies to all wave frequencies When damping zone is chosen HAWASSI VBM1 will automatically set the boundary of the damping zone as a fully transmitted boundary condition PRB with Cr 0 so that all waves will be absorbed or transmitted at the boundary Warning When PRB of HAWASSI VBM1 is used with Cr 0 fully transmitted not all wave frequencies will be transmitted perfectly short waves kh gt z will partly be reflected less than 10 4 1 8 Advanced Settings HAWASSI VBM1 provides an additional GUI so called Advanced Settings to change given standard settings that are used in the simulation The Advanced Settings GUI can be called by clicking the menu Setting gt Advanced Settings see Figure 4 30 and Figure 4 31 Physical parameters that can be changed in Advanced Settings GUI are the value for the gravitational acceleration g m s and the fluid densit
31. dessacnnendidenedseiastzeetedidedentaaanendedenedsniasazentaed 4 22 BAUD A AAD Ne EN ac sets ajaa ag aia a aan as ad aaa aaa anaa a ana e Anaa aa a jarang anga a ol al E 4 22 Fisure 4 2 Domain of computation isosecnesatcasievessdensaostcasevesnsesiexlabcadevesnnes sasatcanesesnses salabcedesesseessansbeaasseeseds 4 23 Figure 4 22 Suggesting for the length Of dX 000 0 anana 4 23 Figure 4273 DB atiy Mery Cy PC sasada aa E aaa a aaa a aa aaa aaa alang Laa aa aaa Na ag aal Ja ae BAEN ad NE AR ag aaa Ng aa Bk ngak aaa 4 23 Fig re 4 24 Setting up Flat WOU Oi accents reas aa aia aa aa Ka gala Ag A ANAN E NAN Ta aa A AA a ga Aga aa a a aa A aaa aga 4 24 Figure 4 25 Setting up Sloping bottom amp Underwater Mountain ccccccceeeeeeeeeseeeeeseeeeeeeeees 4 24 Figure 4 26 Bathymetry profile for the option Underwater Mountain cccccccceeeeeeeeeeseseeeeeeeees 4 24 Figure 4 27 Bathymetry profile for the option Sloping bottom saaaeeeeaaaaeenana nana aana anane nane 4 25 Figure 4 28 How to set boundary condition in the West or left location ccccceeeeeeeeeeeeeeeeeeeeeeeees 4 26 Figure 4 29 How to set boundary condition in the East or right location ccessssssssseseeseeeeeeeees 4 26 Figure 4 30 Accessing Advanced Settings GUI ccccccccccccccccccccesseeesesssseeeeecceeeeeeessaaaaaaassssseseeeeeeeess 4 27 Figure 4 31 GUI of Advanced Settings 0 0
32. e internal flow by specifying time and space information see the left panel in Figure 4 57 Warning Calculation of the internal flow depending on x z and t can create a large amount of data and for fine discretization the calculation time may take long 4 39 Page HAWASSI Plotting Animation Mesh Plot Line Flot tint tend ect A Time Y x F axes z KAE vest Meat cnt cts z m XSpace _ Z bottom top t z axes Z axes Tspace l ce x m xt axes faxes Hz Exim z m lt Internal Flow re tn Interior Flow Calculation Input file file dala aded DATA_TC1_intro_flat t_init t_end X_ west X east bottom dz Z top 50 100 Horizontal Interval 200 1000 Wertical Interval 5 0 1 1 14 Time Interval Post Processing Figure 4015 Animation buwi a i i i i 400 BOO a00 1000 Figure 4 58 Animation of dynamic pressure of internal flow 4 40 Page HAWASSI 5 Getting started test cases HAWASSI VBM1 provides 6 test cases of increasing complexity The first two cases are meant for practicing the software such that the user gets an idea how the software works in handling influx signals boundary conditions different wave models SWE amp VBM etc Test cases 3 up to 6 are cases for which simulation results can be compared with experimental data in a hydrodynamic laboratory the experimental data are ava
33. eeeeeeeeeeeeeeseaaaas 5 55 EE O AAN KAR E E A KR NE E E NE N AN NA O eee er 6 56 6 1 VBM References to basic papers and applications ccccccccececececeeceeseeeesseseseeeeceeeeeeeeeeeaaaas 6 56 oy AJENGE TET EN re eine T E AA EAT AER AA BA AARE E errr 6 57 l 2IPage HAWASSI List of Figures Figure 2 1 At the left a comparison of the dispersion relations of Boussinesq type models by Madsen amp S rensen 1992 Nwogu 1993 and the VBM with parabolic profile by Klopman et al 2010 with the Airy linear theory dispersion relation as function of kh At the right is the comparison of the dispersion relation of VBM with 1 2 3 Airy profiles and a parabolic profile with the Airy linear theory 2 11 Peu al YAS e O a A E ence nt aaa ga aapa pa a a ng aa ag aaa ag NE 4 13 Figure 4 2 The main GUI of HAWASSI VBMI eaaa ee ne nana nana nenen e nana n nana nenen ane nana naen 4 14 Figure 4 3 Choosing Working directory aaa aana Gaga Ng aa KA NGENE A NGAGAGI a ANG ag da kaa ga KG Kaga gaga aii 4 15 Fig re 4 4 Project Name amp USET S 1 Reman Renner Re enn nee Maer NAN eee KAN He nent Gae ada tore E gali aa gaga Na Ag dag ada eee 4 15 Figure 4 5 Choosing model type ian as sa sawana Nga aan ana Ana yaaa aan Ae kaa gansa gana kanan ai IA gak aaa baa anga Aak 4 16 Figure 4 6 Airy user defined value HZ ccccceeeeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 16 Pepe She INGON S
34. eeeeeeeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 5 46 Ficure 5 10 Plotting snapshot at speciile TING 5 aana apa aaa aaa ag a aa dab gae aaa aana aaa dap a a alaga aa aaa gara nani 5 46 Figure 5 11 Load experimental data in the Validation panel cccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 5 47 Pieure 5 12 Sclectiex perimental Cal A aaa sia akan aan aa kn agengan as ba babana aga ea outa oto an bag aaa aaa 5 47 Figure 5 13 Status bar when the experimental data is loaded cece eeeceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 5 48 Figure 5 14 Signal comparison between the experimental data and the simulation 5 48 1 4 Page HAWASSI List of Tables Table 4 1 Data format for user defined initial conditions seseeseeeeeeeenenanana nana nean anana nean eee 4 18 Table 4 2 File format for the user defined ITIK asat kaa aa aa aa anan eaaa aga pana Kana aab aa aa a aa aa 4 20 Table 4 3 File format for User defined bathymetry ccccccccssssssesssseeeeccceeeeeeeeeeeaaaasesseseeseeeeeeeess 4 25 Table 4 4 Output files folders ss sana saa kanan akang ANGEN GA AAN GTA AA AAN EN aiao ENENG BANG NAN Aga na Daka a aaa 4 31 Table 4 5 Data format for the Validation panel ccccccccccccccecceeeeesssssesseeeecceceeeeeeeeaaaaaaasesseseeeeeeesess 4 36 1 5 Page HAWASSI Preamble Waves are fascinating important and challenging The importance can be substantiated from some well k
35. f Boussinesq equations for nearshore wave propagation J Waterw Port Coast Ocean Eng 119 618 638 1993 6 57 Page
36. fined 10 col 4 r f_init f_end nonlinear simulation Filter freq band Hz 0 01 1 7 with the Generation x measurement data Method Area Uni Eroa E TC6_Draupner from MARIN Nonlinear adjustment 4 User s note TCS New Year wave o The Draupner platform was above a t_start dt t_end 30 0 1240 water depth of approximately 70m Project file is loaded Rescale the laboratory data 1 70 to real life data and perform the simulation at real scale What do you observe with respect to the calculation time Calculate the relative computation time the calculation time divided by the length of simulation time o References o A L Latifah amp E van Groesen Coherence and Predictability of Extreme Events in Irregular Waves Nonlin Processes Geophys 19 2012 199 213 ISSN 1023 5809 o Kharif C Phelinovsky E and Slunyaev A Rogue waves in the Ocean Springer Verlag Berlin Heidelberg 2009 O Dysthe K Krogstad H E and Muller P Oceanic rogue waves Annual Review of Fluid Mechanics 40 287 310 20 5 55 Page HAWASSI 6 References 6 1 VBM References to basic papers and applications e D Adytia Simulations of short crested harbour waves with variational Boussinesq modelling In Proceedings JSOPE 2014 2014 912 918 e D Adytia M Woran amp E van Groesen Effect of a possible Anak Krakatau explosion in the Jakarta Bay Proceedings Basic Science International Conference 2012 Malan
37. g Indonesia K1 5 ISBN 978 979 25 6033 6 e 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 e Ivan Lakhturov Optimization of Variational Boussinesq Models PhD Thesis UTwente 9 November 2012 e Didit Adytia Coastal zone simulations with Variational Boussinesq Modelling PhD Thesis UTwente 24 May 2012 e I Lakhturov D Adytia amp E van Groesen Optimized Variational ID Boussineq modelling for broad band waves over flat bottom Wave Motion 49 2012 309 322 e 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 e 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 e G Klopman Variational Boussinesq modelling of surface gravity waves over bathymetry PhD Thesis UTwente 27 May 2010 e 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 e D Adytia amp E van Groesen Variational Boussinesq model for simulations of coastal waves and
38. h the first column is the time vector in s and the second column is the signal data as illustrated in Table 4 2 Table 4 2 File format for the user defined influx s t s t pt After the data file is chosen then there will be a pop up asking for Time interval for soft start see Figure 4 16 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 options Harmonic and Jonswap the interval is set automatically to be 3 times the peak period 4 20 Page HAWASSI Time Signal location User defined INFLUX MARIN Time interval for soft start s if wa Figure 4 16 Time interval for soft start s of the user defined signal Time Signal f init f end Filter freq band Hz bs Filter the influx signal set a frequency band Hz f_init f_end Figure 4 17 Filter time signal In the software it is possible to filter the chosen influx signal by specifying Filter under Time Signal panel HAWASSI VBM1 provides 4 methods for generating waves using embedded influxing see Figure 4 16 Two main types are area and point generation For area generation a confined spatial function is used as a force function in the wave equations As a consequence the corresponding area around the influx location will move vertically with decreasing amplitude away from t
39. he first column first row item can be left blank Note that the header information in the first row should be preceded with Yo so it will not be interpreted as data An illustration of an acceptable data format is shown in Table 4 5 During loading of 4 35 Page 4 HAWASSI the experimental data the location s will be read automatically After the data is loaded the signal and or its spectrum can be compared with option similar as in Plotting panel see Figure 4 50 validation Case File measurement position s Load File SI Load the file of measurement data sia file MARIN 101013_PWG td to be compared with the simulation Figure 4 49 How to load experiment data in the Validation panel Table 4 5 Data format for the Validation panel some header O P t nant nant TING 0 a S a mn t4 nEn next Emt n x xpt n t ne Xm ty pty me ety NEk ne Extn validation Time frame for signal comparison Compare the time signal sampe tinit tend Please provide the time frame 5 p 0 180 53 of the simulation VU En Fi 0 200 0055 5 Amp Select the range of Compare the spectrum spectrum SEN R the signal frequency Hz select amplitude andor energy spectrum b A Energy 3 Camere thearn E de Seci Coefficient of for smoother function forthe spectrum Compare the energy spectrum P default 3 Figure 4 50 How to compare signal
40. he influx point In this generation area the desired wave is not accurate outside the generation area the waves are according to the prescribed signal For point generation a delta Dirac function is used as force function in the wave equations The signal is influxed into the domain from a single point at the influx location the amplitudes of the given input signal are then enlarged to obtain the correct influxed wave Based on the direction of propagation the generation method can be in two directions the bi directional influxing or in one direction the unidirectional influxing The unidirectional influxing is divided into Uni and Uni for the wave direction to the right and to the left respectively Generation Method Area Uni Area Bi direction Area Uni Area Uni Point Bi direction Select generation method of signal influx Figure 4 18 Generation method for embedded influx Advice For influxing of a rather high wave compared to the depth it is advised to use area generation since point generation requires the influxing of even higher waves 4 21 Page HAWASSI When a signal is influxed into a nonlinear wave model the nonlinearity may produces undesirable spurious modes on the generated waves see Liam e a 2012 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 HAWASSI VBMI1 In o
41. igure 4 44 and Figure 4 45 Plotting as Plot the bathymetry Significant Wawe Heig Plot the Significant Wave Height Hs rm A MTA Plot the Maximum Ternporal b Amplitude MITA m Figure 4 44 Option for plotting the bathymetry significant wave height Hs and MTA Plotting Hamitonian Energy gt Plot the Hamiltonian Enerqy as function of time E Wave Disturbances ian Plot the Wave Disturbances Wave disturbances is normalized significant wave height Hs with respect to Hs at influx line Figure 4 45 Option for plotting the wave energy and wave disturbance 4 34 Page HAWASSI Data at a specific time or location can be obtained by selecting the option Buoys and or Snapshot Buoy will extract signal information at a position to be specified in the computation domain such as the time signal and or its spectrum see Figure 4 46 For the spectrum the frequency band option is provided to filter the extracted signal information see Figure 4 47 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 48 Plotting x lacations where Plot signals at certain location s the signalis will be extracted in the darman Buoy s ane NS Dx XL X A C Please provide the location tint tend Tine are of the si signal Plot the signal Mbhuoy location s Signal
42. ilable in the software package For getting started with the software we take the first test case as an example Then we illustrate how to compare simulations with data for the test cases 3 6 After that the test cases are described Acknowledgements We are very grateful to be allowed to use measurement data of MARIN Maritime Research Institute Netherlands Dr T Bunnik Only by testing with realistic data the software can be validated and improved 5 1 Illustrations to use the GUTS 5 1 1 Illustration Test case 1 Steps for getting started the software are listed as follows 1 Open the HAWASSI VBM1 software by clicking the shortcut so called VBM1 in the Desktop or in the Start menu The Main GUI of the software will appear as in Error Reference source not found 2 Select your Working Directory by clicking lia button The user can select any folder that he she wanted In the select folder the software will create Output folder By vemi H Select Directory to Open OT m lt My Documents HAWASSI MBMI ka I Search HAWASSI WBMI p Organize New folder Ji Documentation Output TestCases_VBM1 4 HAWASSLVBML di Documentation gt Output gt TestCases VBM1 Folder HAWASSI MBMI Figure 5 1 Select Working Directory 3 Test cases are stored by default in C Users UserName Documents HAWASSI_VBM1 folder see Figure 5 2 Open test case l by selecting File
43. ing all frames in png jpg format sasana anana anna 4 34 Figure 4 44 Option for plotting the bathymetry significant wave height Hs and MTA eee 4 34 Figure 4 45 Option for plotting the wave energy and wave disturbance cc0cceeeeeeeeeeeeeeeeeeeeeeeeeees 4 34 Figure 4 46 Option for plotting signal s at specific Buoy location s cccceeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 35 Figure 4 47 Option for filtering the extracted signal at specific Buoy location s cceeeeees 4 35 Figure 4 45 Option 40 plota Snapshot sesccecceessascseecerectasiacsaceereeeavanscosebeeeeeuacenacsocteteemcctanacuocbeeeeeecteniensest 4 35 Figure 4 49 How to load experiment data in the Validation panel ccccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 36 Figure 4 50 How to compare signal and or its spectrum in the Validation panel c ceeeeeeeeees 4 36 Figure 4 51 Activating internal flow calculation in the Advanced Settings ccccceeeeeeeeeeeeeeeeeeeees 4 37 Figure 4 52 Calling Post Processing GUI for Internal FIOW 1 0 0 0 aana anana anae 4 37 Figure 4 53 Post Processing GUI for internal flOW cece cece ceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 38 Figure 4 54 Panel for Interior Flow Calculation 0 ccccccceeseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 4 39 Figure 4 55 Notification in the Post Processing GUI for internal
44. 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 clicking the installer and following the instruction in the installation wizard 3 3 Second step Installing HAWASSI VBM1 After the installation of the MCR is done the installation of HAWASSI VBMI1 can be performed by double clicking the installer of HAWASSI VBM1 setup_HAWASSIVBM_v 1 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 VBM1 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 1s 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 VBM1 are provided and placed in My Document HAWASSI_VBMI1 3 12 Page HAWASSI 4 GUI s of HAWASSI VBM1 For ease of operation HAWASSI VB
45. ison of the dispersion relations of Boussinesq type models by Madsen amp S rensen 1992 Nwogu 1993 and the VBM with parabolic profile by Klopman et al 2010 with the Airy linear theory dispersion relation as function of kh At the right is the comparison of the dispersion relation of VBM with 1 2 3 Airy profiles and a parabolic profile with the Airy linear theory 2 4 Units and Computational grid HAWASSI VBM1 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 VBM1 uses a Cartesian coordinate system with a uniform grid in the only horizontal dimension 2 l1 Page HAWASSI 3 Installing HAWASSI VBM1 software HAWASSI VBMI1 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 HAWASSI VBMI1 can be done in two main steps the installation of MCR and the installation of HAWASSI VBM1 itself 3 1 System requirements HAWASSI VBM1 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 VBMI1 package v 1 1 requires MCR
46. l domain can be defined in Xspace as shown in Figure 4 21 by specifying the most left and right boundary Xwest and Xeast respectively and the spatial discretization dx HAWASSI VBM1 uses an equidistance grid with grid size dx The software will automatically estimate a value for dx to be used based on the chosen influx signal and or initial condition If the user input value of dx is 4 22 Page HAWASSI larger than the estimated value dx a pop up warning dialog will notify a suggestion for the dx as illustrated in Figure 4 22 Kwest dx Xeast h Space Domain of the simulation rn Figure 4 21 Domain of computation Figure 4 22 Suggesting for the length of dx Warning Choose an appropriate dx not too small and not too large for the simulation If the dx 15 unnecessarily small the computation will take longer time or the solution may not converge because it cannot satisfy the given tolerance in the odesolver which will make the computation to stop or fail The Bottom profile bathymetry to be used in the simulation can be chosen in the section Bathymetry see Figure 4 23 There are four choices for the bathymetry type For Flat bottom the depth will be constant everywhere and the user will be asked to specify the Depth input see Figure 4 24 Bathymetry Type Flat bottom Flat bottom Select bathymetry type Depth U
47. l influx is created to save the influx signal if it is a signaling problem as a txt ASCH file in the folder Output ProjectName The option is checked by default Time Stepping Time partition 1 Divide the simulationcalculation into several steps a Spit Output Data Puta natural number Split the output simulation data at t_end to be hotstart conditions into the number of time partition in the next simulation _ Output data will be DATA _ Projecth 4 27 Page HAWASSI Time Stepping Time partition 1 Hotstart M Split Output Data Save signal influx Save the influx si Figure 4 32 Time Stepping in Advance Setting HAWASSI VBMI1 can calculate internal flow quantities such as pressure vertical and horizontal velocities as well as accelerations during a certain time interval and vertical discretization To activate this calculation the user should check Internal Flow in the Advance Setting GUI Further details about internal flow will be described in section 0 To save all information that have been configured in the Advanced Settings GUI click Save to cancel all changes click SS 4 1 9 Running Simulation After all parameters and input have been set the user can start the simulation This step is divided into 3 main steps 1 Preparation When the button is pushed the software will extract all input data from the Main GUI and the Advanced Settings GUI
48. lation file DATA TC1 intro flat tL init t end AK west X east bottom dz Z top Time Interval 0 5 100 Horizontal Interval 400 2000 Vertical Interval 14 43 0 778 1 14 Input file loaded Post Processing Velocity Acceleration Horizontal C Vertical Horizontal C Vertical tint t end wedi Kwest Xeast X Z axes bottom top tz axes Figure 4 53 Post Processing GUI for internal flow There are two main panels in the Post Processing GUI for the internal flow i e Interior Flow Calculation and Post Processing In the Interior Flow Calculation panel the user should specify the time interval Time Internal during which the results are desired the Horizontal Interval and the Vertical Interval The procedure to fill in these parameters are illustrated in Figure 4 54 The calculation 4 38 Page HAWASSI process can then be started by pushing Calculate a comment bar in the lower part of the GUI will notify e g Calculating Internal Flow followed by Finished when the calculation is done see Figure 4 55 Z space domain m for the Internal Flow Calculation Should not lower than initial depth sign for under water Interior Flow Calculation Input file uploaded file data loaded v DATA_TC_Harmonicy Simulation data to be processed for Internal Flow t_initt_end X_west X_east z_bottom dz z_top Time In
49. mine the type of vertical potential profile to be used which will determine the dispersive quality of the model as illustrated in Figure 2 1 The user can select the model type in the option Dispersion see Figure 4 5 In the Optimized VBM the software will calculate automatically the optimized wave number s x to be used in the Airy profile s based on the problem to be solved signaling problem and or initial value problem VBM with User defined settings the user can choose frequencies in Hz related to the wave numbers to be used as input in the Airy profile by filling in Airy UserDefined value Hz as shown in Figure 4 6 4 15 Page dH AWASS Suggestion For less experienced users it is advised to perform a trial simulation with the VBM parabolic profile to get an idea how dispersive waves will propagate When the model s dispersion quality 1s too poor for the wave to be simulated the Optimized VBM with 1 profile and after that with 2 profiles can be tried Optimized VBM with 3 profiles is suggested to be used only when simulating very short waves and or when the wave spectrum is very broad Warning 1 The more Airy profiles are being used the more computation time will be needed For rather simple waves simulations with 2 and 3 Airy profiles will approximately cost 1 5 and 2 25 times more CPU time than needed for one Airy profile respectively For complicated waves very steep broad spectrum the calculation time may
50. more important is 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 Warning The nonlinear calculation will take 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 Wave model Dispersion Optimized WBM 1 profile 4 Nonlinearity linear Linear Linear Nonlinear Weakly nonlinear f Figure 4 7 Nonlinearity of VBM 4 1 4 Initial Conditions Initial conditions for the surface elevation 7 and the surface potential 4 have to be specified as shown in Figure 4 7 Choosing option Zero means that the surface elevation 7 and surface potential have 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 three columns the first column specifies the discretization of the horizontal x coordinate and the second and third columns are the data of 7 and respectively If only two columns
51. nderwater moonen Sloping bottom linearly increasing User defined Figure 4 23 Bathymetry type For the bathymetry type Underwater mountain and Sloping bottom linearly increasing the user has to specify additional parameters such as shallowest depth the shallowest depth the of the Underwater mountain or Sloping bottom and Xo Xend the start and end x location of the Underwater mountain or Sloping bottom See Figure 4 25 till Figure 4 27 4 23 Page HAWASSI Bathymetry Figure 4 24 Setting up Flat bottom Bathymetry Tipe Sloping bottom 7 K The shallowest depth rr of Wo wend the underwater mountairsloping bottom Depth rm N Initial depth rr starting and ending x lacations rm of the Underwater mountairnsloping bottom Figure 4 25 Setting up Sloping bottom amp Underwater Mountain shallawest depth still water Shallowest depth Xo Xend Figure 4 26 Bathymetry profile for the option Underwater Mountain 4 24 Page 4 HAWASSI still water Shallowest depth Xo Xend Figure 4 27 Bathymetry profile for the option Sloping bottom For the option User defined bathymetry the user has to specify a file txt dat asc or mat that contains the bathymetry profile through a dialog box The data should consist of two columns the first column is the discretized equidistant x coordinate the second column contains
52. nown 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 and for intercontinental telecommunication through cables e Ocean resources of food and minerals are only at the start of discovery profits from wind parks and harvesting of wave energy in coastal areas is 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 aspec
53. o_flat_InternalFlow_iterl mat 5 980 KE J Output a l HOTSTART_TC1_intro_flat_InternalFlow_t300s tc 17 KB di TC1_intro_flat_InternalFlow B ma INFLUX_TC1_intro_flat_InternalFlow tet 12 KB ates ng El INPUT_TC1_intro_flat_InternalFlow mat 401 KB di an i F LOG_TC1_intro_flat_InternalFlow log 4 KE 2 PSLINTFLOW_TC1_intro_flat_InternalFlow_iterl mat 2 034 KB a 9 items Figure 4 36 Output files 4 30 Page HAWASSI Table 4 4 Output files folders Output Files Folders INPUT_ ProjectName mat Input data for the simulation a result from Preparation Dc DATA_ ProjectName _iter 1 mat Output data from the simulation a result from RUN The iter 1 indicates it is the simulation data of the i th output data i LOG_ ProjectName log A log file of the case ProjectName If the ProjectName is not changed the log file will be added with a log of the new simulation not replacing the log file INFLUX_ ProjectName txt The influx signal data txt or ASCII file that is used in the simulation only available when the user has checked the Save influx signal checkbox in the Advanced Settings GUI 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 has checked the Hotstart checkbox in the Advanced Settings GUI PSL INTFLOW_ ProjectName mat Stored variables that are
54. orizontal accelerations e Time partitioned simulation is possible to reduce computer hardware requirement e Project examples with harmonic focusing and irregular wave above bathymetry The model has been compared with series of experiments in hydrodynamic laboratories see Section 6 2 10 Page 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 the Boussinesq type model that is based on Madsen amp S rensen 1992 adopted in MIKE21 Boussinesq Wave by DHI and Nwogu 1993 adopted in BOUSS2D SMS by Aquaveo Both models have dispersion accuracy up to kh 3 14 which means these are capable to simulate waves with wavelength more than 2 times the depth The dispersion quality of HAWASSI VBM compared with these other two Boussinesq type models 1s illustrated in Figure 2 1 Tk o g 2k 4 o g Airy VBM 1 profile 5 YBM 2 profiles 7 YBM profiles VBM parabolic profile 6 2 4 6 8 10 12 A 5 10 15 20 25 30 kh kh Airy VBM parabolic profile amp Madsen amp Sorensen 92 7 Nwogu 93 a 0 39 Figure 2 1 At the left a compar
55. rder 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 VBM1 provides this facility the user can prescribe the length of a so called nonlinear adjustment zone The length has to be specified related to the peak wavelength see Figure 4 17 Advice For nonlinear simulation the length of Nonlinear adjustment should be at least 2 peak wavelengths For influxing rather high waves with respect to the depth the length should be more than 2 peak wavelengths Nonlinear adjustment aN A For nonlinear simulation provide the length of nonlinear adjustment zone default 2 wavelengths Figure 4 19 Nonlinear adjustment 4 1 6 Simulation Time The user can specify the time interval for the simulation in the box Time see Figure 4 20 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 end time of the simulation Time discretization dt is not the time discretization for calculation of the time integration The HAWASSI VBMI uses automatic internal time stepping in the matlab odesolver tatart dt tend X Time Time of the simulation 5 Figure 4 20 Simulation time 4 1 7 Geometry Bathymetry and Boundary Conditions The computationa
56. roup speed of half the frequency difference Nonlinear effects may lead to drastic changes by nonlinear mode generation By vemi Suggestions File Setting Log Tool Help o Enjoy and study the eee beat pattern for linear C Users DiditiDesktop _VBM1 simulations mee how Wave model __ TT Geo bath amp boundary the composite wave NT Xwest dx Xeast Dispersion Optimized VBM 2 profiles ee travels through the aie ane ii kada lt r Nonlinearity Weakly nonlinear Bathymetry fixed form of the l shallowest depth envelope that travels TYPE Sloping bottom 15 slower casual Depth m 30 7170 5 7470 5 o Imagine what happens Initial conditions if the amplitudes are 2222S not the same BC West Damping investigate with a g Time Signal simulation location o o Nonlinear effects See gg amers y f_int f_end depend on the Fitter freq band Hz 0 02 0 2 amplitude investigate Sa Ki l PEE Sees and be surprised by Method fea Uni arisan andina results of simulations Nonlinear adjustment 3 User s note TC3_Bichr_305007_2prof See the mode generation that takes SS ae 0 1 1800 place low and high frequency generated modes compared to spectra of linear simulations by comparing spectra o Compare result of simulation with the measurement data from MARIN for the given default value of this test case Adytia and Van Groesen 2012 o Contemplate on possible reflections from
57. sa aia aa e a a NG a aaa a NGE a et aa NR A GB D AA Ba a aaa Ga a aa aa Ba a nabi alah bab ah 4 13 4 1 ETE E aaa ap aaa a maa La AE AAT a ag a gek KAG EE aba a aaa ETA a aga baka aa aa aa aha 4 14 4 1 1 Work PUCCIO aaa sasa aine aa aa a agak ak aan aa a akan a aa a Ba Gak baana a a ak aa gha aaa 4 15 4 1 2 Project Name amp User SIN UG aan atan apanane ia aban aja ana ete ati pagaen da aana seeded 4 15 4 1 3 Wave Mode Da andaga san nana aa A no ce ie neces aaa san anana ada sea nana sang andana aaa a 4 15 4 1 4 BCPC ondi OUG 55432 KN ANA E NE KN ANE E TA E KN 4 17 4 1 5 MAAS rn yas cs ANNA E A ER AKAN MEN NA TAK TAN A A TA RA ER EN 4 19 4 1 6 AG LAOK TTNG aane aana aaa aa aaa an nga naa aa aa Ba aaa BA aaa a na aab aaa an ee ag aaa a aaa aa ngak gana aia 4 22 4 1 7 Geometry Bathymetry and Boundary Conditions cccccccccceececceceeeeeeeeesseseeeeeeeeeees 4 22 4 1 8 PN YANG SEHENG ia daa GAN NG a aa aaa a a EEN aa a EN a Na aa a WANGAN eda Ga aa ga 4 20 4 1 9 R nmne SAE sasa kane ajian a pinge dana aaa ban ba e a A gigi dadan ajaga ge kani daa gada napa ag aga ga 4 28 AM MUU TEN aaa aaa an ab aga e a gara a ana a AAN ana Wa aana SA Sa ag a ab abaan TE Gak a ana a aa bah 4 30 4 2 Post Processing GUI for Wave Dynamics cccccccccccecceceeceeeaeeeseseeseeeeeeeeeeeeeseeeaaaaaaeesssees 4 32 4 2 1 EAEan UE TO i AE NE EN NANGEN NG ANENG AA MAEN ANENG LEN NE AE EN RN 4 33 4 2 2 TONEN D sasab eka paia a aaa niaga ef oe a actrees E de ne
58. sinesq character of the codes e The way of approximating K is based on Dirichlet s principle for the boundary value problem in the fluid domain By restricting the set of competing functions in the minimization an approximation of K is obtained The variational derivative 04K y Oy 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 x The choice of these values determines the dispersion 2 9 Page HAWASSI relation Given the spectrum of the influx signal or initial profile of the case under investigation the values of Kn 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 7 a
59. system of elliptic equations has to be solved for the amplitudes of the Airy functions Advice 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 VBMI accounts for the following physical phenomena of waves in 1 HD horizontal direction long crested waves e Wave propagation dispersion and shoaling e Nonlinear wave wave interactions Features of the software include e The quality of dispersion is optimized for the specific wave problem to be simulated which makes it possible to simulate deep ocean waves or very short waves kh 15 or more as well as infragravity waves e Various methods for wave influx using embedded sources and an adaptation zone for influx of highly nonlinear waves e Use of efficient damping zones and partially reflecting walls e Interior Flow Module for the calculation of interior fluid velocities fluid accelerations and the dynamic and total pressure at user defined positions between surface and bottom Facilities of the software include e GUI Graphical User Interface for input of wave characteristics and model parameters e GUI for post processing of the output wave simulation e GUI for post processing of internal flow calculation such as the pressure vertical amp horizontal velocities and vertical amp h
60. terval 50 100 Horizontal Interval 14 54 Vertical Interval 0 40 0 006 0 21 for the Internal Flow Calculation Starting and Ending time s p 3 sl for the Internal Flow Calculation p Should not exceed simulation domain k Calculate the Interior Flow should not exceed the time which data available Figure 4 54 Panel for Interior Flow Calculation Calculating Internal Flow 23 25 14 seconds remaining relative Figure 4 55 Notification in the Post Processing GUI for internal flow When finished the data are stored as INTFLOW ProjectName mat in the Output ProjectName folder and the Post Processing panel will automatically load the data Quantities to be shown can be chosen in the panel such as Dynamic Pressure Velocity amp Acceleration in horizontal and vertical directions see Figure 4 56 Post Processing Simulation data to be Input file uploaded file processed for Post Processing loaded Y Data INTFLOW TO b Pressure Velocity Acceleration Plot the Vertical Velocity EX Dynamic piii Tyertical IN orizontal T wertical Plot the Dynamic Pressure _ pjot the Horizontal Velocity Plot the Horizontal Acceleration Figure 4 56 Panel of Post Processing in the Post Processing GUI for Internal Flow In a Plotting sub panel after specifying ranges results can be plotted in 1D and 2D see the left part of Figure 4 57 In the Animation sub panel an animation can be shown of th
61. the end of wave tank and other possible disturbances that are not simulated 5 52 Page HAWASSI 5 5 Test Case 4 Irregular Wave Propagation above a Sloping Bottom The same lay out as Test Case 2 and3 now for irregular waves Illustration of a hotstart at the initial time of the simulation the basin is already filled with waves that are given the correct velocity to propagate but new waves are influxed from a position in such a way that the combined initial value and influx problem produce a smooth continuation This shows the possibility to time partition one long time influx problem into several shorter time simulations Suggestions By vemi o On the flat part before simi a ana the slope the Working Directory CAWsers DiditiDesktop _VBM1 complicated wave forms of the irregular wave can be observed Wave model i Geo bath amp boundary l Xwest dx Xeast all su g gestions from Dispersion Optimized VBM 2 profiles Xspace 02810400 Test Case 2 apply here Nonlinearity Weakly nonlinear Bathymetry l shallowest depth to investigate observe 7 TYPE Sloping bottom 15 effects related to the rs E sloping bottom MTA s Depth m 30 7170 5 7470 5 Initial conditions can be used for easy investigation User defined o Compare MTA s of a linear and nonlinear Time Signal location simulation a o Compare result of finit f end simulation with the Filter freg band Hz
62. therlands 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 Software Although much care is taken to arrive at trustful results of simulations with HAWASSI Yayasan AB cannot be held responsible for any result of simulations obtained with the software or consequential actions or calculations that are based on the results e g because of possible bugs wrong use of the software or other causes 1 7 Page HAWASSI 1 Introduction This document is the Manual of HAWASSI VBM1 software that serves as a guide for using and running the software HAWASSI VBM1 simulates phase resolved waves in Horizontal Direction 1HD long crested waves as are generated in wave tanks to simulate on scale coastal and oceanic waves above flat and varying bathymetry and with partially reflecting walls and damping zones Section 2 describes the underlying model the peculiarities and the capabilities of the software together with the features of the software it is advised to read this Section before continuing to the rest of the manual Section 3 provides a description for the step by step installation process of the software A condensed description to handle the software regarding GUIs and input output parameters is given in Section 4 Section 5 gives a short tutorial with test cases for
63. to compare simulation results with experimental data After pushing the RUN button plot s will automatically pop up as illustrated in Figure 5 14 BJ Figure 2 Validation TENGGA Meas x 50 Sim Figure 5 14 Signal comparison between the experimental data and the simulation 5 48 Page HAWASSI 5 2 Test Case 1 Wave Propagation above a Flat Bottom This introductory example of waves above a flat bottom is to illustrate the 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 setting is pre programmed the details are given in the GUI plot below E VBM1 File Setting Log Tool Help Working Directory CAUsers DiditiDesktop _ VBMi Wave model Geo bath amp boundary Xwest dx Xeast Dispersion Optimized VBM 1 profile ASpace 400 4 2000 Nonlinearity Linear Bathymetry Type Flat bottom Depth m 30 Initial conditions Boundary BC West Damping Zl Time Signal location Tpi Amp BC East Damping Harmonic g D 121 finit f end Fitter freq band Hz Generation Project Name TC1_intro_flat Method Area Uni User s note Wave propagation above flat bottor Suggestions o Look at the log file to find the calculated wave length How many points will be in one wavelength for
64. to the sloping bottom MTA s can be used _ to investigate such aspects as o Investigate the shoaling amplitude increase above shallower depths o Reflected waves become more pronounced for steeper slopes File Setting Log Tool Help Working Directory C AWsers DiditiDesktop _ VEMI Wave model Dispersion Optimized VBM 1 profile Nonlinearity Li Initial conditions Time Signal Harmonic Filter freq band Hz Generation N Method dhalan ndadi t start dt t_end Time 0 0 5 400 Geo bath amp boundary Kwest dx Xeast MSpace 450 4 2500 Bathymetry shallowest depth TYPE Sloping bottom 15 Xo Xend Depth m 30 1000 1400 Boundary BC West Damping A BC East pamping ai Project Name TC Intro Slope User s note Wave propagation above a sloping RUN STOP 5 511Page da HAWASS 5 4 Test Case 3 Bichromatic Wave Propagation on a Sloping Bottom A linear bichromatic wave is the sum of two monochromatic waves of the same amplitude If the frequencies of the waves are close together a characteristic beat pattern results as a consequence of dispersion the individual waves translate with their own phase speed which is different That causes a wave profile that can be seen as a modulated wave that travels with a speed that is the phase speed of the averaged frequency the envelope travels with a different speed with the smaller speed that is approximately the g
65. ts 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 Broer 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 6 Page HAWASSI 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 Ne
66. y p kg m 4 26 Page HAWASSI Advanced Settings Log Tool Help Parameters K g m s 2 p ka m 3 Adwanced Settings Figure 4 30 Accessing Advanced Settings GUI Time Stepping Figure 4 31 GUI of Advanced Settings Hotstart Split Output Data In Advanced Settings it is possible to 7 Save signal influx modify settings in the Time Stepping such as Time partition Hotstart Save signal influx and Split output data see Figure 4 32 Time partition is an option to Internal Flow Pressure divide the simulation calculation into eee Sie several steps the default is 1 such that at Horizontal Velocity the end of each step the solution is stored at the hard drive of the user s PC this reduces memory RAM requirements and makes in Horizontal Acceleration Vertical Acceleration a long time simulation the simulation data directly available after each iteration in the time partition When the Split output data option is checked the software will split the simulation data into the number of specified Time partition After the simulation the last iteration will automatically be loaded into the Post Processing GUI For saving the last state of the surface elevation 7 and surface potential after each iteration in time partition as a txt ASCII file the option Hotstart should be checked is default The option Save signa

Download Pdf Manuals

Related Search

Related Contents

Philips D 6280 User's Manual  National Instruments NI 78xxR User's Manual  資料4 安全チェックシート  Apricorn EZ Upgrade Universal & SATA - 160GB  MY MCL  PEB-7603 User`s Manual R1.0  Descargar Ficha Técnica  Sharkoon Side Panel Open  （注）GX形に関する記述は含まれていません。 発行年月  TSI MITA 8120 Mask Integrity Test Accessory User Manual  

Copyright © All rights reserved. Failed to retrieve file