In this study,we perform particle-resolved simulations of settling spheroidal particles,considering oblate and prolate spheroids and spheres,and investigate the shape effect on the particle dynamics in suspensions wit...In this study,we perform particle-resolved simulations of settling spheroidal particles,considering oblate and prolate spheroids and spheres,and investigate the shape effect on the particle dynamics in suspensions with volume fraction 1%and 5%.We first examine the single-point statistics of the translational and rotational motion of the settling particles.The horizontal velocity has a symmetrical distribution with standard deviation dependent on the particle shape.The greater horizontal velocity fluctuations of the non-spherical particles,compared to that of spheres,are attributed to the horizontal drift of settling spheroids with oblique orientations induced by the fluid-particle and particle-particle interactions.The fluctuation of particle vertical velocity,instead,is skewed under the effect of wake-induced hydrodynamic interactions.Further,we explore the particle pair statistics,which demonstrate the formation of column-like particle micro-structures for the lowest volume fraction considered.This clustering is more pronounced for spheroidal particles than spheres,due to the stronger attractions among vertically-aligned settling spheroids.Moreover,the particle pair statistics are directly related to the collision rate among the dispersed particles.The local accumulation of oblate/prolate spheroids serves as the major mechanism to promote the particle-particle collisions in dilute suspensions.展开更多
The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRE...The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.展开更多
A discontinuous smoothed particle hydrodynamics(DSPH)method considering block contacts is originally developed to model the cracking,frictional slip and large deformation in rock masses,and is verified by theoretical,...A discontinuous smoothed particle hydrodynamics(DSPH)method considering block contacts is originally developed to model the cracking,frictional slip and large deformation in rock masses,and is verified by theoretical,numerical and/or experimental results.In the DSPH method,cracking is realized by breaking the virtual bonds via a pseudo-spring method based on Mohr–Coulomb failure criteria.The damaged particles are instantaneously replaced by discontinuous particles and the contact bond between the original and discontinuous particles is formed to simulate the frictional slip and separation/contraction between fracture surfaces based on the block contact algorithm.The motion of rock blocks and the contact force of discontinuous particles are determined following Newton's second law.The results indicate that the DSPH method precisely captures the cracking,contact formation and complete failure across six numerical benchmark tests.This single smoothed particle hydrodynamics(SPH)framework could significantly improve computational efficiency and is potentially applicable to broad multi-physical rock engineering problems of different scales.展开更多
The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its wid...The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its widespread application in high Reynolds number flow and shear turbulence.To better simulate these flows,this paper proposes the stability-enhanced VPM based on a Reformulated VPM(RVPM)constrained by conservation of angular momentum,integrating a relaxation scheme to suppress the divergence of the vorticity field,and further coupling the Sub-Grid Scale(SGS)model to account for the turbulence dissipation caused by vortex advection and vortex stretching.The validity of the RVPM is confirmed by simulating an isolated vortex ring's evolution.The results also demonstrate that the relaxation scheme of vorticity enhances the numerical stability of the VPM by mitigating the divergence of the vorticity field.The leapfrogging vortex rings simulation demonstrates that the RVPM with the present SGS model can more precisely feature the leapfrog and fusion of vortex rings and has improved numerical stability in high Reynolds number flows.The round turbulent jet simulation confirms that the stability-enhanced VPM can stably simulate shear turbulence and accurately resolve fluctuating components and Reynolds stresses in the turbulence.展开更多
We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework em...We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework employs SPH for fluid dynamics,addressing large deformations,shocks,and plasma behavior,while FEM resolves electromagnetic fields via Maxwell's equations for magnetic vector and electric scalar potentials,ensuring divergence-free conditions and global current density calculations in conductive region.Operator splitting method couples these modules,enabling real-time integration of magnetic,electric,thermal,and fluid fields.Benchmark tests validate the code against analytical solutions and existing models,including blow-by instability simulations that demonstrate the method's accuracy in capturing fluid-magnetic interactions.Designed for 3D applications,SPMHD offers robust scalability across multiprocessor architectures,establishing it as a versatile tool for plasma physics research.展开更多
An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes(SNN) and local-Delaunay-trian...An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes(SNN) and local-Delaunay-triangulation(LDT) tech-niques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h-adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h-adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.展开更多
Hydrothermal method was used to synthesize nanoscale particles of MnZn ferrites. The crystallites were characterized by XRD, TEM and SEM. The effects of the reaction time, temperature and additives on the product were...Hydrothermal method was used to synthesize nanoscale particles of MnZn ferrites. The crystallites were characterized by XRD, TEM and SEM. The effects of the reaction time, temperature and additives on the product were investigated. Crystallization process would be carried out above 160 ℃ for 5 h or more, higher temperature can reduce the reaction time. Additives were used to remove impurities such as Fe 2O 3, ZnMnO 3.10～15 nm pure slightly agglomerated MnZn ferrite crystallites with a narrow grain size distribution were obtained.展开更多
Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this pa...Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this paper, a particle swarm optimization(PSO) method is introduced to solve and control a symplectic multibody system for the first time. It is first combined with the symplectic method to solve problems in uncontrolled and controlled robotic arm systems. It is shown that the results conserve the energy and keep the constraints of the chaotic motion, which demonstrates the efficiency, accuracy, and time-saving ability of the method. To make the system move along the pre-planned path, which is a functional extremum problem, a double-PSO-based instantaneous optimal control is introduced. Examples are performed to test the effectiveness of the double-PSO-based instantaneous optimal control. The results show that the method has high accuracy, a fast convergence speed, and a wide range of applications.All the above verify the immense potential applications of the PSO method in multibody system dynamics.展开更多
Calcined kaolin/TiO2 composite particle material （CK/TCPM） was prepared with TiO2 coating on the surfaces of calcined kaolin particles by the mechano-chemical method. X-ray diffraction （XRD） and scanning electron ...Calcined kaolin/TiO2 composite particle material （CK/TCPM） was prepared with TiO2 coating on the surfaces of calcined kaolin particles by the mechano-chemical method. X-ray diffraction （XRD） and scanning electron microscope （SEM） were used to investigate the microstructures and morphologies, respectively. The mechanism of the mechano-chemical reaction between calcined kaolin and TiO2 was studied by infrared spectra （IR）. The results show that TiO2 coats evenly on the surfaces of calcined kaolin particles by Si-O-Ti and Al-O-Ti bonds on their interfaces. The hiding power and whiteness of CK/TCPM are 17.12 g/m^2 and 95.7%, respectively, presenting its similarity to TiO2 in pigment properties.展开更多
An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-tri...An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-triangulation（LDT） techniques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h- adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h- adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle d...As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.展开更多
This investigation was undertaken to predict the mass gain （MG） of cobalt electroless deposition （ED） on ceramic SiC particles.Response surface methodology （RSM） based on a full factorial design with three ED pa...This investigation was undertaken to predict the mass gain （MG） of cobalt electroless deposition （ED） on ceramic SiC particles.Response surface methodology （RSM） based on a full factorial design with three ED parameters and 30 runs was used to conduct the experiments and to establish a mathematical model by means of Design-Expert software.Three ED parameters considered were pH,bath temperature and ceramic particle morphology.Analysis of variance was applied to validate the predicted model.The results of confirmation analysis by scanning electron microscopy （SEM） show that the developed models are reasonably accurate.The pH is the most effective parameter for the MG.Also,the highest mass gain is obtained for the lowest pH,highest bath temperatures and heat-treated SiC particles.In addition,the developed model shows that the optimal parameters to get a maximum value of mass gain are pH,bath temperature and ceramic particle state of 8,70 ℃ and heat treatment,respectively.展开更多
Discrete element model was developed to simulate the ellipsoidal particles moving in the moving bed. Multi-element model was used to describe a ellipsoidal particle, the contact detection algorithm of ellipsoidal part...Discrete element model was developed to simulate the ellipsoidal particles moving in the moving bed. Multi-element model was used to describe a ellipsoidal particle, the contact detection algorithm of ellipsoidal particle was developed, and both contact force and gravity force were considered in the models. The simulation results were validated by our experiment. Three algorithms for representing an ellipsoidal particle were compared in macro and micro aspects. The results show that there exists big difference in the microscopic parameters such as kinetic energy, rotational kinetic energy, deformation, contact force and collision number which leads to the difference of macroscopic parameters. The relative error in the discharge rate and tracer particle position is the largest between 3-tangent-element representation and experimental results. The flow pattern is similar for the 5-element and 3-intersection representations. The only difference is the discharge rate of 5-element representation is larger than the experimental value and that of the 3-intersection representation has the contrary result. Finally the 3-intersection- element reoresentation is chosen in the simulation due to less comouting time than that of the 5-element renresentation.展开更多
On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is present...On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.展开更多
Internal erosion occurs when fine particles escape from the soil driven by seepage flow,which is considered to be the crucial factor causing the failure of earth structures filled with gravelly soil.The objective of t...Internal erosion occurs when fine particles escape from the soil driven by seepage flow,which is considered to be the crucial factor causing the failure of earth structures filled with gravelly soil.The objective of this paper is to suggest an appropriate method to assess internal erosion potential of gravelly soil.By analyzing the sensitivity of soil material to internal erosion,the variable(Dc15/df85)max and the content of coarse particles(Pc)are selected as the evaluation indexes(Dc15 and df85 are the diameters of 15%mass passing in the coarse component and 85%mass passing in the fine component,respectively).A series of gravelly soils with different particle size distributions are tested for internal erosion by the self-made permeameter.Based on the test results,an evaluation method for the internal erosion of gravelly soil is proposed.Gravelly soil is prone to internal erosion when 60%≤Pc<95%and(Dc15/df85)max≥9.5.The proposed method shows good accuracy in evaluating the internal erosion of 36 soil samples from other studies,which confirms the reliability of the method.The proposed method makes it possible to accurately assess internal erosion of gravelly soil,and an alternative method is provided for engineers to determine whether there is a risk of internal erosion in earth structures consisting of gravelly soil.展开更多
This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combina...This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combination of the vector mechanics and numerical calculations. It models the analyzed domain composed of finite particles. Newton＇s second law is adopted to describe the motions of all particles. A convected material flame and explicit time integration for the solution procedure is also adopted in this method. By using the FPM, there is no need to solve any nonlinear equations, to calculate the stiffness matrix or equilibrium matrix, which is very helpful in the analysis of kinematically indeterminate structures. The basic formulations for the space bar are derived, following its solution procedures for bar assemblies. Three numerical examples are analyzed using the FPM. Results obtained from both the straight pretension cable and the suspension cable assembly show that the FPM can produce a more accurate analysis result. The motion simulation of the four-bar space assembly demonstrates the capability of this method in the analysis ofkinematically indeterminate structures.展开更多
Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular method...Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.展开更多
A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow t...A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.展开更多
The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO...The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO). Different topologies of a multilayer neural network were studied and the optimum architecture was determined. Property data of 350 compounds were used for training the network. To discriminate different substances the molecular structures defined by the concept of the classical group contribution method were given as input variables. The capabilities of the network were tested with 155 substances not considered in the training step. The study shows that the proposed GCM+ANN+PSO method represent an excellent alternative for the estimation of flash points of organic compounds with acceptable accuracy (AARD = 1.8%; AAE = 6.2 K).展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.92252104,12388101,and 12472224).
文摘In this study,we perform particle-resolved simulations of settling spheroidal particles,considering oblate and prolate spheroids and spheres,and investigate the shape effect on the particle dynamics in suspensions with volume fraction 1%and 5%.We first examine the single-point statistics of the translational and rotational motion of the settling particles.The horizontal velocity has a symmetrical distribution with standard deviation dependent on the particle shape.The greater horizontal velocity fluctuations of the non-spherical particles,compared to that of spheres,are attributed to the horizontal drift of settling spheroids with oblique orientations induced by the fluid-particle and particle-particle interactions.The fluctuation of particle vertical velocity,instead,is skewed under the effect of wake-induced hydrodynamic interactions.Further,we explore the particle pair statistics,which demonstrate the formation of column-like particle micro-structures for the lowest volume fraction considered.This clustering is more pronounced for spheroidal particles than spheres,due to the stronger attractions among vertically-aligned settling spheroids.Moreover,the particle pair statistics are directly related to the collision rate among the dispersed particles.The local accumulation of oblate/prolate spheroids serves as the major mechanism to promote the particle-particle collisions in dilute suspensions.
基金funded by the National Natural Science Foundation of China(No.51978048).
文摘The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.
基金financial support from the National Key Research and Development Program of China(Grant No.2019YFC1509702)the Fundamental Research Funds for the Central Universities in Chinathe National Natural Science Foundation of China(Grant No.42377162).
文摘A discontinuous smoothed particle hydrodynamics(DSPH)method considering block contacts is originally developed to model the cracking,frictional slip and large deformation in rock masses,and is verified by theoretical,numerical and/or experimental results.In the DSPH method,cracking is realized by breaking the virtual bonds via a pseudo-spring method based on Mohr–Coulomb failure criteria.The damaged particles are instantaneously replaced by discontinuous particles and the contact bond between the original and discontinuous particles is formed to simulate the frictional slip and separation/contraction between fracture surfaces based on the block contact algorithm.The motion of rock blocks and the contact force of discontinuous particles are determined following Newton's second law.The results indicate that the DSPH method precisely captures the cracking,contact formation and complete failure across six numerical benchmark tests.This single smoothed particle hydrodynamics(SPH)framework could significantly improve computational efficiency and is potentially applicable to broad multi-physical rock engineering problems of different scales.
基金co-supported by the National Natural Science Foundation of China(No.12402272)the Natural Science Basic Research Program of Shaanxi Province,China(No.2024JC-YBQN-0024)the Fundamental Research Funds for the Central Universities,China(No.D5000240030)。
文摘The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its widespread application in high Reynolds number flow and shear turbulence.To better simulate these flows,this paper proposes the stability-enhanced VPM based on a Reformulated VPM(RVPM)constrained by conservation of angular momentum,integrating a relaxation scheme to suppress the divergence of the vorticity field,and further coupling the Sub-Grid Scale(SGS)model to account for the turbulence dissipation caused by vortex advection and vortex stretching.The validity of the RVPM is confirmed by simulating an isolated vortex ring's evolution.The results also demonstrate that the relaxation scheme of vorticity enhances the numerical stability of the VPM by mitigating the divergence of the vorticity field.The leapfrogging vortex rings simulation demonstrates that the RVPM with the present SGS model can more precisely feature the leapfrog and fusion of vortex rings and has improved numerical stability in high Reynolds number flows.The round turbulent jet simulation confirms that the stability-enhanced VPM can stably simulate shear turbulence and accurately resolve fluctuating components and Reynolds stresses in the turbulence.
基金supported by the Major National Science and Technology Infrastructure(No.2208-000000-04-01249628)the Shanghai Science and Technology Commission(No.21DZ1206500)。
文摘We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework employs SPH for fluid dynamics,addressing large deformations,shocks,and plasma behavior,while FEM resolves electromagnetic fields via Maxwell's equations for magnetic vector and electric scalar potentials,ensuring divergence-free conditions and global current density calculations in conductive region.Operator splitting method couples these modules,enabling real-time integration of magnetic,electric,thermal,and fluid fields.Benchmark tests validate the code against analytical solutions and existing models,including blow-by instability simulations that demonstrate the method's accuracy in capturing fluid-magnetic interactions.Designed for 3D applications,SPMHD offers robust scalability across multiprocessor architectures,establishing it as a versatile tool for plasma physics research.
基金Project supported by the National Natural Science Foundation of China (No.10202018)
文摘An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes(SNN) and local-Delaunay-triangulation(LDT) tech-niques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h-adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h-adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.
文摘Hydrothermal method was used to synthesize nanoscale particles of MnZn ferrites. The crystallites were characterized by XRD, TEM and SEM. The effects of the reaction time, temperature and additives on the product were investigated. Crystallization process would be carried out above 160 ℃ for 5 h or more, higher temperature can reduce the reaction time. Additives were used to remove impurities such as Fe 2O 3, ZnMnO 3.10～15 nm pure slightly agglomerated MnZn ferrite crystallites with a narrow grain size distribution were obtained.
基金Project supported by the National Natural Science Foundation of China(Nos.91648101 and11672233)the Northwestern Polytechnical University(NPU)Foundation for Fundamental Research(No.3102017AX008)the National Training Program of Innovation and Entrepreneurship for Undergraduates(No.S201710699033)
文摘Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this paper, a particle swarm optimization(PSO) method is introduced to solve and control a symplectic multibody system for the first time. It is first combined with the symplectic method to solve problems in uncontrolled and controlled robotic arm systems. It is shown that the results conserve the energy and keep the constraints of the chaotic motion, which demonstrates the efficiency, accuracy, and time-saving ability of the method. To make the system move along the pre-planned path, which is a functional extremum problem, a double-PSO-based instantaneous optimal control is introduced. Examples are performed to test the effectiveness of the double-PSO-based instantaneous optimal control. The results show that the method has high accuracy, a fast convergence speed, and a wide range of applications.All the above verify the immense potential applications of the PSO method in multibody system dynamics.
基金Funded by the National Key Technology R&D Program of China(No.2008BAE60B06)Beijing Municipal Science&Technology Commission (No.Z080003032208015)
文摘Calcined kaolin/TiO2 composite particle material （CK/TCPM） was prepared with TiO2 coating on the surfaces of calcined kaolin particles by the mechano-chemical method. X-ray diffraction （XRD） and scanning electron microscope （SEM） were used to investigate the microstructures and morphologies, respectively. The mechanism of the mechano-chemical reaction between calcined kaolin and TiO2 was studied by infrared spectra （IR）. The results show that TiO2 coats evenly on the surfaces of calcined kaolin particles by Si-O-Ti and Al-O-Ti bonds on their interfaces. The hiding power and whiteness of CK/TCPM are 17.12 g/m^2 and 95.7%, respectively, presenting its similarity to TiO2 in pigment properties.
文摘An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-triangulation（LDT） techniques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h- adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h- adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
基金the financial support provided by the National Key Research and Development Program of China(Grant No.2016YFC0800200)the National Natural Science Foundation of China(Grant Nos.51578494 and 51778568)the Fundamental Research Funds for the Central Universities(Grant No.2019QNA4043).
文摘As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.
文摘This investigation was undertaken to predict the mass gain （MG） of cobalt electroless deposition （ED） on ceramic SiC particles.Response surface methodology （RSM） based on a full factorial design with three ED parameters and 30 runs was used to conduct the experiments and to establish a mathematical model by means of Design-Expert software.Three ED parameters considered were pH,bath temperature and ceramic particle morphology.Analysis of variance was applied to validate the predicted model.The results of confirmation analysis by scanning electron microscopy （SEM） show that the developed models are reasonably accurate.The pH is the most effective parameter for the MG.Also,the highest mass gain is obtained for the lowest pH,highest bath temperatures and heat-treated SiC particles.In addition,the developed model shows that the optimal parameters to get a maximum value of mass gain are pH,bath temperature and ceramic particle state of 8,70 ℃ and heat treatment,respectively.
基金Supported by the Major State Basic Research Development Program of China （2011CB201505）, the National Natural Science Foundation of China （50976025） and the Key Proj ect.of Science and Technology of Henan Province （12B610012）.
文摘Discrete element model was developed to simulate the ellipsoidal particles moving in the moving bed. Multi-element model was used to describe a ellipsoidal particle, the contact detection algorithm of ellipsoidal particle was developed, and both contact force and gravity force were considered in the models. The simulation results were validated by our experiment. Three algorithms for representing an ellipsoidal particle were compared in macro and micro aspects. The results show that there exists big difference in the microscopic parameters such as kinetic energy, rotational kinetic energy, deformation, contact force and collision number which leads to the difference of macroscopic parameters. The relative error in the discharge rate and tracer particle position is the largest between 3-tangent-element representation and experimental results. The flow pattern is similar for the 5-element and 3-intersection representations. The only difference is the discharge rate of 5-element representation is larger than the experimental value and that of the 3-intersection representation has the contrary result. Finally the 3-intersection- element reoresentation is chosen in the simulation due to less comouting time than that of the 5-element renresentation.
基金supported by the National Natural Science Foundation of China (Grant No.10871124)the Innovation Program of Shanghai Municipal Education Commission,China (Grant No.09ZZ99)
文摘On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.
基金financially supported by the National Natural Science Foundation of China(Grant No.41790432)Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA20030301)。
文摘Internal erosion occurs when fine particles escape from the soil driven by seepage flow,which is considered to be the crucial factor causing the failure of earth structures filled with gravelly soil.The objective of this paper is to suggest an appropriate method to assess internal erosion potential of gravelly soil.By analyzing the sensitivity of soil material to internal erosion,the variable(Dc15/df85)max and the content of coarse particles(Pc)are selected as the evaluation indexes(Dc15 and df85 are the diameters of 15%mass passing in the coarse component and 85%mass passing in the fine component,respectively).A series of gravelly soils with different particle size distributions are tested for internal erosion by the self-made permeameter.Based on the test results,an evaluation method for the internal erosion of gravelly soil is proposed.Gravelly soil is prone to internal erosion when 60%≤Pc<95%and(Dc15/df85)max≥9.5.The proposed method shows good accuracy in evaluating the internal erosion of 36 soil samples from other studies,which confirms the reliability of the method.The proposed method makes it possible to accurately assess internal erosion of gravelly soil,and an alternative method is provided for engineers to determine whether there is a risk of internal erosion in earth structures consisting of gravelly soil.
基金supported by the National Natural Science Foundation of China (No. 50638050)the National High-Tech R&D (863) Program (No. 2007AA04Z441), China
文摘This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combination of the vector mechanics and numerical calculations. It models the analyzed domain composed of finite particles. Newton＇s second law is adopted to describe the motions of all particles. A convected material flame and explicit time integration for the solution procedure is also adopted in this method. By using the FPM, there is no need to solve any nonlinear equations, to calculate the stiffness matrix or equilibrium matrix, which is very helpful in the analysis of kinematically indeterminate structures. The basic formulations for the space bar are derived, following its solution procedures for bar assemblies. Three numerical examples are analyzed using the FPM. Results obtained from both the straight pretension cable and the suspension cable assembly show that the FPM can produce a more accurate analysis result. The motion simulation of the four-bar space assembly demonstrates the capability of this method in the analysis ofkinematically indeterminate structures.
基金project of “Natural Science Foundation of China, China (Nos. 5187914, 51679123, 51479095)”
文摘Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.
基金This work was supported by the National Natural Science Foundation of China （No. 50275094）.
文摘A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.
文摘The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO). Different topologies of a multilayer neural network were studied and the optimum architecture was determined. Property data of 350 compounds were used for training the network. To discriminate different substances the molecular structures defined by the concept of the classical group contribution method were given as input variables. The capabilities of the network were tested with 155 substances not considered in the training step. The study shows that the proposed GCM+ANN+PSO method represent an excellent alternative for the estimation of flash points of organic compounds with acceptable accuracy (AARD = 1.8%; AAE = 6.2 K).
