---
title: Coupled Complex Boundary Method (CCBM)
url: https://www.emergentmind.com/topics/coupled-complex-boundary-method-ccbm
type: topic
---

# Coupled Complex Boundary Method (CCBM)

The Coupled Complex Boundary Method (CCBM) is a reformulation strategy for inverse Cauchy problems, free boundary problems, and geometric inverse problems in which overdetermined Dirichlet and Neumann data are encoded into a single complex Robin-type boundary condition for a complex-valued state. First put forward by Cheng et al. (2014) to deal with inverse source problems and subsequently developed in shape optimization and coefficient-identification settings, CCBM replaces boundary overdeterminedness by the requirement that the imaginary part of the complex solution vanish in the domain, and then solves the resulting problem by optimization, adjoint-based sensitivity analysis, and numerical regularization [2205.12620][2302.11828][2403.11819][2508.16898][2204.03416][2602.01656][2509.09533][2605.12202].

## 1. Conceptual principle and historical development

The central idea of CCBM is to transform an overdetermined boundary value problem into a single complex boundary value problem with a complex Robin boundary condition coupling Dirichlet and Neumann data. In the exterior Bernoulli setting, this was introduced as a shape optimization framework in which the free boundary is identified by optimizing a cost function constructed from the imaginary part of the solution in the whole domain [2205.12620]. In inverse Cauchy problems for elliptic equations, the same principle was used to combine measured data on the accessible boundary and missing data on the inaccessible boundary into a single complex formulation, leading to an operator equation suitable for iterative regularization [2204.03416].

In this formulation, one writes the state as \(u=u_1+i u_2\) or \(w=w_r+i\,\Im\{w\}\), depending on the application. The real part carries one component of the physical boundary information and the imaginary part carries the other. If the geometry, source, or coefficient is correct, the imaginary component vanishes in the domain, so reconstruction reduces to minimizing a volumetric misfit built from that imaginary part. In the 2022 Bernoulli paper, this was described as shifting the least-squares criterion from the boundary to the interior [2205.12620].

Subsequent work broadened the method across PDE classes. The same complex-coupling mechanism was used for Stokes free boundary problems and Stokes obstacle identification, for advection–diffusion obstacle recovery from a single pair of Cauchy data, for diffusion coefficient reconstruction, for bioluminescence tomography, and for cavity reconstruction under a homogeneous Robin condition [2302.11828][2403.11819][2508.16898][2602.01656][2509.09533][2605.12202]. This suggests a unifying CCBM template: encode the available Cauchy data into one complex boundary condition, solve a well-posed complex PDE, and identify the unknown by enforcing vanishing imaginary part.

## 2. Canonical mathematical formulations

Across the cited literature, CCBM appears in several PDE settings, but the complex coupling pattern remains stable.

| Setting | Complex boundary formulation | Primary unknown |
|---|---|---|
| Exterior Bernoulli problem [2205.12620] | \(\partial_n u + i u = \lambda\) on the free boundary \(\Sigma\) | Free boundary |
| Stokes free boundary problem [2302.11828] | \(-pn + \alpha u + i (u\cdot n)n = 0\) on \(\Sigma\) | Free boundary |
| Stokes obstacle identification [2403.11819] | \(\sigma(u,p)n + i u = g + i f\) on the accessible boundary \(\Sigma\) | Obstacle shape |
| Advection–diffusion obstacle recovery [2508.16898] | \(\sigma\,\partial_n u + i u = g + i f\) on \(\Sigma\), \(u=0\) on \(\Gamma\) | Inclusion shape |
| Diffusion coefficient recovery [2602.01656] | \(\alpha \partial_n w + i w = g + i f\) on \(\Gamma\) | Diffusion coefficient |
| Bioluminescence tomography [2509.09533] | \(D\partial_\nu u + i\alpha u = g_1 + i\alpha g_2\) on \(\Gamma\) | Source support and intensity |
| Cavity reconstruction with homogeneous Robin condition [2605.12202] | \(\partial_n u + i\rho u = g + i\rho f\) on \(\Sigma\), \(\partial_n u + \alpha u = 0\) on \(\Gamma\) | Unknown interior boundary |

For the advection–diffusion problem, the forward model is
\[
-\operatorname{div}(\sigma \nabla u) + \mathbf{b}\cdot \nabla u = 0
\quad \text{in } \Omega=D\setminus \overline{\omega},
\]
with \(u=0\) on the unknown boundary \(\Gamma\) and \(\sigma\partial_n u + iu = g+if\) on the accessible boundary \(\Sigma\). Decomposition \(u=u_1+i u_2\) yields two coupled real PDEs, and the inverse problem is equivalent to enforcing \(u_2=0\) in \(\Omega\) [2508.16898].

