Gallium melt flow
OpenLB – Open Source Lattice Boltzmann Code › Forums › on Lattice Boltzmann Methods › General Topics › Gallium melt flow
- This topic has 10 replies, 3 voices, and was last updated 2 years, 3 months ago by Alex.
-
AuthorPosts
-
May 18, 2022 at 12:25 pm #6573AlexParticipant
Hello everyone!
I am working on natural convection simulation via the hybrid lattice Boltzmann method. I have solved the benchmark problem of natural convection of air (Prandtl number Pr=0.71) inside the closed differentially heated cavity. In this case, the code works fine even with high Rayleigh numbers. However, when I change the Prandtl number to 0.02 which corresponds to Gallium melt, the code is numerically unstable. I have tried different collision models such as two relaxation time approximation, multiple relaxation time approximation, high order regularization procedure and even LES formulation, but still it has no effect on stability. Have someone ever encountered the same problem?
Thanks,
AlexMay 20, 2022 at 10:40 am #6576mathiasKeymasterHi Alex,
some hints:
– choose larger git sizes, smaller dx
– uMax should stay below 0.12, if choose smaler time step sizes or set smaller uMax (depending on the used unit converter)
– check the source of unstability (stop simulation if uMax>0.3 for the first time in the main loop and write results, output the velocity and locate where the great values orrur.Best
MathiasMay 20, 2022 at 11:50 am #6577AlexParticipantDear Prof. Mathias,
Many thanks for your hints. Grid size variation has no effect on stability. I thought about to capture the moment when crucial error occurs. But what I should do next I don’t know since, for example, in BGK model relaxation time is constant all over the domain. An interesting note that when I apply magnetic field with high Hartmann number, code becomes stable.
kind regards,
AlexMay 20, 2022 at 12:08 pm #6578mathiasKeymasterok, in that case
– check the source of unstability (stop simulation if uMax>0.3 for the first time in the main loop and write results, output the velocity and locate where the great values orrur.
May 21, 2022 at 6:01 am #6579AlexParticipantMany thanks. I have found that the problem of low and high Pr flow simulation via LBM is partially discussed in [https://doi.org/10.3390/fluids6040148]. As Prof. Mathias noted, relaxation time needs to be non-uniformly treated for the grid points.
May 21, 2022 at 8:08 pm #6580mathiasKeymasterYes, this was also done by Max Gaedtke (google for his PhD thesis). There is a link to the code in the appandix. Best Mathias
June 2, 2022 at 6:41 am #6602AlexParticipantUnfortunately, i could not find the thesis of Dr. Gaedtke. But I found the paper of his with the TRT model for gallium and parafin melting. However, the low Pr treatment is not discussed in this work. If someone encounteres the same problem, you can refer to the work of Bawazeer et. al. [https://doi.org/10.1016/j.ijheatmasstransfer.2019.118562]. The authors added a correction term to the force term in order to make the model stable with low Prandtl numbers.
June 2, 2022 at 4:24 pm #6603mathiasKeymasterGaedtke thesis: https://publikationen.bibliothek.kit.edu/1000126904 .
Interesting paper! If you find out more, you could give a short talk in our seminar. That might also be interesting for Fedor, talk with him about it.
Best
MathiasJune 2, 2022 at 4:26 pm #6604jjessbergerParticipantDear Alex,
thanks for your comment. Here, you also find the thesis: https://publikationen.bibliothek.kit.edu/1000126904
Yours, Julius
June 3, 2022 at 4:02 am #6611AlexParticipantDear Prof. Mathias,
Thank you for the link!
kind regards,
AlexJune 3, 2022 at 4:02 am #6612AlexParticipantDear Julius,
Many thanks for the link!
kind regards,
Alex -
AuthorPosts
- You must be logged in to reply to this topic.