The Tsunami model has need to fixed this problem

Question

aoming 2024-4-7，13:16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2103911-the-tsunami-model-has-need-to-fixed-this-problem

评论： aoming 2024-4-8，1:11

hi,

I found that the "https://ww2.mathworks.cn/help/symbolic/solve-a-partial-differential-equation-tsunami-simulation.html#d126e37561", it has the problem that soliton = @(x,t) A.*sech(sqrt(3/4*g*A/H)*(x/c2+t)).^2; ,but it not Conservation of Dimensions . I think that it should be soliton = @(x,t) A.*sech(sqrt(3/4*g*A/H^2)*(x/c2+t)).^2. And I want to confirm that. Thank you!

2 个评论
显示无隐藏无

Torsten 2024-4-7，14:10

在 MATLAB Online 中打开

The unit of

sqrt(3/4*g*A/H^2)

is sqrt(m)/s. I doubt this can be correct.

aoming 2024-4-8，1:09

Thanks! thank you.But I actually have an another problem,

the code :S = fft(soliton(-0.8*TimeScale*c2, linspace(0,TimeScale,2*(Nt/2)-1)))';

S = repmat(S,1,Nx);

I would like to ask where the initial position of the waves is？ In the -0.8*TimeScale*c2？ and why it is minus value？

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

the cyclist 2024-4-7，15:46

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2103911-the-tsunami-model-has-need-to-fixed-this-problem#answer_1437666

Disclaimer: I am not an expert in this!

I was not able to track down the exact (Goring & Reichlan) reference from the documentation, and therefore not the exact formula. But, I did find some related materials.

I'm guessing that you might be being fooled by having seen the formula where kappa is inside sech(), rather than g. Looking over some other reference, it seems to me that

kappa ~ g / H

which is where the "missing" dimension of H is hiding.

Another place that may have fooled you is that in the documentation, A is a dimensionless (divided by H) amplitude.

I hope that helps.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

aoming 2024-4-8，1:11

Thanks，cyclist!

But I actually have an another problem

the code :S = fft(soliton(-0.8*TimeScale*c2, linspace(0,TimeScale,2*(Nt/2)-1)))';

S = repmat(S,1,Nx);

I would like to ask where the initial position of the waves is？ In the -0.8*TimeScale*c2？ and why it is minus value？or the initial position is in the x_max = L2 + 12;？

请先登录，再进行评论。