Message Boards Message Boards

3 Replies
2 Total Likes
View groups...
Share this post:

Inverting dispersion relation

Posted 10 years ago
Hi everyone,

I'm trying to manipulate the dispersion relation for capillary-gravity waves, i.e.
?² = Tanh[kh](gk + ?k³/?)

in order to express k as a function of h. My first pick was to express h as a function of k
h = (1/k) ArcTanh [?² (gk + ?k³/?)^(-1) ]

and then try to invert it with Findroot, typing something like this :
h[x_] := (1/x) ArcTanh[((70 \[Pi])^2)/ (9.81x + (0.0206 x³/950) )]
hmax = 0.01;
hmin = 0.001;
dh = 0.0001;
Sol = Table[{{b -> i}, FindRoot[i == h[x], {x, 0.001}]} // Flatten, {i, hmin, hmax, dh}];

(I replaced some parameters with their experimental values) but this seems to stumble upon some convergence error. I'm not really familiar with Mathematica, and I found nothing that could help me in the doc. Does anyone have an idea on what I should try ?

Thanks in advance for your replies !
POSTED BY: Olivier Mathieu
3 Replies
>>Your Findroot line worked quite nice
Thanks, but it's not mine in fact, I've googled in the past when I had a need for it emoticon

POSTED BY: Ivan Morozov
Posted 10 years ago
Your Findroot line worked quite nice, although my solution seems to become complex much quicker that I expected, but I can figure out a physical explanation for this.

My deepest thanks !
POSTED BY: Olivier Mathieu
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract