The simulation of three-dimensional (3D) non-isothermal, non-Newtonian fluid filling process is an extremely difficult task and remains a challenging problem, which includes polymer melt flow with free surface coupl...The simulation of three-dimensional (3D) non-isothermal, non-Newtonian fluid filling process is an extremely difficult task and remains a challenging problem, which includes polymer melt flow with free surface coupled with transient heat transfer. This paper presents a full 3D non-isothermal two-phase flow model to predict the complex flow in melt filling process, where the Cross-WLF model is applied to characterize the rheological behav- ior of polymer melt. The governing equations are solved using finite volume method with SIMPLEC algorithm on collocated grids and the melt front is accurately captured by a high resolution level set method. A domain exten- sion technique is adopted to deal with the complex cavities, which greatly reduces the computational burden. To verify the validity of the developed 3D approach, the melts filling processes in two thin rectangular cavities (one of them with a cylindrical insert) are simulated. The predicted melt front interfaces are in good agreement with the experiment and commercial software prediction. For a case with a rather complex cavity, the dynamic filling process in a hemispherical shell is successfully simulated. All of the numerical results show that the developed numerical procedure can provide a reasonable orediction for injection molding process.展开更多
Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were de...Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4.展开更多
The National Laboratory of Multiphase Flow in Power Engineering Xi' an Jiaotong University, Xi' an 710049, China Abstract:Calculations of several classic questions in numerical heat transfer with delaying modi...The National Laboratory of Multiphase Flow in Power Engineering Xi' an Jiaotong University, Xi' an 710049, China Abstract:Calculations of several classic questions in numerical heat transfer with delaying modified QUICK scheme show good agreement with benchrnark. Combining volume of fluid model with the well known SIMPLEC method. the reasonanble Tgylor bubble shape has ben obtained by using delaying modified QUICK scheme and donor acceptor scheme of interfacial advection respectively. At the same the the stream field and the liquid film thickness around the Taylor bubble have been well figured out展开更多
基金Supported by the National Basic Research Program of China(2012CB025903)the National Natural Science Foundation of China(91434201,11402210)
文摘The simulation of three-dimensional (3D) non-isothermal, non-Newtonian fluid filling process is an extremely difficult task and remains a challenging problem, which includes polymer melt flow with free surface coupled with transient heat transfer. This paper presents a full 3D non-isothermal two-phase flow model to predict the complex flow in melt filling process, where the Cross-WLF model is applied to characterize the rheological behav- ior of polymer melt. The governing equations are solved using finite volume method with SIMPLEC algorithm on collocated grids and the melt front is accurately captured by a high resolution level set method. A domain exten- sion technique is adopted to deal with the complex cavities, which greatly reduces the computational burden. To verify the validity of the developed 3D approach, the melts filling processes in two thin rectangular cavities (one of them with a cylindrical insert) are simulated. The predicted melt front interfaces are in good agreement with the experiment and commercial software prediction. For a case with a rather complex cavity, the dynamic filling process in a hemispherical shell is successfully simulated. All of the numerical results show that the developed numerical procedure can provide a reasonable orediction for injection molding process.
文摘Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4.
文摘The National Laboratory of Multiphase Flow in Power Engineering Xi' an Jiaotong University, Xi' an 710049, China Abstract:Calculations of several classic questions in numerical heat transfer with delaying modified QUICK scheme show good agreement with benchrnark. Combining volume of fluid model with the well known SIMPLEC method. the reasonanble Tgylor bubble shape has ben obtained by using delaying modified QUICK scheme and donor acceptor scheme of interfacial advection respectively. At the same the the stream field and the liquid film thickness around the Taylor bubble have been well figured out