Particle Contact Research Articles

Ionic Polymerization-Based Synthesis of Bioinspired Adhesive Hydrogel Microparticles with Tunable Morphologies from Microfluidics.

Shape-anisotropic hydrogel microparticles have attracted considerable attention for drug-delivery applications. Particularly, nonspherical hydrogel microcarriers with enhanced adhesive and circulatory abilities have demonstrated value in gastrointestinal drug administration. Herein, inspired by the structures of natural suckers, we demonstrate an ionic polymerization-based production of calcium (Ca)-alginate microparticles with tunable shapes from Janus emulsion for the first time. Monodispersed Janus droplets composed of sodium alginate and nongelable segments were generated using a coflow droplet generator. The interfacial curvatures, sizes, and production frequencies of Janus droplets can be flexibly controlled by varying the flow conditions and surfactant concentrations in the multiphase system. Janus droplets were ionically solidified on a chip, and hydrogel beads of different shapes were obtained. The in vitro and in vivo adhesion abilities of the hydrogel beads to the mouse colon were investigated. The anisotropic beads showed prominent adhesive properties compared with the spherical particles owing to their sticky hydrogel components and unique shapes. Finally, a novel computational fluid dynamics and discrete element method (CFD-DEM) coupling simulation was used to evaluate particle migration and contact forces theoretically. This review presents a simple strategy to synthesize Ca-alginate particles with tunable structures that could be ideal materials for constructing gastrointestinal drug delivery systems.

Discrete Element Method Simulation of Particulate Material Fracture Behavior on a Stretchable Single Filter Fiber with Additional Gas Flow

This study presents a comprehensive discrete element method (DEM) simulation approach for the stretching of a filter fiber with a separated polydisperse particle structure on top. For a realistic interaction between the fiber surface and the particles, the original surface of the polymer fiber was projected onto the surface of the fiber cylinder using surface imaging technologies (atomic force microscopy (AFM) and white-light interferometry). In addition, the adhesive forces between particle–fiber and particle–particle contacts were calibrated in the DEM domain using values from self-conducted AFM measurements. Fiber stretching was implemented by the linear motion of small periodic fiber elements. Discretization problems were resolved through studying the stretching of a fiber segment at the size of 8 mm. A critical fiber element length was discovered to be ≈100 μm for minimizing discretization dependencies during the cracking of the particle structure. The number and density of particle–particle contacts within the particle loading on the fiber were obtained at two different elongation rates. Effects such as densification of the particulate structure and increased detachment due to additional air flow were demonstrated.

Discrete element analysis of the dynamic behavior of ballast track-subgrade system under train dynamic load

To study the dynamic response law of the ballast track-subgrade coupling system under the action of train dynamic load, considering the real shape of ballast particles, a refined discrete element model of the rail-sleeper-ballast bed-subgrade coupling system was established. The established model can realistically simulate the coupling mechanism between irregular ballast particles, between track panel and ballast bed, and between ballast bed and subgrade under train dynamic load. The validity of the established model is verified using on-site measurement results. This study, from a microscopic perspective, focuses on particle contact forces, vibrations, and frictional energy consumption in the rail-sleeper-ballast bed-subgrade coupling system under train dynamic load and investigates the influence of train speed and axle weight. Results show that an increase in train speed and axle weight enhances the number and probability distribution of strong force chains between ballast particles, increasing the risk of particle breakage and deterioration. The influence range and depth in the ballast bed-subgrade system increase with the train's speed and axle weight. Due to the granular nature of the ballast bed, the vibration acceleration distribution of ballast particles under train dynamic loads is relatively random. Overall, the vibration level of the ballast bed increases with train speed and axle weight. The frictional energy consumption of ballast particles rapidly grows as the train's bogie load passes, and the effect of train axle weight on the frictional energy consumption of ballast particles is more significant than the train's speed.

Molecular dynamics simulations of head-on low-velocity collisions between particles.

