Ce travail s’intéresse à un écoulement autour d’un obstacle cylindro-conique chauffé, enrotation à l’intérieur d’une cavité cylindrique. Il s’agit d’étudier l’influence d’uneperturbation de type harmonique sur le comportement dynamique de l’écoulement et sur lemécanisme de l’échange thermique avec l’obstacle. La fréquence de référence del’écoulement périodique a été déterminée par l’analyse de l’évolution du régime vers sonétat stationnaire. L’étude a été menée pour des fréquences variant de 1/3 de la fréquencede référence à 10 fois cette dernière. La périodicité de l’écoulement est assurée par lacondition d’entrée sur le débit supposé périodique. L’étude numérique a été réalisée pourdes valeurs du nombre de Reynolds inférieur à 100 et pour des valeurs de nombre de Taylorinférieur à 500 afin d’éviter la présence des instabilités au sein de l’écoulement. Deuxméthodes numériques ont été utilisées : une méthode de différences-finies combinée avec unschéma ADI et la méthode des éléments-finis. Le fluide considéré est l’air avec despropriétés physiques constantes. Le couplage convection-conduction a été considéré dansune plage de variation du rapport K de la conductivité thermique dusolide à celle du fluide de 10-2 à 103. Les résultats montrent queles composantes de la vitesse varient en phase avec le débit, alors que la pression et lesdifférents paramètres thermiques sont en déphasage avec la vitesse d’entrée. On note aussil’aspect périodique de tous les paramètres dynamiques et thermiques de l’écoulement pourles nombres de Reynolds utilisés. L’influence de la fréquence se manifeste davantage pourses valeurs élevées. Le nombre de Nusselt est proportionnel à la fréquence. Ces résultatssont en bon accord avec ceux données par d’autres études de référence.

By considering the basic stress equations for a unit volume of ice, a set of differential equations describing ice shelf flow is derived. In view of the lack of basal shear stresses at the bottom of ice shelf a model simulation which is restricted to the horizontal dimensions will not imply substantial errors. The model is applied to the Filchner-Ronne Ice Shelf, Antarctica, and model equations are solved in terms of finite differences on a 10 × 10 km grid. Present ice thickness data and boundary conditions, i.e. the balance velocities at the grounding line and strain rates at the ice front are entered as input. Using a non-linear Glen-type flow law (n=3) and a constant depth-averaged flow law parameter, representing an ice temperature of −17°C, a convincing velocity field is derived as a solution of the model equations. The model takes into account restrained flow across ice rumples where sufficient field data are available. A diagnostic run reproducing present velocity magnitudes is followed by two prognostic runs, each representing 2000 years of simulation. Transient ice thickness changes are obtained from solving the mass conservation equation. Two different assumptions concerning basal melting rates demonstrate its importance to ice shelf dynamics. Assumptions are: a) no basal melting, b) basal melting rates (−2m a−1 to +3m a−1) as derived from model results and geophysical field data.