The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/ficti...The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/fictitious domain method used in our DNS is a promising new way for simulation of particulate flows. The settling behaviors of the particle are presented ranging in Reynolds number from 0 to about 700, which showed that our results for low Reynolds numbers agreed well with that reported before. Nevertheless, for higher Reynolds numbers our results were different from theirs. The long-term mean equilibrium positions in our results were all on the centerline, but not at off-center position as reported before. In order to validate our simulation, experiments were also conducted. The results showed that the sedimenting behavior simulated in this paper agreed well with our experiment result.展开更多
The direct-forcing fictitious domain method is extended to simulate the locomotion of a passively pitching foil. Our study focuses on the hysteresis phenomenon that the critical frequency for the reverse of the locomo...The direct-forcing fictitious domain method is extended to simulate the locomotion of a passively pitching foil. Our study focuses on the hysteresis phenomenon that the critical frequency for the reverse of the locomotion direction of the wing in case of decreasing frequency is smaller than that in case of increasing frequency. In our simulations, the hysteresis phenomenon is produced by imposing different initial conditions at a same frequency. Our results indicate that the ratio of the heaving amplitudes of two foil edges is crucial to the direction of the foil＇s horizontal motion, and the amplitude of the leading edge is generally smaller. The critical frequencies for the reverse of the locomotion direction are increased, when the foil-fluid density ratio is decreased or the spring constant is increased. The critical frequencies in the bi-stability regime also depend on the initial velocity imposed, and the hysteresis loop generally becomes larger if the initial velocities are closer to the terminal locomotion velocities of the foil.展开更多
The hemodynamic interactions of red blood cells （RBCs） in a microcapillary flow are investigated in this paper. This kind of interaction is considered as the non-contact mutual interaction of cells, which is importa...The hemodynamic interactions of red blood cells （RBCs） in a microcapillary flow are investigated in this paper. This kind of interaction is considered as the non-contact mutual interaction of cells, which is important in the suspension flow of blood, but not sufficiently understood. The distributed Lagrange multiplier/fictitious domain method in the lattice Boltzmann framework is used to solve the suspension of the RBCs. The modification of the flow due to the cells, the dependence of the cell deformation on the flow and the cell-cell interaction via the fluid are discussed. It is revealed that the long-range hydrodynamic interaction with a long interacting distance, more than about 5 times of the RBC equivalent radius, mainly has effect on the rheology properties of the suspension, such as the mean velocity, and the short-range interaction is sensitive to the shape of the cell in the microcapillary flow. The flow velocity around the cell plays a key role in the cell deformation. In the current configuration of the flow and cells, the cells repel each other along the capillary.展开更多
In this paper, we propose a lattice Boltzmann (LB) method coupled with adirect-forcing fictitious domain (DF/FD) method for the simulation of particle suspensions. This method combines the good features of the LB and ...In this paper, we propose a lattice Boltzmann (LB) method coupled with adirect-forcing fictitious domain (DF/FD) method for the simulation of particle suspensions. This method combines the good features of the LB and the DF/FD methodsby using two unrelated meshes, namely, an Eulerian mesh for the flow domain and aLagrangian mesh for the solid domain, which avoids the re-meshing procedure anddoes not need to calculate the hydrodynamic forces at each time step. The non-slipboundary condition is enforced by introducing a forcing term into the lattice Boltzmann equation, which preserves all remarkable advantages of the LBM in simulatingfluid flows. The present LB-DF/FD method has been validated by comparing its results with analytical results and previous numerical results for a single circular particleand two circular particles settling under gravity. The interaction between particle andwall, the process of drafting-kissing-tumbling (DKT) of two settling particles will bedemonstrated. As a manifestation of the efficiency of the present method, the settlingof a large number (128) of circular particles is simulated in an enclosure.展开更多
A parallel direct-forcing fictitious domain method is applied in fully-resolved numerical simulations of particle-laden turbulent flows in a square duct. The effects of finite-size heavy particles on the mean secondar...A parallel direct-forcing fictitious domain method is applied in fully-resolved numerical simulations of particle-laden turbulent flows in a square duct. The effects of finite-size heavy particles on the mean secondary flow, the mean streamwise velocity, the root-meansquare velocity fluctuation, and the particle concentration distribution are investigated at the friction Reynolds number of 150, the particle volume fraction of 2.36%, the particle diameter of 0.1 duct width, and the Shields number ranging from 1.0 to 0.2 Our results show that the particle sedimentation breaks the up-down symmetry of the mean secondary vortices, and results in a stronger secondary-flow circulation which transports the fluids downward in the bulk center region and upward along the side walls at a low Shields number. This circulation has a significant impact on the distribution of the mean streamwise velocity, whose maximum value occurs in the lower half duct, unlike in the plane channel case. The flow resistance is increased and the turbulence intensity is reduced, as the Shields number is decreased. The particles accumulate preferentially at the face center of the bottom wall, due to the effect of the mean secondary flow. It is observed that the collision model has an important effect on the results, but does not change the results qualitatively.展开更多
A spring model is used to simulate the skeleton structure of the red blood cell (RBC) membrane and to study the red blood cell (RBC) rheology in Poiseuille flow with an immersed boundary method. The lateral migration ...A spring model is used to simulate the skeleton structure of the red blood cell (RBC) membrane and to study the red blood cell (RBC) rheology in Poiseuille flow with an immersed boundary method. The lateral migration properties of many cells in Poiseuille flow have been investigated. The authors also combine the above methodology with a distributed Lagrange multiplier/fictitious domain method to simulate the interaction of cells and neutrally buoyant particles in a microchannel for studying the margination of particles.展开更多
A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immerse...A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.展开更多
In this paper, an absorbing Fictitious Boundary Condition (FBC) is presented to generate an iterative Domain Decomposition Method (DDM) for analyzing waveguide problems.The relaxed algorithm is introduced to improve t...In this paper, an absorbing Fictitious Boundary Condition (FBC) is presented to generate an iterative Domain Decomposition Method (DDM) for analyzing waveguide problems.The relaxed algorithm is introduced to improve the iterative convergence. And the matrix equations are solved using the multifrontal algorithm. The resulting CPU time is greatly reduced.Finally, a number of numerical examples are given to illustrate its accuracy and efficiency.展开更多
文摘The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/fictitious domain method used in our DNS is a promising new way for simulation of particulate flows. The settling behaviors of the particle are presented ranging in Reynolds number from 0 to about 700, which showed that our results for low Reynolds numbers agreed well with that reported before. Nevertheless, for higher Reynolds numbers our results were different from theirs. The long-term mean equilibrium positions in our results were all on the centerline, but not at off-center position as reported before. In order to validate our simulation, experiments were also conducted. The results showed that the sedimenting behavior simulated in this paper agreed well with our experiment result.
基金Project supported by the National Natural Science Foun-dation of China(Grant No.11372275)the Program for New Century Excellent Talents in University
文摘The direct-forcing fictitious domain method is extended to simulate the locomotion of a passively pitching foil. Our study focuses on the hysteresis phenomenon that the critical frequency for the reverse of the locomotion direction of the wing in case of decreasing frequency is smaller than that in case of increasing frequency. In our simulations, the hysteresis phenomenon is produced by imposing different initial conditions at a same frequency. Our results indicate that the ratio of the heaving amplitudes of two foil edges is crucial to the direction of the foil＇s horizontal motion, and the amplitude of the leading edge is generally smaller. The critical frequencies for the reverse of the locomotion direction are increased, when the foil-fluid density ratio is decreased or the spring constant is increased. The critical frequencies in the bi-stability regime also depend on the initial velocity imposed, and the hysteresis loop generally becomes larger if the initial velocities are closer to the terminal locomotion velocities of the foil.
基金supported by the National Natural Science Foundation of China(Grant Nos.11372278,10902098)the Fundamental Research Funds of the Central Universities(Grant No.2010QNA40107)
文摘The hemodynamic interactions of red blood cells （RBCs） in a microcapillary flow are investigated in this paper. This kind of interaction is considered as the non-contact mutual interaction of cells, which is important in the suspension flow of blood, but not sufficiently understood. The distributed Lagrange multiplier/fictitious domain method in the lattice Boltzmann framework is used to solve the suspension of the RBCs. The modification of the flow due to the cells, the dependence of the cell deformation on the flow and the cell-cell interaction via the fluid are discussed. It is revealed that the long-range hydrodynamic interaction with a long interacting distance, more than about 5 times of the RBC equivalent radius, mainly has effect on the rheology properties of the suspension, such as the mean velocity, and the short-range interaction is sensitive to the shape of the cell in the microcapillary flow. The flow velocity around the cell plays a key role in the cell deformation. In the current configuration of the flow and cells, the cells repel each other along the capillary.
基金This work was supported by the Scientific Project of Zhejiang Province of China(No.2008C01024-4).
文摘In this paper, we propose a lattice Boltzmann (LB) method coupled with adirect-forcing fictitious domain (DF/FD) method for the simulation of particle suspensions. This method combines the good features of the LB and the DF/FD methodsby using two unrelated meshes, namely, an Eulerian mesh for the flow domain and aLagrangian mesh for the solid domain, which avoids the re-meshing procedure anddoes not need to calculate the hydrodynamic forces at each time step. The non-slipboundary condition is enforced by introducing a forcing term into the lattice Boltzmann equation, which preserves all remarkable advantages of the LBM in simulatingfluid flows. The present LB-DF/FD method has been validated by comparing its results with analytical results and previous numerical results for a single circular particleand two circular particles settling under gravity. The interaction between particle andwall, the process of drafting-kissing-tumbling (DKT) of two settling particles will bedemonstrated. As a manifestation of the efficiency of the present method, the settlingof a large number (128) of circular particles is simulated in an enclosure.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.11372275,51376162)the Resea-rch Fund for the Doctoral Program of Higher Education of China(Grant No.20130101110035)
文摘A parallel direct-forcing fictitious domain method is applied in fully-resolved numerical simulations of particle-laden turbulent flows in a square duct. The effects of finite-size heavy particles on the mean secondary flow, the mean streamwise velocity, the root-meansquare velocity fluctuation, and the particle concentration distribution are investigated at the friction Reynolds number of 150, the particle volume fraction of 2.36%, the particle diameter of 0.1 duct width, and the Shields number ranging from 1.0 to 0.2 Our results show that the particle sedimentation breaks the up-down symmetry of the mean secondary vortices, and results in a stronger secondary-flow circulation which transports the fluids downward in the bulk center region and upward along the side walls at a low Shields number. This circulation has a significant impact on the distribution of the mean streamwise velocity, whose maximum value occurs in the lower half duct, unlike in the plane channel case. The flow resistance is increased and the turbulence intensity is reduced, as the Shields number is decreased. The particles accumulate preferentially at the face center of the bottom wall, due to the effect of the mean secondary flow. It is observed that the collision model has an important effect on the results, but does not change the results qualitatively.
基金supported by the National Science Foundation of the United States (Nos. ECS-9527123, CTS-9873236, DMS-9973318, CCR-9902035, DMS-0209066, DMS-0443826, DMS-0914788)
文摘A spring model is used to simulate the skeleton structure of the red blood cell (RBC) membrane and to study the red blood cell (RBC) rheology in Poiseuille flow with an immersed boundary method. The lateral migration properties of many cells in Poiseuille flow have been investigated. The authors also combine the above methodology with a distributed Lagrange multiplier/fictitious domain method to simulate the interaction of cells and neutrally buoyant particles in a microchannel for studying the margination of particles.
文摘A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.
文摘In this paper, an absorbing Fictitious Boundary Condition (FBC) is presented to generate an iterative Domain Decomposition Method (DDM) for analyzing waveguide problems.The relaxed algorithm is introduced to improve the iterative convergence. And the matrix equations are solved using the multifrontal algorithm. The resulting CPU time is greatly reduced.Finally, a number of numerical examples are given to illustrate its accuracy and efficiency.
