Free surface flows 117
[i,j]=meshgrid(1:1:nx,1:1:ny);
df1=df(:,3);
df2=reshape(df1,70,40);
df2=df2';
An=ones(size(df2));
idx=find(df2==0);
An(idx)=0;
figure,pcolor(An); hold on;
colormap(bone)
shading flat
% Plot every interval;
% quiver(i,j,u',v','Color','k');
% Plot every second interval;
quiver(i(1:2:end,1:2:end),j(1:2:end,1:2:end),(u(1:2:end,1:2:
end))',(v(1:2:end,1:2:end))','Color','k'); hold on;
PROBLEM 5.4
Solve the same problem by making the following suggested changes while
keeping the rest of the data constant:
1. Change the incoming longshore velocity (left-hand boundary), so that
it linearly varies from 0 at the coastline to 0.4 m/s at the open-sea
boundary. Run the simulation and discuss the results as compared to
those of the constant longshore velocity.
2. Change linearly the water depth from 4 metres at the coastline to 8
metres at the open-sea boundary. Run the simulation and discuss the
results as compared to those of constant depth domain.
3. Remove the entrance jetty and conduct the simulation with only the
L-shaped breakwater. Compare the simulation results with those
from the original breakwater configuration.
4. Change the Smagorinsky parameter to 0.1 and 1.0. Run the simulations, compare and discuss the results.
5. Shorten the horizontal arm of the breakwater from 100 m to 50 m,
while keeping intact the rest of the configuration. Run the program
and compare the simulation results with those from the original
breakwater configuration.
Example 5.5
This exercise investigates the wind-generated flow pattern in the
Thermaikos Gulf, near the City of Thessaloniki, Greece (Figure 5.15).
The topography is given in the text file ThermD.txt. The data are provided for a rectangular domain of 44 km × 46 km, on a square grid of
2 km × 2 km. The inland ‘dry’ areas are defined by a zero elevation,
and the gulf water depths by a positive number for every cell center
(i,j). The radiation condition is applied to the open-sea boundary, and
[i,j]=meshgrid(1:1:nx,1:1:ny);
df1=df(:,3);
df2=reshape(df1,70,40);
df2=df2';
An=ones(size(df2));
idx=find(df2==0);
An(idx)=0;
figure,pcolor(An); hold on;
colormap(bone)
shading flat
% Plot every interval;
% quiver(i,j,u',v','Color','k');
% Plot every second interval;
quiver(i(1:2:end,1:2:end),j(1:2:end,1:2:end),(u(1:2:end,1:2:
end))',(v(1:2:end,1:2:end))','Color','k'); hold on;
PROBLEM 5.4
Solve the same problem by making the following suggested changes while
keeping the rest of the data constant:
1. Change the incoming longshore velocity (left-hand boundary), so that
it linearly varies from 0 at the coastline to 0.4 m/s at the open-sea
boundary. Run the simulation and discuss the results as compared to
those of the constant longshore velocity.
2. Change linearly the water depth from 4 metres at the coastline to 8
metres at the open-sea boundary. Run the simulation and discuss the
results as compared to those of constant depth domain.
3. Remove the entrance jetty and conduct the simulation with only the
L-shaped breakwater. Compare the simulation results with those
from the original breakwater configuration.
4. Change the Smagorinsky parameter to 0.1 and 1.0. Run the simulations, compare and discuss the results.
5. Shorten the horizontal arm of the breakwater from 100 m to 50 m,
while keeping intact the rest of the configuration. Run the program
and compare the simulation results with those from the original
breakwater configuration.
Example 5.5
This exercise investigates the wind-generated flow pattern in the
Thermaikos Gulf, near the City of Thessaloniki, Greece (Figure 5.15).
The topography is given in the text file ThermD.txt. The data are provided for a rectangular domain of 44 km × 46 km, on a square grid of
2 km × 2 km. The inland ‘dry’ areas are defined by a zero elevation,
and the gulf water depths by a positive number for every cell center
(i,j). The radiation condition is applied to the open-sea boundary, and