The particle contact model is important for powder simulations. Although several contact models have been proposed, their validity has not yet been well established. Therefore, we perform molecular dynamics (MD) simulations to clarify the particle interaction. We simulate head-on collisions of two particles with impact velocities less than a few percent of the sound velocity to investigate the dependence of the interparticle force and the coefficient of restitution on the impact velocity and particle radius. In this study, we treat particles with a radius of 10-100nm and perform simulations with up to 0.2 billion atoms. We find that the interparticle force exhibits hysteresis between the loading and unloading phases. Larger impact velocities result in strong hysteresis and plastic deformation. For all impact velocities and particle radii, the coefficient of restitution is smaller than that given by the Johnson-Kendall-Robert theory, which is a contact model that gives the force between elastic spherical particles. A contact model of inelastic particles cannot reproduce our MD simulations. In particular, the coefficient of restitution is significantly reduced when the impact velocity exceeds a certain value. This significant energy dissipation cannot be explained even by the contact models including plastic deformation. We also find that the coefficient of restitution increases with increasing particle radius. We also find that the previous contact models including plastic deformation cannot explain the strong energy dissipation obtained in our MD simulations, although they agree with the MD results for very low impact velocities. Accordingly, we have constructed a new dissipative contact model in which the dissipative force increases with the stress generated by collisions. The new stress-dependent model successfully reproduces our MD results over a wider range of impact velocities than the conventional models do. In addition, we proposed another, simpler, dissipative contact model that can also reproduce the MD results.

Improving the thermal-mechanical performance of bio-treated backfill materials by addition of magnetic iron oxide nanoparticles (nano-Fe3O4)

The thermal conductivity of backfill materials directly affects the heat transfer efficiency between energy geo-structures and the surrounding stratum. Microbially induced carbonate precipitation (MICP) possesses great potential for improving the thermal conductivity of backfill materials. Magnetic iron oxide nanoparticles (i.e., nano-Fe3O4) have been proven to enhance bacterial biochemical activity by altering the permeability of bacterial biofilms, thus potentially improving the MICP process. It was supposed to enhance the thermal conductivity of backfill materials, allowing for applying energy geo-structures in arid environments. In this study, MICP in a solution environment was conducted to analyze bacterial urease activity and morphology of precipitation at varying nano-Fe3O4 contents. Additionally, sand columns treated with MICP and different nano-Fe3O4 contents were performed to obtain the thermal conductivity and unconfined compressive strength (UCS) through the transient plane source (TPS) method and uniaxial compression (UC) experiment. The mineral type, precipitation morphology, and microstructure were identified using scanning electron microscopy (SEM) and X-ray diffraction (XRD). The mechanism of the effect of nano-Fe3O4 on bacterial urease activity and thermal-mechanical behaviors was also discussed. The results indicated that the nano-Fe3O4 could enhance bacterial urease activity and promote vaterite precipitation in the solution environment. Conversely, when applied to MICP-treated sand, nano-Fe3O4 could facilitate calcite formation. Increasing the nano-Fe3O4 content showed a positive correlation with increased thermal conductivity and UCS. Specifically, the optimal values of thermal conductivity and UCS increased by 2.42 times and 2.39 times, respectively, compared to MICP-treated specimens without nano-Fe3O4. Microstructure analysis revealed that calcite precipitation at the particle contact served a dual function: cementing particles, thereby improving the mechanical strength and simultaneously acting as a "thermal bridge" to enhance the thermal conductivity. Furthermore, this study provides a new perspective on utilizing magnetized bacteria to reinforce specific locations within rocks and soils in the presence of an external magnetic field.

Removal of heavy metals from industrial paint effluent using coconut shell-derived activated carbon

This study evaluated the adsorption of heavy metals from paint industry effluent onto activated carbon produced from Coconut shells. The paint effluent samples were characterized before and after treatment via Atomic Absorption Spectroscopy (AAS) to identify the heavy metal ions present in the samples. The adsorption efficiency of the activated carbon and the effect of varying adsorbent particle size, impregnation ratio, and contact time were determined. The activated carbon from coconut shells was characterized via Scanning Electron Microscopy (SEM) and Fourier Transform Infrared Spectroscopy (FTIR) to determine the surface morphology and surface-active functional group. The SEM revealed development in pore size after adding a chemical activator to create functional adsorption sites. The finding concluded that the smaller particle size of adsorbents, higher surface area, higher contact time, and higher impregnation ratio increase the efficiency of the adsorption process.

