A new finite-difference time-domain scheme is presented for the propagation of VLF–LF radio waves in the Earth–ionosphere waveguide. The new scheme relies on the implicit solution of the auxiliary equation that governs the current density in the ionosphere. The advantages and drawbacks of the new scheme are discussed. Its main advantage is its stability condition, which is the same as that of the FDTD method in a vacuum. This permits the time step of the calculation to be increased and then the overall computational time to be reduced. Numerical experiments demonstrate the accuracy of the new scheme and the reduction of the computational time.

Cet article décrit un nouveau schéma aux différences finies pour la propagation VLF–LF dans le guide d'onde Terre–ionosphère. Ce schéma repose sur une solution implicite de l'équation auxiliaire qui gouverne la densité de courant dans l'ionosphère. Les avantages et inconvénients du nouveau schéma sont discutés. Son principal avantage tient dans sa condition de stabilité qui est identique à celle de la méthode FDTD dans le vide. Cela permet d'augmenter le pas temporel de la résolution FDTD et, en conséquence, de réduire le temps de calcul. Des résultats de calcul illustrent la bonne précision du nouveau schéma et les réductions de temps calcul qu'il permet d'obtenir.