An immersed-boundary method is used to investigate the flapping wings with different aspect ratios ranging from 1 to 5.The numerical results on wake structures and the performance of the propulsion are given.Unlike th...An immersed-boundary method is used to investigate the flapping wings with different aspect ratios ranging from 1 to 5.The numerical results on wake structures and the performance of the propulsion are given.Unlike the case of the two-dimensional flapping foil,the wing-tip vortices appear for the flow past a three-dimensional flapping wing,which makes the wake vortex structures much different.The results show that the leading edge vortex merges into the trailing edge vortex,connects with the wing tip vortices and then sheds from the wing.A vortex ring forms in the wake,and exhibits different patterns for different foil aspect ratios.Analysis of hydrodynamic performances shows that both thrust coefficient and efficiency of the flapping wing increase with increasing aspect ratio.展开更多
Two-dimensional numerical simulations are performed to study the propulsive performance of fish-like swimming foils using the immersed-boundary method. A single fish as well as two fishes in tandem arrangement are stu...Two-dimensional numerical simulations are performed to study the propulsive performance of fish-like swimming foils using the immersed-boundary method. A single fish as well as two fishes in tandem arrangement are studied. First, the effect of the phase speed on the propulsive performance of a single fish is analyzed. The wake structures and pressure distribution near the wavy fish are also examined. The results show good correlation with those by previous researchers. Second, two tandem fishes with the same phase speed and amplitude are studied. The results show that the fish situated directly behind another one endure a higher thrust than that of a single one.展开更多
The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive comp...The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.展开更多
Nowadays,the design of fixed packed bed reactors still relies on empirical correlations,which,especially for small tube to particle diameter ratios,are mostly too inaccurate because of the presence of wall effects.The...Nowadays,the design of fixed packed bed reactors still relies on empirical correlations,which,especially for small tube to particle diameter ratios,are mostly too inaccurate because of the presence of wall effects.Therefore,the simulation of fixed packed bed reactors plays an important role to predict and control the flow and process parameters in such,nowadays and in the future.Because of its straightforward applicability to non-uniform packings with particles of arbitrary shapes,the immersed boundary method(IBM)has advantages over other numerical methods and is used more and more frequently.This paper compares two approaches of IBMs for the simulation of fixed bed reactors with spherical shaped particles.The classic,smooth approach is compared to the straightforward to implement blocked-off method for velocity fields above the fixed bed for particle Reynolds numbers of 300 and 500.Results from experimental inline PIV-measurements of the reactor to be simulated serve as a basis for comparison.Very good agreement with the experiment is found for both simulation methodologies with higher resolutions,considering the more stable flow at a particle Reynolds number of 300.Differences in the different IBM approaches occurred for the more unsteady flow at a particle Reynolds number of 500.Compared to the blocked-off method,the smooth IBM reflects the formation of additional jets and recirculation zones better right above the bed,though increasing the fluid mesh resolution improves the accuracy of the blocked-off method.Overall,a more diffusive behaviour is found for the blocked-off simulations due to the stairstep representation,which is avoided by using interpolation stencils as in the smooth IBM.With higher mesh refinement in the blocked-off IBM this effect can be reduced,but this also increases the computational effort.展开更多
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.展开更多
The present paper presents an experimental and numerical investigation of the dispersion of the gaseous jet flow and co-flow for the simple unit cell(SUC)and body-centred cubic(BCC)configuration of particles in packed...The present paper presents an experimental and numerical investigation of the dispersion of the gaseous jet flow and co-flow for the simple unit cell(SUC)and body-centred cubic(BCC)configuration of particles in packed beds.The experimental setup is built in such a way that suitable and simplified boundary conditions are imposed for the corresponding numerical framework,so the simulations can be done under very similar conditions as the experiments.Accordingly,a porous plate is used for the co-flow to achieve the uniform velocity and the fully developed flow is ensured for the jet flow.The SUC and BCC particle beds consist of 3D-printed spheres,and the non-isotropy near the walls is mostly eliminated by placing half-spheres at the channel walls.The flow velocities are analysed directly at the exit of the particle bed for both beds over 36 pores for the SUC configuration and 60 pores for the BCC configuration,for particle Reynolds numbers of 200,300,and 400.Stereo particle image velocimetry is experimentally arranged in such a way that the velocities over the entire region at the exit of the packed bed are obtained instantaneously.The numerical method consists of a state-of-the-art immersed boundary method with adaptive mesh refinement.The paper presents the pore jet structure and velocity field exiting from each pore for the SUC and BCC packed particle beds.The numerical and experimental studies show a good agreement for the SUC configuration for all flow velocities.For the BCC configuration,some differences can be observed in the pore jet flow structure between the simulations and the experiments,but the general flow velocity distribution shows a good overall agreement.The axial velocity is generally higher for the pores located near the centre of the packed bed than for the pores near the wall.In addition,the axial velocities are observed to increase near the peripheral pores of the packed bed.This behaviour is predominant for the BCC configuration as compared to the SUC configuration.The velocities near the per展开更多
Simulating unsteady turbulent flow in turbomachines is still challenging due to the complexity of blade geometry and relative motion between rotor and stator.This study presents an Immersed Boundary Method(IBM)for hig...Simulating unsteady turbulent flow in turbomachines is still challenging due to the complexity of blade geometry and relative motion between rotor and stator.This study presents an Immersed Boundary Method(IBM)for high-Reynolds turbomachinery internal flows,and shows the advantage of the automatic grid generation techniques and flexible moving boundary treatments.The wall functions are used in the present method to alleviate the wall resolution restriction of turbulence simulation.The Two-Dimensional(2-D)IBM solver,which was previously developed and tested for a low-speed compressor,is further validated for a well-documented Low-Pressure Turbine(LPT)cascade.Both the blade loading and the total pressure losses in the wake are well captured by the present 2-D solver.The complex Three-Dimensional(3-D)effects in turbomachines motivate the further development of an extended 3-D IBM solver by using a curvilinear-coordinate system that facilitates the hub and casing boundary treatment.The good performance of the 3-D solver is demonstrated through comparison with CFX solver solutions for the rotor configuration of Advanced Noise Control Fan(ANCF).Further effects of the grid resolution on capturing the blade wake are discussed.The results indicate that the present 3-D solver is capable of reproducing the evolution of the blade wake with suitable computational grid.展开更多
This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid mome...This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid momentum equations are solved by using the fractional step method,while the energy equation is solved by employing the second-order Adams-Bashforth method.For efficient load balancing between the CPU cores,the coupling between particles and fluid is obtained by applying the body force in the fluid equations,which depends on the solid volume fraction of particles contained in each grid cell,and then by linearly interpolating the particle temperature and velocity on the fluid grid cell(in place of the delta function commonly used in literature).Several test cases from the literature are studied,and good agreement is observed between the simulation results and the literature.Finally,a scaling study on multiple core machines is performed,demonstrating the proposed IBM's capabilities for a significant reduction in processing time.展开更多
A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface i...A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.展开更多
A resolved CFD-DEM method is proposed to simulate the fluid-particle interaction for large complex granules.The airflow in a vertical sinter fixed bed is numerically studied using this method.The multi-sphere clumped ...A resolved CFD-DEM method is proposed to simulate the fluid-particle interaction for large complex granules.The airflow in a vertical sinter fixed bed is numerically studied using this method.The multi-sphere clumped method is used to create irregular sinter particles in DEM.The immersed boundary method and dynamic cell refinement are applied to describe the fluid flow around particles with higher resolution,by which the fluid-particle interaction can be simulated more accurately.The simulation results presented the packing voidage distributions and the airflow fields in the sinter beds of different single and mixed particle size ranges.The bed pressure drops were simulated and the results were compared with the corresponding experimental ones.The good agreement indicated that the proposed resolved CFD-DEM method is an effective tool to model the fluid-particle interaction for irregular large granules in the gas-solid multi-phase systems.展开更多
The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural cha...The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural changes of such a network induced by diabetes.A cubic region representing local skin tissue is selected as the computational domain,which in turn includes two intravascular and extravascular sub-domains.To save computational resources,the capillary network is reduced to a 1D pipeline model and embedded into the extravascular region.On the basis of the immersed boundary method(IBM)strategy,fluid and heat fluxes across a capillary wall are distributed to the surrounding tissue nodes by a delta function.We consider both steady and periodic blood pressure conditions at the entrances of the capillary network.Under steady blood pressure conditions,both the interstitial fluid pressure and tissue temperature around the capillary network are larger than those in other places.When the periodic blood pressure condition is considered,tissue temperature tends to fluctuate with the same frequency of the forcing,but the related waveform displays a smaller amplitude and a certain time(phase)delay.When the connectivity of capillary network is diminished,the capacity of blood redistribution through the capillary network becomes weaker and a subset of the vessel branches lose blood flow,which further aggravates the amplitude attenuation and time delay of the skin temperature fluctuation.展开更多
A numerical study of linear wave scattering over a floating platform has been simulated by an efficient numericalmodel in this letter.The non-hydrostatic model is used to simulate the free surface and the uneven botto...A numerical study of linear wave scattering over a floating platform has been simulated by an efficient numericalmodel in this letter.The non-hydrostatic model is used to simulate the free surface and the uneven bottom.For thesolid body modelling,the immersed boundary method(IBM)is implemented by introducing a virtual boundaryforce into the momentum equations to emulate the boundary conditions.This implementation enhances theability of the model to simulate interactions between waves and floating structures.A numerical case involvingwave interactions with a floating platform is studied to validate the numerical model.By simulating the wavepropagation,the numerical model captures the variation of the wave scattering very well,which verifies theperformance of the numerical model and the robust strategy of the IBM.展开更多
The interaction between free surface flow and structure is investigated using a new level set immersed boundary method.The incorporation of an improved immersed boundary method with a free surface capture scheme imple...The interaction between free surface flow and structure is investigated using a new level set immersed boundary method.The incorporation of an improved immersed boundary method with a free surface capture scheme implemented in a Navier-Stokes solver allows the interaction between fluid flow with free surface and moving body/bodies of almost arbitrary shape to be modelled.A new algorithm is proposed to locate exact forcing points near solid boundaries,which provides an accurate numerical solution.The discretized linear system of the Poisson pressure equation is solved using the Generalized Minimum Residual(GMRES)method with incomplete LU preconditioning.Uniform flow past a cylinder at Reynolds number Re=100 is modelled using the present model and results agree well with the experiment and numerical data in the literature.Water exit and entry of a cylinder at the prescribed velocity is also investigated.The predicted slamming coefficient is in good agreement with experimental data and previous numerical simulations using a ComFlow model.The vertical slamming force and pressure distribution for the free falling wedge is also studied by the present model and comparisons with available theoretical solutions and experimental data are made.展开更多
Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls usi...Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.展开更多
An efficient immersed boundary-lattice Boltzmann method(IB-LBM)is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows.Stress LBM based on Giesekus and Oldroyd-B constitutive equa...An efficient immersed boundary-lattice Boltzmann method(IB-LBM)is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows.Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor.A boundary thickening-based direct forcing IB method is adopted to solve the particle-fluid interactions with high accuracy for non-slip boundary conditions.A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computa-tional parameters.An asynchronous calculation strategy is adopted to further improve the computing efficiency.The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods.The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed.Subsequently,504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh-Taylor instabilities was observed.At last,simulations of a dense flow involving 11001 particles,the largest number of particles to date,were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid.The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time.These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.展开更多
基金supported by the Doctoral Research Foundation of Chinese Universities (Grant No. 20070335066)the National Natural Science Foundation of China (Grant Nos. 50735004, 10802075)
文摘An immersed-boundary method is used to investigate the flapping wings with different aspect ratios ranging from 1 to 5.The numerical results on wake structures and the performance of the propulsion are given.Unlike the case of the two-dimensional flapping foil,the wing-tip vortices appear for the flow past a three-dimensional flapping wing,which makes the wake vortex structures much different.The results show that the leading edge vortex merges into the trailing edge vortex,connects with the wing tip vortices and then sheds from the wing.A vortex ring forms in the wake,and exhibits different patterns for different foil aspect ratios.Analysis of hydrodynamic performances shows that both thrust coefficient and efficiency of the flapping wing increase with increasing aspect ratio.
基金Project supported by the National Natural Science Foundation of China (Grant No: 10472104) and National Laboratory of Hydrodynamics of China.
文摘Two-dimensional numerical simulations are performed to study the propulsive performance of fish-like swimming foils using the immersed-boundary method. A single fish as well as two fishes in tandem arrangement are studied. First, the effect of the phase speed on the propulsive performance of a single fish is analyzed. The wake structures and pressure distribution near the wavy fish are also examined. The results show good correlation with those by previous researchers. Second, two tandem fishes with the same phase speed and amplitude are studied. The results show that the fish situated directly behind another one endure a higher thrust than that of a single one.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.12072217 and 42077254)the Natural Science Foundation of Hunan Province,China(Grant No.2022JJ30567).
文摘The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.
基金funded by the Deutsche Forschungsgemeinschaft(DFG,German Research Foundation)-Project-ID 422037413-TRR 287Gefordert durch die Deutsche Forschungsgemeinschaft(DFG)-Projektnummer 422037413-TRR 287.
文摘Nowadays,the design of fixed packed bed reactors still relies on empirical correlations,which,especially for small tube to particle diameter ratios,are mostly too inaccurate because of the presence of wall effects.Therefore,the simulation of fixed packed bed reactors plays an important role to predict and control the flow and process parameters in such,nowadays and in the future.Because of its straightforward applicability to non-uniform packings with particles of arbitrary shapes,the immersed boundary method(IBM)has advantages over other numerical methods and is used more and more frequently.This paper compares two approaches of IBMs for the simulation of fixed bed reactors with spherical shaped particles.The classic,smooth approach is compared to the straightforward to implement blocked-off method for velocity fields above the fixed bed for particle Reynolds numbers of 300 and 500.Results from experimental inline PIV-measurements of the reactor to be simulated serve as a basis for comparison.Very good agreement with the experiment is found for both simulation methodologies with higher resolutions,considering the more stable flow at a particle Reynolds number of 300.Differences in the different IBM approaches occurred for the more unsteady flow at a particle Reynolds number of 500.Compared to the blocked-off method,the smooth IBM reflects the formation of additional jets and recirculation zones better right above the bed,though increasing the fluid mesh resolution improves the accuracy of the blocked-off method.Overall,a more diffusive behaviour is found for the blocked-off simulations due to the stairstep representation,which is avoided by using interpolation stencils as in the smooth IBM.With higher mesh refinement in the blocked-off IBM this effect can be reduced,but this also increases the computational effort.
基金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.
基金funded by the Deutsche Forschungsgemeinschaft(DFG,German Research Foundation)-Project-ID 422037413-TRR 287Gefördert durch die Deutsche Forschungsgemeinschaft(DFG)-Projektnummer 422037413-TRR 287.
文摘The present paper presents an experimental and numerical investigation of the dispersion of the gaseous jet flow and co-flow for the simple unit cell(SUC)and body-centred cubic(BCC)configuration of particles in packed beds.The experimental setup is built in such a way that suitable and simplified boundary conditions are imposed for the corresponding numerical framework,so the simulations can be done under very similar conditions as the experiments.Accordingly,a porous plate is used for the co-flow to achieve the uniform velocity and the fully developed flow is ensured for the jet flow.The SUC and BCC particle beds consist of 3D-printed spheres,and the non-isotropy near the walls is mostly eliminated by placing half-spheres at the channel walls.The flow velocities are analysed directly at the exit of the particle bed for both beds over 36 pores for the SUC configuration and 60 pores for the BCC configuration,for particle Reynolds numbers of 200,300,and 400.Stereo particle image velocimetry is experimentally arranged in such a way that the velocities over the entire region at the exit of the packed bed are obtained instantaneously.The numerical method consists of a state-of-the-art immersed boundary method with adaptive mesh refinement.The paper presents the pore jet structure and velocity field exiting from each pore for the SUC and BCC packed particle beds.The numerical and experimental studies show a good agreement for the SUC configuration for all flow velocities.For the BCC configuration,some differences can be observed in the pore jet flow structure between the simulations and the experiments,but the general flow velocity distribution shows a good overall agreement.The axial velocity is generally higher for the pores located near the centre of the packed bed than for the pores near the wall.In addition,the axial velocities are observed to increase near the peripheral pores of the packed bed.This behaviour is predominant for the BCC configuration as compared to the SUC configuration.The velocities near the per
基金co-supported by the National Natural Science Foundation of China(No.52022009)the Science Center for Gas Turbine Project of China(No.P2022-A-II-003-001)the Key Laboratory Foundation,China(No.2021-JCJQ-LB-062-0102).
文摘Simulating unsteady turbulent flow in turbomachines is still challenging due to the complexity of blade geometry and relative motion between rotor and stator.This study presents an Immersed Boundary Method(IBM)for high-Reynolds turbomachinery internal flows,and shows the advantage of the automatic grid generation techniques and flexible moving boundary treatments.The wall functions are used in the present method to alleviate the wall resolution restriction of turbulence simulation.The Two-Dimensional(2-D)IBM solver,which was previously developed and tested for a low-speed compressor,is further validated for a well-documented Low-Pressure Turbine(LPT)cascade.Both the blade loading and the total pressure losses in the wake are well captured by the present 2-D solver.The complex Three-Dimensional(3-D)effects in turbomachines motivate the further development of an extended 3-D IBM solver by using a curvilinear-coordinate system that facilitates the hub and casing boundary treatment.The good performance of the 3-D solver is demonstrated through comparison with CFX solver solutions for the rotor configuration of Advanced Noise Control Fan(ANCF).Further effects of the grid resolution on capturing the blade wake are discussed.The results indicate that the present 3-D solver is capable of reproducing the evolution of the blade wake with suitable computational grid.
文摘This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid momentum equations are solved by using the fractional step method,while the energy equation is solved by employing the second-order Adams-Bashforth method.For efficient load balancing between the CPU cores,the coupling between particles and fluid is obtained by applying the body force in the fluid equations,which depends on the solid volume fraction of particles contained in each grid cell,and then by linearly interpolating the particle temperature and velocity on the fluid grid cell(in place of the delta function commonly used in literature).Several test cases from the literature are studied,and good agreement is observed between the simulation results and the literature.Finally,a scaling study on multiple core machines is performed,demonstrating the proposed IBM's capabilities for a significant reduction in processing time.
基金supported by the Open Project of Key Laboratory of Aerospace EDLA,CASC(No.EDL19092208)。
文摘A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.
基金the financial support for this work from the National Natural Science Foundation of China(grant No.52104340)China Postdoctoral Science Foundation(grant No.2020M672425)+1 种基金The Key Research and Development Program of Hubei Province(grant No.2022BCA058)Natural Science Foundation of Hubei Province(grant No.2020CFB133).
文摘A resolved CFD-DEM method is proposed to simulate the fluid-particle interaction for large complex granules.The airflow in a vertical sinter fixed bed is numerically studied using this method.The multi-sphere clumped method is used to create irregular sinter particles in DEM.The immersed boundary method and dynamic cell refinement are applied to describe the fluid flow around particles with higher resolution,by which the fluid-particle interaction can be simulated more accurately.The simulation results presented the packing voidage distributions and the airflow fields in the sinter beds of different single and mixed particle size ranges.The bed pressure drops were simulated and the results were compared with the corresponding experimental ones.The good agreement indicated that the proposed resolved CFD-DEM method is an effective tool to model the fluid-particle interaction for irregular large granules in the gas-solid multi-phase systems.
基金This study was supported by National Natural Science Foundation of China(NSFC No.51576033)Dalian Innovative Funding of Science and Technology(2018J12SN076)NSFC No 11602053.
文摘The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural changes of such a network induced by diabetes.A cubic region representing local skin tissue is selected as the computational domain,which in turn includes two intravascular and extravascular sub-domains.To save computational resources,the capillary network is reduced to a 1D pipeline model and embedded into the extravascular region.On the basis of the immersed boundary method(IBM)strategy,fluid and heat fluxes across a capillary wall are distributed to the surrounding tissue nodes by a delta function.We consider both steady and periodic blood pressure conditions at the entrances of the capillary network.Under steady blood pressure conditions,both the interstitial fluid pressure and tissue temperature around the capillary network are larger than those in other places.When the periodic blood pressure condition is considered,tissue temperature tends to fluctuate with the same frequency of the forcing,but the related waveform displays a smaller amplitude and a certain time(phase)delay.When the connectivity of capillary network is diminished,the capacity of blood redistribution through the capillary network becomes weaker and a subset of the vessel branches lose blood flow,which further aggravates the amplitude attenuation and time delay of the skin temperature fluctuation.
基金supported by Shanghai 2021“Science and Technology Innovation Action Plan”:Scientific and Technological Projects for Social Development(Grant No.21DZ1202701).
文摘A numerical study of linear wave scattering over a floating platform has been simulated by an efficient numericalmodel in this letter.The non-hydrostatic model is used to simulate the free surface and the uneven bottom.For thesolid body modelling,the immersed boundary method(IBM)is implemented by introducing a virtual boundaryforce into the momentum equations to emulate the boundary conditions.This implementation enhances theability of the model to simulate interactions between waves and floating structures.A numerical case involvingwave interactions with a floating platform is studied to validate the numerical model.By simulating the wavepropagation,the numerical model captures the variation of the wave scattering very well,which verifies theperformance of the numerical model and the robust strategy of the IBM.
基金the support of the South West of England Regional Development Agency in funding this project through Peninsular Research Institute for Marine Renewable Energy(http://www.primare.org/).
文摘The interaction between free surface flow and structure is investigated using a new level set immersed boundary method.The incorporation of an improved immersed boundary method with a free surface capture scheme implemented in a Navier-Stokes solver allows the interaction between fluid flow with free surface and moving body/bodies of almost arbitrary shape to be modelled.A new algorithm is proposed to locate exact forcing points near solid boundaries,which provides an accurate numerical solution.The discretized linear system of the Poisson pressure equation is solved using the Generalized Minimum Residual(GMRES)method with incomplete LU preconditioning.Uniform flow past a cylinder at Reynolds number Re=100 is modelled using the present model and results agree well with the experiment and numerical data in the literature.Water exit and entry of a cylinder at the prescribed velocity is also investigated.The predicted slamming coefficient is in good agreement with experimental data and previous numerical simulations using a ComFlow model.The vertical slamming force and pressure distribution for the free falling wedge is also studied by the present model and comparisons with available theoretical solutions and experimental data are made.
基金supported by the National Natural Science Foundation of China (NSFC) Basic Science Center Program for “Multiscale Problems in Nonlinear Mechanics” (Grant No. 11988102)the NSFC Program (Grant No. 11772337)+3 种基金the Science Challenge Program (Grant No. TZ2016001)the Strategic Priority Research Program,Chinese Academy of Sciences (CAS) (Grant No. XDB22040104)the Key Research Program of Frontier Sciences, CAS (Grant No. QYZDJ-SSW-SYS002)the CAS Center for Excellence in Complex System Mechanics
文摘Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.
基金supported by the National Natural Science Foundation of China(Grant Nos.51876075and 51876076).
文摘An efficient immersed boundary-lattice Boltzmann method(IB-LBM)is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows.Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor.A boundary thickening-based direct forcing IB method is adopted to solve the particle-fluid interactions with high accuracy for non-slip boundary conditions.A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computa-tional parameters.An asynchronous calculation strategy is adopted to further improve the computing efficiency.The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods.The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed.Subsequently,504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh-Taylor instabilities was observed.At last,simulations of a dense flow involving 11001 particles,the largest number of particles to date,were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid.The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time.These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.
