---
title: 'LCRM: Hybrid Interface Tracking for Multiphase Flow'
url: https://www.emergentmind.com/topics/level-contour-reconstruction-method-lcrm
type: topic
---

# LCRM: Hybrid Interface Tracking for Multiphase Flow

Searching arXiv for Level Contour Reconstruction Method and related papers to ground the article.
Found relevant papers:
- 2509.24722 — "Direct numerical simulation of two-phase flows with surfactant-induced surface viscous effects"
- 1702.02478 — "An Interface-Tracking Technique for Multiphase Flow with Soluble Surfactant"
- 0901.1773 — likely Shin et al. IJNMF 2009 LCRM-related preprint/result

The Level Contour Reconstruction Method (LCRM) is an interface–tracking scheme for three-dimensional multiphase flow that combines an explicit Lagrangian interface with an Eulerian distance function used for periodic interface reconstruction. In the cited formulations, LCRM is presented as a hybrid front-tracking/level-set approach that retains sharp interfaces, accurate curvature, and good mass conservation while handling breakup and coalescence without explicit connectivity bookkeeping. It has been used for surface-tension-dominated multiphase flow, insoluble and soluble surfactant transport, Marangoni stresses, and surfactant-induced surface viscous stresses within the Boussinesq–Scriven constitutive framework [1702.02478, 2509.24722].

## 1. Hybrid formulation and defining characteristics

LCRM was originally developed by Shin & Juric as a hybrid scheme that keeps an explicit Lagrangian interface but periodically reconstructs that interface from a distance or level-set field defined on the Eulerian grid. In the descriptions given for both soluble-surfactant and surface-viscous flow formulations, the method is neither a classical front-tracking algorithm with persistent global mesh connectivity nor a pure level-set method in which a signed-distance field is advected and reinitialized. Instead, the interface is advected explicitly as discrete Lagrangian triangular elements, and, when needed, a signed distance function is reconstructed from that interface and used to rebuild the interface as the zero level contour of the distance field in each cell [1702.02478, 2509.24722].

This hybridization addresses two persistent difficulties in three-dimensional interfacial CFD. Classical front tracking gives sharp interfaces and exact geometric control, but mesh-quality operations, remeshing, and connectivity management become difficult under large deformation and especially under breakup or coalescence. Pure level-set methods handle topology change elegantly and provide smooth normals and curvature, but they are susceptible to mass loss and represent the interface over several cells. LCRM is designed to avoid the connectivity management and mesh-quality operations of classical front tracking in 3D while also avoiding the mass loss and curvature errors inherent to pure level-set advection [1702.02478].

A common misconception is to identify LCRM with a standard level-set method. The cited formulations explicitly distinguish the two: the distance function is computed directly from the Lagrangian interface, and, in the surface-viscous implementation, the level-set role is geometric and topological rather than kinematic. The interface kinematics remain purely Lagrangian, with marker motion governed by
\[
\frac{d\mathbf x_f}{dt}=\mathbf u_f=\mathbf u(\mathbf x=\mathbf x_f),
\]
integrated with a second-order Runge–Kutta scheme [2509.24722].

## 2. Geometric representation, reconstruction, and surface operators

In three dimensions, the interface is represented by a triangulated Lagrangian surface. Each triangular element \(e\) is defined by three vertices, with area and normal given by
\[
A_e=\frac{1}{2}\,|\mathbf x_{02}\times \mathbf x_{03}|,\qquad
\mathbf n_f=\frac{\mathbf x_{02}\times \mathbf x_{03}}{|\mathbf x_{02}\times \mathbf x_{03}|}.
\]
In the soluble-surfactant formulation, the interfacial surfactant concentration \(\Gamma\) is stored at the center of each element; in the surface-viscous formulation, the surface carries \(\Gamma\), \(\sigma(\Gamma)\), \(\mu_s(\Gamma)\), \(\mu_d(\Gamma)\), and related interfacial stresses [1702.02478, 2509.24722].

