À quoi sert ce calculateur
Cet outil détermine le s-ième zéro positif des fonctions de Bessel de première espèce \(J_{v}(x)\) et de seconde espèce \(Y_{v}(x)\) (aussi appelée fonction de Neumann ou de Weber), pour un ordre réel \(v\). Notés \(j_{v,s}\) et \(y_{v,s}\), ces zéros interviennent partout en physique et en ingénierie : vibrations des membranes circulaires (peaux de tambour), conduction de la chaleur dans les cylindres, guides d'ondes électromagnétiques ou encore séries de Fourier-Bessel. Il s'agit d'un outil de mathématiques pures dédié aux fonctions spéciales : il est universel et ne dépend d'aucune région ni d'aucune unité.
Comment l'utiliser
Saisissez l'ordre v (un nombre réel, généralement compris entre 0 et 200) et l'indice du zéro s (un entier positif : 1, 2, 3, …). Le calculateur renvoie \(j_{v,s}\), la s-ième racine positive de \(J_{v}(x)\), ainsi que \(y_{v,s}\), la s-ième racine positive de \(Y_{v}(x)\). Par exemple, \(v = 0\) et \(s = 1\) correspondent au mode fondamental de vibration d'une peau de tambour.
La formule et la méthode
Ces deux fonctions sont solutions de l'équation de Bessel \(x^2 y'' + x y' + (x^2 - v^2)y = 0\). \(J_{v}\) est évaluée à partir de sa série entière convergente ; \(Y_{v}\) repose sur la relation $$Y_{v} = \frac{J_{v}\cos(v\pi) - J_{-v}}{\sin(v\pi)},$$ le cas des ordres entiers étant traité comme une limite numérique. Pour localiser le s-ième zéro, on part de l'estimation asymptotique de McMahon, on encadre la racine, puis on l'affine par dichotomie jusqu'à convergence. Comme \(Y_{v}(x)\) tend vers \(-\infty\) lorsque \(x \to 0\), la recherche débute en un petit \(x\) positif afin d'éviter la singularité logarithmique.
Exemple détaillé
Pour \(v = 0\) et \(s = 1\) : le premier zéro positif de \(J_{0}(x)\) vaut \(2{,}4048255577\), et le premier zéro positif de \(Y_{0}(x)\) vaut \(0{,}8935769663\). Pour \(v = 1\) et \(s = 1\) : \(j_{1,1} = 3{,}8317059702\) et \(y_{1,1} = 2{,}1971413260\).
Que signifient les zéros dans les applications
Les zéros des fonctions de Bessel ne sont pas une simple curiosité abstraite — ce sont les valeurs propres discrètes qui émergent chaque fois qu'un problème d'ondes, de chaleur ou de potentiel est posé sur un domaine circulaire ou cylindrique. La condition aux limites force la partie radiale de la solution à s'annuler (ou sa dérivée à s'annuler) à la frontière, et cette condition n'est satisfaite que aux zéros \(j_{v,s}\).
Membrane circulaire vibrante (peau de tambour)
Pour un tambour circulaire idéal de rayon \(a\) avec un bord encastré, le déplacement se sépare en modes étiquetés par \((v,s)\), où \(v\) compte les diamètres nodaux angulaires et \(s\) compte les cercles nodaux radiaux. Les fréquences propres autorisées sont \(f_{v,s}=\frac{c}{2\pi a}\,j_{v,s}\), où \(c\) est la vitesse des ondes. Le ton fondamental utilise \(j_{0,1}=2.4048256\) ; les valeurs \(s\) plus élevées et les valeurs \(v\) plus élevées élèvent toutes deux la hauteur, et parce que les \(j_{v,s}\) ne sont pas des multiples entiers les uns des autres, les harmoniques d'un tambour sont inharmoniques.
Conduction thermique cylindrique
Lors de la résolution de l'équation de la chaleur dans un long cylindre avec une surface à température fixe, les fonctions propres radiales sont \(J_0(\lambda_s r/a)\) avec \(\lambda_s=j_{0,s}\). Chaque mode décroît dans le temps comme \(\exp\!\left(-\alpha (j_{0,s}/a)^2 t\right)\), donc le plus petit zéro \(j_{0,1}\) gouverne le profil de température qui décroît le plus lentement et persiste le plus longtemps. Un \(s\) plus grand donne des valeurs propres plus grandes et donc une décroissance plus rapide.
Fréquences de coupure des guides d'ondes
Dans un guide d'ondes métallique circulaire creux de rayon \(a\), les modes transverses magnétiques (TM) se coupent à des fréquences déterminées par \(j_{v,s}\) et les modes transverses électriques (TE) par les zéros de la dérivée \(J_v'\). Pour les modes TM, la coupure est \(f_{c}=\frac{c\,j_{v,s}}{2\pi a}\) ; c'est au-dessus de cette fréquence seulement que le mode se propage. Le mode TM le plus bas (TM\(_{01}\)) utilise à nouveau \(j_{0,1}\).
Série de Fourier–Bessel
Toute fonction raisonnable sur un disque peut être développée comme \(f(r)=\sum_{s=1}^{\infty} c_s\,J_v\!\left(j_{v,s}\,r/a\right)\). Les zéros mis à l'échelle \(j_{v,s}/a\) agissent exactement comme les nombres d'ondes \(n\pi/L\) d'une série de sinus de Fourier ordinaire, et l'orthogonalité des \(J_v(j_{v,s}r/a)\) sur le disque (avec poids \(r\)) permet de calculer les coefficients par intégration, \(c_s=\frac{2}{a^2 J_{v+1}^2(j_{v,s})}\int_0^a f(r)\,J_v\!\left(j_{v,s}r/a\right) r\,dr\).
Comment les zéros se déplacent
Pour un ordre fixe, augmenter l'indice \(s\) augmente \(j_{v,s}\) par étapes qui se rapprochent de \(\pi\) (l'oscillation devient quasi périodique loin de l'origine). Pour un indice fixe, augmenter l'ordre \(v\) repousse le premier zéro vers l'extérieur à peu près comme \(j_{v,1}\approx v + 1.8557\,v^{1/3}+\cdots\) pour grand \(v\), reflétant comment le terme centrifuge \(v^2/x^2\) dans l'équation de Bessel retarde le début de l'oscillation. Le point clé pratique : une complexité angulaire plus élevée (\(v\)) et plus de nœuds radiaux (\(s\)) correspondent tous deux à des valeurs propres plus grandes et donc à des fréquences plus élevées ou à une décroissance plus rapide.
Pour un exemple concret de tambour avec \(c=100\ \text{m/s}\) et \(a=0.30\ \text{m}\), la fréquence fondamentale est \(f_{0,1}=\frac{100}{2\pi(0.30)}\,(2.4048256)\approx 127.6\ \text{Hz}\). La première harmonique utilise \(j_{1,1}=3.8317060\), donnant \(\approx 203\ \text{Hz}\), donc le rapport harmonique-fondamental est \(j_{1,1}/j_{0,1}\approx 1.593\) — perceptiblement pas une octave ou une quinte, ce qui explique pourquoi les tambours sonnent non mélodiques comparés aux cordes.
FAQ
Pourquoi le premier zéro de \(Y_{v}\) est-il inférieur à 1 ? \(Y_{0}\) diverge vers \(-\infty\) à l'origine et traverse l'axe des abscisses vers \(x \approx 0{,}894\) avant d'atteindre son premier maximum ; son premier zéro est donc bien plus petit que celui de \(J_{0}\).
v peut-il être non entier ? Oui. La formule définissant \(Y_{v}\) s'applique directement à tout \(v\) non entier, et les ordres entiers sont traités comme une limite continue.
Quelle est la précision des résultats ? Les calculs utilisent l'arithmétique en double précision, ce qui offre environ 10 chiffres significatifs pour des valeurs modérées de \(v\) et \(s\).
