La modélisation des champs sonores en acoustique architecturale est d'autant plus complexe que la géométrie d'étude est elle-même complexe. En particulier, le domaine de propagation peut présenter de nombreux volumes couplés de tailles très différentes, connectés par de faibles surfaces de couplage. Dans ces conditions, les formulations analytiques classiques, basées par exemple sur la théorie de la réverbération, ou les méthodes numériques, basées par exemple sur des méthodes de tracé de rayons ou de faisceaux sonores, sont difficilement applicables, voire impossibles. Ces dernières années, une solution alternative, fondée sur l'utilisation d'une équation de diffusion pour la densité d'énergie sonore, a été développée et validée. Dans ce présent travail, nous présentons les développements permettant ainsi de modéliser les champs sonores en acoustique des salles, en prenant en compte les différents phénomènes propagatifs mis en jeu, tel que l'atténuation atmosphérique, l'absorption par les parois, la transmission acoustique à travers les cloisons et la réflexion mixte au niveau des surfaces. Dans un second temps, nous présentons la méthode retenue pour la résolution numérique des équations de diffusion mises en jeu. A cet effet, une interface logicielle a été spécialement développée, permettant de manipuler de manière fonctionnelle les données géométriques et acoustiques du problème, avant l'exécution d'un logiciel de calcul multi-physique, en s'affranchissant ainsi du formalisme spécifique associé au modèle de diffusion. Cette interface permet ensuite le post-traitement de l'ensemble des résultats, notamment afin de déterminer et de représenter différents paramètres classiques d'acoustique des salles (niveaux sonores, temps de réverbération…).