The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S）...The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.展开更多
To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced v...To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.展开更多
We numerically study the dynamics of particle crystals in annular microchannels by the immersed-boundary（IB）lattice Boltzmann（LB） coupled model, analyze the fluid-particle interactions during the migration of part...We numerically study the dynamics of particle crystals in annular microchannels by the immersed-boundary（IB）lattice Boltzmann（LB） coupled model, analyze the fluid-particle interactions during the migration of particles,and reveal the underlying mechanism of a particle focusing on the presence of fluid flows. The results show that the Reynolds and Dean numbers are key factors influencing the hydrodynamics of particles. The particles migrate onto their equilibrium tracks by adjusting the Reynolds and Dean numbers. Elliptical tracks of particles during hydrodynamic focusing can be predicted by the IB-LB model. Both the small Dean number and the small particle can lead to a small size of the focusing track. This work would possibly facilitate the utilization of annular microchannel flows to obtain microfluidic flowing crystals for advanced applications in biomedicine and materials synthesis.展开更多
In this work,the deformation of free interface during water entry and exit of a circular cylinder is investigated numerically by using the two-dimensional(2D)immersed boundary-multiphase lattice Boltzmann flux solver(...In this work,the deformation of free interface during water entry and exit of a circular cylinder is investigated numerically by using the two-dimensional(2D)immersed boundary-multiphase lattice Boltzmann flux solver(IB-MLBFS).The fluid domain is discretized by finite volume discretization,and the flux on the grid interface is evaluated by lattice Boltzmann equations.Both the implicit velocity correction and the surface flux correction are implemented by using the immersed boundary-method to consider the fluid-structure interaction and the contact interface between the multiphase fluids and the structure.First,the water entry of a circular cylinder is simulated and the results are compared with the experiment,which considered the length-diameter ratio of the circular cylinder.The reliability of 2D simulation is verified and the deformation of the free interface is well investigated.Afterward,the water exit of a circular cylinder with constant velocity is simulated,which is less researched.In addition,the results show the advantage of present IB-MLBFS to some extent.Finally,the water exit and re-entry of a circular cylinder are presented,and the results present the complex deformation of the free interface and the dynamic response of the moving structure.Based on the numerical results,the free interface of the multiphase fluids is well captured,and the contact interface on the boundary of the moving structure is accurately presented by the IB-MLBFS.展开更多
The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle...The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.展开更多
基金Supported by the Aeronautical Science Foundation of China(20111453012)the National Defense Pre-Research Foundation of China(9140A13040111HK0329)~~
文摘The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.
基金Supported by the National Natural Science Foundation of China(52201350,52201394,and 52271301)the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.
基金Supported by the National Natural Science Foundation of China under Grant Nos 51728601 and 51771118the Natural Science Foundation of Jiangsu Province under Grant No BK20150600
文摘We numerically study the dynamics of particle crystals in annular microchannels by the immersed-boundary（IB）lattice Boltzmann（LB） coupled model, analyze the fluid-particle interactions during the migration of particles,and reveal the underlying mechanism of a particle focusing on the presence of fluid flows. The results show that the Reynolds and Dean numbers are key factors influencing the hydrodynamics of particles. The particles migrate onto their equilibrium tracks by adjusting the Reynolds and Dean numbers. Elliptical tracks of particles during hydrodynamic focusing can be predicted by the IB-LB model. Both the small Dean number and the small particle can lead to a small size of the focusing track. This work would possibly facilitate the utilization of annular microchannel flows to obtain microfluidic flowing crystals for advanced applications in biomedicine and materials synthesis.
基金Supported by the National Natural Science Foundation of China(52061135107)the Fundamental Research Fund for the Central Universities(DUT20TD108,DUT20LAB308)the Liao Ning Revitalization Talents Program(XLYC1908027),and Dalian Innovation Research Team in Key Areas(2020RT03).
文摘In this work,the deformation of free interface during water entry and exit of a circular cylinder is investigated numerically by using the two-dimensional(2D)immersed boundary-multiphase lattice Boltzmann flux solver(IB-MLBFS).The fluid domain is discretized by finite volume discretization,and the flux on the grid interface is evaluated by lattice Boltzmann equations.Both the implicit velocity correction and the surface flux correction are implemented by using the immersed boundary-method to consider the fluid-structure interaction and the contact interface between the multiphase fluids and the structure.First,the water entry of a circular cylinder is simulated and the results are compared with the experiment,which considered the length-diameter ratio of the circular cylinder.The reliability of 2D simulation is verified and the deformation of the free interface is well investigated.Afterward,the water exit of a circular cylinder with constant velocity is simulated,which is less researched.In addition,the results show the advantage of present IB-MLBFS to some extent.Finally,the water exit and re-entry of a circular cylinder are presented,and the results present the complex deformation of the free interface and the dynamic response of the moving structure.Based on the numerical results,the free interface of the multiphase fluids is well captured,and the contact interface on the boundary of the moving structure is accurately presented by the IB-MLBFS.
文摘The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.
