In this example we compare the numerical predictions of Fluidity with a well-established two-dimensional cartesian geometry benchmark result for Stokes flow (Blankenbach et al., 1989). Fluidity’s results from this, and similar cases, are published in Davies et al. (2011), where solution strategies employed for solving the Stokes system are also discussed.
The specific example considered is steady-state isoviscous convection at a Rayleigh number (Ra) of 105, in a two-dimensional square domain of unit dimensions. Boundary conditions for temperature are T = 0 at the surface, T = 1 at the base, with insulating (homogeneous Neumann) side-walls. For velocity, free-slip and no normal flow boundary conditions are specified at all boundaries.
The results presented in the figures above. The final steady-state solution consists of one convective cell, with hot upwelling flow confined to the left hand side of the domain and cold down-welling flow on the right hand side. Note that, excluding the location of upwelling flow at x = 0 or x = 1, results are insensitive to the initial condition. Quantitatively, results show excellent agreement with the benchmark predictions of Blankenbach et al. (1989), with solution accuracy improving with increased resolution, as expected. A higher resolution case (at a resolution of 96 × 96 elements), achieves results that are almost identical to the benchmark solution.
Further details about Fluidity's Stokes flow model can be found in the following paper: