The discovery of vast quantities of hydrocarbons during exploration of the Mackenzie Delta has led to the detection of numerous anomalous deep holes in various delta channels in the past four decades. The proposed pipeline crossing from the Mackenzie Delta to northern Alberta has renewed interest in studying the stability of these deep holes. The main goal of this research project was to study ice-covered flow velocity and bed shear stress distributions using a previously calibrated computational fluid dynamics (CID) model for a 30 m deep hole in the East Channel of the Mackenzie Delta to provide an initial assessment of its stability. In this study, a previously developed 3D CFD model using the FLUENT code with the renormalization group theory (RNG) k-s turbulence closure model with nonequilibrium wall functions was adapted to represent ice-covered flow conditions for the study reach. The numerical model was applied to simulate four flow conditions (500, 720, and 1000 [m.sup.3]/s along with a moving ice cover during a flow rate of 720 [m.sup.3]/s). Results indicated that a single vortex was formed near the inner bend above the hole and deposition may be a key process occurring during ice-covered flow conditions. Key words: computational fluid dynamics, finite element, ice-covered flow, bed shear stress, deep hole. L'exploration du delta du Mackenzie a la suite de la decouverte de vastes quantites d'hydrocarbures a mene a la detection de plusieurs fosses d'affouillement profondes anormales dans les divers canaux du delta an cours des quarante dernieres annees. Le franchissement propose du pipeline du delta du Mackenzie an nord de ''Alberta a renouvele l'interet de verifier la stabilite de ces fosses profondes. Le principal objectif de ce projet de recherche etait d'etudier la vitesse d'ecoulement sous convert de glace et les distributions des contraintes de cisaillement du lit en utilisant un modele numerique de la dynamique des fluides etalonne pour une fosse profonde de 30 m dans le canal Est du delta du Mackenzie afin de fournir une evaluation initiale de sa stabilite. Dans cette etude, un modele numerique tridimensionnel de la dynamique des fluides, developpe anterieurement, utilisant le code FLUENT et le modele de turbulence k-s ' RNG ' muni de fonctions de non equilibre des parois, a ete adapte de maniere a representer les conditions d'ecoulement sous convert de glace pour la portee de l'etude. Le modele numerique a ensuite ete applique pour simuler quatre conditions d'ecoulement (500, 720, 1000 [m.sup.3]/s de meme qu'un convert de glace se deplargant avec un ecoulement de 720 [m.sup.3]/s). Les resultats indiquent qu'un seul vortex s'est forme pres de la portion interne du coude du fleuve au-dessus de la fosse d'affouillement et que la deposition pent etre un processus important survenant durant les conditions d'ecoulement sous convert de glace. Mots-cles : dynamique des fluides numerique, elements fmis, ecoulement sous convert de glace, contraintes de cisaillement du lit, fosse d'affouillement profonde. [Traduit par la Redaction], Introduction The discovery of vast quantities of hydrocarbon reserves in the Mackenzie Delta-Beaufort Sea region during the late 1960s prompted the oil and gas industry to become actively interested in [...]