206 Computational Modelling in Hydraulic and Coastal Engineering
index=k
end
[i,j]=meshgrid(1:1:nx,1:1:ny);
df1=fDep(:,3);
df2=reshape(df1,70,40);
df2=df2';
An=ones(size(df2));
idx=find(df2== 0);
An(idx)=0;
figure;
pcolor(An);
colormap(jet);
% colormap(gray(3))
shading flat;
hold on;
i=1:nx;
j=1:ny;
Cplot=cmax(i,j);
% hold on;
contour(i,j,Cplot')
% colorbar('vert')
PROBLEM 8.1
Solve the same problem by making the suggested modifications while keeping the rest of the data constant:
1. Place the contaminant point source at location x = 225 m, y = 190 m.
2. Remove the horizontal arm of the breakwater.
3. Replace the point source by a line source extending from point (50 m,
25 m) to point (50 m, 75 m).
4. Remove the left vertical section of the breakwater.
5. Change the decay coefficient to 0.01 s –1 (decay) and –0.01 s –1 (growth).
Run the simulations, compare the results and comment on the changes
observed regarding the spreading of the contaminant plume.
Example 8.2
This application is based on the bathymetry and hydrodynamic data
from Example 5.5 and investigates the advection and diffusion of a
pollutant point source in the Thermaikos Gulf. The results are indicative of the contaminant transport under the forcing of a northwest wind
of constant speed of 9.94 m/s. Other data used for the simulation are
Concentration (strength) of point source = 1 Kg/m 3
Decay coefficient = 10 –5 s –1
Location of the point source: x = 24 km; y = 24 km