The Eulerian component of LCRM is a signed distance field, denoted \(\varphi(\mathbf x)\) or \(\phi(\mathbf x)\), constructed directly from the piecewise planar Lagrangian interface rather than by solving a transport PDE. Reconstruction is performed cell-by-cell. Each Eulerian cell is subdivided into tetrahedra, and the zero isocontour inside each tetrahedron is intersected to create local triangles. The cited descriptions emphasize a tetra-marching procedure that ensures a unique, non-ambiguous reconstruction in each cell. In the surface-viscous implementation, reconstruction is performed typically every 25 time steps, with the frequency chosen such that the displacement of markers per step does not exceed the smallest grid spacing. This regular rebuilding keeps elements well-shaped, avoids stretching and clustering, and permits topology change to be inherited from the distance field [1702.02478, 2509.24722].

Although curvature is computed from the Lagrangian surface, normals for several local operations are often obtained from the distance function,
\[
\mathbf n(\mathbf x)=\frac{\nabla\varphi(\mathbf x)}{|\nabla\varphi(\mathbf x)|}.
\]
At edge centers, this normal is used to define the binormal direction \(\mathbf p\) and tangential direction \(\mathbf t\), which enter the computation of surface gradients of \(\Gamma\), \(\sigma(\Gamma)\), and, in the Boussinesq–Scriven implementation, velocity-derived surface tensors. This geometric split—Lagrangian curvature and area, Eulerian normals and local orientation—is one of the characteristic features of LCRM [1702.02478].

## 3. One-fluid formulation, surfactant transport, and Marangoni coupling

The flow is formulated on a fixed Eulerian grid using a one-fluid or single-field incompressible Navier–Stokes system with discontinuous properties. In the surface-viscous formulation,
\[
\rho(\mathbf x,t)=\rho_1\mathcal H(\mathbf x,t)+\rho_2\big(1-\mathcal H(\mathbf x,t)\big),\qquad
\mu(\mathbf x,t)=\mu_1\mathcal H(\mathbf x,t)+\mu_2\big(1-\mathcal H(\mathbf x,t)\big),
\]
while the soluble-surfactant formulation uses an indicator \(I(\mathbf x,t)\) in an equivalent single-field representation. Incompressibility is enforced through
\[
\nabla\cdot\mathbf u=0.
\]
Interfacial forces enter the momentum equation through sharply localized force densities distributed from the Lagrangian interface to the Eulerian grid by a discrete Peskin-type kernel [1702.02478, 2509.24722].

For surface tension and Marangoni stresses, the interfacial force is written as the sum of a normal curvature term and a tangential term,
\[
\mathbf F(\mathbf x)=\int_\Gamma \sigma\,\kappa\,\mathbf n\,\delta(\mathbf x-\mathbf x_f)\,dA
+\int_\Gamma \nabla_s\sigma\,\delta(\mathbf x-\mathbf x_f)\,dA.
\]
The normal part is implemented through a hybrid curvature formulation that aligns the discrete curvature force with the pressure gradient operator and reduces parasitic currents. The tangential part is assembled edge-by-edge after decomposing \(\nabla_s\sigma\) into binormal and tangential components [1702.02478].

In the insoluble-surfactant formulation, the surface concentration satisfies
\[
\frac{D_s\Gamma}{Dt}=D_\Gamma \nabla_s^2\Gamma,
\]
with the surface material derivative defined by
\[
\frac{D_s(\cdot)}{Dt}=\frac{\partial(\cdot)}{\partial t}+\nabla_s\big((\cdot)\mathbf u_s\big).
\]
In the soluble-surfactant extension, the interfacial conservation law is written as
\[
\frac{\partial \Gamma}{\partial t}+\nabla_s\cdot(\Gamma \mathbf u_s)=D_s\nabla_s^2\Gamma+S_\Gamma,
\]
where the adsorption/desorption source is
\[
S_\Gamma=k_a C_s(\Gamma_\infty-\Gamma)-k_d\Gamma.
\]
Bulk concentration in one phase evolves through
\[
\frac{\partial C}{\partial t}+\mathbf u\cdot\nabla C=\nabla\cdot(D_b\nabla C),
\]
and is coupled to the interface by the sharp Neumann condition
\[
\mathbf n\cdot\nabla C\big|_\Gamma=-\frac{S_\Gamma}{D_{b2}}.
\]
In practice this boundary condition is enforced through a ghost-cell method adapted from the “sharp energy” method of Shin & Choi [1702.02478].

