Optimal isoparametric finite elements and error estimates for domains involving curved boundaries
SIAM Journal on Numerical Analysis
Mixed and hybrid finite element methods
Mixed and hybrid finite element methods
SIAM Journal on Scientific Computing
Structural boundary design via level set and immersed interface methods
Journal of Computational Physics
Shapes and geometries: analysis, differential calculus, and optimization
Shapes and geometries: analysis, differential calculus, and optimization
Journal of Computational Physics
Algebraic Mesh Quality Metrics
SIAM Journal on Scientific Computing
Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator
FCRC '96/WACG '96 Selected papers from the Workshop on Applied Computational Geormetry, Towards Geometric Engineering
Structural optimization using sensitivity analysis and a level-set method
Journal of Computational Physics
Locally Conservative Fluxes for the Continuous Galerkin Method
SIAM Journal on Numerical Analysis
Journal of Computational Physics
Peristaltic transport of non-Newtonian fluid in a diverging tube with different wave forms
Mathematical and Computer Modelling: An International Journal
Peristaltic transport in a tapered tube
Mathematical and Computer Modelling: An International Journal
Applying a second-kind boundary integral equation for surface tractions in Stokes flow
Journal of Computational Physics
Journal of Computational and Applied Mathematics
Hi-index | 31.45 |
Transport is a fundamental aspect of biology and peristaltic pumping is a fundamental mechanism to accomplish this; it is also important to many industrial processes. We present a variational method for optimizing the wave shape of a peristaltic pump. Specifically, we optimize the wave profile of a two dimensional channel containing a Navier-Stokes fluid with no assumption on the wave profile other than it is a traveling wave (e.g. we do not assume it is the graph of a function). Hence, this is an infinite-dimensional optimization problem. The optimization criteria consists of minimizing the input fluid power (due to the peristaltic wave) subject to constraints on the average flux of fluid and area of the channel. Sensitivities of the cost and constraints are computed variationally via shape differential calculus and we use a sequential quadratic programming (SQP) method to find a solution of the first order KKT conditions. We also use a merit-function based line search in order to balance between decreasing the cost and keeping the constraints satisfied when updating the channel shape. Our numerical implementation uses a finite element method for computing a solution of the Navier-Stokes equations, adjoint equations, as well as for the SQP method when computing perturbations of the channel shape. The walls of the channel are deformed by an explicit front-tracking approach. In computing functional sensitivities with respect to shape, we use L^2-type projections for computing boundary stresses and for geometric quantities such as the tangent field on the channel walls and the curvature; we show error estimates for the boundary stress and tangent field approximations. As a result, we find optimized shapes that are not obvious and have not been previously reported in the peristaltic pumping literature. Specifically, we see highly asymmetric wave shapes that are far from being sine waves. Many examples are shown for a range of fluxes and Reynolds numbers up to Re=500 which illustrate the capabilities of our method.