Overview
A second-order diffusive-interface immersed boundary method has been developed for simulating incompressible flow, specifically addressing scenarios involving phase change and moving interfaces. This method introduces a smooth extension strategy to restore formal second-order spatial accuracy, particularly in regions near phase-changing boundaries where traditional methods often exhibit reduced accuracy.
Research Context
Accurate resolution of interfacial gradients is a critical requirement for simulating two-phase flows, especially in applications that involve phase transitions or active matter. Diffuse-interface immersed boundary methods (IBMs) are known for their efficiency in addressing such problems. However, a common limitation of these traditional methods is a reduction to first-order accuracy when applied near phase-changing boundaries. This reduction in accuracy is attributed primarily to local derivative discontinuities.
The problem of achieving higher-order accuracy in the presence of boundaries or interfaces has been a focus in computational fluid dynamics. The efficiency of IBMs makes them attractive, but the accuracy degradation near interfaces, particularly those undergoing phase change, presents a computational challenge that impacts the reliability of simulations for complex interfacial phenomena.
Approach
The proposed method introduces a smooth extension strategy designed to restore formal second-order spatial accuracy. The core of this strategy involves extrapolating the scalar field across the interface. This extrapolation is structurally designed to ensure continuity of derivatives. The researchers clarified that the root cause of the accuracy reduction in traditional IBMs near phase-changing boundaries is the presence of local derivative discontinuities.
To maintain the divergence-free condition, which is fundamental in incompressible fluid solvers, this smooth extension protocol is applied exclusively to the scalar transport equations. The velocity field within the framework maintains the standard treatment associated with diffuse-interface methods. This differential application ensures that the modifications for scalar accuracy do not inadvertently compromise the fundamental properties of the incompressible flow solver.
Findings
- The main reason for the reduction to first-order accuracy in traditional diffuse-interface immersed boundary methods near phase-changing boundaries was identified as local derivative discontinuities.
- A smooth extension strategy was successfully developed and implemented to restore formal second-order spatial accuracy by extrapolating the scalar field across the interface.
- This smooth extension strategy structurally ensures derivative continuity.
- The smooth extension was applied exclusively to the scalar transport equations to preserve the divergence-free condition in incompressible fluid solvers.
- The velocity field in the proposed method retains the standard diffuse-interface treatment.
- The framework was systematically validated against classical phase-change benchmarks, specifically one-dimensional evaporation and boiling problems.
- The method was also applied to the spontaneous autophoretic motion of isotropic particles.
- Numerical results confirmed the capability of the developed method in resolving complex multi-physics boundary couplings.
Why This Matters
The ability to accurately resolve interfacial gradients with higher-order spatial accuracy in two-phase flows is crucial for realistic simulations of processes involving phase transitions or active matter. The developed method addresses a identified limitation of traditional techniques, offering enhanced fidelity for a range of complex physical phenomena where interfacial dynamics are key.