Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data imm...Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.展开更多
The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic perf...The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic performance of a butterfly like model is studied using a solver based on immersed boundary method and adaptive mesh.For aerodynamic performance,adding a"forward-backward-forward"sweeping motion to the forewing makes it more suitable for fast cruising flight,and compared to the situation without forewing sweeping,the drag is reduced by 46%and the lift to drag ratio is increased by 45%.On the contrary,adding a"backward-for ward-backward"sweeping motion to the forewing increases lift and makes it more suitable for climb flight.For downstroke and middle to late upstroke,the forewing sweeping affects the Leading-Edge Vortex(LEV)through two factors:sweeping velocity and forward sweeping angle,and their effects are coupled.A large forward sweeping velocity can enhance the strength of LEV,while a large forward sweeping angle can weaken it.For early upstroke,the forewing sweeping can affect the wake capture mechanism,sweeping backward can enhance it while sweeping forward can weaken it.The findings in this work provide insight into the design of butterfly like Micro AirVehicles(MAVs).展开更多
The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mes...The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.展开更多
The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lat...The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lateral distance are numerically simulated using a penalty-immersed boundary method.The effects of the flapping phase and lateral distance on the propulsive performance of two fish meetings are analyzed.Results show that,when two plates meet,if their leading edges diverge laterally,the individual plate can efficiently and rapidly move apart from the other horizontally.If their leading edges converge laterally,the plate motion can be retarded,leading to high energy consumption.Moreover,an increasing lateral distance between two plates significantly weakens the fluid-structure interactions,resulting in an exponential decline in mean cruising speed.A quantitative force analysis based on vortex dynamic theory is performed to gain physics insight into the hydrodynamic interaction mechanism.It is found that lateral separation between the two leading edges enhances the vorticity generation and boundary vorticity flux on the surface of the plate,subsequently reinforcing the thrust effect and increasing horizontal velocity.This study offers insight into the hydro-dynamic mechanisms of the fluid-structure interactions among fish moving toward each other and suggests potential strategies for enhancing the maneuverability of robotic fish in complex environment.展开更多
To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced v...To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.展开更多
The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle...The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.展开更多
A facile ammonium-dichromate solution immersion method was introduced to synthesize the copperwettable Cr3C2 coating on and inside the carbon-carbon （C/C） preform. The formation mechanism and the microstructures of ...A facile ammonium-dichromate solution immersion method was introduced to synthesize the copperwettable Cr3C2 coating on and inside the carbon-carbon （C/C） preform. The formation mechanism and the microstructures of the Cr3C2 coatings were studied. The contact angle between molten copper and the C/C decreased from 140°to 60°, demonstrating the significant improvement in the wettability. The Cr3C2- coated C/C-Cu composite with only 4.2% porosity and 3.69 gcm^-3 density was manufactured through copper infiltration. As a result, the thermal and electrical conductivity of the modified C/C-Cu increased significantly due to the infiltrated copper. Also the mechanical properties of the composites including both the flexural and compressive strengths were enhanced by over 100%. The modified C/C-Cu composite exhibited lower friction coefficients and wear rates for different load levels than those of the commercial C/Cu composite. These results demonstrate the potential of the modified C/C-Cu material for use in electrical contacts.展开更多
In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body...In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.展开更多
An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in ...An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.展开更多
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 hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simu...A hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simulate the incompressible flow field.The sound field is predicted by discretizing acoustic perturbation equations(APEs)with dispersionrelation-preserving space scheme and low-dispersion and low-dissipation Runge-Kutta time integration.A sharp-interface IBM based on ghost-cell is implemented for present two-step DNS-APE approach to deal with complex moving bodies with Cartesian grids.The present method is validated through simulations of sound generation caused by flow past a rotating cylinder,an oscillating cylinder,and tandem oscillating and stationary cylinders.The sound generated by typical kinds of complicated bio-inspired locomotions,i.e.,flapping flight by wings of varied shapes and collective undulatory swimming in tandem,are investigated using present method.The results demonstrate potential of the hybrid approach in addressing flowinduced sound generation and propagation with complex moving boundaries in a fluid medium,especially for the sound characteristics of bio-mimetic flows,which might shed lights on investigations on bio-acoustics,ethology of complex animal system,and related bio-mimetic design for quietness.展开更多
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.展开更多
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.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S）...The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.展开更多
[Objective] The aim was to explore release characteristics of vinyl chlo- ride-vinyl acetate copolymer controlled-release N fertilizer and the effects on minerat nitrogen in soils. [Method] Vinyl chloride-vinyl acetat...[Objective] The aim was to explore release characteristics of vinyl chlo- ride-vinyl acetate copolymer controlled-release N fertilizer and the effects on minerat nitrogen in soils. [Method] Vinyl chloride-vinyl acetate copolymer and hydroxyl-modi- fied VCNAc were taken as coating materials to prepare slow release fertilizer. Nutri- ent release characteristics of VC/VAc slow release fertilizer was evaluated by water immersion method and the effects of VC/VAc slow release fertilizer on mineral ni- trogen were researched by pot experiment. [Result] The release periods of VC-VAc controlled-release urea and hydroxyl-modified VC/VAc coated urea were 60 and 50 d, respectively. Furthermore, the content of ammonium nitrogen reached the peak on the 30th d and the content of nitrate nitrogen reached the peak on the 60th d in soils in treatments with VCNAc and hydroxyl-modified VC/VAc; the content of nitrate nitrogen rose again on the 120th d in the treatment with VC/VAc. In terms of wheat yield, different treatments showed insignificant differences and rice yield in the treatment with VCNAc was significantly higher than that in the treatment with hy- droxyl-modified VCNAc （P〈0.05）. [Conclusion] The release days of slow controlled- release fertilizer vary upon pot experiment method and water immersion method. Slow controlled-release fertilizer is not suitable for monoculture, due to long fertilizer efficiency, but multiple cropping would be optimal for its role to be fully exploited.展开更多
This paper introduces MultiPHydro,an in-house computational solver developed for simulating hydrodynamic and multiphase fluid—body interaction problems,with a specialized focus on multiphase flow dynamics.The solver ...This paper introduces MultiPHydro,an in-house computational solver developed for simulating hydrodynamic and multiphase fluid—body interaction problems,with a specialized focus on multiphase flow dynamics.The solver employs the boundary data immersion method(BDIM)as its core numerical framework for handling fluid—solid interfaces.We briefly outline the governing equations and physical models integrated within MultiPHydro,including weakly-compressible flows,cavitation modeling,and the volume of fluid(VOF)method with piecewise-linear interface reconstruction.The solver’s accuracy and versatility are demonstrated through several numerical benchmarks:single-phase flow past a cylinder shows less than 10%error in vortex shedding frequency and under 4%error in hydrodynamic resistance;cavitating flows around a hydrofoil yield errors below 7%in maximum cavity length;water-entry cases exhibit under 5%error in displacement and velocity;and water-exit simulations predict cavity length within 7.2%deviation.These results confirm the solver’s capability to reliably model complex fluid-body interactions across various regimes.Future developments will focus on refining mathematical models,improving the modeling of phase-interaction mechanisms,and implementing GPU-accelerated parallel algorithms to enhance compatibility with domestically-developed operating systems and deep computing units(DCUs).展开更多
A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the ...A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.展开更多
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.展开更多
The flapping motion of a flexible propulsor near the ground was simulated using the immersed boundary method. The hydrodynamic benefits of the propulsor near the ground were explored by varying the heaving frequency (...The flapping motion of a flexible propulsor near the ground was simulated using the immersed boundary method. The hydrodynamic benefits of the propulsor near the ground were explored by varying the heaving frequency (St) of the leading edge of the flexible propulsor. Propulsion near the ground had some advantages in generating thrust and propelling faster than propulsion away from the ground. The mode analysis and flapping amplitude along the Lagrangian coordinate were examined to analyze the kinematics as a function of the ground proximity (d) and St. The trailing edge amplitude (\(a_\mathrm{tail}\)) and the net thrust (\(\overline{{F}}_x\)) were influenced by St of the flexible propulsor. The vortical structures in the wake were analyzed for different flapping conditions.展开更多
基金Supported by the National Key R&D Program of China(2022YFB3303500).
文摘Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.
基金supported by the National Natural Science Foundation of China(No.52175277)The computations were conducted in Computing Center in Xi'an.The authors thankfully acknowledge these institutions.
文摘The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic performance of a butterfly like model is studied using a solver based on immersed boundary method and adaptive mesh.For aerodynamic performance,adding a"forward-backward-forward"sweeping motion to the forewing makes it more suitable for fast cruising flight,and compared to the situation without forewing sweeping,the drag is reduced by 46%and the lift to drag ratio is increased by 45%.On the contrary,adding a"backward-for ward-backward"sweeping motion to the forewing increases lift and makes it more suitable for climb flight.For downstroke and middle to late upstroke,the forewing sweeping affects the Leading-Edge Vortex(LEV)through two factors:sweeping velocity and forward sweeping angle,and their effects are coupled.A large forward sweeping velocity can enhance the strength of LEV,while a large forward sweeping angle can weaken it.For early upstroke,the forewing sweeping can affect the wake capture mechanism,sweeping backward can enhance it while sweeping forward can weaken it.The findings in this work provide insight into the design of butterfly like Micro AirVehicles(MAVs).
基金partially supported by the National Natural Science Foundation of China(Nos.92271103,12202191)。
文摘The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.
基金joined PI of Westlake University(Grant Nos.041030150118 and 103110556022101)Scientific Research Funding Project of Westlake University(Grant No.2021WUFP017).
文摘The role of hydrodynamic effect in the meeting of multiple fish is a fascinating topic.The interactions of two self-propelled flexi-ble plates swimming in opposite directions horizontally and maintaining a certain lateral distance are numerically simulated using a penalty-immersed boundary method.The effects of the flapping phase and lateral distance on the propulsive performance of two fish meetings are analyzed.Results show that,when two plates meet,if their leading edges diverge laterally,the individual plate can efficiently and rapidly move apart from the other horizontally.If their leading edges converge laterally,the plate motion can be retarded,leading to high energy consumption.Moreover,an increasing lateral distance between two plates significantly weakens the fluid-structure interactions,resulting in an exponential decline in mean cruising speed.A quantitative force analysis based on vortex dynamic theory is performed to gain physics insight into the hydrodynamic interaction mechanism.It is found that lateral separation between the two leading edges enhances the vorticity generation and boundary vorticity flux on the surface of the plate,subsequently reinforcing the thrust effect and increasing horizontal velocity.This study offers insight into the hydro-dynamic mechanisms of the fluid-structure interactions among fish moving toward each other and suggests potential strategies for enhancing the maneuverability of robotic fish in complex environment.
基金Supported by the National Natural Science Foundation of China(52201350,52201394,and 52271301)the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.
文摘The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.
基金the financial support from of the National Basic Research Program of China (Nos. 2012CB619600 and 2011CB012803)
文摘A facile ammonium-dichromate solution immersion method was introduced to synthesize the copperwettable Cr3C2 coating on and inside the carbon-carbon （C/C） preform. The formation mechanism and the microstructures of the Cr3C2 coatings were studied. The contact angle between molten copper and the C/C decreased from 140°to 60°, demonstrating the significant improvement in the wettability. The Cr3C2- coated C/C-Cu composite with only 4.2% porosity and 3.69 gcm^-3 density was manufactured through copper infiltration. As a result, the thermal and electrical conductivity of the modified C/C-Cu increased significantly due to the infiltrated copper. Also the mechanical properties of the composites including both the flexural and compressive strengths were enhanced by over 100%. The modified C/C-Cu composite exhibited lower friction coefficients and wear rates for different load levels than those of the commercial C/Cu composite. These results demonstrate the potential of the modified C/C-Cu material for use in electrical contacts.
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.
基金supported by the National Natural Science Foundation of China(No.51174236)the National Basic Research Program of China(973 Program)(No.2011CB606306)the Opening Project of State Key Laboratory of Porous Metal Materials(No.PMM-SKL-4-2012)
文摘An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.52301380 and 11602277).
文摘A hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simulate the incompressible flow field.The sound field is predicted by discretizing acoustic perturbation equations(APEs)with dispersionrelation-preserving space scheme and low-dispersion and low-dissipation Runge-Kutta time integration.A sharp-interface IBM based on ghost-cell is implemented for present two-step DNS-APE approach to deal with complex moving bodies with Cartesian grids.The present method is validated through simulations of sound generation caused by flow past a rotating cylinder,an oscillating cylinder,and tandem oscillating and stationary cylinders.The sound generated by typical kinds of complicated bio-inspired locomotions,i.e.,flapping flight by wings of varied shapes and collective undulatory swimming in tandem,are investigated using present method.The results demonstrate potential of the hybrid approach in addressing flowinduced sound generation and propagation with complex moving boundaries in a fluid medium,especially for the sound characteristics of bio-mimetic flows,which might shed lights on investigations on bio-acoustics,ethology of complex animal system,and related bio-mimetic design for quietness.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
基金Supported by the Aeronautical Science Foundation of China(20111453012)the National Defense Pre-Research Foundation of China(9140A13040111HK0329)~~
文摘The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.
基金Supported by National Department Public Benefit Research Foundation(201203013)Modern Agricultural Industry Technology System(CARS-11-B-15)+2 种基金IPNI Project(JIANGSU-10)Special Fund for Agro-scientific Research in the Public Interest(201003014-1-2)Jiangsu Agriculture S&T Self-Innovation Project[CX(12)3037]~~
文摘[Objective] The aim was to explore release characteristics of vinyl chlo- ride-vinyl acetate copolymer controlled-release N fertilizer and the effects on minerat nitrogen in soils. [Method] Vinyl chloride-vinyl acetate copolymer and hydroxyl-modi- fied VCNAc were taken as coating materials to prepare slow release fertilizer. Nutri- ent release characteristics of VC/VAc slow release fertilizer was evaluated by water immersion method and the effects of VC/VAc slow release fertilizer on mineral ni- trogen were researched by pot experiment. [Result] The release periods of VC-VAc controlled-release urea and hydroxyl-modified VC/VAc coated urea were 60 and 50 d, respectively. Furthermore, the content of ammonium nitrogen reached the peak on the 30th d and the content of nitrate nitrogen reached the peak on the 60th d in soils in treatments with VCNAc and hydroxyl-modified VC/VAc; the content of nitrate nitrogen rose again on the 120th d in the treatment with VC/VAc. In terms of wheat yield, different treatments showed insignificant differences and rice yield in the treatment with VCNAc was significantly higher than that in the treatment with hy- droxyl-modified VCNAc （P〈0.05）. [Conclusion] The release days of slow controlled- release fertilizer vary upon pot experiment method and water immersion method. Slow controlled-release fertilizer is not suitable for monoculture, due to long fertilizer efficiency, but multiple cropping would be optimal for its role to be fully exploited.
文摘This paper introduces MultiPHydro,an in-house computational solver developed for simulating hydrodynamic and multiphase fluid—body interaction problems,with a specialized focus on multiphase flow dynamics.The solver employs the boundary data immersion method(BDIM)as its core numerical framework for handling fluid—solid interfaces.We briefly outline the governing equations and physical models integrated within MultiPHydro,including weakly-compressible flows,cavitation modeling,and the volume of fluid(VOF)method with piecewise-linear interface reconstruction.The solver’s accuracy and versatility are demonstrated through several numerical benchmarks:single-phase flow past a cylinder shows less than 10%error in vortex shedding frequency and under 4%error in hydrodynamic resistance;cavitating flows around a hydrofoil yield errors below 7%in maximum cavity length;water-entry cases exhibit under 5%error in displacement and velocity;and water-exit simulations predict cavity length within 7.2%deviation.These results confirm the solver’s capability to reliably model complex fluid-body interactions across various regimes.Future developments will focus on refining mathematical models,improving the modeling of phase-interaction mechanisms,and implementing GPU-accelerated parallel algorithms to enhance compatibility with domestically-developed operating systems and deep computing units(DCUs).
基金supported by the National Natural Science Foundation of China (11302238, 11232011. and 11572331)support from the Strategic Priority Research Program (XDB22040104)+1 种基金the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences (QYZDJ-SSW-SYS002)the National Basic Research Program of China (973 Program 2013CB834100: Nonlinear science)
文摘A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.
基金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 the Creative Research Initiatives (Grant 2016-004749) program of the National Research Foundation of Korea (MSIP)
文摘The flapping motion of a flexible propulsor near the ground was simulated using the immersed boundary method. The hydrodynamic benefits of the propulsor near the ground were explored by varying the heaving frequency (St) of the leading edge of the flexible propulsor. Propulsion near the ground had some advantages in generating thrust and propelling faster than propulsion away from the ground. The mode analysis and flapping amplitude along the Lagrangian coordinate were examined to analyze the kinematics as a function of the ground proximity (d) and St. The trailing edge amplitude (\(a_\mathrm{tail}\)) and the net thrust (\(\overline{{F}}_x\)) were influenced by St of the flexible propulsor. The vortical structures in the wake were analyzed for different flapping conditions.
