HLBM: particleparticle interaction and fluidparticle momentum exchange
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › HLBM: particleparticle interaction and fluidparticle momentum exchange
 This topic has 2 replies, 2 voices, and was last updated 1 month, 1 week ago by Junwei Guo.

AuthorPosts

October 7, 2020 at 12:07 am #5205Junwei GuoParticipant
Dear OpenLB developers:
I am pretty interested in your HLBM method for resolved fluidparticle flow. Let start the discussion with your example code dkt2d and your publication [1].
1. Hydrodynamic force from fluid on particle. If I am correct, the force acting on particle is computed via “void ParticleDynamics3D<T, DESCRIPTOR>::computeBoundaryForce()” in the source file “olb1.3r1/src/particles/hlbm/hlbmDynamics3D.hh”. I check the code “BlockLatticePorousMomentumLossForce3D<T, DESCRIPTOR>::BlockLatticePorousMomentumLossForce3D()” in source file “olb1.3r1/src/functors/lattice/blockLatticeLocalF3D.hh”. There is a loop over a particle covered lattice cell to compute the momentum exchange:
// compute momentum exchange force on particle
T tmpForce[3] = {0.,0.,0.};
tmpForce[0] += this>_blockLattice.get(iX, iY, iZ).template getFieldPointer<descriptors::VELOCITY_NUMERATOR>()[0];
tmpForce[1] += this>_blockLattice.get(iX, iY, iZ).template getFieldPointer<descriptors::VELOCITY_NUMERATOR>()[1];
tmpForce[2] += this>_blockLattice.get(iX, iY, iZ).template getFieldPointer<descriptors::VELOCITY_NUMERATOR>()[2];I know the code is summing the numerator (Ref. [1], equation below Eq. 6) to recover the hydrodynamic force. From Ref. [1], the hydrodynamic force is computed from Eq. 7. The problem is, how does the numerator related to the bounceback equation (equation above Eq. 7 from Ref. [1]), how it is implemented in OpenLB code?
2. Particleparticle interaction. I talked to Dr. Krause. The code now is limited to hard sphere collision. I am fine with that. But how is it implement in HLBM? I did a test. When two particle is in contact. The numerator is contacted zone is very high, so the particle can separate. My question is how is the numerator is related to particle collision?
Thanks a lot. I am developing a particleflow model that may discuss latter.
Regards,
Junwei Guo (PhD candidate)[1] Krause, M. J.; Klemens, F.; Henn, T.; Trunk, R. & Nirschl, H.
Particle flow simulations with homogenised lattice Boltzmann methods
Particuology, Elsevier, 2017, 34, 1–13October 12, 2020 at 9:28 am #5224NicolasParticipantDear Junwei Guo,
after intruducing HLBM to OpenLB in 2017 there have been various minor, but also some major changes, including the momentum exchange algorithm. A lot of them are documented in a more recent publication on HLBM and the OpenLB user guide. You are strongly advised to contemplate both. Another more uptodate publication is being worked on at the moment, however no specific date for it can be given for now. Here, some further changes to the algorithm will be discussed, which will eventually be part of our next release, which is scheduled for the end of this very month. In case you are looking for the greatest consistency between publications and the implementation, you might need to await the upcoming release and our next HLBM related publication. Another option is of course to attend our next OpenLB spring school, which will give you the opportunity to talk with the developers directly.
I hope, this is of any help for now.Best,
NicolasOctober 18, 2020 at 11:27 pm #5234Junwei GuoParticipantDear Nicolas,
I am looking forward to attending the next OpenLB spring school. That would happen next year. Can you please contact the developers to solve this problem?
Thanks again.
Regards,
Junwei GUo 
AuthorPosts
 You must be logged in to reply to this topic.