Characterization and Performance of Peanut Shells in Caffeine and Triclosan Removal in Batch and Fixed-Bed Column Tests.

Peanut shells' adsorption performance in caffeine and triclosan removal was studied. Peanut shells were analyzed for their chemical composition, morphology, and surface functional groups. Batch adsorption and fixed-bed column experiments were carried out with solutions containing 30 mg/L of caffeine and triclosan. The parameters examined included peanut shell particle size (120-150, 300-600, and 800-2000 µm), adsorbent dose (0.02-60 g/L), contact time (up to 180 min), bed height (4-8 cm), and hydraulic loading rate (2.0 and 4.0 m3/m2-day). After determining the optimal adsorption conditions, kinetics, isotherm, and breakthrough curve models were applied to analyze the experimental data. Peanut shells showed an irregular surface and consisted mainly of polysaccharides (around 70% lignin, cellulose, and hemicellulose), with a specific surface area of 1.7 m2/g and a pore volume of 0.005 cm3/g. The highest removal efficiencies for caffeine (85.6 ± 1.4%) and triclosan (89.3 ± 1.5%) were achieved using the smallest particles and 10.0 and 0.1 g/L doses over 180 and 45 min, respectively. Triclosan showed easier removal compared to caffeine due to its higher lipophilic character. The pseudo-second-order kinetics model provided the best fit with the experimental data, suggesting a chemisorption process between caffeine/triclosan and the adsorbent. Equilibrium data were well-described by the Sips model, with maximum adsorption capacities of 3.3 mg/g and 289.3 mg/g for caffeine and triclosan, respectively. In fixed-bed column adsorption tests, particle size significantly influenced efficiency and hydraulic behavior, with 120-150 µm particles exhibiting the highest adsorption capacity for caffeine (0.72 mg/g) and triclosan (143.44 mg/g), albeit with clogging issues. The experimental data also showed good agreement with the Bohart-Adams, Thomas, and Yoon-Nelson models. Therefore, the findings of this study highlight not only the effective capability of peanut shells to remove caffeine and triclosan but also their versatility as a promising option for water treatment and sanitation applications in different contexts.

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles.

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely, dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data are relevant to another type of non-DLVO interaction, namely, the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. Our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

A numerical study on the effect of particle surface wettability on homogeneous and heterogeneous condensation processes in supersonic flows

Supersonic separators are an emerging device for processing natural gas impurities. Natural gas contains solid particulate impurities in addition to gaseous impurities, which can trigger heterogeneous condensation within the supersonic separator. There is a lack of research on the effect of particle wettability on supersonic condensation flow, although there have been studies on the effect of particle surface wettability on heterogeneous condensation in this area. Therefore, in this paper, the supersonic condensation flow is modeled considering the contact angle, and the effects of particle contact angle on heterogeneous condensation and homogeneous condensation are investigated separately. The results show that homogeneous condensation is inhibited and heterogeneous condensation is promoted by decreasing the particle contact angle. The reason is the smaller contact angle leads to the reduction of the free energy barrier for heterogeneous nucleation, which makes the conditions for heterogeneous nucleation lower and condensation more likely to occur.

Contact stress decomposition in large amplitude oscillatory shear of concentrated noncolloidal suspensions