The surface tension is surfactant-dependent. The cited formulations use a Langmuir–Szyszkowski equation of state,
\[
\sigma(\Gamma)=\sigma_0+RT\Gamma_\infty\ln\left(1-\frac{\Gamma}{\Gamma_\infty}\right)
=\sigma_0\left[1+\beta_s\ln\left(1-\frac{\Gamma}{\Gamma_\infty}\right)\right],
\qquad
\beta_s=\frac{RT\Gamma_\infty}{\sigma_0},
\]
and, in some soluble-surfactant benchmarks, a linear equation of state,
\[
\sigma=\sigma_s\left(1-\beta_s\frac{\Gamma}{\Gamma_\infty}\right).
\]
Marangoni stresses then arise from
\[
\mathbf f_M=\nabla_s \sigma(\Gamma).
\]
A notable numerical feature is that surface gradients of \(\Gamma\) and \(\sigma\) are computed with the same probing technique: the distance function is used to place edge-centered probe points along \(\mathbf p\) and \(\mathbf t\), and two-point finite differences provide directional derivatives while remaining on the reconstructed interface [1702.02478, 2509.24722].

## 4. Boussinesq–Scriven surface rheology within LCRM

The 2025 formulation extends LCRM to interfacial flows with surfactant-induced surface viscous stresses. From an extrinsic viewpoint, the interface is treated as a two-dimensional manifold embedded in three-dimensional space. The surface identity tensor is
\[
\mathbf I_s=\mathbf I-\mathbf n\mathbf n,
\]
and, for an inviscid interface, the surface-excess pressure tensor is
\[
\mathbf P_0^{\mathbf S}=\sigma \mathbf I_s.
\]
Its surface divergence yields
\[
\mathbf p_0^{\mathbf S}=\nabla_s\cdot(\sigma \mathbf I_s)=\nabla_s\sigma+\kappa \sigma \mathbf n,
\]
which recovers the normal capillary force and Marangoni tangential force [2509.24722].

