fcmubore.m

波浪数值模拟

dx = 0.1;
nx = 1001;

xu = [0:nx-1]'*dx;
xeta = [1:nx-1]*dx + dx/2';

h=1;
eta0 = 0.1;
U0 = 0.0976;

%eta0 = 0.2;
%U0 = 0.1909; % eta0= 0.2;

%eta0 = 0.3;
%U0 =  0.2804; % eta = 0.3;

a = 5;

Ucr = 0.5*U0*(1 - tanh(xu/a));
Ue = 0.5*U0*(1 - tanh(xeta/a));
Ecr = Ue + 0.25*Ue.^2;

h = 1;

subplot(2,1,1)
plot(xeta,Ecr)
ylabel('\eta (m)')

subplot(2,1,2)
plot(xu,Ucr)
ylabel('U (m/s)')