La simulation numérique de la propagation d'ondes de choc acoustiques nécessite de prendre en compte des échelles temporelles et spatiales très différentes (plusieurs ordres de grandeur). Les temps de calcul élevés rendent les approches directes dans des configurations réalistes (équations d'Euler ou de Navier-Stokes) prohibitives pour le moment. Ceci justifie le recours à des équations modèles. Le modèle le plus courant est basé sur l'approximation parabolique standard, dont le domaine de validité reste néanmoins limité aux faisceaux collimatés. Pour aller au-delà et s'affranchir de ces limitations, une nouvelle approche nommée HOWARD ("Heterogenous One Way Approximation for the Resolution of Diffraction") a été mise au point. Elle permet de simuler la propagation d'ondes de choc acoustiques en milieux hétérogènes sans limite de validité angulaire, contrairement aux approximations paraboliques utilisées habituellement. Seul le champ rétro-diffusé est négligé (approximation de milieu faiblement hétérogène). Dans un premier temps, cette méthode sera exposée en détail et comparée aux méthodes existantes. Puis, différentes configurations seront étudiées pour en illustrer les potentialités.