The concentrated noncolloidal suspensions show complex rheological behavior, which is related to the existence of contact stress. However, determining the contact stress in time-varying flow like oscillatory shear is challenging. Herein, we propose a contact stress decomposition method to decompose the total stress directly into contact stress and hydrodynamic stress in large amplitude oscillatory shear (LAOS). The results of hydrodynamic stress and contact stress are consistent with those determined by the shear reversal experiment. The contact stress decomposition also explains the failure of the Cox–Merz rule in noncolloidal suspensions because the particle contacts exist in steady shear but are absent in small amplitude oscillatory shear. The intracycle and intercycle of contact stress are further analyzed through the general geometric average method. The intracycle behaviors exhibit strain hardening, strain softening, and shear thickening. The intercycle behaviors show bifurcations in stress-strain and stress-strain rate relations, where the transition strains at different concentrations define the state boundaries between the discrete particle contacts, the growing of particle contacts, and the saturated contacts. We also established a phenomenological constitutive model using a structural parameter to describe the shear effect on the buildup and breakdown of particle contacts. The contact stress of noncolloidal suspensions with wide ranges of particle concentrations and strain amplitudes under LAOS can be well described by the model.

Formation dynamics of branching structure in the slippery DLCA model.

We numerically investigated the aggregation dynamics and resulting network structures of colloidal gels using the slippery diffusion-limited cluster aggregation (DLCA) model. In this model, bonds are irreversibly formed upon the particle contacts, but the angles among them are not fixed, unlike the conventional DLCA. This allows clusters to be deformed in the process of aggregation. By characterizing the aggregation dynamics and using a reduced network scheme, our simulation revealed two distinct branching structure formation routes depending on the particle volume fraction ϕ. In lower volume fraction systems (ϕ ≤ 8%), the deformations of small-size clusters proceed prior to the percolation. When the Maxwell criterion is satisfied and the clusters become mechanically stable, the formation of the branching structure is nearly completed. After forming the branching structures, they aggregate and form a larger percolating network. Then, the aggregation proceeds through the elongation and straightening of the chain parts of the network. In higher volume fraction systems (ϕ > 8%), on the other hand, the clusters percolate, and a fine and homogeneous branching structure is formed at the early stage of the aggregation. In the aging stage, it collapses into a denser and more heterogeneous structure and becomes more stable. Our quantitative analyses of the branching structure will shed light on a new strategy for describing the network formation and elasticity of colloidal gels.

Investigation of the Bonding Mechanism of Al Powder Particles through Pulse Current Sintering Technology

Compared with traditional powder metallurgy, pulse current sintering is an advanced powder-forming technology, but its bonding mechanism is still an open topic for debate. In this paper, pulse current sintering is used as the connection technology and millimeter-sized Al particles are used as the research object. In the whole sintering process, no pressure was loaded; the function of the pulse current was the only source of heat with which to achieve the bonding of Al particles. The bonding mechanism of pulse current sintering was investigated from the perspective of material connection behavior. The results show that the pulse current density of the particle surface reaches 3.48 × 105 A/m2 instantly, while the current density of the particle center is only 8187 A/m2 at the initial stage, which is the main difference between pulse current sintering and traditional powder metallurgy sintering. With the densification process, the current density and temperature distribution in the contact region as well as the center of Al particles contact region tend to be more consistent. Finally, dense interfacial bonding was obtained, and the contact region of Al particles also demonstrated a high hardness value of 0.6385 GPa and yield strength value of 212.83 MPa. The whole process can be considered as a comprehensive action of melting (evaporation), diffusion, and plastic deformation. Based on the above results, a new technology, named high-frequency pulse current sintering, was proposed.

La0.8Sr0.2Ga0.8Mg0.2O3−δ − molten carbonate membrane for high − temperature CO2 separation

A perovskite structure La0.8Sr0.2Ga0.8Mg0.2O3−δ (LSGM), is a promising support material for dual − phase membrane for high − temperature CO2 separation owing to its high oxygen ion conductivity. The microstructure plays an important role in CO2 separation performance of the membrane. The result shows that effective particle − to − particle contact is a crucial guarantee for improving CO2 permeability. The uniformly distributed small pores of the support can greatly reduce the leakage. The Al2O3 coating on the porous support can improve the wettability between LSGM and molten carbonate, simultaneously increasing CO2 permeation flux while reducing leakage. The Al2O3 coated LSGM − based dual − phase membrane achieves a CO2 permeation flux of 0.34 mL min−1 cm−2 at 750 °C, with no detectable leakage. Above 750 °C, a reaction occurs between LSGM and molten carbonate, leading to a decline in membrane performance. For LSGM − carbonate dual − phase membranes, it is recommended to operate at temperatures not exceeding 750 °C.

