Numerical simulation of unsteady flows in Czochralski crystal growth by lattice Boltzmann methods

详细信息    查看全文

• 作者：Haibo Huang^a ; ^{huanghb@ustc.edu.cn" class="auth_mail} ; Xi-yun Lu^a ; Manfred Krafczyk^b
• 关键词：Lattice Boltzmann ; Thermal ; Axisymmetric ; Crystal growth ; Convection&ndash ; diffusion
• 刊名：International Journal of Heat and Mass Transfer
• 出版年：July 2014
• 年：2014
• 卷：74
• 期：Complete
• 页码：156-163
• 全文大小：753 K

文摘

This study presents model extensions for a lattice Boltzmann (LB) approach to thermal axisymmetric flow including swirl or rotation. An incompressible axisymmetric lattice Boltzmann D2Q9 model was applied to solve the axial and radial velocities through inserting source terms into the two-dimensional lattice Boltzmann equation. The equations governing azimuthal (or swirling) velocity and the temperature were also solved by the LBM. It is found that this scheme is much more stable and consistent compared to previous hybrid schemes. It provides a significant advantage in simulation of melt flows with high Reynolds number and high Grashof number. The present scheme was validated by comparing the LB results with benchmark solutions for melt flow in Czochralski crystal growth. Unsteady flows with high Grashof numbers were studied in detail. The critical Grashof number for the onset of the oscillation is found to be about b=MathURL&_method=retrieve&_eid=1-s2.0-S0017931014002087&_mathId=si53.gif&_user=111111111&_pii=S0017931014002087&_rdoc=1&_issn=00179310&md5=0b06519c4a07517bae481ff8709a56bf" title="Click to view the MathML source">2.5×10⁶$2.5\times {10}^6$. The oscillation amplitude b=MathURL&_method=retrieve&_eid=1-s2.0-S0017931014002087&_mathId=si54.gif&_user=111111111&_pii=S0017931014002087&_rdoc=1&_issn=00179310&md5=f6c3b87ef4714c6d173454989ec637aa" title="Click to view the MathML source">蠄b>maxb> is proportional to b=MathURL&_method=retrieve&_eid=1-s2.0-S0017931014002087&_mathId=si55.gif&_user=111111111&_pii=S0017931014002087&_rdoc=1&_issn=00179310&md5=37f4c77fc2a3be4b1b345c744fa728b8" title="Click to view the MathML source">(Gr-Grb>cb>)^0.5 for b=MathURL&_method=retrieve&_eid=1-s2.0-S0017931014002087&_mathId=si56.gif&_user=111111111&_pii=S0017931014002087&_rdoc=1&_issn=00179310&md5=7e16163d5744c12d97fa45d24a3b57a0" title="Click to view the MathML source">2.5×10⁶<Gr<6×10⁶$2.5\times {10}^6<\mathit{Gr}<6\times {10}^6$. The frequencies and flow patterns of the unsteady flows are also analyzed. The distributions of the mean quantities of the temperature and rms of temperature at Grashof number as high as b=MathURL&_method=retrieve&_eid=1-s2.0-S0017931014002087&_mathId=si57.gif&_user=111111111&_pii=S0017931014002087&_rdoc=1&_issn=00179310&md5=896567e1a7fefdd49c585af947b541bb" title="Click to view the MathML source">6×10⁷$6\times {10}^7$ is found to be similar to those obtained by 3D simulations.