The Boussinesq–Scriven constitutive law introduces surface shear and dilatational viscosities. The surface viscous stress is written as
\[
\mathbf P_\mu^s=
\left(\mu_d^{\mathbf S}-\mu_s^{\mathbf S}\right)(\nabla_s\cdot \mathbf u_f)\,\mathbf I_s
+2\mu_s^{\mathbf S}\mathbf D_s,
\]
where the adopted form of the surface rate-of-deformation tensor is the extrinsic–Lagrangian Scriven form,
\[
2\mathbf D_s=\nabla_s\mathbf u_f\cdot \mathbf I_s+\mathbf I_s\cdot(\nabla_s\mathbf u_f)^T.
\]
The total surface lineal force due to surface viscosity becomes
\[
\mathbf F_v=\mathbf F_v^d+\mathbf F_v^s
=\int_{A'} \nabla_s \cdot \left[(\mu_d-\mu_s)(\nabla_s \cdot \mathbf u)\mathbf I_s\right]\delta(\mathbf x-\mathbf x_f)\,dA'
+\int_{A'} \nabla_s \cdot \big(2\mu_s \mathbf D_s\big)\delta(\mathbf x-\mathbf x_f)\,dA' .
\]
This formulation is implemented directly on the Lagrangian interface and then transferred to the Eulerian grid through the same discrete delta machinery used for capillary and Marangoni forces [2509.24722].

For the dilatational contribution, the scalar
\[
\sigma_{vis}=(\mu_d-\mu_s)(\nabla_s\cdot \mathbf u)
\]
is defined as a surface-viscous tension. The corresponding force decomposes into a normal term and a tangential term,
\[
\mathbf F_v^d=
\underbrace{\int_{A_e}\sigma_{vis}\,\kappa\,\mathbf n\,\delta(\mathbf x-\mathbf x_f)\,dA_e}_{\mathbf F_n}
+
\underbrace{\int_{A_e}\nabla_s\sigma_{vis}\,\delta(\mathbf x-\mathbf x_f)\,dA_e}_{\mathbf F_t}.
\]
The normal component is cast in the same hybrid form used for capillary forces, while the tangential component uses probing to compute \(\nabla_s \sigma_{vis}\) edge-by-edge. The paper explicitly notes that the dilatational part is implemented using exactly the same machinery as capillary forces in LCRM, with \(\sigma\) replaced by \(\sigma_{vis}\) [2509.24722].

For the shear contribution, the force is converted to an edge integral through the surface Gauss theorem,
\[
\mathbf F_v^s
=
\int_C \left[ \mu_s \left( \nabla_s\mathbf u \cdot \mathbf I_s + \mathbf I_s\cdot(\nabla_s\mathbf u)^T \right) \right]\cdot \mathbf p \,\delta(\mathbf x-\mathbf x_f)\,dl,
\]
so accurate evaluation of \(\nabla_s\mathbf u\) at edge midpoints is required. The implementation exploits the fact that the full Eulerian gradient \(\nabla \mathbf u\) is already available on the grid. At an edge midpoint,
\[
\mathbf A=(\nabla_s\mathbf u)
=(\nabla\mathbf u)-\mathbf n\big(\mathbf n\cdot(\nabla\mathbf u)\big),
\]
with \(\mathbf n\) obtained from the distance function. Then
\[
\mathbf B=\mathbf A\cdot \mathbf I_s,\qquad
2\mathbf D_s=\mathbf B+\mathbf B^T,
\]
and the surface divergence required in \(\sigma_{vis}\) is obtained on the fly as
\[
\nabla_s\cdot\mathbf u=\mathrm{tr}(\mathbf A).
\]
This is one of the clearest examples of why the Lagrangian–Eulerian division of labor is useful: the interface remains sharp and explicitly represented, while the distance field supplies robust normals and the Eulerian grid supplies the full velocity gradient needed for surface rheology [2509.24722].

## 5. Numerical realization, force spreading, and parallel structure

The cited LCRM implementations use a staggered MAC grid and a projection method for incompressible Navier–Stokes. Convective terms are discretized with a second-order ENO scheme, and the pressure Poisson equation is solved with a parallel hybrid multigrid/GMRES solver that is reported as robust up to density ratios \(O(10^4)\). Interface points are advanced by second-order Runge–Kutta, and interfacial quantities are updated on the moving triangulated surface [1702.02478, 2509.24722].

Transfer between the Lagrangian surface and the Eulerian grid is handled through a compact-support discrete delta kernel. In the surface-viscous formulation, the kernel is written as
\[
D_{ijk}(\mathbf x_f)=
\frac{\delta(x_f/h_x-i)\,\delta(y_f/h_y-j)\,\delta(z_f/h_z-k)}{h_x h_y h_z},
\]
and any scalar or vector quantity \(X_f\) defined on the interface is distributed by
\[
M_{sijk}=\sum_f X_f D_{ijk}(\mathbf x_f)\Delta A_f.
\]
The soluble-surfactant formulation uses the same immersed-boundary style transfer to construct Eulerian fields \(M_s\) and \(Q\), after which
\[
\Gamma(\mathbf x_{ijk})=\frac{M_{s,ijk}}{Q_{ijk}}.
\]
This Eulerian distribution preserves surfactant mass, facilitates coupling to the bulk concentration equation, and provides a field from which \(\sigma(\Gamma)\) can be evaluated and re-interpolated after interface reconstruction [1702.02478, 2509.24722].

The probing technique is central to the numerical treatment of surface gradients. At an edge center, the normal from the distance function defines \(\mathbf p\) and \(\mathbf t\), and probe points are placed at a distance \(d_l\approx h\) in appropriate directions while remaining on the interface through the \(\varphi=0\) constraint. Directional derivatives such as
\[
(\nabla_s \Gamma)_{\mathbf p}\approx \frac{\Gamma_{\text{out}}-\Gamma_{\text{in}}}{2d_l},
\qquad
(\nabla_s \Gamma)_{\mathbf t}\approx \frac{\Gamma_{\text{right}}-\Gamma_{\text{left}}}{2d_l}
\]
are then used for diffusion, Marangoni stress, and, by direct analogy, gradients of surface-viscous tension. Because the method avoids connectivity-based stencils, these calculations remain local and are naturally adapted to distributed-memory parallelism [1702.02478, 2509.24722].

The code base discussed in the cited papers is BLUE. The surface-viscous study reports that BLUE is massively parallel and tested up to 131,072 cores, while the soluble-surfactant study emphasizes that interface operations such as distance computation, delta distribution, and probing are local plus halo exchanges. This locality is presented as a major advantage over classical front-tracking schemes that rely on global mesh adjacency, which are described as cumbersome to parallelize [1702.02478, 2509.24722].

## 6. Verification, applications, limitations, and broader significance

The published LCRM formulations are supported by a broad benchmark suite. For soluble surfactant transport, the uniformly expanding sphere test shows that numerical \(\Gamma(t)\) matches the analytical solution
\[
\Gamma(t)=\Gamma_0\left(\frac{R_0}{R(t)}\right)^2
\]
extremely well on a coarse \(64^3\) grid with 12 cells across the diameter and 1000 reconstruction steps, indicating very good total surfactant mass conservation and no degradation due to repeated interface reconstruction. For a stationary droplet with surface diffusion, one-dimensional profiles of \(\Gamma(\theta,t)\) match the exact solution closely even on \(64^3\), and grid refinement from \(64^3\) to \(256^3\) yields between first- and second-order convergence in \(L^1\) norm. Marangoni migration velocities converge toward the YGB prediction, and simple-shear drop deformation compares well with experiments and prior simulations; the reported rotation angle is 0.26 rad with surfactant versus 0.27 rad experimentally, and 0.65 rad without surfactant versus 0.61 rad experimentally. Total surfactant mass is conserved to within a few tenths of a percent on refined grids. The same formulation also handles adsorption, desorption, breakup into main and satellite drops, and large-scale annular film flow with insoluble and soluble surfactant [1702.02478].

For surface-viscous interfacial rheology, LCRM is validated against several canonical problems. In simple shear, mesh refinement with \(R/\Delta x=4,8,16\) shows convergence of the major and minor axes \(L,B\) and deformation \(D=(L-B)/(L+B)\), with errors less than a few percent relative to Flumerfelt’s small-deformation theory as the grid is refined. Temporal and steady deformations compare closely with Pozrikidis’ BEM results, Flumerfelt’s analytical expression in the Stokes limit, and previous front-tracking simulations. Increasing surface shear viscosity slows deformation and reduces the steady deformation parameter, while large \(Bq_s\) produces nearly rigid-like surface circulation and reduced tangential velocity gradients. For a rising surface-viscous drop in quiescent fluid, terminal velocity decreases with increasing \(Bq_d\), and fine-grid errors are reported as \(\lesssim 3\)–\(4\%\). In vertically oscillated two-layer systems, threshold acceleration \(F_c\) for Faraday instability increases with the total Boussinesq number and matches two-dimensional ALE results within about \(4\%\) across cases; separate effects of \(Bq_s\) and \(Bq_d\) are resolved, with slightly larger errors for shear viscosity. In Faraday-wave-induced atomization, the method reproduces deep craters, jet-like ligaments, and pinch-off events in the clean case, while added surface viscosity damps ligament formation, delays pinch-off, and changes droplet sizes and topology [2509.24722].

These results delimit the method’s current strengths and its acknowledged challenges. The cited comparisons state that LCRM offers excellent mass conservation, reported as less than \(0.1\%\) error, improved curvature estimation, and reduced parasitic currents, all of which are especially important when interfacial forces depend sensitively on curvature and surface gradients. At the same time, the surface-viscous study notes that accuracy of surface shear viscosity appears slightly more sensitive than dilatational viscosity, with Faraday-wave thresholds deviating by up to about \(4\%\), and that high resolution is required to reduce underestimation or overestimation in demanding DNS settings. The soluble-surfactant formulation also notes that explicit treatment of diffusion and surfactant can impose time-step constraints, and its bulk diffusion examples are limited to one phase, although extension to both phases is described as straightforward [1702.02478, 2509.24722].

Within the cited research trajectory, LCRM functions as a general three-dimensional interface framework in which surface PDEs, Marangoni forcing, and full Boussinesq–Scriven interfacial rheology can be implemented in a single Navier–Stokes solver. This suggests a methodological niche between implicit-capture methods and connectivity-heavy front-tracking schemes: surface fields remain naturally defined on a Lagrangian manifold, while topology change, probing, and geometric reconstruction are delegated to a locally reconstructed Eulerian distance field. A plausible implication is that this division of labor is the central reason LCRM has been extended from capillary flow to soluble surfactant transport and then to surfactant-induced surface viscosity without changing its core interface philosophy [1702.02478, 2509.24722].

Source: https://www.emergentmind.com/topics/level-contour-reconstruction-method-lcrm