Micromechanical analysis of contact erosion under cyclic loads using the coupled CFD‒DEM method

Different layers of soil often have distinct particle sizes. When exposed to the natural environment, soil is easily affected by natural rainfall, rising groundwater levels, and human activities, leading to particle contact erosion, which reduces the safety and service performance of the soil structure. In this paper, a coupled computational fluid dynamics–discrete element method (CFD–DEM) model was employed to investigate the particle migration phenomena, mechanical response of contact interfaces, variations in flow fields, and macroscopic deformation during the contact erosion process under cyclic loads at different frequencies and amplitudes. The conclusions are presented as follows: (1) Within one cycle of cyclic loading, both compression during loading and stress relaxation during unloading are the main factors triggering the migration of fine particles. (2) The migration and loss of fine particles mainly occur in the early stages of cyclic loading, where strong contact force chains are formed within the fine particle layer, leading to significant plastic deformation of the soil at the macroscopic level. (3) Under cyclic loading, changes in the soil pore structure cause an upwards hydraulic gradient in the initial quiescent water flow field. This hydraulic gradient can rupture weak contact force chains and cause particle pumping. (4) Increasing the frequency and amplitude of cyclic loading intensifies the erosion of fine particles, causing greater axial deformation of the soil. Compared to cyclic loading frequency, the amplitude of cyclic loading has a greater impact on contact erosion.

Applying SPH-DEM theory to the motion of particulate solid-liquid flow on continuous screening

Challenges in current research on fine mesh wet screening (≥50 mesh) include representing fine wire boundaries, contentious porous media handling, and low computational efficiency. To address these challenges, this study employs an averaged SPH-DEM method for modeling solid-liquid flow on screen. It integrates the Symplectic scheme and quaternion method to enhance computational accuracy and scalability. The fine wire boundaries are discretized into spheres to facilitate Hertz contact of DEM particles, and a novel screen aperture-fluid resistance model is introduced to calculate the resistance exerted on SPH particles by screen surface. The simulation results of dam break and periodic screening exhibit robust consistency with experimental data and other models, effectively validating the feasibility and accuracy of this methodology. Moreover, simulations of continuous screening using a 0.2 mm wire demonstrate that liquid enhances particle passage, further substantiating the practical viability of employing this model for simulating fine mesh wet screening.

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

An electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y2O3. Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

Discrete element modelling of the elastic-plastic and viscoelastic properties of a lithium-ion battery electrode layer

Mechanical degradation mechanisms are one of the leading causes of charge capacity loss in lithium-ion batteries. This study further develops a discrete element method (DEM) simulation framework, which investigates how the local contact behaviour affect the global mechanical properties of the active layer. The local microstructure consists of active particles held together by a binder domain, making up a granular medium. This study investigates the impact of the layer's global properties from the type of particle contact model. Experiments were also performed to measure size distribution and the material properties of the active material. The time dependency of the active layer, stemming from the viscoelastic binder domain, was studied in relaxation simulations, which were based on experimental measurements.

Investigation of the Shear Mechanism at Sand-Concrete Interface under the Influence of the Concave Groove Angle of the Contact Surface

