In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along wi...In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil. The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.展开更多
Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （M...Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.展开更多
To predict the flow evolution of fish swimming problems,a flow solver based on the immersed boundary lattice Boltzmann method is developed.A flexible iterative algorithm based on the framework of implicit boundary for...To predict the flow evolution of fish swimming problems,a flow solver based on the immersed boundary lattice Boltzmann method is developed.A flexible iterative algorithm based on the framework of implicit boundary force correction is used to save the computational cost and memory,and the momentum forcing is described by a simple direct force formula without complicated integral calculation when the velocity correction at the boundary node is determined.With the presented flow solver,the hydrodynamic interaction between the fish-induced dynamic stall vortices and the incoming vortices in unsteady flow is analyzed.Numerical simulation results unveil the mechanism of fish exploiting vortices to enhance their own hydrodynamic performances.The superior swimming performances originate from the relative movement between the“merged vortex”and the locomotion of the fishtail,which is controlled by the phase difference.Formation conditions of the“merged vortex”become the key factor for fish to exploit vortices to improve their swimming performance.We further discuss the effect of the principal components of locomotion.From the results,we conclude that lateral translation plays a crucial role in propulsion while body undulation in tandem with rotation and head motion reduce the locomotor cost.展开更多
Dynamic fluid-solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method （IBM） is a convenient approac...Dynamic fluid-solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method （IBM） is a convenient approach to handle fluid-solid interactions with complex geometries. In this work, Uhlmann＇s direct-forcing IBM is improved and implemented on a supercomputer with CPU-GPU hybrid architec- ture. The direct-forcing IBM is modified as follows： the Poisson＇s equation for pressure is solved before evaluation of the body force, and the force is only distributed to the Cartesian grids inside the immersed boundary. A multidirect forcing scheme is used to evaluate the body force. These modifications result in a divergence-free flow field in the fluid domain and the no-slip boundary condition at the immersed boundary simultaneously. This method is implemented in an explicit finite-difference fractional-step scheme, and validated by 2D simulations of lid-driven cavity flow, Couette flow between two concentric cylinders and flow over a circular cylinder. Finally, the method is used to simulate the sedimentation of two circular particles in a channel. The results agree very well with previous experimental and numerical data, and are more accurate than the conventional direct-forcing method, especially in the vicinity of a moving boundary.展开更多
Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Latti...Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.展开更多
A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier...A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier-Stokes equations using the projection-based fractional step method in a fixed staggered Cartesian grid system.A volume of fluid(VOF)method with second-order accuracy is employed to trace the free surface.To represent the intricate surface geometry,the structure is discretized using the unstructured triangle mesh.Additionally,a ray tracing method is employed to classify fluid and solid points.A high-order stable scheme has been introduced to reconstruct the local velocity at interface points.Three FSI problems,including wave evolution around a breakwater,interaction between a periodic wave train and a moving float,and a 3-D moving object interacting with the free surface,were investigated to validate the accuracy and stability of the proposed model.The numerical results are in good agreement with the experimental data.Additionally,we evaluated the computational performance of the proposed GPU-based model.The GPU-based model achieved a 42.29 times speedup compared with the single-core CPU-based model in the three-dimension test.Additionally,the results regarding the time cost of each code section indicate that achieving more significant acceleration is associated with solving the turbulence,advection,and diffusion terms,while solving the pressure Poisson equation(PPE)saves the most time.Furthermore,the impact of grid number on computational efficiency indicates that as Fluid-solid interaction(FSI)immersed boundary method(IBM)graphical processing unit(GPU)two-phase flow moving rigid bodythe number of grids increases,the GPU-based model outperforms the multi-core CPU-based model.展开更多
In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as mad...In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as made of fluid and introduce two virtual forcing terms.First,a virtual fluid force distributed only on the solid region is appended to the momentum equation to make the region behave like a real solid body and satisfy the prescribed velocity.Second,a virtual heat source located inside the solid region near the boundary is added to the energy transport equation to impose the thermal boundary condition on the solid boundary.We take the implicit second-order backward differentiation to discretize the time variable and employ the Choi-Moin projection scheme to split the coupled system.As for spatial discretization,second-order centered differences over a staggered Cartesian grid are used on the entire domain.The advantages of this method are its conceptual simplicity and ease of implementation.Numerical experiments are performed to demonstrate the high performance of the proposed method.Convergence tests show that the spatial convergence rates of all unknowns seem to be super-linear in the 1-norm and 2-norm while less than linear in the maximum norm.展开更多
文摘In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil. The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.
基金Project(2012CB932202)supported by the National Basic Research Program of ChinaProjects(50890174,50971088)supported by the National Natural Science Foundation of China
文摘Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.
基金supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)。
文摘To predict the flow evolution of fish swimming problems,a flow solver based on the immersed boundary lattice Boltzmann method is developed.A flexible iterative algorithm based on the framework of implicit boundary force correction is used to save the computational cost and memory,and the momentum forcing is described by a simple direct force formula without complicated integral calculation when the velocity correction at the boundary node is determined.With the presented flow solver,the hydrodynamic interaction between the fish-induced dynamic stall vortices and the incoming vortices in unsteady flow is analyzed.Numerical simulation results unveil the mechanism of fish exploiting vortices to enhance their own hydrodynamic performances.The superior swimming performances originate from the relative movement between the“merged vortex”and the locomotion of the fishtail,which is controlled by the phase difference.Formation conditions of the“merged vortex”become the key factor for fish to exploit vortices to improve their swimming performance.We further discuss the effect of the principal components of locomotion.From the results,we conclude that lateral translation plays a crucial role in propulsion while body undulation in tandem with rotation and head motion reduce the locomotor cost.
基金supported by the National Natural Science Foundation of China(NSFC) under Grant Nos.21225628,51106168 and 11272312the "Strategic Priority Research Program" of Chinese Academy of Sciences(CAS) under Grant No.XDA07080102
文摘Dynamic fluid-solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method （IBM） is a convenient approach to handle fluid-solid interactions with complex geometries. In this work, Uhlmann＇s direct-forcing IBM is improved and implemented on a supercomputer with CPU-GPU hybrid architec- ture. The direct-forcing IBM is modified as follows： the Poisson＇s equation for pressure is solved before evaluation of the body force, and the force is only distributed to the Cartesian grids inside the immersed boundary. A multidirect forcing scheme is used to evaluate the body force. These modifications result in a divergence-free flow field in the fluid domain and the no-slip boundary condition at the immersed boundary simultaneously. This method is implemented in an explicit finite-difference fractional-step scheme, and validated by 2D simulations of lid-driven cavity flow, Couette flow between two concentric cylinders and flow over a circular cylinder. Finally, the method is used to simulate the sedimentation of two circular particles in a channel. The results agree very well with previous experimental and numerical data, and are more accurate than the conventional direct-forcing method, especially in the vicinity of a moving boundary.
基金WJD,JYZ,CLC,ZX,and ZGY were supported by the National Natural Science Foundation of China(Grant Number 51705143)the Education Department of Hunan Province(Grant Number 22B0464)the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant Number QL20230249).
文摘Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.
基金supported by the Key Research and Development Program of Yunnan Province(Grant No.202203AA080009)the Open Fund of State Key Laboratory of Hydraulics and Mountain River Engineering,Sichuan University(Grant No.SKHL2208).
文摘A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier-Stokes equations using the projection-based fractional step method in a fixed staggered Cartesian grid system.A volume of fluid(VOF)method with second-order accuracy is employed to trace the free surface.To represent the intricate surface geometry,the structure is discretized using the unstructured triangle mesh.Additionally,a ray tracing method is employed to classify fluid and solid points.A high-order stable scheme has been introduced to reconstruct the local velocity at interface points.Three FSI problems,including wave evolution around a breakwater,interaction between a periodic wave train and a moving float,and a 3-D moving object interacting with the free surface,were investigated to validate the accuracy and stability of the proposed model.The numerical results are in good agreement with the experimental data.Additionally,we evaluated the computational performance of the proposed GPU-based model.The GPU-based model achieved a 42.29 times speedup compared with the single-core CPU-based model in the three-dimension test.Additionally,the results regarding the time cost of each code section indicate that achieving more significant acceleration is associated with solving the turbulence,advection,and diffusion terms,while solving the pressure Poisson equation(PPE)saves the most time.Furthermore,the impact of grid number on computational efficiency indicates that as Fluid-solid interaction(FSI)immersed boundary method(IBM)graphical processing unit(GPU)two-phase flow moving rigid bodythe number of grids increases,the GPU-based model outperforms the multi-core CPU-based model.
基金the Ministry of Science and Technology of Taiwan under grants MOST 107-2115-M-035-007-MY2(C.-S.You)MOST 106-2115-M-005-005-MY2(P.-W.Hsieh)MOST 106-2115-M-008-014-MY2(S.-Y.Yang)。
文摘In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as made of fluid and introduce two virtual forcing terms.First,a virtual fluid force distributed only on the solid region is appended to the momentum equation to make the region behave like a real solid body and satisfy the prescribed velocity.Second,a virtual heat source located inside the solid region near the boundary is added to the energy transport equation to impose the thermal boundary condition on the solid boundary.We take the implicit second-order backward differentiation to discretize the time variable and employ the Choi-Moin projection scheme to split the coupled system.As for spatial discretization,second-order centered differences over a staggered Cartesian grid are used on the entire domain.The advantages of this method are its conceptual simplicity and ease of implementation.Numerical experiments are performed to demonstrate the high performance of the proposed method.Convergence tests show that the spatial convergence rates of all unknowns seem to be super-linear in the 1-norm and 2-norm while less than linear in the maximum norm.
基金Projects support by National Nature Science Foundation of China(11172069)undergraduate key reform project Curriculum system of theories and applications of modern continuum mechanics issued by Shanghai Municipal Education Commission in 2011
