Write matlab codes to determine the wave characteristics
Course:- MATLAB Programming
Reference No.:- EM131057630

Assignment Help
Assignment Help >> MATLAB Programming

Question 1: 

A wave with the deep-water characteristics of (T = 6.8 s and H = 5.0 m) is propagating toward a coast, where the water depth is 6.2 meters. At the coast, the surface elevation can be described as

η(t) = (H/2) cos(ω/t), and

the horizontal and vertical water particle velocities of the wave can be calculated with

u(t) = (ωH/2)(cosh(k(z+d))/sinh(kd))cos(ωt)

v(t) = -(ωH/2)(sinh(k+(z+d))/sinh(kd))sin(ωt)

where g is the gravity, H is the wave amplitude, k is the wave number, d is the mean water depth, ω is the angular frequency, z is the location below the mean water level (MWL).


Write MATLAB codes to-

(i) Determine the wave characteristics (i.e. wave height, wave length, wave number, wave celerity) in deep water and at the coast (use the Dispersion relationship and fzero function in MATLAB);

(ii) Calculate the horizontal and vertical water particle velocities for t varies for 2 wave periods and z varies from 0 to -d (with an interval of 1 m).

(iii) Generate figures to present the time series (t varies from 0 to 2T) of the water particle velocities at the locations of the trough, z = -H/2 m and at the bottom. (iii) present the vertical profile of the magnitude of horizontal water particle velocity from the MWL to the sea bottom.

(iv) generate a similar graph as below for one wave period and at the mean  water level, to show the relationship between the surface elevation, the horizontal and vertical  water particle velocities.


Question 2-

Consider a 5km stretch of coast oriented in the north-south direction with the ocean to the east. The predominant wave direction is from the east-south-east. At the southern end the typical breaker height is 1.4 m and the breaker angle is 10o. At the northern end, the breaker height is 1.45m and the breaker angle is 12o. The breaker parameter γb = 0.8.

The beach profiles along the section are similar with slopes near the break point of 1/40. The sand is made of quartz (s = 2.63, p = 0.28) with a median grain size of 0.22mm and measureable seasonal bed level changes are restricted to depths less than 6 metres. The berm height is 3 m AHD.


The average shorenormal sediment transport rates (Qx) for 1993-2015 were saved in the data file "sediment.xlsm". There are no sinks and sources for sediment transport in the control domain. The erosion rate (metres of shoreline retreat rate) can be calculated using (see details in Coastal Process module lecture notes)

(1-p)(hc + B) ∂xs/∂t = -Qx + ∂Qy/∂y + Qsink - Qsource

Qy = (K/16(s-1)√γ) √gHb5/2 sin2θb


p is the sediment porosity,

xs is the shoreline coordinate

Qsource is a sediment input (e.g. river discharge, beach nourishment)

Qsink is a sediment loss (e.g. dredging)

K ≈ 0.77 is an empirical coefficient which has a weak dependence on grain size.

s is the specific weight

Hb is the breaker height

γ is the breaker index

θb is the wave crest angle at the break point

all other variables are as defined in the figure.

The aerial photographs of the stretch of the coast from 1993 to 2015 indicate the shoreline location, xs which were saved in the file "shoreline.xlsm".

Use Matlab to develop a numerical model to calculate the location of the shoreline from 1993 to 2015, i.e. xs. (assume xs = 0 in 1993)  Evaluate the accuracy of the model.

Attachment:- Assignment.rar

Put your comment

Ask Question & Get Answers from Experts
Browse some more (MATLAB Programming) Materials
Set up the quasi-steady equations of motion for a pitching and plunging 2D flat plate. Calculate natural frequency and damping ratio. Plot the variation of the natural frequen
Write a MATLAB code to find the stator core dimension D and L and number of slots and winding turns for a delta connected squirrel cage induction motor then start the specif
Compute and plot separate correlation functions of white noise data and of sinusoidal data. Do this using both the biased and unbiased autocorrelation functions. Are the resul
Prepare a program in Matlab to add two numbers from two arrays using try catch construct. Initialize two arrays of different size, and then prompt the user to input indexes
Performs curve-fits of all the data - Determines what type of loading has been applied based on which curve-fit is the best. 6.Demonstrates that the chosen curve-fit is the
A building of 10 floors is used for parking cars. The area of each floor can be used for parking 100 cars. Assume that the building spaces as elements of a matrix, write a Mat
you will develop a policy statement and the processes and procedures that can be used to drive the design, implementation, and maintenance of a wireless network in a SOHO (
Write a program to compute A. Test it with P = $55,000 and interest rate of 6.6% (i = 0.066). Compute results for n = 1, 2, 3, 4, and 5 and display the results as a table with