For the Stokes obstacle problem, CCBM uses a complex-valued velocity and pressure, with the complex Robin-like condition \(\sigma(u,p)n + i u = g + i f\) on the accessible boundary and \(u=0\) on the obstacle boundary. The obstacle is recovered by minimizing a cost functional based on the imaginary parts of the velocity and pressure over the whole domain [2403.11819]. For free boundary Stokes flow, the overdetermined stress-free and no-penetration conditions are merged into \(-pn + \alpha u + i (u\cdot n)n = 0\), again with recovery achieved by forcing the imaginary components to vanish [2302.11828].

In BLT, the complex coupling is parameter-dependent:
\[
-\nabla\cdot(D\nabla u)+\mu_a u=\phi(x)\chi_{\Omega_0},
\qquad
D\partial_\nu u + i\alpha u = g_1 + i\alpha g_2,
\]
and the support \(\Omega_0\) of the source is sought together with its intensity \(\phi\) [2509.09533]. In diffusion coefficient recovery, the unknown is no longer a shape but the spatially varying coefficient \(\alpha(x)\), yet the same principle is retained: solve a complexified PDE and reconstruct \(\alpha\) by driving \(\Im\{w(\alpha)\}\) toward zero [2602.01656].

## 3. Cost functionals, adjoints, and shape calculus

The defining CCBM objective is a domain integral built from the imaginary part of the complex state. In the exterior Bernoulli problem and in advection–diffusion obstacle recovery, the standard choice is
\[
J(\Omega)=\frac12\|u_2\|_{L^2(\Omega)}^2
=\frac12\int_\Omega |u_2(x)|^2\,dx,
\]
and minimizing \(J\) over admissible shapes yields the sought domain [2205.12620][2508.16898]. In the Stokes free boundary problem, the analogue is
\[
J(\Omega)=\frac12\|i\|_{L^2(\Omega)^d}^2+\frac12\|i_i\|_{L^2(\Omega)}^2,
\]
where \(i\) and \(i_i\) are the imaginary parts of velocity and pressure [2302.11828]. In obstacle identification for Stokes flow, the cost is
\[
J(\Omega)=\frac12\left(\|u_i\|_{L^2(\Omega)^d}^2+\|p_i\|_{L^2(\Omega)}^2\right),
\]
and \(J=0\) precisely when the current shape matches the true obstacle [2403.11819].

These functionals are paired with adjoint equations and boundary representations of the shape derivative. For advection–diffusion, the first-order shape derivative in direction \(\theta\) is
\[
dJ(\Omega)[\theta]=\int_\Gamma G(u,p)\,n\cdot\theta\,ds,
\qquad
G(u,p)=\sigma\left(\partial_n u_1\partial_n p_2-\partial_n u_2\partial_n p_1\right),
\]
where the adjoint \(p=p_1+i p_2\) solves
\[
\operatorname{div}(\sigma\nabla p)+\mathbf{b}\cdot\nabla p+(\operatorname{div}\mathbf{b})p=u_2
\quad\text{in }\Omega,
\]
with \(p=0\) on \(\Gamma\) and \(\sigma\partial_n p + p\,\mathbf{b}\cdot n - ip =0\) on \(\Sigma\) [2508.16898].

For the Stokes obstacle problem, the rearrangement method gives the boundary form
\[
dJ(\Omega)[V]
=
\int_\Gamma
\left(
\Im\left\{\overline{\sigma(v,q)n}\cdot u\right\}
+\frac12|u_i|^2
\right)(n\cdot V)\,ds,
\]
with \(v,q\) the adjoint state [2403.11819]. For cavity reconstruction with homogeneous Robin condition, the derivative is written as
\[
dJ(\Omega)[\theta]
=
\int_\Gamma
\left[
\frac12|u_2|^2+\Psi_\Gamma(u,p)
\right](\theta\cdot n)\,ds,
\]
where
\[
\Psi_\Gamma(u,p)
=
(\nabla u_1\cdot \nabla p_2-\nabla u_2\cdot \nabla p_1)
+
(\alpha^2-\alpha\kappa)(p_1u_2-p_2u_1),
\]
and \(\kappa\) is the mean curvature [2605.12202].

