In this paper,the relaxation algorithm and two Uzawa type algorithms for solving discretized variational inequalities arising from the two-phase Stefan type problem are proposed.An analysis of their convergence is pre...In this paper,the relaxation algorithm and two Uzawa type algorithms for solving discretized variational inequalities arising from the two-phase Stefan type problem are proposed.An analysis of their convergence is presented and the upper bounds of the convergence rates are derived.Some numerical experiments are shown to demonstrate that for the second Uzawa algorithm which is an improved version of the first Uzawa algorithm,the convergence rate is uniformly bounded away from 1 if τh^-2 is kept bounded,where τ is the time step size and h the space mesh size.展开更多
Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to so...Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to solve the gas-phase governing equ-ations and the characteristic method to follow the tracks of particles, we then obtainedthe full coupled numerical method of two-phase.transonic, turbulent flow. Here, par- ticle size may be grouped, the subsonic boundary condition at entry of nozzle is ireatedby quasi-characteristic method in reference plane and the algebraic model is used forturbulent flow. These methods are applied in viscous two-phase flow. calculation of ro-cket nozzle and in the prediciton of thrust and specific impulse for solid propellant ro-cket motor. The calculation results are in good agreement with the measurerment va-lues. Moreover, the influences of different particle radius, different particle mass frac-tion and particle size grouped on flow field have been discussed, and the influences of particle two-dimensional radial velosity component and viscosity on specific impulse ofrocket motor have been analysed.The method of this paper possesses the advantage of saving computer time. More important, the effect is more obvious for the calculation of particle size being grouped.展开更多
In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows...In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows. As the application of linear elements in SUPG/PSPG schemes gives rise to inconsistency in stabilization terms due to the inability to regenerate the diffusive term from viscous stresses, the numerical accuracy would deteriorate dramatically. To address this issue, projections of convection and pressure gradient terms are constructed and incorporated into the stabilization formulation in our method. This would substantially recover the consistency and free the practitioner from burdensome computations of most items in the residual. Moreover, the XFEM is employed to consider in a convenient way the fluid properties that have interfacial jumps leading to discontinuities in the velocity and pressure fields as well as the projections. A number of numerical examples are analyzed to demonstrate the complete recovery of consistency, the reproduction of interfacial discontinuities and the ability of the proposed projection-based SUPG/PSPG XFEM to model two-phase flows with open and closed interfaces.展开更多
The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the d...The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the development of a fully implicit discretization method.The main advantage of the fully implicit method is its unconditional stability.Newton's scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations.However,the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers.In this paper,a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations.The drawbacks of Newton's method are avoided in the developed model.Different limiter functions are considered,and the stabilized method is developed under steady and transient conditions.The results obtained by the numerical modeling are in good agreement with the experimental data.As expected,the results prove that the developed model is not restricted by any stability limit.展开更多
基金supported by the National Natural Science Foundation (10871179) of China
文摘In this paper,the relaxation algorithm and two Uzawa type algorithms for solving discretized variational inequalities arising from the two-phase Stefan type problem are proposed.An analysis of their convergence is presented and the upper bounds of the convergence rates are derived.Some numerical experiments are shown to demonstrate that for the second Uzawa algorithm which is an improved version of the first Uzawa algorithm,the convergence rate is uniformly bounded away from 1 if τh^-2 is kept bounded,where τ is the time step size and h the space mesh size.
文摘Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to solve the gas-phase governing equ-ations and the characteristic method to follow the tracks of particles, we then obtainedthe full coupled numerical method of two-phase.transonic, turbulent flow. Here, par- ticle size may be grouped, the subsonic boundary condition at entry of nozzle is ireatedby quasi-characteristic method in reference plane and the algebraic model is used forturbulent flow. These methods are applied in viscous two-phase flow. calculation of ro-cket nozzle and in the prediciton of thrust and specific impulse for solid propellant ro-cket motor. The calculation results are in good agreement with the measurerment va-lues. Moreover, the influences of different particle radius, different particle mass frac-tion and particle size grouped on flow field have been discussed, and the influences of particle two-dimensional radial velosity component and viscosity on specific impulse ofrocket motor have been analysed.The method of this paper possesses the advantage of saving computer time. More important, the effect is more obvious for the calculation of particle size being grouped.
文摘In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows. As the application of linear elements in SUPG/PSPG schemes gives rise to inconsistency in stabilization terms due to the inability to regenerate the diffusive term from viscous stresses, the numerical accuracy would deteriorate dramatically. To address this issue, projections of convection and pressure gradient terms are constructed and incorporated into the stabilization formulation in our method. This would substantially recover the consistency and free the practitioner from burdensome computations of most items in the residual. Moreover, the XFEM is employed to consider in a convenient way the fluid properties that have interfacial jumps leading to discontinuities in the velocity and pressure fields as well as the projections. A number of numerical examples are analyzed to demonstrate the complete recovery of consistency, the reproduction of interfacial discontinuities and the ability of the proposed projection-based SUPG/PSPG XFEM to model two-phase flows with open and closed interfaces.
文摘The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the development of a fully implicit discretization method.The main advantage of the fully implicit method is its unconditional stability.Newton's scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations.However,the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers.In this paper,a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations.The drawbacks of Newton's method are avoided in the developed model.Different limiter functions are considered,and the stabilized method is developed under steady and transient conditions.The results obtained by the numerical modeling are in good agreement with the experimental data.As expected,the results prove that the developed model is not restricted by any stability limit.
