No label defined (Q4947176)

A numerical algorithm is considered for the plain problem on a convective viscous fluid flow in the domain \(G=\{0\leq x\leq L;\;0\leq y\leq f(x)\}\). Boussinesq equations with the variables \(\psi\), \(\omega\), \(\theta\) are used for describing the flow, where \(\psi\) is the stream function, \(\omega\) is the function of a vortex of a velocity, and \(\theta\) is the temperature. An explicit form of the boundary conditions for \(\psi\) and \(\omega\) on the free boundary \(y=f(x)\) are obtained. The Thom formula is used for finding the vortex on the right-hand walls \(x=0\), \(x=L\), and \(y=0\). Using the transformation \(x=\xi\), \(y=f(\xi)\eta\), the author maps the domain \(G\) with the moving boundary \(y=f(x)\) onto the rectangle \(0\leq\xi\leq L\), \(0\leq\eta\leq 1\) in which the problem is solved numerically by an iterative scheme of stabilizing correction. The results obtained are compared with the results previously published by \textit{S.~Chippada}, \textit{T.~C.~Jue}, and \textit{B.~Ramaswamy} [Int. J. Numer. Methods Eng. 38, No. 2, 335-351 (1995; Zbl 0823.76046)]. The Stephan phase transition problem is also considered. The domain \(G\) is divided by a surface \(S(t)\) into two subdomains \(G_1\) and \(G_2\) which are occupied by the liquid and rigid phases of the substance, respectively. In the whole domain \(G\), the problem is solved by the fictitious domains method. The above-described algorithm is applied to a domain occupied by the liquid phase. The shock-capturing method allows us to determine the surface \(S(t)\) automatically from the value of the temperature. Results of calculations of melting and crystallization processes are represented.