CH1 - Calculating the offshore wave climate



MATLAB Assignment - Calculating

the offshore wave climate

For many coastal engineering studies and projects, a time-series of wind measurements is used to calculate a time-series of wave heights.

For tutorial question 16, nonograms (graphs from which wave height and wave period are estimated) were used. Obviously though, the method used in 16 is impractical for large series of measurements.

Your task is to write a program to hindcast a time-series of wind measurements, using MATLAB.

Instructions:

You have been given a time-series of hourly wind measurements (velocity and direction) from Falsterbo Peninsula in the south of Sweden. Assuming the wind field over the fetch is homogenous and described by the measured winds at Falsterbo, and ignoring accumulation of wave energy, calculate the wave climate at Skanör harbour. Also assume that you are hindcasting for deep water.

To do this, write a generalized program in Matlab to calculate hourly values of Hs and Ts generated by a measured time-series of winds at a location for various fetches, F. By generalized, it is meant that the program may be used for any input wind measurement time-series.

Preparation:

Read pages 3-24 to 3-32 of the Shore Protection Manual. These pages discuss adjusting the measured wind speed to a standard 10 m wind speed and how to calculate the wind stress factor, UA, in which the wave hindcasting equations are presented.

Read pages 3-44 to 3-66 (it’s not that much! Most is diagrams). These pages discuss the formulae for hindcasting of waves from the wind stress factor, UA.

Do tutorial question 16. You could also start looking at steps 1 to 4 for ‘starting out

hints’ if you want to get a head-start.

[pic]

Input Data is for Falsterbo Pensinsula, Sweden. You should calculate the wave heights for deep water outside of Skanör, which is a coastline oriented facing approximately WNW. Fetches were calculated to be:

[pic]

Calculating the offshore wave climate: Writing a

Matlab routine.

This sheet should provide you with a starting point for the project:

1. Write down what you want your program to output. In what format do you want the output? (scalar, vector, matrix? Will you graph the output?)

2. What is the input data required? In what format is the input data (scalar, vector, or matrix?)

3. Write down in detail the steps you will have to go through to hindcast the waves (e.g. formatting input data, adjusting the wind speed, determining the fetch which is dependent on direction, calculating H and T, determining if it is fetch or duration limited etc. ). You should also write down the equations you will use for each step. Show this to Caroline, and she will let you know if you missed anything prior to starting programming.

4. You are now ready to start programming!

Some helpful hints for programming:

1. Set up your program by following the steps outlined above, especially what you wrote down for (3).

2. Use commenting throughout the program. To insert a comment (a line not to be read by the program), use % and then follow with your comment. This not only makes it easier for me to mark your work, it makes it easier for you to find mistakes in the program.

3. Using ‘;’ at the end of a line prevents the output of that line being displayed as the program runs.

4. Use the debug function within Matlab (use run in debug mode). If the program will not run, step through the program and make sure it is calculating each line as you expect it to.

5. Besides doing normal calculations with scalars and vectors (don’t forget Matlab differentiates between dot and cross products when working with vectors!), you will probably need to use “for…end”, and “if….elseif…..end” structures.

Calculating the offshore wave climate: Data

You are given a time-series of wind measurements at Falsterbo peninsula for 1996. The wind measurements are taken at a height of 12 m above sea level, over sea, and the temperature differences are not known. Each wind measurement is taken for a duration of 3 hours and the wind is assumed to blow steadily at the measured velocity and direction for those 3 hours.

Download the wind data file from the course web site.

To check that your program is giving the correct answers you can double-check the values for the dates in the table below:

|Row |Date |Time |H |T |UA |

|1 |1/01/1996 |1:00:00 |0 |0 |16.4 |

|481 |1/03/1996 |1:00:00 |0.24 |2.36 |4.25 |

|969 |1/05/1996 |1:00:00 |0 |0 |2.99 |

|1457 |1/07/1996 |1:00:00 |0.85 |3.88 |11.53 |

|1953 |1/09/1996 |1:00:00 |0.71 |3.61 |9.98 |

|2441 |1/11/1996 |1:00:00 |0 |0 |11.53 |

|2446 |1/11/1996 |16:00:00 |1.33 |4.63 |16.43 |

|2549 |14/11/1996 |13:00:00 |0 |0 |0 |

|2608 |21/11/1996 |22:00:00 |2.04 |5.34 |25.22 |

Report:

• Hand in a copy of your program as a text file.

• Make a similar table as the one above showing your own results for H, T, and UA at the dates/times listed above.

• Also, what was the maximum wave height and period for 1996. Write down which date and time at which the maximum wave height occurred.

• Finally, write a paragraph describing what the results are, and what assumptions might affect the results.

................
................

In order to avoid copyright disputes, this page is only a partial summary.

Google Online Preview   Download