A fast adjoint-free shape gradient for p-Laplacian compliance: algorithmic speedup and engineering benchmarks
Shape optimization of nonlinear membranes and non-Newtonian flow domains governed by the $p$-Laplace operator is hampered by the cost of an auxiliary linearised adjoint solve at every iteration. We show that this solve is unnecessary: a self-adjoint identity reduces the adjoint state to a multiple of the primal state...