Decoupled two‐grid finite element method for the time‐dependent natural convection problem I: Spatial discretization

In this article, a decoupled two grid finite element method (FEM) is proposed and analyzed for the nonsteady natural convection problem using the coarse grid numerical solutions to decouple the nonlinear coupled terms, and the corresponding optimal error estimates are derived. Compared with the standard Galerkin FEM and the usual two‐grid FEM, our algorithm not only keeps good accuracy but also saves a lot of computational cost. Some numerical examples are provided to verify the performances of the decoupled two‐grid FEM. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the decoupled two‐grid FEM for the nonsteady natural convection problem. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2135–2168, 2015     

Formation of longitudinal aggregation of inclusions in bulk sapphire and yttrium‐aluminum garnet grown by horizontal directed crystallization method

It is shown that the formation of longitudinal aggregation of inclusions in bulk sapphire and yttrium‐aluminum garnet (YAG) grown by the method of HDC is caused by local accumulation of impurities, disturbance of morphological stability of the crystallization front and capture of inclusions and impurities in the nodal region of the melt two‐vortex convection. Studied is the influence of thermal and geometrical parameters of the melt and the shape of the crystallization front on the conditions of the formation of the capture of inclusions.     

竖壁传热对方腔内Benard对流的影响

本文采用具有QUICK差分格式的SIMPLE算法对边界竖壁传热的方腔内空气Benard对流进行了数值计算,根据计算结果探讨了竖壁传热对Benard对流的影响.计算表明,在所考虑的几何和物理条件下,所有竖壁绝热时,腔内流体形成平行与短轴的多个涡卷;竖壁存在传热时,腔内流体形成平行于长轴的两个涡卷,并且平行于长轴竖壁的传热热流方向相反时,涡卷的旋转方向也相反;垂直于长轴蛏壁的传热对近壁附近的流动有一定影响.     

Finite element modelling of natural-convection-controlled change of phase

The problem of phase change in the presence of natural convection has been investigated. A model has been proposed based on the treatment of the release/absorption of latent heat as a heat source/sink in combination with the standard Galerkin finite element method with a primitive variable formulation on a fixed grid. To demonstrate the capabilities of the model, three cases of phase change of an aluminium alloy in the presence of natural convection arc considered, i.e. solidification, melting and combined solidification and melting. The solidification of water in a square cavity is modelled as another example, taking into account the density extremum, and the results are compared with a previously published work.     

Energy-based decomposition of friction drag in turbulent square-duct flows

The generation of friction drag in turbulent duct flows has direct connection with statistical quantities and corresponding turbulence dynamics in the duct cross-section. In this study, we generalize the RD identity (Renard and Deck, 2016) to a ‘two-dimensional’ form which we exploit to decompose the mean friction drag in turbulent square-duct flows into contributions associated with viscosity, turbulence and cross-stream convection. The friction Reynolds number of the duct flows ranges from 220 to 2000. The scaling, spatial distribution and local normalization of the contributions to friction are investigated and compared with those in pipe and channel flows. As in other canonical flows, we find logarithmic growth of the turbulent contribution in contrast to the viscous one, the former thus becoming dominant at high enough Reynolds numbers. Whereas cross-stream convection has no net effect on friction, its contribution may be locally comparable to the other two, hence may be responsible for redistribution of friction along the duct perimeter.     

环形浅液池内热毛细对流的热力学特性

为了了解径向温度梯度作用下环形浅液池内硅熔体热毛细对流的热力学特性,利用有限差分法进行了非稳态三维数值模拟。液池外半径r0=50 mm,内半径ri=15 mm,深度为d=3 mm。结果表明,当温度梯度较小时,流动为稳定轴对称流动,系统总熵产较小;随着温度梯度的增加,流动将失去其稳定性,首先转化为径向脉动波,此时系统总熵产呈周期性变化;温度梯度再增加时,流动转化为热流体波,系统总熵产较大,但不再随时间变化。     

统一的对流扩散型可压缩流体力学方程与解法

流体力学的动量方程、能量方程、湍动能方程和耗散方程都具有对流扩散方程的形式,但连续方程却不是对流扩散型的。对于可压缩问题,本文通过合理的数学推导,不作任何近似、假定与简化,得到一个全新的连续方程形式．该连续方程以压力为未知变量,并具有对流扩散型形式,使得所有的流体动力学方程组都具有完全统一的方程形式,给出了这种三维对流扩散方程组的有限精确差分计算格式。对流体力学的进一步发展具有一定意义．     

热声波数值模拟的虚假振荡研究

采用可压缩流动的SIMPLE算法对一维封闭空腔内由边界突然加热所引起的非稳态热声波进行了数值模拟,对流-扩散项采用了中心差分、一阶迎风差分、QUICK、及MIJSCL等不同格式。计算表明各种格式均存在不同程度的虚假振荡现象,其大小与热声波的强度及离散格式的形式等多种因素有关。这些结果对热声波的进一步研究及高效可靠的对流差分格式的开发具有重要意义。     

Three-dimensional mixed convection flow with variable thermal conductivity and frictional heating

In this article,three-dimensional mixed convection flow over an exponentially stretching sheet is investigated.Energy equation is modelled in the presence of viscous dissipation and variable thermal conductivity.Temperature of the sheet is varying exponentially and is chosen in a form that facilitates the similarity transformations to obtain self-similar equations.Resulting nonlinear ordinary differential equations are solved numerically employing the Runge-Kutta shooting method.In order to check the accuracy of the method,these equations are also solved using bvp4c built-in routine in Matlab.Both solutions are in excellent agreement.The effects of physical parameters on the dimensionless velocity field and temperature are demonstrated through various graphs.The novelty of this analysis is the self-similar solution of the threedimensional boundary layer flow in the presence of mixed convection,viscous dissipation and variable thermal conductivity.     

Optical Trapping with Focused Surface Waves

Near-field optical trapping can be realized with focused evanescent waves that are excited at the water–glass interface due to the total internal reflection, or with focused plasmonic waves excited on the water–gold interface. Herein, the performance of these two kinds of near-field optical trapping techniques is compared using the same optical microscope configuration. Experimental results show that only a single-micron polystyrene bead can be trapped by the focused evanescent waves, whereas many beads are simultaneously attracted to the center of the excited region by focused plasmonic waves. This difference in trapping behavior is analyzed from the electric field intensity distributions of these two kinds of focused surface waves and the difference in trapping behavior is attributed to photothermal effects due to the light absorption by the gold film.