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 11 months, 4 weeks ago by Alex.
May 18, 2022 at 12:25 pm #6573
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?
AlexMay 20, 2022 at 10:40 am #6576
– 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.
MathiasMay 20, 2022 at 11:50 am #6577
Dear 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.
AlexMay 20, 2022 at 12:08 pm #6578
ok, 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 #6579
Many 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 #6580
Yes, this was also done by Max Gaedtke (google for his PhD thesis). There is a link to the code in the appandix. Best MathiasJune 2, 2022 at 6:41 am #6602
Unfortunately, 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 #6603
Gaedtke 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.
MathiasJune 2, 2022 at 4:26 pm #6604jjessbergerParticipant
thanks for your comment. Here, you also find the thesis: https://publikationen.bibliothek.kit.edu/1000126904
Yours, JuliusJune 3, 2022 at 4:02 am #6611
Dear Prof. Mathias,
Thank you for the link!
AlexJune 3, 2022 at 4:02 am #6612
Many thanks for the link!
- You must be logged in to reply to this topic.