An effective implementation of the bidirectional buffer: Towards laminar and turbulent open-boundary flows

To effectively handle flows characterized by strong backflow and multiple open boundaries within particle-based frameworks, this study introduces three enhancements to improve the consistency, independence, and accuracy of the buffer-based open boundary condition in SPHinXsys. First, to improve the buffer consistency, the hysteresis-based particle generation scheme is introduced to prevent the excessive particle addition induced by strong backflow. Secondly, the independence of the bidirectional buffer is enhanced through region-constrained and independent labeling schemes, which effectively eliminate buffer interference and erroneous particle deletion in complex open-boundary flows. Thirdly, the original zeroth-order consistent pressure boundary condition is upgraded to first-order consistency by introducing a mirror boundary treatment for the correction matrix near the open boundary. The three improvements are applicable to both laminar and turbulent flow simulations. The turbulent implementation is based on the rigorously validated weakly compressible smoothed particle hydrodynamics coupled with the Reynolds-averaged Navier-Stokes (WCSPH-RANS) method. Four test cases, including straight and U-shaped channel flows, a 2D plane jet, and the flow in a 3D self-rotational micro-mixer, are conducted to comprehensively validate the proposed improvements. Among these cases, the turbulent plane jet is successfully simulated at a moderate resolution within a very compact computational domain involving strong backflow, a condition that is usually challenging for mesh-based methods. The three improvements require only minor modifications to the code framework, yet they yield significant performance gains.