Skip to content

Effect On Particle Due to Fluid Density Change

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Effect On Particle Due to Fluid Density Change

Viewing 6 posts - 1 through 6 (of 6 total)
  • Author
    Posts
  • #8111
    Jijo
    Participant

    Hello OpenLB community,

    I am currently working on thermally heated resolved particles.I have been investigating whether these particles are influenced by changes in fluid density, particularly in relation to the Boussinesq approximation. I experimented with such test case but I can not notice any effect on the particle’s velocity when including this effect. Any idea on such scenario?

    Regards,

    #8113
    jan
    Participant

    Dear Jijo,

    gravity is applied as a (constant) external acceleration, see accExt and apply_gravity in examples/particles/dkt2d. So the changes of the fluid density don’t affect the particle.

    Also, currently the external acceleration is defined for the whole particle system. Updating it within the ParticleManager would therefore update it for all particles. In your particular case, I think it would be best to manually apply the gravity at each timestep. For example, you could evaluate the density at the current particle position if an averaged density is sufficient. If you need different forces on the surfaces, then a more sophisticated approach would be necessary.

    Best regards,
    Jan

    #8117
    Jijo
    Participant

    Dear Jan,

    Thank you for your reply. I appreciate your constant support and help. The problem I am having is that the fluid itself is not flowing upward due to natural convection and I think such issue is due to the force field.

    It seems that the porousparticleBGK dynamics does not contain a force field and I believe that buoyancy forces are coupled using Boussinesq forces by the Guo forcing scheme. What do you think I should do in this case?

    Regards,

    #8119
    jan
    Participant

    Dear Jijo,

    if you solely need a force field for the coupling, then perhaps try PorousParticleKupershtokhForcedBGKdynamics and see if that works.

    Best regards,
    Jan

    #9327
    Jijo
    Participant

    Dear Jan,

    Thank you for your response. The PorousParticleKupershtokhForcedBGKdynamics scheme works perfectly; however, I am encountering discrepancies with other forcing schemes. Specifically, when applying PorousParticleGuoForcedBGKdynamics to the isothermal DKT problem, the particle settling behavior changes significantly compared to the original scheme. Do you have any idea into what might be causing this issue?

    For reference, I am testing with the DKT2D example code.

    I appreciate any guidance you can provide.

    Best regards,
    Anas

    • This reply was modified 4 months, 1 week ago by Jijo.
    #9397
    Christoph
    Participant

    Dear Jijo,

    the error in results of various forcing schemes was investigated by Trunk et al. with a case of a single settling sphere in https://www.mdpi.com/2079-3197/9/2/11 . In this paper they found an error of 5 to 8% in the cases they studied when using the Guo forcing scheme.

    The DKT2d case is more suszeptible to differences in the forcing scheme, as the settling spheres interfere with each other. This likely causes the large discrepancies you found.

    Best regards,
    Christoph

Viewing 6 posts - 1 through 6 (of 6 total)
  • You must be logged in to reply to this topic.