Periodic Outlet Boundary Conditions ( nonequilibrium extrapolation method)
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Periodic Outlet Boundary Conditions ( nonequilibrium extrapolation method)
- This topic has 35 replies, 5 voices, and was last updated 12 months ago by jan.
-
AuthorPosts
-
March 29, 2023 at 2:44 pm #7324Henderson_CParticipant
Hello everyone,
I was wondering how OpenLB tackles the migration of certain object (droplet, bubbles, ….) through the domain outlet. I read some papers, and many has recommended the use of nonequilibrium extrapolation method. Is there an implementation of such method in OpenLB?
Regards,
March 29, 2023 at 4:07 pm #7325FBukreevKeymasterHello,
we have particle counter by the Lagrangian particles and the mass counter by Eulerian particles at the outlet. Which one do you want to use?
Greetings
FedorMarch 29, 2023 at 4:09 pm #7326Henderson_CParticipantDear Fedor,
Thank you for replying. Can you elaborate on the difference between them?
Thanks
March 29, 2023 at 4:16 pm #7327FBukreevKeymasterIn first case you track the separate points moving with Newtonian equation system and in the Eulerian case you regard particles as a concentration cloud.
March 29, 2023 at 4:19 pm #7328Henderson_CParticipantDear Fedor,
Can you show me both implementations as I am now indecisive of which method to use.
March 29, 2023 at 4:20 pm #7329FBukreevKeymasterLook in the examples/particles/bifurcation/eulerLagrange and eulerEuler. There you can see it clear.
March 29, 2023 at 8:13 pm #7332Henderson_CParticipantDear Fedor,
I want the particle to move to the flow inlet when it passes through the flow domain outlet. I have set the periodicity to true, but when the particle touches the outlet, it sticks to the outlet wall.
regards,
- This reply was modified 1 year, 8 months ago by Henderson_C.
March 29, 2023 at 8:33 pm #7334FBukreevKeymasterWhich particles do you use? Lagrangian or Eulerian? If you have set the periodicity on the particle Lattice there must be no outlet and inlet at all.
March 29, 2023 at 8:49 pm #7335Henderson_CParticipantDear Fedor,
I am not using the super particle system. I am simply using the single particle simulation. The particle is flowing through a channel with an inlet and outlet. To achieve fully hydrodynamic condition, I set the periodicity in the x-direction to be true (cuboidGeometry.setPeriodicity(false, true)). I thought once I enable this feature, any flow interruption due to the particle motion (Brownian motion) beside the outlet will be shown in the inlet of the channel. But that is not the case.
March 30, 2023 at 10:10 am #7336FBukreevKeymasterHello,
dou you simulate the resolved moving particle?
March 30, 2023 at 12:39 pm #7337Henderson_CParticipantDear Fedor,
Yes, similar to one used in dkt2d example.
regards,
- This reply was modified 1 year, 8 months ago by Henderson_C.
March 30, 2023 at 1:54 pm #7339FBukreevKeymasterYou can add the periodicity vector (false,true) at the end of the ParticleManager, look also at the ParticleManger code in src.
March 30, 2023 at 2:09 pm #7340Henderson_CParticipantDear Fedor,
By that you mean modify the source code (src/particles/particlemanager) with the periodicity vector?
regards,
March 30, 2023 at 2:24 pm #7341FBukreevKeymasterI mean that you can look at the parameters of the particlemanager in the src and than put the periodicity vector as argument on the right place in the particlemanager declaration in your simulation setup.
March 30, 2023 at 8:20 pm #7342Henderson_CParticipantDear Fedor,
I am not developer not I am expert in c++ but to understand correctly, I need to modify the ParticleManager constructor to accept the periodicity vector as an argument, and initialize the periodicity member variable with it. Is that correct? When you mention the parameters, what do you mean exactly? I am really confused on how to solve my issue.
Your continuous support is really appreciated.
Regards, -
AuthorPosts
- You must be logged in to reply to this topic.