Dans cet exposé, nous ferons un panorama de méthodes et simulations numériques pour les fluides à seuil, basées sur des méthodes de dualité. Dans un premier temps, nous présenterons le problème des équations de type Bingham dans un canal en expansion-contraction qui permet d'obtenir des couches limites viscoplastiques. Nous revisiterons la théorie asymptotique d'Oldroyd (1947) dans le cas où les nombres caractéristiques sont modérés. Cette étude mélange simulations HPC et allers-retours avec des expériences physiques d'IRSTEA. Une seconde partie traitera ensuite d'un modèle original de Saint-Venant-Bingham pour ces fluides viscoplastiques, en lien avec des applications géophysiques. Nous proposons un nouveau schéma volumes-finis qui couple dualité et techniques équilibrées. Ses propriétés sont illustrées sur un prototype d'avalanche de neige dense dans le couloir de Taconnaz (massif du Mont-Blanc).