Topology optimization of Stokes flow with traction boundary conditions using low-order finite elements

We consider topology optimization of Stokes flow with traction boundary conditions using finite elements with low-order velocity-approximation and an element-wise constant hydrostatic pressure. The finite element formulation is stabilized using a penalty on the jump in pressure between adjacent elements. Convergence of solutions to the finite element-discretized topology optimization problem is shown, and several optimization problems are solved using a preconditioned conjugate gradient solver for the finite element matrix problem. Stable convergence to high-quality designs without an excessive number of linear solver iterations is observed, and it is seen that the finite element formulation is not particularly sensitive to the choice of the pressure jump penalty parameter, thus making it a practically useful method. (C) 2021 The Author(s). Published by Elsevier B.V.