An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

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.

Simulation of nanofluid natural convection based on single-particle hydrodynamics in energy-conserving dissipative particle dynamics(eDPD)

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.

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

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.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

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.

Smoothed-Particle Hydrodynamics Simulation of Ship Motion and Tank Sloshing under the Effect of Regular Waves

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.

Updated Lagrangian Particle Hydrodynamics (ULPH)Modeling of Natural Convection Problems

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.

On the Vacuum Hydrodynamics of Moving Bodies—The Theory of General Singularity

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.

A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

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.

Application of Smoothed Particle Hydrodynamics(SPH)for the Simulation of Flow-Like Landslides on 3D Terrains

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.

An efficient non-iterative smoothed particle hydrodynamics fluid simulation method with variable smoothing length

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.

Numerical and experimental study on the particle erosion and gas–particle hydrodynamics in an integral multi-jet swirling spout-fluidized bed

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.

Simulation of Natural Convection Flow with Magneto-Hydrodynamics in a Wavy Top Enclosure with a Semi-Circular Heater

Natural convection flow in enclosure has different applications such as room ventilation, heat exchangers, the cooling system of a building etc. The Finite-Element method based on the Galerkin weighted residual approach is used to solve two-dimensional governing mass, momentum and energy-equations for natural convection flow in the presence of a magnetic field on a roof top with semi-circular heater. In the enclosure the horizontal lower wall was heated, the vertical two walls were adiabatic, inside the semi-circular heater, the wavy top wall cooled. The parameters Rayleigh number, Hartmann number and Prandtl number are considered. The effects of the Hartmann number and Rayleigh number on the streamlines, isotherms, velocity profiles and average Nusselt number are examined graphically. The local Nusselt number and the average Nusselt number of the heated portion of the enclosure with the semi-circular heater are presented in this paper. Finally, for the validation of the existing work, the current results are compared with published results and the auspicious agreement is achieved.

Simulation of mould filling process using smoothed particle hydrodynamics

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.

Smoothed particle hydrodynamics modeling and simulation of foundry 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 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.

三维Electro-Hydrodynamics方程组光滑解的延拓准则

研究了三维Electro-Hydrodynamics方程组局部光滑解的延拓准则,利用负指标Besov空间中的对数Sobolev不等式和Kato-Ponce交换子估计方法,在Besov空间中获得了只涉及速度场水平分量的最优延拓准则,改进和推广了已有关于Electro-Hydrodynamics方程组正则性准则的部分结论,是Electro-Hydrodynamics方程组正则性理论的有益补充。

Leaching hydrodynamics of weathered elution-deposited rare earth ore

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). [

Application of CFD Modeling to Hydrodynamics of CycloBio Fluidized Sand Bed in Recirculating Aquaculture Systems

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.

Learning from Fish: Kinematics and Experimental Hydrodynamics for Roboticists

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.

Hydrodynamics and Mass Transfer in a Modified Three-phase Airlift Loop Reactor

A modified internal-loop airlif reactor （MIALR） with a continuous slurry phase was studied to investigate the local hydrodynamic characteristics, including gas holdup, bubble size, bubble rise velocity and local mass transfer properties. Based on the analysis of geometrical construction and fluid properties of gas and slurry, MIALR was divided into six flow regions. In these flow regions, the local hydrodynamic characteristics were investigated over a wide range of operating variables. Furthermore, a new method was developed to measure the dissolved oxygen concentration. The volumetric mass-transfer coefficient in six flow regions was also calculated for comparison.

Experimental and numerical investigations of scale-up effects on the hydrodynamics of slurry bubble columns

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.