Contact angle in ShanChen multi component model
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Contact angle in ShanChen multi component model
 This topic has 4 replies, 2 voices, and was last updated 2 years, 10 months ago by savis.

AuthorPosts

August 30, 2019 at 5:04 am #4532zshi6193Participant
Dear OpenLB,
In ShanChen multicomponent model, I use the following methods to set up the contact angle:
BounceBack<T, DESCRIPTOR> bounceBackRhoPhase1( .2);
BounceBack<T, DESCRIPTOR> bounceBackRhoPhase2( .8);The sum of the above two fictitious density is 1. To change them accordingly, different contact angles can be achieved. But the question is that simulation is hard to converge. For example, if I set the fictitious density for phase 2 is 0.8. Initially, it can form a droplet shape on the surface, but it continues to wet the surface until becoming flat.
I am just wondering whether I miss something when setting up the contact angle?
Thanks
Best regards,
SimonAugust 30, 2019 at 7:03 am #4533zshi6193ParticipantBased on the paper (Proposed approximation for contact angles in ShanandChentype multicomponent multiphase lattice Boltzmann models, 2007), is it possible to add the adhesion force between fluid and solid for the consideration of contact angle by using “sLatticephase1.addLatticeCoupling”?
Thanks
Best regards,
SimonSeptember 3, 2019 at 6:25 pm #4534savisParticipantDear Simon,
I have taken a look at the paper and I think I have some idea of how this can be implemented. Unfortunately it isn’t particularly easy.
There are several ShanChen lattice couplings (such as ShanChenForcedGenerator2D, for which the force is calculated here), however these all only include the fluidfluid interactions and they will need to be edited to allow fluidsolid interactions.
Firstly, one (or two) new G interaction parameters need to be provided for the solid interaction.
Secondly, the extra contributions to the force need to be evaluated (similar to lines 133/134 in the linked code).
Finally, the extra terms would need to be included in the calculation of the force (lines 148/150).For the extra contribution to the force an indicator for the solid boundary is needed. However I am not sure how this could be accessed within the lattice coupling (I believe you need the superGeometry). If anyone else has any idea about how to do this that would be appreciated.
Best wishes,
SamSeptember 23, 2019 at 8:58 am #4551zshi6193ParticipantDear Sam,
Thanks for your time to look into the problem. Appreciate.
I propose a new way for this problem. If the ShanChen model allows three components, one of components could be the boundary. A new G can be introduced to define the interaction between the boundary and the fluid. In this way, the contact angle can be achieved. Compared to define the contact angle via virtual wall density, the new proposed way can yield numerical stable results. Would you think that it can work in OpenLB?Many thanks
Best regards,
Simon.September 24, 2019 at 12:53 pm #4571savisParticipantDear Simon,
To my knowledge, there is no three component implementation of the ShanChen model in OpenLB. So this approach would likely involve much more work than simply adding a solid interaction to the existing method.
Best wishes,
Sam 
AuthorPosts
 You must be logged in to reply to this topic.