A method that series perturbations approximate solutions to N-S equations with boundary conditions was discussed and adopted. Then the method was proved in which the asymptotic solutions of viscous fluid flow past a s...A method that series perturbations approximate solutions to N-S equations with boundary conditions was discussed and adopted. Then the method was proved in which the asymptotic solutions of viscous fluid flow past a sphere were deducted. By the ameliorative asymptotic expansion matched method, the matched functions, are determined easily and the ameliorative curve of drag coefficient is coincident well with measured data in the case that Reynolds number is less than or equal to 40 000.展开更多
In this paper we consider numerical simulation of incompressible viscous flow in a infinite slip channel. Local artificial boundary conditions at an artificial boundary are derived by the continuity of velocity and no...In this paper we consider numerical simulation of incompressible viscous flow in a infinite slip channel. Local artificial boundary conditions at an artificial boundary are derived by the continuity of velocity and normal stress at the segment artificial boundary. The the original problem is reduced to a boundary value problem on a bounded computational domain. Numerical example shows that our artificial boundary conditions are very effective.展开更多
For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a s...For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a structure-preserving method based on a discrete variational derivative method.Furthermore,to prevent the vertex concentration that may lead to numerical instability,we discretely introduce Deckelnick’s tangential velocity.Here,a modification term is introduced in the process of adding tangential velocity.This modified term enables the method to reproduce the equations’properties while preventing vertex concentration.Numerical experiments demonstrate that the proposed approach captures the equations’properties with high accuracy and avoids the concentration of vertices.展开更多
The steady bifurcation flows in a spherical gap (gap ratio =0.18) with rotating inner and stationary outer spheres are simulated numerically for Reci<Re<1500 1500 by solving steady axisymmetric incompressible Na...The steady bifurcation flows in a spherical gap (gap ratio =0.18) with rotating inner and stationary outer spheres are simulated numerically for Reci<Re<1500 1500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775<Re<1220 and three steady stable flows with 0, 1, or 2 vortices for 1 220<Re<l500. The formation of different flows at the same Reynolds number is related with different initial conditions which can be generated by different accelerations of the inner sphere. Generation of asro-or two-vortex flow depends mainly on the acceleration, but that of one-vortex flow also depends on the perturbation breaking the equatorial symmetry. The mechanism of development of a saddle point in the meridional plane at higher Re number and its role in the formation of two-vortex flow are analy2Ed.展开更多
文摘A method that series perturbations approximate solutions to N-S equations with boundary conditions was discussed and adopted. Then the method was proved in which the asymptotic solutions of viscous fluid flow past a sphere were deducted. By the ameliorative asymptotic expansion matched method, the matched functions, are determined easily and the ameliorative curve of drag coefficient is coincident well with measured data in the case that Reynolds number is less than or equal to 40 000.
文摘In this paper we consider numerical simulation of incompressible viscous flow in a infinite slip channel. Local artificial boundary conditions at an artificial boundary are derived by the continuity of velocity and normal stress at the segment artificial boundary. The the original problem is reduced to a boundary value problem on a bounded computational domain. Numerical example shows that our artificial boundary conditions are very effective.
基金This work was supported by JSPS KAKENHI Grant Nos.19K14590,21K18301,Japan.
文摘For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a structure-preserving method based on a discrete variational derivative method.Furthermore,to prevent the vertex concentration that may lead to numerical instability,we discretely introduce Deckelnick’s tangential velocity.Here,a modification term is introduced in the process of adding tangential velocity.This modified term enables the method to reproduce the equations’properties while preventing vertex concentration.Numerical experiments demonstrate that the proposed approach captures the equations’properties with high accuracy and avoids the concentration of vertices.
基金Project supported by the National Natural Science Foundation of China.
文摘The steady bifurcation flows in a spherical gap (gap ratio =0.18) with rotating inner and stationary outer spheres are simulated numerically for Reci<Re<1500 1500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775<Re<1220 and three steady stable flows with 0, 1, or 2 vortices for 1 220<Re<l500. The formation of different flows at the same Reynolds number is related with different initial conditions which can be generated by different accelerations of the inner sphere. Generation of asro-or two-vortex flow depends mainly on the acceleration, but that of one-vortex flow also depends on the perturbation breaking the equatorial symmetry. The mechanism of development of a saddle point in the meridional plane at higher Re number and its role in the formation of two-vortex flow are analy2Ed.