Shape optimization of peristaltic pumping

Shawn W. Walker, Michael Shelley

Research output: Contribution to journalArticle

Abstract

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 L2-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.

Original languageEnglish (US)
Pages (from-to)1260-1291
Number of pages32
JournalJournal of Computational Physics
Volume229
Issue number4
DOIs
StatePublished - Feb 20 2010

Fingerprint

shape optimization
Shape optimization
pumping
quadratic programming
Quadratic programming
tangents
Fluids
fluid power
Differentiation (calculus)
Fluxes
differential calculus
costs
optimization
fluids
sensitivity
sine waves
profiles
biology
traveling waves
Navier-Stokes equation

Keywords

  • Finite element method
  • PDE-constrained optimization
  • Peristalsis
  • Pumping
  • Shape optimization
  • Variational method

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

Shape optimization of peristaltic pumping. / Walker, Shawn W.; Shelley, Michael.

In: Journal of Computational Physics, Vol. 229, No. 4, 20.02.2010, p. 1260-1291.

Research output: Contribution to journalArticle

Walker, Shawn W. ; Shelley, Michael. / Shape optimization of peristaltic pumping. In: Journal of Computational Physics. 2010 ; Vol. 229, No. 4. pp. 1260-1291.
@article{1a7f10a358b146bcbee731b9541a6303,
title = "Shape optimization of peristaltic pumping",
abstract = "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 L2-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.",
keywords = "Finite element method, PDE-constrained optimization, Peristalsis, Pumping, Shape optimization, Variational method",
author = "Walker, {Shawn W.} and Michael Shelley",
year = "2010",
month = "2",
day = "20",
doi = "10.1016/j.jcp.2009.10.030",
language = "English (US)",
volume = "229",
pages = "1260--1291",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "4",

}

TY - JOUR

T1 - Shape optimization of peristaltic pumping

AU - Walker, Shawn W.

AU - Shelley, Michael

PY - 2010/2/20

Y1 - 2010/2/20

N2 - 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 L2-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.

AB - 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 L2-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.

KW - Finite element method

KW - PDE-constrained optimization

KW - Peristalsis

KW - Pumping

KW - Shape optimization

KW - Variational method

UR - http://www.scopus.com/inward/record.url?scp=72049098545&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=72049098545&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2009.10.030

DO - 10.1016/j.jcp.2009.10.030

M3 - Article

AN - SCOPUS:72049098545

VL - 229

SP - 1260

EP - 1291

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 4

ER -