Nous allons maintenant considérer le deuxième cas important qui va intervenir lorsqu'on s'intéresse à la propagation en régime non-linéaire, c'est celui où on a, non pas un ensemble discret de fréquences dans le milieu, mais un continuum de fréquences, c'est-à-dire un spectre qui pourra avoir une certaine largeur, et qui va correspondre au cas d'une impulsion brève ou d'une impulsion ultra-brève. Et on va résoudre ce problème dans ce qu'on appelle l'approximation de l'onde lentement variable, qui n'est autre que l'approximation paraxiale appliquée au cas de la propagation non-linéaire. Donc on va à nouveau partir évidemment de l'équation de propagation avec dans le membre de droite le terme non-linéaire, et puis dans le membre de gauche l'équation de propagation linéaire habituelle, et on va effectuer un certain nombre d'approximations. Tout d'abord, on va supposer que le milieu est isotrope, ce qui nous permettra de négliger la double réfraction, et donc ce terme en gradient divergence E, ensuite, on va faire l'approximation scalaire, c'est-à-dire qu'on va considérer que on aura des faisceaux faiblement divergents, et donc qui vont tous se propager, enfin que notre faisceau va se propager selon l'axe z, et donc que le champ électrique, qui est transverse, sera dans le plan x y, et par exemple, polarisé dans une direction x, et c'est ce qui nous permettra de supposer uniquement la composante de ce champ selon, par exemple, l'axe x. Donc dans la suite, j'aurai uniquement une grandeur scalaire : le champ réel E. Ensuite on va évidemment utiliser le champ complexe qu'on avait déjà introduit dès la deuxième semaine, en considérant uniquement les composants de fréquences positives du champ réel E de t, et donc tout ça nous permet d'arriver à cette équation de propagation, qui est à peu près la même qu'en haut, simplement on a plus le terme de gradient divergence E, et puis on travaille avec des grandeurs complexes. Alors comme d'habitude, on a ici dans le membre de gauche l'équation de propagation linéaire, mais où on a rajouté un terme source correspondant à la réponse non-linéaire du milieu. Pour la partie linéaire, on va supposer que la réponse du milieu est caractérisée par son vecteur d'onde, k de oméga qui sera, comme d'habitude, égal à racine carrée de un plus khi un de oméga, donc qu'on écrira encore n de oméga fois oméga sur c. Donc, on a cette équation de propagation, qu'on va écrire dans l'espace des fréquences, donc en passant de l'espace direct où le champ dépend de t, à l'espace de Fourier où le champ va dépendre de oméga, par contre je ne touche pas à la variable r, et donc je vais directement avoir ici dans l'espace de Fourier, le laplacien de E de r et de oméga, ensuite, la dérivée seconde du champ par rapport au temps ou la dérivée seconde de la polarisation par rapport au temps, donc ça c'est quelque chose qu'on a vu déjà un certain nombre de fois, vous savez donc qu'on va faire apparaître moins oméga deux sur c deux avec la susceptibilité, et on aura finalement, en utilisant le vecteur d'onde du milieu on pourra écrire cela sous la forme plus k de oméga au carré multiplié par E de r et de oméga égal toujours à l'espace de Fourier ici la dérivée seconde de p par rapport au temps va me donner un moins oméga deux, donc j'aurai moins oméga deux sur epsilon zéro c deux fois la polarisation non-linéaire p n l de r et de oméga. Donc on retrouve évidemment ce qu'on a vu dans la vidéo précédente, à gauche l'équation de Helmholtz, et puis ici, le terme source, la seule différence c'est qu'au lieu d'avoir à faire à des fréquences discrètes, et bien on a un ensemble continu de fréquences. Bien, donc on arrive à cette équation de propagation, avec dans le membre de gauche l'équation de Helmholtz, et dans le membre de droite le terme source lié à la polarisation non-linéaire, et on va procéder comme dans la vidéo précédente, en faisant l'approximation paraxiale, c'est-à-dire qu'on va supposer qu'on a un faisceau lumineux qui se propage selon l'axe z, avec une faible diverge, et on va introduire l'enveloppe u de r du faisceau multiplié par une porteuse en exponentiel i k de oméga z. Donc c'est aussi comme ce qu'on avait fait la semaine dernière dans le cas de l'approximation paraxiale pour la propagation d'un faisceau lumineux monochromatique, simplement ici, c'est pas un faisceau monochromatique, on a une dépendance en oméga, ce qui fait que la porteuse ici va dépendre de la fréquence oméga, mais c'est la seule différence. Donc comme précédemment, on pourra écrire le lapalcien de E en dérivant cette expression-là, on va évidemment avoir des termes supplémentaires pour la dérivée par rapport à z, donc des termes que vous connaissez bien, et l'approximation paraxiale, comme vous le savez, consiste à négliger la dérivée seconde de u par rapport à z, et donc on va négliger ce terme non seulement par rapport à la dérivée première, mais également par rapport au laplacien transverse de l'enveloppe u. Donc, si je remplace dans l'équation de propagation, j'obtiens une équation qui ressemble beaucoup à ce qu'on avait vu dans la vidéo précédente pour des ondes monochromatiques, simplement, au lieu d'avoir un ensemble discret d'équations différentielles, et bien j'ai un continuum d'équations différentielles pour toutes les valeurs possibles de oméga. L'intérêt évidemment, c'est qu'on a ici une équation différentielle du premier ordre en z, et ce qu'on va faire, c'est qu'on va faire passer, enfin, on va diviser l'équation par deux i k de oméga pour avoir directement une équation qui porte sur d u sur d z, donc on aura l'équation de propagation avec donc la dérivée première de u par rapport à z, un terme qui fait intervenir le laplacien transverse mais également la dépendance en fréquence ici, et puis à droite le terme source de l'équation de propagation. On aboutit donc à cette équation de propagation pour l'enveloppe u de r et de oméga associée au champ électrique. C'est une équation qui est satisfaisante puisqu'on a une dérivée première par rapport à z, c'était l'intérêt de l'approximation paraxiale. L'inconvénient c'est que dans la définition de u, on a ici une porteuse qui va dépendre de oméga, ce qui va rendre l'exploitation de cette équation un peu compliquée. Donc pour palier à ce problème, on va revenir à la variable originale, à la variable E de r et de oméga, et pour ça il suffit simplement de remplacer ici u par le produit E exponentiel moins i k z dans l'équation. Donc je vais devoir calculer la dérivée première de u par rapport à z, donc ça c'est très simple puisque u c'est un produit entre E de x y z oméga et exponentiel moins i k z, donc quand je dérive ce produit j'ai un premier terme qui est la dérivée de E par rapport à z, et puis un deuxième terme qui va être la dérivée de exponentiel moins i k de oméga z qui va s'écrire moins i k de oméga multiplié par E avec évidemment en facteur exponentiel moins i k de oméga z. Donc je n'ai plus qu'à remplacer dans cette équation d u sur d z par sa valeur pour obtenir l'équation différentielle suivante : donc d u sur d z, j'aurai les deux termes d sur d z moins i k de oméga, le laplacien je le retrouve inchangé puisque je n'ai pas changé la dépendance en x et en y entre ces deux expressions, j'aurai évidemment exponentiel moins i k de oméga z en facteur mais que je peux simplifier avec le terme qui est ici, et donc finalement on retrouve ici le terme source sans le facteur exponentiel moins i k de oméga z. Donc voilà l'équation différentielle qu'on obtient finalement sur le champ électrique E de r et de oméga dans l'espace des fréquences. Alors cette équation différentielle, c'est l'équation qu'on appelle l'équation de propagation non-linéaire dans le cadre de l'approximation de l'onde lentement variable, c'est ce qu'on appelle en anglais Slowly-Evolving Wave Approximation, c'est une équation qui a été établie assez récemment en 1997, tout simplement parce que contrairement à ce qu'on appelle l'approximation de l'enveloppe lentement variable, elle ne fait pas d'hypothèse sur la durée de l'impulsion, elle peut s'appliquer pour des impulsions très très brèves qui sont précisément les impulsions qu'on a été capable de produire à la fin du vingtième siècle. Donc cette équation elle va être évidemment très utile et, en particulier, comme c'est une équation différentielle du premier ordre en z, elle se prêtera bien à une résolution numérique où on pourra prendre en compte la valeur exacte de k de oméga selon le matériau considéré, et on pourra aussi prendre en compte les effets spatiaux avec le laplacien transverse.
