Adjoint-based Control of Three Dimensional Stokes Droplets

TL;DR

This paper introduces an adjoint-based optimization method for controlling the deformation of neutrally buoyant droplets in Stokes flow, leveraging shape calculus and boundary integral methods for efficient gradient computation.

Contribution

It presents a novel continuous adjoint formulation combined with high-order boundary integral discretization for droplet shape control in Stokes flow.

Findings

Efficient adjoint evolution equation expressed via scalar field

High accuracy and stability demonstrated on control problems

Effective gradient computation for shape optimization

Abstract

We develop a continuous adjoint formulation and implementation for controlling the deformation of clean, neutrally buoyant droplets in Stokes flow through farfield velocity boundary conditions. The focus is on dynamics where surface tension plays an important role through the Young-Laplace law. To perform the optimization, we require access to first-order gradient information, which we obtain from the linearized sensitivity equations and their corresponding adjoint by applying shape calculus to the space-time tube formed by the interface evolution. We show that the adjoint evolution equation can be efficiently expressed through a scalar adjoint transverse field. The optimal control problem is discretized by high-order boundary integral methods using Quadrature by Expansion coupled with a spherical harmonic representation of the droplet surface geometry. We show the accuracy and stability of the scheme on several tracking-type control problems.