In the present study,the nanofliud natural convection is investigated by the energy-conserving dissipative particle dynamics(eDPD)method,where the nanoparticles are considered at the single-particle level.The thermal ...In the present study,the nanofliud natural convection is investigated by the energy-conserving dissipative particle dynamics(eDPD)method,where the nanoparticles are considered at the single-particle level.The thermal expansion coefficientβand the viscosityμof the simulated system containing nanoparticles are calculated and found to be in close alignment with the previous simulation results.The single-particle hydrodynamics in e DPD enables simulations of nanofluid natural convection with higher Rayleigh numbers and greater nanoparticle volume fractions.Additionally,this approach is utilized to simulate the nanoparticle distribution during the enhanced heat transfer process in the nanofluid natural convection.The localized aggregation of nanoparticles enhances the heat transfer performance of the nanofluid under specific Rayleigh numbers and nanoparticles volume fractions.展开更多
Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strat...Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.展开更多
In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SP...In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.展开更多
Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the ...Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.展开更多
A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction pro...A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.展开更多
Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to...Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.展开更多
Natural convection is a heat transfer mechanism driven by temperature or density differences,leading to fluid motion without external influence.It occurs in various natural and engineering phenomena,influencing heat t...Natural convection is a heat transfer mechanism driven by temperature or density differences,leading to fluid motion without external influence.It occurs in various natural and engineering phenomena,influencing heat transfer,climate,and fluid mixing in industrial processes.This work aims to use the Updated Lagrangian Particle Hydrodynamics(ULPH)theory to address natural convection problems.The Navier-Stokes equation is discretized using second-order nonlocal differential operators,allowing a direct solution of the Laplace operator for temperature in the energy equation.Various numerical simulations,including cases such as natural convection in square cavities and two concentric cylinders,were conducted to validate the reliability of the model.The results demonstrate that the proposed model exhibits excellent accuracy and performance,providing a promising and effective numerical approach for natural convection problems.展开更多
The Theory of General Singularity is presented, unifying quantum field theory, general relativity, and the standard model. This theory posits phonons as fundamental excitations in a quantum vacuum, modeled as a Bose-E...The Theory of General Singularity is presented, unifying quantum field theory, general relativity, and the standard model. This theory posits phonons as fundamental excitations in a quantum vacuum, modeled as a Bose-Einstein condensate. Through key equations, the role of phonons as intermediaries between matter, energy, and spacetime geometry is demonstrated. The theory expands Einsteins field equations to differentiate between visible and dark matter, and revises the standard model by incorporating phonons. It addresses dark matter, dark energy, gravity, and phase transitions, while making testable predictions. The theory proposes that singularities, the essence of particles and black holes, are quantum entities ubiquitous in nature, constituting the very essence of elementary particles, seen as micro black holes or quantum fractal structures of spacetime. As the theory is refined with increasing mathematical rigor, it builds upon the foundation of initial physical intuition, connecting the spacetime continuum of general relativity with the hydrodynamics of the quantum vacuum. Inspired by the insights of Tesla and Majorana, who believed that physical intuition justifies the infringement of mathematical rigor in the early stages of theory development, this work aims to advance the understanding of the fundamental laws of the universe and the perception of reality.展开更多
To adapt to the change of aquaculture workshop site,optimize the shape of aquaculture tanks and improve the utilization rate of breeding space,it is necessary to determine the appropriate length width ratio parameters...To adapt to the change of aquaculture workshop site,optimize the shape of aquaculture tanks and improve the utilization rate of breeding space,it is necessary to determine the appropriate length width ratio parameters of aquaculture tanks.In this paper,computational fluid dynamics(CFD)technology is adopted to study the flow field performance of aquaculture tanks with different L/B ratios(L:the length;B:the width,of aquaculture tank)and different jet direction conditions(lengthways jet and widthways jet).A three-dimensional numerical calculation model of turbulence in rounded rectangle aquaculture tanks in dual-diagonal-inlet layout was established.Jet directions are arranged lengthways and widthways,and the water flow velocity,resistance coefficient change,vorticity,etc.are analyzed under two working conditions.Results show that the flow field performance in aquaculture tank decreases with the increase of the L/B ratio.The flow field performed well when L/B was 1.0-1.3,sharply dropped at 1.4-1.6,and poor at 1.7-1.9.The results provided a theoretical basis for the design and optimization in flow field performance of the industrialized circulating aquaculture tanks.展开更多
In this paper,using the computational fluid dynamics based on Euler Lagrange and the commercial software Barracuda VR,the gas-particle hydrodynamics and the erosion of particles on the inner wall and internal componen...In this paper,using the computational fluid dynamics based on Euler Lagrange and the commercial software Barracuda VR,the gas-particle hydrodynamics and the erosion of particles on the inner wall and internal components of the spouted bed in the integrated multi-jet swirling spout-fluidized bed(IMSSFB)are studied.Erosion experiments have obtained the characterization of particle erosion on internal components and verified the relevant numerical models.The results show that:the particle distribution within the IMSSFB is uneven due to the cyclonic effect of the axial swirl vane(ASV),resulting in particle erosion for the ASV being concentrated on one side;when the gas reaches the top,too high an erosion gas velocity leads to gas backflow.As the filling height increases,there is a tendency for the erosion position of the particles on the ASV to expand upwards.However,the effect of increasing gas velocity on the erosion position is insignificant.展开更多
Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for ...Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for risk assessment and mitigation design.The smoothed particle hydrodynamics(SPH)method has been successfully applied to the simulation of two-dimensional(2D)and three-dimensional(3D)flow-like landslides.However,the influence of boundary resistance on the whole process of landslide failure is rarely discussed.In this study,a boundary condition considering friction is proposed and integrated into the SPH method,and its accuracy is verified.Moreover,the Navier-Stokes equation combined with the non-Newtonian fluid rheologymodel was utilized to solve the dynamic behavior of the flow-like landslide.To verify its performance,the Shuicheng landslide event,which occurred in Guizhou,China,was taken as a case study.In the 2D simulation,a sensitivity analysis was conducted,and the results showed that the shearing strength parameters have more influence on the computation accuracy than the coefficient of viscosity.Afterwards,the dynamic characteristics of the landslide,such as the velocity and the impact area,were analyzed in the 3D simulation.The simulation results are in good agreement with the field investigations.The simulation results demonstrate that the SPH method performs well in reproducing the landslide process,and facilitates the analysis of landslide characteristics as well as the affected areas,which provides a scientific basis for conducting the risk assessment and disaster mitigation design.展开更多
In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise a...In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise accuracy in fluid regions such as splashes and surfaces.Attempts to address this problem used variable smoothing lengths.Yet the existing methods are computationally complex and non-efficient,because the smoothing length is typically calculated using iterative optimization.Here,we propose an efficient non-iterative SPH fluid simulation method with variable smoothing length(VSLSPH).VSLSPH correlates the smoothing length to the density change,and adaptively adjusts the smoothing length of particles with high accuracy and low computational cost,enabling large time steps.Our experimental results demonstrate the advantages of the VSLSPH approach in terms of its simulation accuracy and efficiency.展开更多
The implementation of high pressure die casting (HPDC) filling process modeling based on smoothed particle hydrodynamics (SPH) was discussed. A new treatment of inlet boundary was established by discriminating flu...The implementation of high pressure die casting (HPDC) filling process modeling based on smoothed particle hydrodynamics (SPH) was discussed. A new treatment of inlet boundary was established by discriminating fluid particles from inlet particles. The roles of artificial viscosity and moving least squares method in the present model were compared in the handling pressure oscillation. The final model was substantiated by simulating filling process in HPDC in both two and three dimensions. The simulated results from SPH and finite difference method (FDM) were compared with the experiments. The results show the former is in a better agreement with experiments. It demonstrates the efficiency and precision of this SPH model in describing flow pattern in filling process.展开更多
A numerical model of foundry filling process was established based on the smoothed particle hydrodynamics(SPH)method.To mimic the constraints that the solid mold prescribes on the filling fluid,a composite treatment...A numerical model of foundry filling process was established based on the smoothed particle hydrodynamics(SPH)method.To mimic the constraints that the solid mold prescribes on the filling fluid,a composite treatment to the solid boundaries is elaborately designed.On solid boundary surfaces,boundary particles were set,which exert Lennard-Jones force on approaching fluid particles;inside the solid mold,ghost particles were arranged to complete the compact domain of near-boundary fluid particles.Water analog experiments were conducted in parallel with the model simulations.Very good agreement between experimental and simulation results demonstrates the success of model development.展开更多
Both porosity ( φ ) and permeability ( k ) of the weathered elution deposited rare earth ores are basic hydrodynamic parameters for RE leaching. The relationship between k and φ of two typical rare earth ores of Sou...Both porosity ( φ ) and permeability ( k ) of the weathered elution deposited rare earth ores are basic hydrodynamic parameters for RE leaching. The relationship between k and φ of two typical rare earth ores of South China in the packed bed was investigated by measuring the flow ( Q ) under various leaching pressure difference (Δ p ). The experimental results show that the relationship between k and φ is unique, moreover the relationship between Q and Δ p is in accord with the Darcy’s law. The effects of the type of ores, the leaching reagents and its concentration, the granule ore size on the leaching permeability have also been investigated. It is demonstrated that k H (for heavy RE ore, k H=35.98?mm 2)> k M H (for middle heavy RE ore, k M H =28.50?mm 2), whereas k (NH 4NO 3)> k (NH 4Cl)> k [(NH 4) 2SO 4], and the k value increases with increasing leaching reagents concentration and granule ore size( k 0.60~0.75?mm =99.96?mm 2, k 0.125~0.60?mm =11.83?mm 2, k 0.074~0.125?mm =0.84?mm 2). [展开更多
To improve the efficiency of a CycloBio fluidized sand bed(CB FSB) in removal of dissolved wastes in recirculating aquaculture systems, the hydrodynamics of solid-liquid flow was investigated using computational fluid...To improve the efficiency of a CycloBio fluidized sand bed(CB FSB) in removal of dissolved wastes in recirculating aquaculture systems, the hydrodynamics of solid-liquid flow was investigated using computational fluid dynamics(CFD) modeling tools. The dynamic characteristics of silica sand within the CB FSB were determined using three-dimensional, unsteady-state simulations with the granular Eulerian multiphase approach and the RNG k-ε turbulence model, and the simulation results were validated using available lab-scale measurements. The bed expansion of CB FSB increased with the increase in water inflow rate in numerical simulations. Upon validation, the simulation involving 0.55 mm particles, the Gidaspow correlation for drag coefficient model and the Syamlal-O'Brien correlation for kinetic granular viscosity showed the closest match to the experimental results. The volume fraction of numerical simulations peaked as the wall was approached. The hydrodynamics of a pilot-scale CB FSB was simulated in order to predict the range of water flow to avoid the silica sand overflowing. The numerical simulations were in agreement with the experimental results qualitatively and quantitatively, and thus can be used to study the hydrodynamics of solid-liquid multiphase flow in CB FSB, which is of importance to the design, optimization, and amplification of CB FSBs.展开更多
Radius of ceramic cone can largely contribute into final solution of analytic models of penetration into ceramic/metal targets.In the present research,a modified model based on radius of ceramic cone was presented for...Radius of ceramic cone can largely contribute into final solution of analytic models of penetration into ceramic/metal targets.In the present research,a modified model based on radius of ceramic cone was presented for ceramic/aluminum targets.In order to investigate and evaluate accuracy of the presented analytic model,obtained results were compared against the results of the Florence’s analytic model and also against numerical modeling results.The phenomenon of impact onto ceramic/aluminum composites were modeled using smoothed particle hydrodynamics(SPH)implemented utilizing ABAQUS Software.Results indicated that,with increasing initial velocity and ceramic thickness and decreasing support layer thickness,the radius of ceramic cone decreases;this ends up increasing residual velocity of the projectile and penetration time and extending the area across which the pressure is distributed.These findings indicate enhanced levels of target energy absorption and the required energy for bending and tensioning the target.As such,it can be observed that,at the same thickness and areal density,the ceramic target has its efficiency enhanced with increasing ceramic thickness and decreasing the support layer thickness.Finally,the results revealed that the associated data with SPH confirm the modified analytic model at higher accuracy than the Florence’s analytic model.展开更多
Over the past 20 years, experimental analyses of the biomechanics of locomotion in fishes have generated a number of key findings that are relevant to the construction of biomimetic fish robots. In this paper, we pres...Over the past 20 years, experimental analyses of the biomechanics of locomotion in fishes have generated a number of key findings that are relevant to the construction of biomimetic fish robots. In this paper, we present 16 results from recent experimental research on the mechanics, kinematics, fluid dynamics, and control of fish locomotion that summarize recent work on fish biomechanics. The findings and principles that have emerged from biomechanical studies of fish locomotion provide important insights into the functional design of fishes and suggest specific design features relevant to construction of robotic fish-inspired vehicles that underlie the high locomotor performance exhibited by fishes.展开更多
Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid con...Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid concentrations(0–30 vol%). Radial profiles of time-averaged gas holdup, axial liquid velocity, and turbulent kinetic energy were measured by using in-house developed conductivity probes and Pavlov tubes. Effects of column diameter, superficial gas velocity, and solid concentration were investigated in a wide range of operating conditions. Experimental results indicated that the average gas holdup remarkably increases with superficial gas velocity, and the radial profiles of investigated flow properties become steeper at high superficial gas velocities. The axial liquid velocities significantly increase with the growth of the column size, whereas the gas holdup was slightly affected. The presence of solid in bubble columns would inhibit the breakage of bubbles, which results in an increase in bubble rise velocity and a decrease in gas holdup, but time-averaged axial liquid velocities remain almost the same as that of the hollow column. Furthermore, a 2-D axisymmetric k–ε model was used to simulate heterogeneous bubbly flow using commercial code FLUENT 6.2. The lateral lift force and the turbulent diffusion force were introduced for the determination of gas holdup profiles and the effects of solid concentration were considered as the variation of average bubble diameter in the model. Results predicted by the CFD simulation showed good agreement with experimental data.展开更多
Numerical simulation of the morphology of a droplet deposited on a solid surface requires an efficient description of the three-phase contact line. In this study, a simple method of implementing the contact angle is p...Numerical simulation of the morphology of a droplet deposited on a solid surface requires an efficient description of the three-phase contact line. In this study, a simple method of implementing the contact angle is proposed, combined with a robust smoothed particle hydrodynamics multiphase algorithm (Zhang 2015). The first step of the method is the creation of the virtual liquid-gas interface across the solid surface by means of dummy particles, thus the calculated surface tension near the triple point serves to automatically modulate the dynarnic contact line towards the equilibrium state. We simulate the evolution process of initially square liquid lumps on fiat and curved surfaces. The predictions of droplet profiles are in good agreement with the analytical solutions provided that the macroscopic contact angle is accurately implemented. Compared to the normal correction method, the present method is straightforward without the need to manually alter the normal vectors. This study presents a robust algorithm capable of capturing the physics of the static welling. It may hold great potentials in bio-inspired superhydrophobic surfaces, oil displacement, microfluidics, ore floatation, etc.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11872283 and 2002212)the Sailing Program of Shanghai,China(No.20YF1432800)。
文摘In the present study,the nanofliud natural convection is investigated by the energy-conserving dissipative particle dynamics(eDPD)method,where the nanoparticles are considered at the single-particle level.The thermal expansion coefficientβand the viscosityμof the simulated system containing nanoparticles are calculated and found to be in close alignment with the previous simulation results.The single-particle hydrodynamics in e DPD enables simulations of nanofluid natural convection with higher Rayleigh numbers and greater nanoparticle volume fractions.Additionally,this approach is utilized to simulate the nanoparticle distribution during the enhanced heat transfer process in the nanofluid natural convection.The localized aggregation of nanoparticles enhances the heat transfer performance of the nanofluid under specific Rayleigh numbers and nanoparticles volume fractions.
基金funded by the National Key Research and Development Program of China(Grant No.2023YFC3008300,Grant No.2019YFC1509702)the National Natural Science Foundation of China(Grant No.42172296).
文摘Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.
基金supported by the National Natural Science Foundation of China(Grant No.52201323).
文摘In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.
基金Project supported by the National Natural Science Foundation of China(No.52109068)the Water Conservancy Technology Project of Jiangsu Province of China(No.2022060)。
文摘Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.
基金open foundation of the Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanicsthe Open Foundation of Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment.
文摘A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.
基金the National Natural Science Foundation of China(No.52271316)the Guangdong Basic and Applied Basic Research Foundation(No.2023A1515030262).
文摘Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.
基金support from the National Natural Science Foundations of China(Nos.11972267 and 11802214)the Open Foundation of the Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics and the Open Foundation of Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment.
文摘Natural convection is a heat transfer mechanism driven by temperature or density differences,leading to fluid motion without external influence.It occurs in various natural and engineering phenomena,influencing heat transfer,climate,and fluid mixing in industrial processes.This work aims to use the Updated Lagrangian Particle Hydrodynamics(ULPH)theory to address natural convection problems.The Navier-Stokes equation is discretized using second-order nonlocal differential operators,allowing a direct solution of the Laplace operator for temperature in the energy equation.Various numerical simulations,including cases such as natural convection in square cavities and two concentric cylinders,were conducted to validate the reliability of the model.The results demonstrate that the proposed model exhibits excellent accuracy and performance,providing a promising and effective numerical approach for natural convection problems.
文摘The Theory of General Singularity is presented, unifying quantum field theory, general relativity, and the standard model. This theory posits phonons as fundamental excitations in a quantum vacuum, modeled as a Bose-Einstein condensate. Through key equations, the role of phonons as intermediaries between matter, energy, and spacetime geometry is demonstrated. The theory expands Einsteins field equations to differentiate between visible and dark matter, and revises the standard model by incorporating phonons. It addresses dark matter, dark energy, gravity, and phase transitions, while making testable predictions. The theory proposes that singularities, the essence of particles and black holes, are quantum entities ubiquitous in nature, constituting the very essence of elementary particles, seen as micro black holes or quantum fractal structures of spacetime. As the theory is refined with increasing mathematical rigor, it builds upon the foundation of initial physical intuition, connecting the spacetime continuum of general relativity with the hydrodynamics of the quantum vacuum. Inspired by the insights of Tesla and Majorana, who believed that physical intuition justifies the infringement of mathematical rigor in the early stages of theory development, this work aims to advance the understanding of the fundamental laws of the universe and the perception of reality.
基金Supported by the National Natural Science Foundation of China(No.31872609)the Innovation Support Program for High-level Talents of Dalian City(No.2019RD12)+1 种基金the Key Research Project of Liaoning Provincial Department of Education in 2022(No.LJKZZ 20220091)the earmarked fund for CARS-49(CARS-49)。
文摘To adapt to the change of aquaculture workshop site,optimize the shape of aquaculture tanks and improve the utilization rate of breeding space,it is necessary to determine the appropriate length width ratio parameters of aquaculture tanks.In this paper,computational fluid dynamics(CFD)technology is adopted to study the flow field performance of aquaculture tanks with different L/B ratios(L:the length;B:the width,of aquaculture tank)and different jet direction conditions(lengthways jet and widthways jet).A three-dimensional numerical calculation model of turbulence in rounded rectangle aquaculture tanks in dual-diagonal-inlet layout was established.Jet directions are arranged lengthways and widthways,and the water flow velocity,resistance coefficient change,vorticity,etc.are analyzed under two working conditions.Results show that the flow field performance in aquaculture tank decreases with the increase of the L/B ratio.The flow field performed well when L/B was 1.0-1.3,sharply dropped at 1.4-1.6,and poor at 1.7-1.9.The results provided a theoretical basis for the design and optimization in flow field performance of the industrialized circulating aquaculture tanks.
基金supported by the National Natural Science Foundation of China(22178286)Shaanxi Qin Chuangyuan“scientist and engineer”team construction project(2022KXJ-041)。
文摘In this paper,using the computational fluid dynamics based on Euler Lagrange and the commercial software Barracuda VR,the gas-particle hydrodynamics and the erosion of particles on the inner wall and internal components of the spouted bed in the integrated multi-jet swirling spout-fluidized bed(IMSSFB)are studied.Erosion experiments have obtained the characterization of particle erosion on internal components and verified the relevant numerical models.The results show that:the particle distribution within the IMSSFB is uneven due to the cyclonic effect of the axial swirl vane(ASV),resulting in particle erosion for the ASV being concentrated on one side;when the gas reaches the top,too high an erosion gas velocity leads to gas backflow.As the filling height increases,there is a tendency for the erosion position of the particles on the ASV to expand upwards.However,the effect of increasing gas velocity on the erosion position is insignificant.
文摘Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for risk assessment and mitigation design.The smoothed particle hydrodynamics(SPH)method has been successfully applied to the simulation of two-dimensional(2D)and three-dimensional(3D)flow-like landslides.However,the influence of boundary resistance on the whole process of landslide failure is rarely discussed.In this study,a boundary condition considering friction is proposed and integrated into the SPH method,and its accuracy is verified.Moreover,the Navier-Stokes equation combined with the non-Newtonian fluid rheologymodel was utilized to solve the dynamic behavior of the flow-like landslide.To verify its performance,the Shuicheng landslide event,which occurred in Guizhou,China,was taken as a case study.In the 2D simulation,a sensitivity analysis was conducted,and the results showed that the shearing strength parameters have more influence on the computation accuracy than the coefficient of viscosity.Afterwards,the dynamic characteristics of the landslide,such as the velocity and the impact area,were analyzed in the 3D simulation.The simulation results are in good agreement with the field investigations.The simulation results demonstrate that the SPH method performs well in reproducing the landslide process,and facilitates the analysis of landslide characteristics as well as the affected areas,which provides a scientific basis for conducting the risk assessment and disaster mitigation design.
基金the Key Program of National Natural Science Foundation of China,No.62237001National Natural Science Foundation for Excellent Young Scholars,No.6212200101+2 种基金National Natural Science Foundation for General Program,Nos.62176066 and 61976052Guangdong Provincial Science and Technology Innovation Strategy Fund,No.2019B121203012and Guangzhou Science and Technology Plan,No.202007040005.
文摘In classical smoothed particle hydrodynamics(SPH)fluid simulation approaches,the smoothing length of Lagrangian particles is typically constant.One major disadvantage is the lack of adaptiveness,which may compromise accuracy in fluid regions such as splashes and surfaces.Attempts to address this problem used variable smoothing lengths.Yet the existing methods are computationally complex and non-efficient,because the smoothing length is typically calculated using iterative optimization.Here,we propose an efficient non-iterative SPH fluid simulation method with variable smoothing length(VSLSPH).VSLSPH correlates the smoothing length to the density change,and adaptively adjusts the smoothing length of particles with high accuracy and low computational cost,enabling large time steps.Our experimental results demonstrate the advantages of the VSLSPH approach in terms of its simulation accuracy and efficiency.
基金Project (2009Z001) supported by the Important Item in Guangdong-Hong Kong Key Project, ChinaProject (2010B090400297) supported by the Cooperation Project in Industry, Education and Research of Guangdong Province and Ministry of Education of China
文摘The implementation of high pressure die casting (HPDC) filling process modeling based on smoothed particle hydrodynamics (SPH) was discussed. A new treatment of inlet boundary was established by discriminating fluid particles from inlet particles. The roles of artificial viscosity and moving least squares method in the present model were compared in the handling pressure oscillation. The final model was substantiated by simulating filling process in HPDC in both two and three dimensions. The simulated results from SPH and finite difference method (FDM) were compared with the experiments. The results show the former is in a better agreement with experiments. It demonstrates the efficiency and precision of this SPH model in describing flow pattern in filling process.
基金Project(2011006B)supported by the Open Project of National Engineering Research Center of Near-Shape Forming for Metallic Materials,ChinaProject(FJ)supported by the CAS"100 talents"Plan
文摘A numerical model of foundry filling process was established based on the smoothed particle hydrodynamics(SPH)method.To mimic the constraints that the solid mold prescribes on the filling fluid,a composite treatment to the solid boundaries is elaborately designed.On solid boundary surfaces,boundary particles were set,which exert Lennard-Jones force on approaching fluid particles;inside the solid mold,ghost particles were arranged to complete the compact domain of near-boundary fluid particles.Water analog experiments were conducted in parallel with the model simulations.Very good agreement between experimental and simulation results demonstrates the success of model development.
文摘Both porosity ( φ ) and permeability ( k ) of the weathered elution deposited rare earth ores are basic hydrodynamic parameters for RE leaching. The relationship between k and φ of two typical rare earth ores of South China in the packed bed was investigated by measuring the flow ( Q ) under various leaching pressure difference (Δ p ). The experimental results show that the relationship between k and φ is unique, moreover the relationship between Q and Δ p is in accord with the Darcy’s law. The effects of the type of ores, the leaching reagents and its concentration, the granule ore size on the leaching permeability have also been investigated. It is demonstrated that k H (for heavy RE ore, k H=35.98?mm 2)> k M H (for middle heavy RE ore, k M H =28.50?mm 2), whereas k (NH 4NO 3)> k (NH 4Cl)> k [(NH 4) 2SO 4], and the k value increases with increasing leaching reagents concentration and granule ore size( k 0.60~0.75?mm =99.96?mm 2, k 0.125~0.60?mm =11.83?mm 2, k 0.074~0.125?mm =0.84?mm 2). [
基金supported by the 12th Five-year National Technology Support Project(2011BAD13B04)
文摘To improve the efficiency of a CycloBio fluidized sand bed(CB FSB) in removal of dissolved wastes in recirculating aquaculture systems, the hydrodynamics of solid-liquid flow was investigated using computational fluid dynamics(CFD) modeling tools. The dynamic characteristics of silica sand within the CB FSB were determined using three-dimensional, unsteady-state simulations with the granular Eulerian multiphase approach and the RNG k-ε turbulence model, and the simulation results were validated using available lab-scale measurements. The bed expansion of CB FSB increased with the increase in water inflow rate in numerical simulations. Upon validation, the simulation involving 0.55 mm particles, the Gidaspow correlation for drag coefficient model and the Syamlal-O'Brien correlation for kinetic granular viscosity showed the closest match to the experimental results. The volume fraction of numerical simulations peaked as the wall was approached. The hydrodynamics of a pilot-scale CB FSB was simulated in order to predict the range of water flow to avoid the silica sand overflowing. The numerical simulations were in agreement with the experimental results qualitatively and quantitatively, and thus can be used to study the hydrodynamics of solid-liquid multiphase flow in CB FSB, which is of importance to the design, optimization, and amplification of CB FSBs.
文摘Radius of ceramic cone can largely contribute into final solution of analytic models of penetration into ceramic/metal targets.In the present research,a modified model based on radius of ceramic cone was presented for ceramic/aluminum targets.In order to investigate and evaluate accuracy of the presented analytic model,obtained results were compared against the results of the Florence’s analytic model and also against numerical modeling results.The phenomenon of impact onto ceramic/aluminum composites were modeled using smoothed particle hydrodynamics(SPH)implemented utilizing ABAQUS Software.Results indicated that,with increasing initial velocity and ceramic thickness and decreasing support layer thickness,the radius of ceramic cone decreases;this ends up increasing residual velocity of the projectile and penetration time and extending the area across which the pressure is distributed.These findings indicate enhanced levels of target energy absorption and the required energy for bending and tensioning the target.As such,it can be observed that,at the same thickness and areal density,the ceramic target has its efficiency enhanced with increasing ceramic thickness and decreasing the support layer thickness.Finally,the results revealed that the associated data with SPH confirm the modified analytic model at higher accuracy than the Florence’s analytic model.
基金This work was supported in part by an Office of Naval Research (USA) MURI grantgrants from the National Science Foundation to George Lauder
文摘Over the past 20 years, experimental analyses of the biomechanics of locomotion in fishes have generated a number of key findings that are relevant to the construction of biomimetic fish robots. In this paper, we present 16 results from recent experimental research on the mechanics, kinematics, fluid dynamics, and control of fish locomotion that summarize recent work on fish biomechanics. The findings and principles that have emerged from biomechanical studies of fish locomotion provide important insights into the functional design of fishes and suggest specific design features relevant to construction of robotic fish-inspired vehicles that underlie the high locomotor performance exhibited by fishes.
基金Supported by the National High Technology Research and Development Program of China(2011AA05A205)the National Natural Science Foundation of China(U1162125,U1361112)
文摘Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid concentrations(0–30 vol%). Radial profiles of time-averaged gas holdup, axial liquid velocity, and turbulent kinetic energy were measured by using in-house developed conductivity probes and Pavlov tubes. Effects of column diameter, superficial gas velocity, and solid concentration were investigated in a wide range of operating conditions. Experimental results indicated that the average gas holdup remarkably increases with superficial gas velocity, and the radial profiles of investigated flow properties become steeper at high superficial gas velocities. The axial liquid velocities significantly increase with the growth of the column size, whereas the gas holdup was slightly affected. The presence of solid in bubble columns would inhibit the breakage of bubbles, which results in an increase in bubble rise velocity and a decrease in gas holdup, but time-averaged axial liquid velocities remain almost the same as that of the hollow column. Furthermore, a 2-D axisymmetric k–ε model was used to simulate heterogeneous bubbly flow using commercial code FLUENT 6.2. The lateral lift force and the turbulent diffusion force were introduced for the determination of gas holdup profiles and the effects of solid concentration were considered as the variation of average bubble diameter in the model. Results predicted by the CFD simulation showed good agreement with experimental data.
基金by the National NaturalScience Foundation of China (Grants 11672335 and 11611530541)China Postdoctoral Science Foundation (Grant 2017M622307)+2 种基金ShandongNatural Science Foundation (Grant ZR201709210320)FundamentalResearch Funds for the Central Universities (Grant18CX0213A)the Endeavour Australia Cheung Kong ResearchFellowship Scholarship from the Australian government.
文摘Numerical simulation of the morphology of a droplet deposited on a solid surface requires an efficient description of the three-phase contact line. In this study, a simple method of implementing the contact angle is proposed, combined with a robust smoothed particle hydrodynamics multiphase algorithm (Zhang 2015). The first step of the method is the creation of the virtual liquid-gas interface across the solid surface by means of dummy particles, thus the calculated surface tension near the triple point serves to automatically modulate the dynarnic contact line towards the equilibrium state. We simulate the evolution process of initially square liquid lumps on fiat and curved surfaces. The predictions of droplet profiles are in good agreement with the analytical solutions provided that the macroscopic contact angle is accurately implemented. Compared to the normal correction method, the present method is straightforward without the need to manually alter the normal vectors. This study presents a robust algorithm capable of capturing the physics of the static welling. It may hold great potentials in bio-inspired superhydrophobic surfaces, oil displacement, microfluidics, ore floatation, etc.