Nonsmooth Exact Penalization Second-order Methods for incompressible Bingham flows

TL;DR

This paper introduces a novel second-order method for solving nonsmooth optimization problems arising from the exact penalization of incompressibility in Bingham fluid flows, avoiding the need for divergence-free approximation schemes.

Contribution

It proposes a new algorithm leveraging generalized second-order information for nonsmooth penalized problems, specifically for incompressible Bingham flows, improving enforcement of divergence-free conditions.

Findings

The method effectively enforces divergence-free velocity fields.

It outperforms traditional inexact penalization approaches.

The approach is applicable to nonsmooth optimization problems in fluid dynamics.

Abstract

We consider the exact penalization of the incompressibility condition $div(u)=0$ for the velocity field of a Bingham fluid in terms of the $L^1$-norm. This penalization procedure results in a nonsmooth optimization problem for which we propose an algorithm using generalized second-order information. Our method solves the resulting nonsmooth problem by considering the steepest descent direction and extra generalized second-order information associated to the nonsmooth term. This method has the advantage that the divergence-free property is enforced by the descent direction proposed by the method without the need of build-in divergence-free approximation schemes. The inexact penalization approach, given by the $L^2$-norm, is also considered in our discussion and comparison.