A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the n...A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).展开更多
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.展开更多
Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and ...Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.展开更多
A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders...A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations(LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution.The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.展开更多
Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit ap...Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.展开更多
The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics. The effect of the regularized delta function on the accuracy is a...The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics. The effect of the regularized delta function on the accuracy is an important subject in the property study. A method of manufactured solutions is used in the research. The computational code is first verified to be mistake-free by using smooth manufactured solutions. Then, a jump in the manufactured solution for pressure is introduced to study the accuracy of the immersed boundary method. Four kinds of regularized delta functions are used to test the effect on the accuracy analysis. By analyzing the discretization errors, the accuracy of the immersed boundary method is proved to be first-order. The results show that the regularized delta function cannot improve the accuracy, but it can change the discretization errors in the entire computational domain.展开更多
The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evalua...The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evaluate the capability of the curvilinear immersed boundary(CURVIB)method in predicting near-wall velocity and pressure fluctuations in turbulent channel flows.Simulation results show that quantities including the time-averaged streamwise velocity,the rms(root-mean-square)of velocity fluctuations,the rms of vorticity fluctuations,the shear stresses,and the correlation coefficients of u'and v"computed from the CURVIB simulations are in good agreement with those from the body-fitted simulations.More importantly,it is found that the time-averaged pressure,the rms and wavenumber-frequency spectra of pressure fluctuations computed using the CURVIB method agree well with the body-fitted results.展开更多
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.展开更多
The performance of a single or the collection of microswimmers strongly depends on the hydrodynamic coupling among their constituents and themselves.We present a numerical study for a single and a pair of microswimmer...The performance of a single or the collection of microswimmers strongly depends on the hydrodynamic coupling among their constituents and themselves.We present a numerical study for a single and a pair of microswimmers based on lattice Boltzmann method(LBM)simulations.Our numerical algorithm consists of two separable parts.Lagrange polynomials provide a discretization of the microswimmers and the lattice Boltzmann method captures the dynamics of the surrounding fluid.The two components couple via an immersed boundary method.We present data for a single swimmer system and our data also show the onset of collective effects and,in particular,an overall velocity increment of clusters of swimmers.展开更多
The immersed boundary method has emerged as an efficient approach for the simulation of finite-sized solid particles in complex fluid flows.However,one of the well known shortcomings of the method is the limited suppo...The immersed boundary method has emerged as an efficient approach for the simulation of finite-sized solid particles in complex fluid flows.However,one of the well known shortcomings of the method is the limited support for the simulation of light particles,i.e.particles with a density lower than that of the surrounding fluid,both in terms of accuracy and numerical stability.Although a broad literature exists,with several authors reporting different approaches for improving the stability of the method,most of these attempts introduce extra complexities and are very costly from a computational point of view.In this work,we introduce an effective force stabilizing technique,allowing to extend the stability range of the method by filtering spurious oscillations arising when dealing with light-particles,pushing down the particle-to-fluid density ratio as low as 0.04.We thoroughly validate the method comparing with both experimental and numerical data available in literature.展开更多
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.展开更多
In this paper,we analyze the stability of the Immersed Boundary Methodapplied to a membrane-fluid system with a plasma membrane immersed in an incompressibleviscous fluid.We show that for small deformations,the planar...In this paper,we analyze the stability of the Immersed Boundary Methodapplied to a membrane-fluid system with a plasma membrane immersed in an incompressibleviscous fluid.We show that for small deformations,the planar rest state isstable for a membrane with bending rigidity.The smoothed version,using a standardregularization technique for the singular force,is also shown to be stable.Furthermore,we show that the coupled fluid-membrane system is stiff and smoothing helpsto reduce the stiffness.Compared to the system of elastic fibers immersed in an incompressiblefluid,membrane with bending rigidity consist of a wider range of decayrates.Therefore numerical instability could occur more easily for an explicit methodwhen the time step size is not sufficiently small,even though the continuous problemis stable.展开更多
One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.Hi...One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.High resolution is required to capture the flow characteristics near the fluid-structure interface to enhance accuracy of the solutions within proximity of the thin deformable body.Hence,in this work,an algorithm is developed to simulate fluid-structure interactions of moving deformable structures with very thin thicknesses.In this algorithm,adaptive mesh refinement(AMR)is integrated with immersed boundary finite element method(IBFEM)with two-stage pressure-velocity corrections.Despite successive interpolation of the flow field by IBM,the governing equations were solved using a fixed structured mesh,which significantly reduces the computational time associated with mesh reconstruction.The cut-cell IBM is used to predict the body forces while FEM is used to predict deformation of the thin elastic structure in order to integrate the motions of the fluid and solid at the interface.AMR is used to discretize the governing equations and obtain solutions that efficiently capture the thin boundary layer at the fluid-solid interface.The AMR-IBFEM algorithm is first verified by comparing the drag coefficient,lift coefficient,and Strouhal number for a benchmark case(laminar flow past a circular cylinder at Re=100)and the results showed good agreement with those of other researchers.The algorithm is then used to simulate 2-D laminar flows past stationary and moving thin structures positioned perpendicular to the freestream direction.The results also showed good agreement with those obtained from the arbitrary Lagrangian-Eulerian(ALE)algorithm for elastic thin boundaries.It is concluded that the AMR-IBFEM algorithm is capable of predicting the characteristics of laminar flow past an elastic structure with acceptable accuracy(error of-0.02%)with only-1%of the computational time for simulations with full mesh refinement.展开更多
A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.Thi...A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.This is used to represent actin biopolymers,which are constituent elements of the cytoskeleton,a complex network-like structure that plays a fundamental role in shape morphology.An extension of the traditional immersed boundary method to include a stochastic stress tensor is also proposed in order to model the thermal fluctuations in the fluid at smaller scales.By way of validation,the response of a single,massless,inextensible semiflexible filament immersed in a thermally fluctuating fluid is obtained using the suggested numerical scheme and the resulting time-averaged contraction of the filament is compared to the theoretical value obtained from the worm-like chain model.展开更多
This paper aims to study the numerical features of a coupling scheme between the immersed boundary(IB)method and the lattice Boltzmann BGK(LBGK)model by four typical test problems:the relaxation of a circular membrane...This paper aims to study the numerical features of a coupling scheme between the immersed boundary(IB)method and the lattice Boltzmann BGK(LBGK)model by four typical test problems:the relaxation of a circular membrane,the shearing flow induced by a moving fiber in the middle of a channel,the shearing flow near a non-slip rigid wall,and the circular Couette flow between two inversely rotating cylinders.The accuracy and robustness of the IB-LBGK coupling scheme,the performances of different discrete Dirac delta functions,the effect of iteration on the coupling scheme,the importance of the external forcing term treatment,the sensitivity of the coupling scheme to flow and boundary parameters,the velocity slip near non-slip rigid wall,and the origination of numerical instabilities are investigated in detail via the four test cases.It is found that the iteration in the coupling cycle can effectively improve stability,the introduction of a second-order forcing term in LBGK model is crucial,the discrete fiber segment length and the orientation of the fiber boundary obviously affect accuracy and stability,and the emergence of both temporal and spatial fluctuations of boundary parameters seems to be the indication of numerical instability.These elaborate results shed light on the nature of the coupling scheme and may benefit those who wish to use or improve the method.展开更多
In this paper,we present an immersed boundary method for simulating moving contact lines with surfactant.The governing equations are the incompressible Navier-Stokes equations with the usual mixture of Eulerian fluid ...In this paper,we present an immersed boundary method for simulating moving contact lines with surfactant.The governing equations are the incompressible Navier-Stokes equations with the usual mixture of Eulerian fluid variables and Lagrangian interfacial markers.The immersed boundary force has two components:one from the nonhomogeneous surface tension determined by the distribution of surfactant along the fluid interface,and the other from unbalanced Young’s force at the moving contact lines.An artificial tangential velocity has been added to the Lagrangian markers to ensure that the markers are uniformly distributed at all times.The corresponding modified surfactant equation is solved in a way such that the total surfactant mass is conserved.Numerical experiments including convergence analysis are carefully conducted.The effect of the surfactant on the motion of hydrophilic and hydrophobic drops are investigated in detail.展开更多
A coupling framework that leverages the advantages of the diffuse and sharp interface immersed boundary(IB)methods is presented for handling the interaction among particles and particles with the static complex geomet...A coupling framework that leverages the advantages of the diffuse and sharp interface immersed boundary(IB)methods is presented for handling the interaction among particles and particles with the static complex geometries of the environment.In the proposed coupling approach,the curvilinear IB method is employed to represent the static complex geometries,a variant of the direct forcing IB method is proposed for simulating particles,and the discrete element method is employed for particle-particle and particle-wall collisions.The proposed approach is validated using several classical benchmark problems,which include flow around a sphere,sedimentation of a sphere,collision of two sedimenting spheres,and collision between a particle and a flat wall,with the present predictions showing an overall good agreement with the results reported in the literature.The capability of the proposed framework is further demonstrated by simulating the interaction between multiple particles and a wall-mounted cylinder,and the particle-laden turbulent flow over periodic hills.The proposed method provides an efficient way to simulate particle-laden turbulent flows in environments with complex boundaries.展开更多
We extend the immersed boundary(IB)method to simulate the dynamics of a 2D dry foam by including the topological changes of the bubble network.In the article[Y.Kim,M.-C.Lai,and C.S.Peskin,J.Comput.Phys.229:5194-5207,2...We extend the immersed boundary(IB)method to simulate the dynamics of a 2D dry foam by including the topological changes of the bubble network.In the article[Y.Kim,M.-C.Lai,and C.S.Peskin,J.Comput.Phys.229:5194-5207,2010],we implemented an IB method for the foam problem in the two-dimensional case,and tested it by verifying the von Neumann relation which governs the coarsening of a two-dimensional dry foam.However,the method implemented in that article had an important limitation;we did not allow for the resolution of quadruple or higher order junctions into triple junctions.A total shrinkage of a bubble with more than four edges generates a quadruple or higher order junction.In reality,a higher order junction is unstable and resolves itself into triple junctions.We here extend the methodology previously introduced by allowing topological changes,and we illustrate the significance of such topological changes by comparing the behaviors of foams in which topological changes are allowed to those in which they are not.展开更多
The paper presents a novel pressure-corrected formulation of the immersed boundary method(IBM)for the simulation of fully compressible non-Boussinesq natural convection flows.The formulation incorporated into the pres...The paper presents a novel pressure-corrected formulation of the immersed boundary method(IBM)for the simulation of fully compressible non-Boussinesq natural convection flows.The formulation incorporated into the pressure-based fractional step approach facilitates simulation of the flows in the presence of an immersed body characterized by a complex geometry.Here,we first present extensive grid independence and verification studies addressing incompressible pressure-driven flow in an extended channel and non-Boussinesq natural convection flow in a differentially heated cavity.Next,the steady-state non-Boussinesq natural convection flow developing in the presence of hot cylinders of various diameters placed within a cold square cavity is thoroughly investigated.The obtained results are presented and analyzed in terms of the spatial distribution of path lines and temperature fields and of heat flux values typical of the hot cylinder and the cold cavity surfaces.Flow characteristics of multiple steady-state solutions discovered for several configurations are presented and discussed in detail.展开更多
The immersed boundary(IB)method is an approach to problems of fluid-structure interaction in which an elastic structure is immersed in a viscous incompressible fluid.The IB formulation of such problems uses a Lagrangi...The immersed boundary(IB)method is an approach to problems of fluid-structure interaction in which an elastic structure is immersed in a viscous incompressible fluid.The IB formulation of such problems uses a Lagrangian description of the structure and an Eulerian description of the fluid.It is well known that some versions of the IB method can suffer from poor volume conservation.Methods have been introduced to improve the volume-conservation properties of the IB method,but they either have been fairly specialized,or have used complex,nonstandard Eulerian finite-difference discretizations.In this paper,we use quasi-static and dynamic benchmark problems to investigate the effect of the choice of Eulerian discretization on the volume-conservation properties of a formally second-order accurate IB method.We consider both collocated and staggered-grid discretization methods.For the tests considered herein,the staggered-grid IB scheme generally yields at least a modest improvement in volume conservation when compared to cell-centered methods,and in many cases considered in this work,the spurious volume changes exhibited by the staggered-grid IB method are more than an order of magnitude smaller than those of the collocated schemes.We also compare the performance of cell-centered schemes that use either exact or approximate projection methods.We find that the volumeconservation properties of approximate projection IB methods depend strongly on the formulation of the projection method.When used with the IB method,we find that pressure-free approximate projection methods can yield extremely poor volume conservation,whereas pressure-increment approximate projection methods yield volume conservation that is nearly identical to that of a cell-centered exact projection method.展开更多
基金financially supported by the Fundamental Research Funds for the Central Universities(Grant Nos.2018B56414 and2019B12014)the National Natural Science Foundation of China(Grant No.51609077)
文摘A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).
基金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 National Natural Science Foundation of China(Grant No.51405375)National Key Basic Research and Development Program of China(973 Program,Grant No.2011CB706606)
文摘Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.
文摘A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations(LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution.The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.
基金supported by the US National Science Foundation (DMS-0713718)
文摘Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.
基金Project supported by the National Natural Science Foundation of China (No. 11102108)the Shanghai Leading Academic Discipline Project (No. B206)
文摘The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics. The effect of the regularized delta function on the accuracy is an important subject in the property study. A method of manufactured solutions is used in the research. The computational code is first verified to be mistake-free by using smooth manufactured solutions. Then, a jump in the manufactured solution for pressure is introduced to study the accuracy of the immersed boundary method. Four kinds of regularized delta functions are used to test the effect on the accuracy analysis. By analyzing the discretization errors, the accuracy of the immersed boundary method is proved to be first-order. The results show that the regularized delta function cannot improve the accuracy, but it can change the discretization errors in the entire computational domain.
基金This work was supported by the National Natural Science Foundation of China(NSFC)Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”(No.11988102)the Strategic Priority Research Program,Chinese Academy of Sciences(CAS)(No.XDB22040104).
文摘The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evaluate the capability of the curvilinear immersed boundary(CURVIB)method in predicting near-wall velocity and pressure fluctuations in turbulent channel flows.Simulation results show that quantities including the time-averaged streamwise velocity,the rms(root-mean-square)of velocity fluctuations,the rms of vorticity fluctuations,the shear stresses,and the correlation coefficients of u'and v"computed from the CURVIB simulations are in good agreement with those from the body-fitted simulations.More importantly,it is found that the time-averaged pressure,the rms and wavenumber-frequency spectra of pressure fluctuations computed using the CURVIB method agree well with the body-fitted results.
基金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.
基金The DFG Priority Programme SPP 1726“Microswimmers—From Single Particle Motion to Collective Behaviour”(HA 4382/5-1)and SFB 1411(Project-ID 416229255)supported this work.
文摘The performance of a single or the collection of microswimmers strongly depends on the hydrodynamic coupling among their constituents and themselves.We present a numerical study for a single and a pair of microswimmers based on lattice Boltzmann method(LBM)simulations.Our numerical algorithm consists of two separable parts.Lagrange polynomials provide a discretization of the microswimmers and the lattice Boltzmann method captures the dynamics of the surrounding fluid.The two components couple via an immersed boundary method.We present data for a single swimmer system and our data also show the onset of collective effects and,in particular,an overall velocity increment of clusters of swimmers.
文摘The immersed boundary method has emerged as an efficient approach for the simulation of finite-sized solid particles in complex fluid flows.However,one of the well known shortcomings of the method is the limited support for the simulation of light particles,i.e.particles with a density lower than that of the surrounding fluid,both in terms of accuracy and numerical stability.Although a broad literature exists,with several authors reporting different approaches for improving the stability of the method,most of these attempts introduce extra complexities and are very costly from a computational point of view.In this work,we introduce an effective force stabilizing technique,allowing to extend the stability range of the method by filtering spurious oscillations arising when dealing with light-particles,pushing down the particle-to-fluid density ratio as low as 0.04.We thoroughly validate the method comparing with both experimental and numerical data available in literature.
基金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.
基金the National Natural Science Foundation of China(No.10472070)。
文摘In this paper,we analyze the stability of the Immersed Boundary Methodapplied to a membrane-fluid system with a plasma membrane immersed in an incompressibleviscous fluid.We show that for small deformations,the planar rest state isstable for a membrane with bending rigidity.The smoothed version,using a standardregularization technique for the singular force,is also shown to be stable.Furthermore,we show that the coupled fluid-membrane system is stiff and smoothing helpsto reduce the stiffness.Compared to the system of elastic fibers immersed in an incompressiblefluid,membrane with bending rigidity consist of a wider range of decayrates.Therefore numerical instability could occur more easily for an explicit methodwhen the time step size is not sufficiently small,even though the continuous problemis stable.
文摘One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.High resolution is required to capture the flow characteristics near the fluid-structure interface to enhance accuracy of the solutions within proximity of the thin deformable body.Hence,in this work,an algorithm is developed to simulate fluid-structure interactions of moving deformable structures with very thin thicknesses.In this algorithm,adaptive mesh refinement(AMR)is integrated with immersed boundary finite element method(IBFEM)with two-stage pressure-velocity corrections.Despite successive interpolation of the flow field by IBM,the governing equations were solved using a fixed structured mesh,which significantly reduces the computational time associated with mesh reconstruction.The cut-cell IBM is used to predict the body forces while FEM is used to predict deformation of the thin elastic structure in order to integrate the motions of the fluid and solid at the interface.AMR is used to discretize the governing equations and obtain solutions that efficiently capture the thin boundary layer at the fluid-solid interface.The AMR-IBFEM algorithm is first verified by comparing the drag coefficient,lift coefficient,and Strouhal number for a benchmark case(laminar flow past a circular cylinder at Re=100)and the results showed good agreement with those of other researchers.The algorithm is then used to simulate 2-D laminar flows past stationary and moving thin structures positioned perpendicular to the freestream direction.The results also showed good agreement with those obtained from the arbitrary Lagrangian-Eulerian(ALE)algorithm for elastic thin boundaries.It is concluded that the AMR-IBFEM algorithm is capable of predicting the characteristics of laminar flow past an elastic structure with acceptable accuracy(error of-0.02%)with only-1%of the computational time for simulations with full mesh refinement.
文摘A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.This is used to represent actin biopolymers,which are constituent elements of the cytoskeleton,a complex network-like structure that plays a fundamental role in shape morphology.An extension of the traditional immersed boundary method to include a stochastic stress tensor is also proposed in order to model the thermal fluctuations in the fluid at smaller scales.By way of validation,the response of a single,massless,inextensible semiflexible filament immersed in a thermally fluctuating fluid is obtained using the suggested numerical scheme and the resulting time-averaged contraction of the filament is compared to the theoretical value obtained from the worm-like chain model.
基金the National Natural Science Foundation of China(NSFC,Grant numbers 10572106,10872153 and 11172219)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20130141110013)。
文摘This paper aims to study the numerical features of a coupling scheme between the immersed boundary(IB)method and the lattice Boltzmann BGK(LBGK)model by four typical test problems:the relaxation of a circular membrane,the shearing flow induced by a moving fiber in the middle of a channel,the shearing flow near a non-slip rigid wall,and the circular Couette flow between two inversely rotating cylinders.The accuracy and robustness of the IB-LBGK coupling scheme,the performances of different discrete Dirac delta functions,the effect of iteration on the coupling scheme,the importance of the external forcing term treatment,the sensitivity of the coupling scheme to flow and boundary parameters,the velocity slip near non-slip rigid wall,and the origination of numerical instabilities are investigated in detail via the four test cases.It is found that the iteration in the coupling cycle can effectively improve stability,the introduction of a second-order forcing term in LBGK model is crucial,the discrete fiber segment length and the orientation of the fiber boundary obviously affect accuracy and stability,and the emergence of both temporal and spatial fluctuations of boundary parameters seems to be the indication of numerical instability.These elaborate results shed light on the nature of the coupling scheme and may benefit those who wish to use or improve the method.
基金supported in part by National Science Council of Taiwan under research grant NSC-97-2628-M-009-007-MY3 and MoE-ATU projectsupported by grants from the Natural Science and Engineering Research Council(NSERC)of Canada and the Mathematics of Information Technology and Complex Systems(MITACS)of Canada.
文摘In this paper,we present an immersed boundary method for simulating moving contact lines with surfactant.The governing equations are the incompressible Navier-Stokes equations with the usual mixture of Eulerian fluid variables and Lagrangian interfacial markers.The immersed boundary force has two components:one from the nonhomogeneous surface tension determined by the distribution of surfactant along the fluid interface,and the other from unbalanced Young’s force at the moving contact lines.An artificial tangential velocity has been added to the Lagrangian markers to ensure that the markers are uniformly distributed at all times.The corresponding modified surfactant equation is solved in a way such that the total surfactant mass is conserved.Numerical experiments including convergence analysis are carefully conducted.The effect of the surfactant on the motion of hydrophilic and hydrophobic drops are investigated in detail.
基金supported by NSFC Basic Science Center Program for"Multiscale Problems in Nonlinear Mechanics"(Grant No.11988102)National Natural Science Foundation of China(Grant No.12172360)Institute of Mechanics CAS,and Chinese Academy of Sciences.J.Qin also acknowledges the support from China Postdoctoral Science Foundation(Grant No.2021M693241).
文摘A coupling framework that leverages the advantages of the diffuse and sharp interface immersed boundary(IB)methods is presented for handling the interaction among particles and particles with the static complex geometries of the environment.In the proposed coupling approach,the curvilinear IB method is employed to represent the static complex geometries,a variant of the direct forcing IB method is proposed for simulating particles,and the discrete element method is employed for particle-particle and particle-wall collisions.The proposed approach is validated using several classical benchmark problems,which include flow around a sphere,sedimentation of a sphere,collision of two sedimenting spheres,and collision between a particle and a flat wall,with the present predictions showing an overall good agreement with the results reported in the literature.The capability of the proposed framework is further demonstrated by simulating the interaction between multiple particles and a wall-mounted cylinder,and the particle-laden turbulent flow over periodic hills.The proposed method provides an efficient way to simulate particle-laden turbulent flows in environments with complex boundaries.
基金supported by National Research Foundation of Korea Grant funded by the Korean Government(2010-0006165)The second author was supported by the Chung-Ang University Research Scholarship Grant in 2010The third author is supported in part by National Science Council of Taiwan under research grant NSC-97-2628-M-009-007-MY3,NSC-98-2115-M-009-014-MY3,and the support of NCTS in Taiwan.
文摘We extend the immersed boundary(IB)method to simulate the dynamics of a 2D dry foam by including the topological changes of the bubble network.In the article[Y.Kim,M.-C.Lai,and C.S.Peskin,J.Comput.Phys.229:5194-5207,2010],we implemented an IB method for the foam problem in the two-dimensional case,and tested it by verifying the von Neumann relation which governs the coarsening of a two-dimensional dry foam.However,the method implemented in that article had an important limitation;we did not allow for the resolution of quadruple or higher order junctions into triple junctions.A total shrinkage of a bubble with more than four edges generates a quadruple or higher order junction.In reality,a higher order junction is unstable and resolves itself into triple junctions.We here extend the methodology previously introduced by allowing topological changes,and we illustrate the significance of such topological changes by comparing the behaviors of foams in which topological changes are allowed to those in which they are not.
基金financial support for this work(grant 218-11-038).
文摘The paper presents a novel pressure-corrected formulation of the immersed boundary method(IBM)for the simulation of fully compressible non-Boussinesq natural convection flows.The formulation incorporated into the pressure-based fractional step approach facilitates simulation of the flows in the presence of an immersed body characterized by a complex geometry.Here,we first present extensive grid independence and verification studies addressing incompressible pressure-driven flow in an extended channel and non-Boussinesq natural convection flow in a differentially heated cavity.Next,the steady-state non-Boussinesq natural convection flow developing in the presence of hot cylinders of various diameters placed within a cold square cavity is thoroughly investigated.The obtained results are presented and analyzed in terms of the spatial distribution of path lines and temperature fields and of heat flux values typical of the hot cylinder and the cold cavity surfaces.Flow characteristics of multiple steady-state solutions discovered for several configurations are presented and discussed in detail.
基金support from American Heart Association grant 10SDG4320049National Science Foundation grants DMS 1016554 and OCI 1047734.
文摘The immersed boundary(IB)method is an approach to problems of fluid-structure interaction in which an elastic structure is immersed in a viscous incompressible fluid.The IB formulation of such problems uses a Lagrangian description of the structure and an Eulerian description of the fluid.It is well known that some versions of the IB method can suffer from poor volume conservation.Methods have been introduced to improve the volume-conservation properties of the IB method,but they either have been fairly specialized,or have used complex,nonstandard Eulerian finite-difference discretizations.In this paper,we use quasi-static and dynamic benchmark problems to investigate the effect of the choice of Eulerian discretization on the volume-conservation properties of a formally second-order accurate IB method.We consider both collocated and staggered-grid discretization methods.For the tests considered herein,the staggered-grid IB scheme generally yields at least a modest improvement in volume conservation when compared to cell-centered methods,and in many cases considered in this work,the spurious volume changes exhibited by the staggered-grid IB method are more than an order of magnitude smaller than those of the collocated schemes.We also compare the performance of cell-centered schemes that use either exact or approximate projection methods.We find that the volumeconservation properties of approximate projection IB methods depend strongly on the formulation of the projection method.When used with the IB method,we find that pressure-free approximate projection methods can yield extremely poor volume conservation,whereas pressure-increment approximate projection methods yield volume conservation that is nearly identical to that of a cell-centered exact projection method.