The method of fundamental solutions for Brinkman flows. Part II. Interior domains

In part I, we considered the application of the method of fundamental solutions (MFS) for solving numerically the Brinkman fluid flow in the unbounded porous medium outside obstacles of known or unknown shapes. In this companion paper we consider the corresponding interior problem for the Brinkman flow in a bounded porous medium which contains an unknown rigid inclusion D ⊂ Ω . The inclusion D is to be identified by a pair of Cauchy data represented by the fluid velocity and traction on the boundary ∂ Ω . The fluid velocity and pressure of the incompressible viscous flow in the porous medium Ω \ D ¯ are approximated by linear combinations of fundamentals solutions for the Brinkman system with sources ( singularities ) placed outside the closure of the solution domain, i.e. in D ∪ ( R 2 \ Ω ¯ ) , assuming, for simplicity, that we analyse planar domains. By further assuming that the unknown obstacle D is star-shaped (with respect to the origin), the inverse problem recasts as the minimization of the nonlinear Tikhonov’s regularization functional with respect to the MFS expansion coefficients and the discretized polar radii defining D . This minimization subject to simple bounds on the variables is solved numerically using the MATLAB optimization toolbox routine lsqnonlin.