A “random-type” sand–concrete interface shear test was developed based on the sand cone method, with a focus on the most commonly encountered triangular contact surface morphology. A “regular-type” triangular interface, matched in roughness to the “random-type”, was meticulously designed. This “regular-type” interface features five distinct triangular groove inclinations: 18°, 33°, 50°, 70°, and 90°. A series of sand–concrete interface direct shear tests were conducted under consistent compaction conditions to investigate the impact of varying compaction densities and triangular groove inclinations on the shear strength at the interface. Particle flow simulations were utilized to examine the morphology of the shear band and the characteristics of particle migration influenced by the triangular contact surface. This analysis is aimed at elucidating the influence of the inclination of the triangular groove on the shear failure mechanism at the sand–concrete interface. The findings indicate that: (1) The morphology of the interface significantly impacts the shear strength of the sand–concrete interface, while the shape of the stress-displacement curve experiences minimal alteration. (2) At smaller inclination angles, particle contact forces are arranged in a wave-like configuration around the sawtooth tip, resulting in a non-uniform stress distribution along the sawtooth slope. However, as the inclination angle grows, the stress concentration at the sawtooth tip diminishes, and the stress distribution across the sawtooth slope becomes more consistent. (3) Particle migration is significantly influenced by the sawtooth’s inclination angle. At lower angles, particles climb the structure’s tip through sliding and rolling. As the angle increases, particle motion shifts to shear, accompanied by a transition in friction from surface friction to internal shear friction. This leads to the formation of a wider shear band and an increase in shear strength.

Effects of scaling criteria on modelling of multi-phase flow in the packed bed using coarse grain CFD-DEM

The numerical modelling with computational fluid dynamics coupled with discrete element method (CFD-DEM) could be of high fidelity but rather time consuming, which leads to the development of coarse grain DEM (CFD-CGDEM). This paper describes the effects of scaling criteria for CFD-CGDEM on the modelling of multi-phase flow in packed bed. Four representative scaling criteria for particles contact interaction are discussed to analyze the packing structure and pressure drop and compare the relative deviation between CFD-CGDEM and CFD-DEM. It is found that when the spring constant in CFD-CGDEM is enlarged by α3 or α2 times or remains the same as that in original CFD-DEM, as α is defined as coarsen degree which is referred as the ratio between parcel size and particle size, a relatively large spring constant should be chosen to maintain the similar porosity and pressure drop. However, when the spring constant in CFD-CGDEM was scaled moderately by α times, the modelling results of CFD-CGDEM are comparable to data from CFD-DEM regardless of the value of spring constant. Furthermore, a quantitative analysis shows that the relative deviation in porosity and pressure drop between CFD-DEM and CFD-CGDEM will be enhanced as coarsen degree increases but the variation is insignificant with domain size for periodic boundary. The increase of parcel size could affect the simulation of porosity in packed bed and deviate the predicted hydrodynamics away from the CFD-DEM. In addition, the comparison of computational efficiency illustrates the speed up of CFD-CGDEM increases with the coarsen degree while the relative deviation also increases. This study provides instructive guidance for the application of CFD-CGDEM to model the multiphase flow in packed beds.

Triaxial Test and Discrete Element Numerical Simulation of Geogrid-Reinforced Clay Soil

Indoor triaxial tests on geogrid-reinforced clay elucidate the macroscopic changes in soil strength indices post-reinforcement, yet the underlying mechanisms of strength enhancement require further investigation. By conducting indoor triaxial tests and establishing a corresponding discrete element numerical model, we can delve into the fine-scale mechanisms of geogrid-reinforced soil. This includes analyzing changes in fine-scale parameters such as porosity, the coordination number, and contact stress between soil particles. The findings suggest that an increase in the number of geogrid reinforcement layers leads to a more pronounced improvement in peak strength and cohesion, albeit with minimal impact on the internal friction angle of the specimens. Furthermore, analysis of the triaxial test curves of reinforced soils indicates that the stress–strain relationship adheres to the Duncan–Chang model. Parameters derived from this model have been validated against experimental data, confirming their accuracy. The discrete element model was used to analyze the variations in fine-scale parameters such as porosity and coordination number. It revealed that reinforcement reduces the fluctuation amplitude of porosity and significantly increases the number of particle contacts, resulting in a denser soil structure. Further analysis of the change in contact stress between particles in the discrete element model revealed that the contact force between particles increased significantly after reinforcement and that the reinforcement played a role in restraining the soil particles and dispersing the reinforcement stress, which explains the increase in the strength of the mesh-reinforced clays from another perspective. This further elucidates the strength enhancement mechanism in geogrid-reinforced clay, offering a new perspective on the mechanical behavior and strength development of such materials.