Free surface flows 113
Computer code 5.4
% Example 5.4 Longshore Current Past a Marina Breakwater
% ho = Water depth [m];
% u = Velocity in the x-direction [m/s];
% v = Velocity in the y-direction [m/s];
% cv = Longshore current velocity [m/s];
% ev = Eddy viscosity [m^2/s];
% sp = Smagorinsky parameter for horizontal eddy coefficient
[dimensionless];
% cf = Coriolis coefficient [1/s];
% fb = bed friction [dimensionless];
% Dd = Spatial step in both x and y directions [m];
% Dt = Time step [s];
% nx = Number of computational steps in the x-direction;
% ny = Number of computational steps in the y-direction;
% nt = Number of time steps;
clc; clear all; close all;
% Input data;
g=9.81;
ho=4;
cv=0.2;
sp=0.3;
40
35
30
25
20
15
10
5
10
20
30
40
50
60
70
Figure 5.14 Circulation pattern around the marina breakwater.