The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed o...The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed on improving the accuracy and efficiency of numerical simulation of similarproblems of fluid flow.In particular,the importance of suitable specification of the distortion func-tion to enforcing the far field boundarv conditions is emphasized.展开更多
A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced ...A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain.展开更多
Two strategies extended the single-cascade methods from a compressible three-dimensional inverse method for radial and mixed flow turbomachines to two three-dimensional multi-cascade co-design methods for single-stage...Two strategies extended the single-cascade methods from a compressible three-dimensional inverse method for radial and mixed flow turbomachines to two three-dimensional multi-cascade co-design methods for single-stage centrifugal compressors.These two three-dimensional methods and a typical quasi-threedimensional streamline curvature through-flow inverse method were employed to design the same subsonic high-speed single-stage centrifugal compressors.The compressor performances were simulated by a commercial Reynolds averaged Navier-Stokes(RANS) equations solver.The studies show that two three-dimensional codesign methods are reasonable and feasible.It was found that : firstly the blade camber angle designed by the three-dimensional methods was larger than that designed by the quasi-three-dimensional method;and secondly with regard to two three-dimensional methods with different boundary conditions,the co-design result differences between the diffusers were small,but those between the deswirlers were relatively large.展开更多
Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fre...Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.展开更多
A magnetohydrodynamic (MHD) formulation is derived to investigate and compare the mitigation effects of both the sheared axial flow and finite Larmor radius (FLR) on the Rayleigh-Taylor (RT) instability in Z-pinch imp...A magnetohydrodynamic (MHD) formulation is derived to investigate and compare the mitigation effects of both the sheared axial flow and finite Larmor radius (FLR) on the Rayleigh-Taylor (RT) instability in Z-pinch implosions. The sheared axial flow is introduced into MHD equations in a conventional way and the FLR effect into the equations via /t → -i(w+ik⊥2pi2Ωi,), as proposed in our previous paper [Chin. Phys. Lett. 2002, 19:217] , where k⊥2 pi2 is referred to FLR effect from the general kinetic theory of magnetized plasma. Therefore the linearized continuity and momentum equations for the perturbed mass-density and velocity include both the sheared axial flow and the FLR effect. It is found that the effect of sheared axial flow with a lower peak velocity can mitigate RT instability in the whole wavenumber region and the effect of sheared axial flow with a higher one can mitigate RT instability only in the large wavenumber region (for normalized wavenumber k】2.4); The effect of FLR can mitigate RT instability in the whole wavenumber region and the mitigation effect is stronger than that of the sheared axial flow with a lower peak velocity in the almost whole wavenumber region.展开更多
The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fit...The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.展开更多
This paper presents a generalized framework of stochastic modeling for particle kinetics in wall-bounded flow.We modified a reflected Brownian motion process and straightforwardly obtained a Kramers equation for parti...This paper presents a generalized framework of stochastic modeling for particle kinetics in wall-bounded flow.We modified a reflected Brownian motion process and straightforwardly obtained a Kramers equation for particle probability density function(PDF).After the wall effects were accounted for as a drift from zero in the mean displacement and suppression in the diffusivity of a particle,an analytical solution was worked out for PDF.Three distinguishable mechanisms were identified to affect the profile of particle probability distribution:external forces,turbophoresis effect,and wall-drift effect.The proposed formulation covers the Huang et al.(2009)model of a wall that produces electrostatic repulsion force and van der Waals force,as well as Monte-Carlo solutions for the Peter and Barenbrug(2002)model under a variety of relaxation times.Moreover,it successfully reproduces the two patterns of particle concentration profiles observed in experiments of sediment-laden open-channel flows.The strength of the wall-drift effect was found to be connected with the interaction frequency between particle and wall.Further exploration of the relationship among flow turbulence,particle inertia,and particle concentration is worthwhile.展开更多
基金Supported by the National Natural Science Foundation of China.
文摘The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed on improving the accuracy and efficiency of numerical simulation of similarproblems of fluid flow.In particular,the importance of suitable specification of the distortion func-tion to enforcing the far field boundarv conditions is emphasized.
文摘A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain.
基金Programme of Introducing Talents of Discipline to Universities(B08009)
文摘Two strategies extended the single-cascade methods from a compressible three-dimensional inverse method for radial and mixed flow turbomachines to two three-dimensional multi-cascade co-design methods for single-stage centrifugal compressors.These two three-dimensional methods and a typical quasi-threedimensional streamline curvature through-flow inverse method were employed to design the same subsonic high-speed single-stage centrifugal compressors.The compressor performances were simulated by a commercial Reynolds averaged Navier-Stokes(RANS) equations solver.The studies show that two three-dimensional codesign methods are reasonable and feasible.It was found that : firstly the blade camber angle designed by the three-dimensional methods was larger than that designed by the quasi-three-dimensional method;and secondly with regard to two three-dimensional methods with different boundary conditions,the co-design result differences between the diffusers were small,but those between the deswirlers were relatively large.
文摘Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.
基金This work was supported by the National Natural Science Foundation of China No.10035020.
文摘A magnetohydrodynamic (MHD) formulation is derived to investigate and compare the mitigation effects of both the sheared axial flow and finite Larmor radius (FLR) on the Rayleigh-Taylor (RT) instability in Z-pinch implosions. The sheared axial flow is introduced into MHD equations in a conventional way and the FLR effect into the equations via /t → -i(w+ik⊥2pi2Ωi,), as proposed in our previous paper [Chin. Phys. Lett. 2002, 19:217] , where k⊥2 pi2 is referred to FLR effect from the general kinetic theory of magnetized plasma. Therefore the linearized continuity and momentum equations for the perturbed mass-density and velocity include both the sheared axial flow and the FLR effect. It is found that the effect of sheared axial flow with a lower peak velocity can mitigate RT instability in the whole wavenumber region and the effect of sheared axial flow with a higher one can mitigate RT instability only in the large wavenumber region (for normalized wavenumber k】2.4); The effect of FLR can mitigate RT instability in the whole wavenumber region and the mitigation effect is stronger than that of the sheared axial flow with a lower peak velocity in the almost whole wavenumber region.
文摘The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.
基金supported by the National Natural Science Foundation of China(Grant Nos.51379100 and 51039003)
文摘This paper presents a generalized framework of stochastic modeling for particle kinetics in wall-bounded flow.We modified a reflected Brownian motion process and straightforwardly obtained a Kramers equation for particle probability density function(PDF).After the wall effects were accounted for as a drift from zero in the mean displacement and suppression in the diffusivity of a particle,an analytical solution was worked out for PDF.Three distinguishable mechanisms were identified to affect the profile of particle probability distribution:external forces,turbophoresis effect,and wall-drift effect.The proposed formulation covers the Huang et al.(2009)model of a wall that produces electrostatic repulsion force and van der Waals force,as well as Monte-Carlo solutions for the Peter and Barenbrug(2002)model under a variety of relaxation times.Moreover,it successfully reproduces the two patterns of particle concentration profiles observed in experiments of sediment-laden open-channel flows.The strength of the wall-drift effect was found to be connected with the interaction frequency between particle and wall.Further exploration of the relationship among flow turbulence,particle inertia,and particle concentration is worthwhile.