The Bernoulli analysis goes further: it proves the existence of the shape derivative of the complex state, computes the shape gradient of the cost functional, characterizes its shape Hessian at the optimal domain under a strong and then a mild regularity assumption, and proves compactness of the latter expression [2205.12620]. The associated Riesz operator is compact in \(H^{1/2}(\Sigma^*)\), not coercive in \(H^1(\Sigma^*)\), which the paper interprets as instability typical of free boundary inverse problems [2205.12620].

## 4. Numerical algorithms and regularization mechanisms

A recurrent numerical component in CCBM is Sobolev-gradient descent. Direct use of the boundary shape gradient can damage mesh quality, so several papers replace the raw boundary density by an \(H^1\)-type Riesz representative. In advection–diffusion, the regularized deformation field \(\xi\) is obtained from
\[
c_b\int_\Omega \nabla \xi:\nabla \varphi
+
(1-c_b)\int_\Gamma \nabla_\Gamma \xi:\nabla_\Gamma \varphi
=
-\int_\Gamma G\,n\cdot\varphi,
\]
and the domain is updated by \(x\mapsto x+t\xi(x)\) [2508.16898]. Closely related \(H^1\)-type systems appear in the Bernoulli, Stokes, and Robin-cavity papers, where finite elements are used for the state, adjoint, and gradient equations, and the stepsize is chosen by Armijo or adaptive line search [2205.12620][2302.11828][2403.11819][2605.12202].

Two later developments introduce explicit constraint handling through ADMM. In advection–diffusion shape recovery, an auxiliary variable \(v\) enforces bounds \(a\leq u_1\leq b\) on the real part of the state, and the augmented Lagrangian
\[
\mathcal{L}_\beta(\omega,v;\lambda)
=
J(\Omega)+I_\mathcal{K}(v)+\frac{\beta}{2}\|u_1-v\|^2+\lambda(u_1-v)
\]
leads to alternating updates of the shape, the projected auxiliary variable, and the multiplier [2508.16898]. In cavity reconstruction with a homogeneous Robin condition, the same strategy is applied to \(\Re\{u\}\) through the constraint \(a\leq \Re\{u(x)\}\leq b\), with a projection step
\[
v^{k+1}
=
\operatorname{proj}_{\mathcal K}
\left(\Re\{u^{k+1}\}+\frac{\lambda^k}{\gamma}\right)
\]
followed by the multiplier update \(\lambda^{k+1}=\lambda^k+\gamma(\Re\{u^{k+1}\}-v^{k+1})\) [2605.12202]. The advection–diffusion paper also compares exact and partial shape gradients and reports that partial gradient selection can improve efficiency [2508.16898].

CCBM has also been used outside shape descent. In inverse Cauchy problems, the complex reformulation yields an operator equation \(K\phi=f^\delta\), with \(K\phi=\tilde u_2\), and this is solved by Golub–Kahan bidiagonalization with Givens rotation under a discrepancy-principle stopping rule. The resulting CCBM-GKB method is analyzed as an iterative regularization scheme [2204.03416].

For coefficient recovery, a modified CCBM introduces a gradient-weighted \(H^1\)-type term in the misfit,
\[
J(\alpha)
=
\frac12\left(
w_0\|\Im\{w(\alpha)\}\|_{L^2(\Omega)}^2
+
w_1\|\nabla \Im\{w(\alpha)\}\|_{L^2(\Omega)}^2
\right),
\]
together with Tikhonov regularization and Sobolev-gradient descent. A projection-based extension and a pick-a-point strategy are then used for piecewise-constant coefficients [2602.01656]. In BLT, the method is embedded in a perimeter- and volume-regularized level-set framework, and first-order optimality makes the source intensity explicit through \(\phi=-(1/\varepsilon)w_1|_{\Omega_0}\), reducing the algorithm to a shape-only optimization problem [2509.09533].

## 5. Applications and reported empirical behavior

The Bernoulli study reports successful reconstructions in two and three spatial dimensions. In 2D, concentric circles, L-shape, and ribbon domains were examined; in 3D, the tests included a perturbed sphere, torus, and four disjoint spheres. The paper states that CCBM converges robustly to exact solutions, sometimes outperforming the Kohn–Vogelius method for coarse meshes and in computational time per iteration, and that the moving-mesh strategy avoided remeshing [2205.12620].

For Stokes free boundary problems, the 2D experiments showed that CCBM and classical tracking Dirichlet data converge to essentially the same optimal shape and display similar cost and gradient norm decay. In 3D, CCBM yielded smoother and more stable optimal surfaces, better preserved topology, and showed greater robustness and faster convergence for the same computational cost [2302.11828]. For immersed-obstacle detection in Stokes flow, the method reconstructed convex, nonconvex, sharply concave, and cornered obstacles in both two and three dimensions, remained robust under Gaussian noise up to 30%, and did so without perimeter or volume functional penalization [2403.11819].

