I am developing code to simulate a multi-component Shan-Chen-based flow in a channel. My code works perfectly for single-component flows; However, whenever I switch to multi-component flows, corner nodes trouble a lot. Please help in this regard. How to calculate Shan-Chen forces at boundary nodes? How to deal with corner nodes? Where are the corner nodes for the velocity inlet? Where are the corner nodes for the top and bottom walls? For the outlet domain, what nodes should be considered?