A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
In this study,a large eddy simulation(LES)for fully-developed turbulent flows through a duct of regular-polygon cross-section using the immersed boundary(IB)method is performed.In case of the turbulent flow through th...In this study,a large eddy simulation(LES)for fully-developed turbulent flows through a duct of regular-polygon cross-section using the immersed boundary(IB)method is performed.In case of the turbulent flow through the square duct,though there are some disagreements of the mean quantities related with the streamwise velocity among the present LES,the previous direct numerical simulation(DNS)and the LES without the IB method,and the present LES can reproduce the secondary flow of the DNS and LES.The LES result for ten types of regular-polygon duct shows that the secondary-flow speed decreases as the number of sides of the regular polygon n increases and that the secondary flow in case of the regular icosagon duct disappears like the turbulent pipe flow.In case of low n,the behavior of the turbulent structures near the side center is different from that near the vertex.展开更多
Vortex-induced vibration(VIV)of cylindrical structures is a critical fluid-structure interaction(FSI)phenomenon in ocean engineering.Simulating VIV accurately can be computationally expensive.This study presents a gra...Vortex-induced vibration(VIV)of cylindrical structures is a critical fluid-structure interaction(FSI)phenomenon in ocean engineering.Simulating VIV accurately can be computationally expensive.This study presents a graphics processing unit(GPU)-accelerated simulation model for VIV utilizing the immersed boundary lattice Boltzmann method(IB-LBM),aiming to reduce computational costs while preserving accuracy.The program is developed using machine learning library JAX,which enables parallelism on GPU and multi-GPU platforms.The model incorporates multi-GPU parallelization and multi-block grid refinement strategies to enhance computational efficiency.Validation against existing high-fidelity simulation data demonstrates good agreement.Performance tests show significant speed-ups with GPU acceleration compared to traditional CPU-based approaches.These results underscore the potential of the developed simulator as an efficient and reliable tool for in-depth parametric studies and practical engineering analysis of VIV,facilitating more rapid design iterations and risk assessments for offshore structures.展开更多
Based on the immersed boundary method, a numerical simulation for an oscillating airfoil is established and a preliminary analysis of the oscillating airfoil is presented with an emphasis on the physical understanding...Based on the immersed boundary method, a numerical simulation for an oscillating airfoil is established and a preliminary analysis of the oscillating airfoil is presented with an emphasis on the physical understanding of fluid-structure interaction. In order to validate the method, two simulation cases： oscillating circular cylinder at low K-C number and two degrees of freedom oscillating cylinder are carded out first and the results are in good agreement with the previous re：searches. In the oscillating airfoil simulation, it is found that the reduced velocity U^＊. is a very sensitive factor and especially U^＊-2.8 is the critical stable boundary in the present work. The method shows the predominance of time saving in computational process for such a complicated fluid-structure interac- tion problem.展开更多
In this paper, a method that combines the characteristic-based split finite element method（CBS-FEM） and the direct forcing immersed boundary（IB） method is proposed for the simulation of incompressible viscous flow...In this paper, a method that combines the characteristic-based split finite element method（CBS-FEM） and the direct forcing immersed boundary（IB） method is proposed for the simulation of incompressible viscous flows. The structured triangular meshes without regarding the location of the physical boundary of the body is adopted to solve the flow, and the no-slip boundary condition is imposed on the interface. In order to improve the computational efficiency, a grid stretching strategy for the background structured triangular meshes is adopted. The obtained results agree very well with the previous numerical and experimental data. The order of the numerical accuracy is shown to be between 1 and 2. Moreover, the accuracy control by adjusting the number density of the mark points purely at certain stages is explored, and a second power law is obtained. The numerical experiments for the flow around a cylinder behind a backward-facing step show that the location of the cylinder can affect the sizes and the shapes of the corner eddy and the main recirculation region. The proposed method can be applied further to the fluid dynamics with complex geometries, moving boundaries, fluid-structure interactions, etc..展开更多
Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numeri...Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numerical method is applied to the simulation of the flow passing an oscillating circular cylinder which is forced to move in X direction under prescribed motions in water at rest at low Keulegan-Carpenter numbers.Then vor-tex-induced vibration of a cylinder with two degrees of freedom which oscillates in in-line direction and transverse direction is simulated using this method.The results are in good agreement with the previous research.Then the method is extended to the oscillating cascade simulation of making various comparisons.It is found that the IBPA(inter blade phase angle) will change as the time goes on,because of the non-uniformity of the flow in the circumferential direction,until the oscillating cascade goes to a stable situation.The reduced velocity and the number of blades are chosen to investigate the effects of them on IBPA.The results indicate that both the reduced velocity and the number of blades are the main factors which influence IBPA.It is worth noting that the coupling process is not necessary to generate any body-fitting grids,which makes it much faster in computational process for such a complicated fluid-structure interaction problem.展开更多
In this study,fish behavior and fish injury at different operating conditions are investigated via numerical simulation to evaluate the fish-friendliness of an axial pump that comprises an inlet pipe,a rotor with six ...In this study,fish behavior and fish injury at different operating conditions are investigated via numerical simulation to evaluate the fish-friendliness of an axial pump that comprises an inlet pipe,a rotor with six blades,a stator with eight vanes,and an outlet pipe.To precisely obtain the flow field when the fish passes through the axial pump,a hybrid large eddy simulation and immersed boundary method is adopted with the full consideration of the fluid-structure interaction comprehensively.The results indicate that the collision between the fish and the wall of flow components in the axial pump is concentrated near the inlet of the rotor,which results in the complexity of the fish trajectory,especially under the large flow rate condition.It is noted that the fish is likely to move in the reverse direction of the main flow after the impact with the rotor blade if the flow rate coefficient is too large,which increases the possibility of collision between the fish and the rotor blade.It is also indicated that the primary factor affecting the strike injury on the fish when it passes through the axial pump is the strike between the fish and the leading edge of the rotor blade.In addition,the strike injury becomes more significant as the flow rate coefficient increases.Furthermore,the results demonstrate that the fish may simultaneously suffer from strike,pressure,and shear stress injuries,once the collision between the fish and the wall of flow components occurs in the axial pump,thus aggravating the combined damage on the fish.Based on these results,it is recommended that hydraulic machinery should not be operated at large flow rates during fish migration from the view of fish-friendliness.展开更多
Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results s...Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results show that fishlike forked configurations have better locomotion performance compared with unforked plates. Based on our results, the caudal fin in carangi- form mode has greater thrust, and the lunate tail fin in thtmniform mode has higher efficiency. These findings are qualitatively con- sistent with biological observations of fish swimming. Analysis of wake topology shows that the wake of the forked plate consists of a chain of alternating reverse horseshoe-like vortical structures. These structures induce a backward jet and generate a positive thrust. Moreover, this backward jet has a more favorable direction compared with that behind an unforked plate.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
文摘In this study,a large eddy simulation(LES)for fully-developed turbulent flows through a duct of regular-polygon cross-section using the immersed boundary(IB)method is performed.In case of the turbulent flow through the square duct,though there are some disagreements of the mean quantities related with the streamwise velocity among the present LES,the previous direct numerical simulation(DNS)and the LES without the IB method,and the present LES can reproduce the secondary flow of the DNS and LES.The LES result for ten types of regular-polygon duct shows that the secondary-flow speed decreases as the number of sides of the regular polygon n increases and that the secondary flow in case of the regular icosagon duct disappears like the turbulent pipe flow.In case of low n,the behavior of the turbulent structures near the side center is different from that near the vertex.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFC2806300)the National Natural Science Foundation of China(Grant No.52401338)+1 种基金supported by the State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University(Grant No.GKZD010089)the Development and Application Project of Ship CAE Software,and the Seed Foundation of Tianjin University.
文摘Vortex-induced vibration(VIV)of cylindrical structures is a critical fluid-structure interaction(FSI)phenomenon in ocean engineering.Simulating VIV accurately can be computationally expensive.This study presents a graphics processing unit(GPU)-accelerated simulation model for VIV utilizing the immersed boundary lattice Boltzmann method(IB-LBM),aiming to reduce computational costs while preserving accuracy.The program is developed using machine learning library JAX,which enables parallelism on GPU and multi-GPU platforms.The model incorporates multi-GPU parallelization and multi-block grid refinement strategies to enhance computational efficiency.Validation against existing high-fidelity simulation data demonstrates good agreement.Performance tests show significant speed-ups with GPU acceleration compared to traditional CPU-based approaches.These results underscore the potential of the developed simulator as an efficient and reliable tool for in-depth parametric studies and practical engineering analysis of VIV,facilitating more rapid design iterations and risk assessments for offshore structures.
基金supported by NSFC under grants 50736007 and 50136010
文摘Based on the immersed boundary method, a numerical simulation for an oscillating airfoil is established and a preliminary analysis of the oscillating airfoil is presented with an emphasis on the physical understanding of fluid-structure interaction. In order to validate the method, two simulation cases： oscillating circular cylinder at low K-C number and two degrees of freedom oscillating cylinder are carded out first and the results are in good agreement with the previous re：searches. In the oscillating airfoil simulation, it is found that the reduced velocity U^＊. is a very sensitive factor and especially U^＊-2.8 is the critical stable boundary in the present work. The method shows the predominance of time saving in computational process for such a complicated fluid-structure interac- tion problem.
基金Project supported by the National High Technology Re-search and Development Program of China(863 Program,Grant No.2012AA011803)the National Natural Scientific Foundation of China(Grant No.11172241)the University Foundation for Fundamental Research of NPU(Grant No.JCY-20130121)
文摘In this paper, a method that combines the characteristic-based split finite element method（CBS-FEM） and the direct forcing immersed boundary（IB） method is proposed for the simulation of incompressible viscous flows. The structured triangular meshes without regarding the location of the physical boundary of the body is adopted to solve the flow, and the no-slip boundary condition is imposed on the interface. In order to improve the computational efficiency, a grid stretching strategy for the background structured triangular meshes is adopted. The obtained results agree very well with the previous numerical and experimental data. The order of the numerical accuracy is shown to be between 1 and 2. Moreover, the accuracy control by adjusting the number density of the mark points purely at certain stages is explored, and a second power law is obtained. The numerical experiments for the flow around a cylinder behind a backward-facing step show that the location of the cylinder can affect the sizes and the shapes of the corner eddy and the main recirculation region. The proposed method can be applied further to the fluid dynamics with complex geometries, moving boundaries, fluid-structure interactions, etc..
文摘Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numerical method is applied to the simulation of the flow passing an oscillating circular cylinder which is forced to move in X direction under prescribed motions in water at rest at low Keulegan-Carpenter numbers.Then vor-tex-induced vibration of a cylinder with two degrees of freedom which oscillates in in-line direction and transverse direction is simulated using this method.The results are in good agreement with the previous research.Then the method is extended to the oscillating cascade simulation of making various comparisons.It is found that the IBPA(inter blade phase angle) will change as the time goes on,because of the non-uniformity of the flow in the circumferential direction,until the oscillating cascade goes to a stable situation.The reduced velocity and the number of blades are chosen to investigate the effects of them on IBPA.The results indicate that both the reduced velocity and the number of blades are the main factors which influence IBPA.It is worth noting that the coupling process is not necessary to generate any body-fitting grids,which makes it much faster in computational process for such a complicated fluid-structure interaction problem.
基金supported by the National Natural Science Foundation of China(Grant Nos.51776102 and 91852103)the Institute for Guo Qiang,Tsinghua University(Grant No.2019GQG1019)+1 种基金the Tsinghua National Laboratory for Information Science and Technologythe China Scholarship Council for sponsoring her visit to the University of Minnesota。
文摘In this study,fish behavior and fish injury at different operating conditions are investigated via numerical simulation to evaluate the fish-friendliness of an axial pump that comprises an inlet pipe,a rotor with six blades,a stator with eight vanes,and an outlet pipe.To precisely obtain the flow field when the fish passes through the axial pump,a hybrid large eddy simulation and immersed boundary method is adopted with the full consideration of the fluid-structure interaction comprehensively.The results indicate that the collision between the fish and the wall of flow components in the axial pump is concentrated near the inlet of the rotor,which results in the complexity of the fish trajectory,especially under the large flow rate condition.It is noted that the fish is likely to move in the reverse direction of the main flow after the impact with the rotor blade if the flow rate coefficient is too large,which increases the possibility of collision between the fish and the rotor blade.It is also indicated that the primary factor affecting the strike injury on the fish when it passes through the axial pump is the strike between the fish and the leading edge of the rotor blade.In addition,the strike injury becomes more significant as the flow rate coefficient increases.Furthermore,the results demonstrate that the fish may simultaneously suffer from strike,pressure,and shear stress injuries,once the collision between the fish and the wall of flow components occurs in the axial pump,thus aggravating the combined damage on the fish.Based on these results,it is recommended that hydraulic machinery should not be operated at large flow rates during fish migration from the view of fish-friendliness.
基金the National Natural Science Foundation of China (Grant No. 10832010)the Innovation Project of the Chinese Academy of Sciences (Grant No. KJCX2-YW-L05)the 111 Project (Grant No. B07033)
文摘Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results show that fishlike forked configurations have better locomotion performance compared with unforked plates. Based on our results, the caudal fin in carangi- form mode has greater thrust, and the lunate tail fin in thtmniform mode has higher efficiency. These findings are qualitatively con- sistent with biological observations of fish swimming. Analysis of wake topology shows that the wake of the forked plate consists of a chain of alternating reverse horseshoe-like vortical structures. These structures induce a backward jet and generate a positive thrust. Moreover, this backward jet has a more favorable direction compared with that behind an unforked plate.