In advection–diffusion shape recovery, standard CCBM was reported to faithfully recover convex and mildly non-convex shapes in noise-free scenarios, while struggling with concavities under noisy data and under the combined effects of advection and diffusion. The ADMM-modified CCBM substantially improved reconstruction for non-convex and concave shapes, worked without explicit perimeter regularization, and was described as robust to moderate noise up to 5–7% Gaussian noise in the reported tests. Cost functional histories showed monotonic decrease and stabilization, with gradient norms confirming convergence [2508.16898].

For diffusion coefficient recovery, the modified CCBM with the weighted \(H^1\)-type misfit was observed to yield more stable reconstructions and to reduce certain high-frequency artifacts when the \(H^1\)-weight was sufficiently large but not excessive. Across the reported numerical scenarios, it often showed favorable stability and robustness relative to several classical boundary-based formulations, and the projection-based extension supported stable recovery of piecewise-constant diffusion coefficients when all subdomains shared a portion of the boundary [2602.01656].

In BLT, the parameter-dependent CCBM was combined with perimeter and volume regularizations and a level-set representation. The reported 2D experiments included smooth, non-convex, polygonal, closely situated, and nested sources. The method reconstructed both support and intensity, handled multiple and nested sources, and achieved lower area and intensity errors than the comparison method labeled “G.Z.”, while remaining robust under noisy boundary data [2509.09533]. In cavity reconstruction with a homogeneous Robin condition, unconstrained CCBM gave good recoveries for smooth or convex cavities and large \(\alpha\), whereas the ADMM-constrained version reduced sensitivity to initialization and noise and improved reconstructions for challenging concave geometries [2605.12202].

For inverse Cauchy completion, the CCBM-GKB algorithm was proved to be a regularization method and was described numerically as much faster than the classic Landweber method. The numerical discussion states that the required number of iterations for a desired error is an order of magnitude smaller for CCBM-GKB than for CCBM-Landweber [2204.03416].

## 6. Limitations, interpretive issues, and related terminology

The published CCBM literature does not claim that the complex reformulation removes ill-posedness. In the Bernoulli analysis, compactness of the shape Hessian is proved and linked directly to instability typical of free boundary inverse problems [2205.12620]. In cavity reconstruction with a homogeneous Robin condition, it is stated explicitly that a single measurement may correspond to infinitely many admissible domains, even though shape optimization can still yield reasonable reconstructions [2605.12202]. A plausible implication is that CCBM regularizes the computational treatment of such problems without converting them into uniquely solvable inverse maps.

Several limitations recur across applications. In the Stokes obstacle paper, reconstructions degrade for very small obstacles far from the boundary or for multiple distant inclusions, and the Lagrangian boundary-deformation approach cannot change topology unless one introduces tools such as topological gradients [2403.11819]. By contrast, the BLT framework uses a level-set representation and therefore naturally accommodates topological changes such as splitting, merging, and nesting [2509.09533]. In diffusion coefficient recovery, performance is explicitly problem- and parameter-dependent, and excessive \(H^1\)-weights may oversmooth or slow convergence [2602.01656]. In the Robin-cavity and BLT settings, the efficacy of the method depends on algorithmic parameter tuning and on the richness of the boundary data; BLT further links the coupling parameter and Tikhonov parameter through the condition \(\alpha=\mathcal{O}(\sqrt{\varepsilon})\) to obtain uniform boundedness of the reconstructed intensity [2605.12202][2509.09533].

The Stokes free boundary paper notes additional practical issues: CCBM requires a PDE solver that can handle complex-valued systems efficiently, and the relation between the imaginary part and the physical quantities is indirect, which can make diagnostics less intuitive [2302.11828]. These are method-intrinsic consequences of replacing an overdetermined real problem by a single complex one.

A nomenclature issue also arises. In the inverse-problem literature discussed here, CCBM denotes the **Coupled Complex Boundary Method**. This should not be conflated with the distinct “complex-boundary treatment method” for cut-Cartesian finite-volume schemes [2304.04216], nor with the “fully coupled-channel complex scaling method” used for the \(K^-pp\) system in few-body quantum physics [1702.08002]. The shared use of “complex boundary” or “coupled complex” language does not indicate methodological identity.

Source: https://www.emergentmind.com/topics/coupled-complex-boundary-method-ccbm