Using cemented paste backfill to tackle the phosphogypsum stockpile in China:A down-to-earth technology with new vitalities in pollutant retention and CO_(2) abatement

Phosphogypsum(PG),a hard-to-dissipate by-product of the phosphorus fertilizer production industry,places strain on the biogeochemical cycles and ecosystem functions of storage sites.This pervasive problem is already widespread worldwide and requires careful stewardship.In this study,we review the presence of potentially toxic elements(PTEs)in PG and describe their associations with soil properties,anthropogenic activities,and surrounding organisms.Then,we review different ex-/in-situ solutions for promoting the sustainable management of PG,with an emphasis on in-situ cemented paste backfill,which offers a cost-effective and highly scalable opportunity to advance the value-added recovery of PG.However,concerns related to the PTEs'retention capacity and long-term effectiveness limit the implementation of this strategy.Furthermore,given that the large-scale demand for ordinary Portland cement from this conventional option has resulted in significant CO_(2) emissions,the technology has recently undergone additional scrutiny to meet the climate mitigation ambition of the Paris Agreement and China's Carbon Neutrality Economy.Therefore,we discuss the ways by which we can integrate innovative strategies,including supplementary cementitious materials,alternative binder solutions,CO_(2) mineralization,CO_(2) curing,and optimization of the supply chain for the profitability and sustainability of PG remediation.However,to maximize the co-benefits in environmental,social,and economic,future research must bridge the gap between the feasibility of expanding these advanced pathways and the multidisciplinary needs.

Anisotropic shearing mechanism of Kangding slate:Experimental investigation and numerical analysis

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.

Microstructure evolution and strengthening mechanism of high -performance powder metallurgy TA15 titanium alloy by hot rolling

Hot deformation of sintered billets by powder metallurgy(PM)is an effective preparation technique for titanium alloys,which is more significant for high-alloying alloys.In this study,Ti–6.5Al–2Zr–Mo–V(TA15)titanium alloy plates were prepared by cold press-ing sintering combined with high-temperature hot rolling.The microstructure and mechanical properties under different process paramet-ers were investigated.Optical microscope,electron backscatter diffraction,and others were applied to characterize the microstructure evolution and mechanical properties strengthening mechanism.The results showed that the chemical compositions were uniformly dif-fused without segregation during sintering,and the closing of the matrix craters was accelerated by increasing the sintering temperature.The block was hot rolled at 1200℃ with an 80%reduction under only two passes without annealing.The strength and elongation of the plate at 20–25℃ after solution and aging were 1247 MPa and 14.0%,respectively,which were increased by 24.5%and 40.0%,respect-ively,compared with the as-sintered alloy at 1300℃.The microstructure was significantly refined by continuous dynamic recrystalliza-tion,which was completed by the rotation and dislocation absorption of the substructure surrounded by low-angle grain boundaries.After hot rolling combined with heat treatment,the strength and plasticity of PM-TA15 were significantly improved,which resulted from the dense,uniform,and fine recrystallization structure and the synergistic effect of multiple slip systems.

Genesis of the Nuri Cu-W-Mo Deposit,Tibet,China:Constraints from in situ Trace Elements and Sr Isotopic Analysis of Scheelite

The Nuri deposit is the only Cu-W-Mo polymetallic deposit with large-scale WO3 resources in the eastern section of the Gangdese metallogenic belt,Tibet,China.However,the genetic type of this deposit has been controversial since its discovery.Based on a study of the geological characteristics of the deposit,this study presents mineralization stages,focusing on the oxide stage and the quartz-sulfide stage where scheelite is mainly formed,referred to as Sch-A and Sch-B,respectively.Through LA-ICP-MS trace element and Sr isotope analyses,the origin,evolutionary process of the oreforming fluid and genesis of the ore deposit are investigated.Scanning Electron Microscope-Cathodoluminescence(SEMCL)observations reveal that Sch-A consists of three generations,with dark gray homogenous Sch-A1 being replaced by relatively lighter and homogeneous Sch-A2 and Sch-A3,with Sch-A2 displaying a gray CL image color with vague and uneven growth bands and Sch-A3 has a light gray CL image color with hardly any growth band.In contrast,Sch-B exhibits a‘core-rim’structure,with the core part(Sch-B1)being dark gray and displaying a uniform growth band,while the rim part(Sch-B2)is light gray and homogeneous.The normalized distribution pattern of rare earth elements in scheelite and Sr isotope data suggest that the early ore-forming fluid in the Nuri deposit originated from granodiorite porphyry and,later on,some country rock material was mixed in,due to strong water-rock interaction.Combining the O-H isotope data further indicates that the ore-forming fluid in the Nuri deposit originated from magmatic-hydrothermal sources,with contributions from metamorphic water caused by water-rock interaction during the mineralization process,as well as later meteoric water.The intense water-rock interaction likely played a crucial role in the precipitation of scheelite,leading to varying Eu anomalies in different generations of scheelite from the oxide stage to the quartz-sulfide stage,while also causing a gradual decrease in oxygen fugacity(fO2)and a slow rise in pH value.Additionally,the high Mo and low Sr contents in the scheelite are consistent with typical characteristics of magmatic-hydrothermal scheelite.Therefore,considering the geological features of the deposit,the geochemical characteristics of scheelite and the O-H isotope data published previously,it can be concluded that the genesis of the Nuri deposit belongs to porphyry-skarn deposit.

NiCoZn/C@melamine sponge-derived carbon composites with highperformance electromagnetic wave absorption

NiMZn/C@melamine sponge-derived carbon(MSDC)composites(M=Co,Fe,and Mn)were prepared by a vacuum pump-ing solution method followed by carbonization.A large number of carbon nanotubes(CNTs)homogeneously attached to the surfaces of the three-dimensional cross-linked of the sponge-derived carbon in the NiCoZn/C@MSDC composite,and CNTs were detected in the NiFeZn/C@MSDC and NiMnZn/C@MSDC composites.Ni_(3)ZnC_(0.7),Ni_(3)Fe,and MnO in-situ formed in the NiFeZn/C@MSDC and NiMnZn/C@MSDC composites.The CNTs in the NiCoZn/C@MSDC composite efficiently modulated its complex permittivity.Thus,the composite exhibited the best performance among the composites,with the minimum reflection loss(RL_(min))of-33.1 dB at 18 GHz and thickness of 1.4 mm.The bandwidth for RL of≤-10 dB was up to 5.04 GHz at the thickness of 1.7 mm and loading of 25wt%.The op-timized impedance matching,enhanced interfacial and dipole polarization,remarkable conduction loss,and multiple reflections and scat-tering of the incident microwaves improved the microwave absorption performance.The effects of Co,Ni,and Fe on the phase and mor-phology provided an alternative way for developing highly efficient and broadband microwave absorbers.

Numerical modeling of fracture propagation of supercritical CO_(2)compound fracturing

The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.

Multi-layer perceptron-based data-driven multiscale modelling of granular materials with a novel Frobenius norm-based internal variable

One objective of developing machine learning(ML)-based material models is to integrate them with well-established numerical methods to solve boundary value problems(BVPs).In the family of ML models,recurrent neural networks(RNNs)have been extensively applied to capture history-dependent constitutive responses of granular materials,but these multiple-step-based neural networks are neither sufficiently efficient nor aligned with the standard finite element method(FEM).Single-step-based neural networks like the multi-layer perceptron(MLP)are an alternative to bypass the above issues but have to introduce some internal variables to encode complex loading histories.In this work,one novel Frobenius norm-based internal variable,together with the Fourier layer and residual architectureenhanced MLP model,is crafted to replicate the history-dependent constitutive features of representative volume element(RVE)for granular materials.The obtained ML models are then seamlessly embedded into the FEM to solve the BVP of a biaxial compression case and a rigid strip footing case.The obtained solutions are comparable to results from the FEM-DEM multiscale modelling but achieve significantly improved efficiency.The results demonstrate the applicability of the proposed internal variable in enabling MLP to capture highly nonlinear constitutive responses of granular materials.

Biomedical rare-earth magnesium alloy:Current status and future prospects

Biomedical magnesium(Mg)alloys have garnered significant attention because of their unique biodegradability,favorable biocompatibility,and suitable mechanical properties.The incorporation of rare earth(RE)elements,with their distinct physical and chemical properties,has greatly contributed to enhancing the mechanical performance,degradation behavior,and biological performance of biomedical Mg alloys.Currently,a series of RE-Mg alloys are being designed and investigated for orthopedic implants and cardiovascular stents,achieving substantial and encouraging research progress.In this work,a comprehensive summary of the state-of-the-art in biomedical RE-Mg alloys is provided.The physiological effects and design standards of RE elements in biomedical Mg alloys are discussed.Particularly,the degradation behavior and mechanical properties,including their underlying action are studied in-depth.Furthermore,the preparation techniques and current application status of RE-Mg alloys are reviewed.Finally,we address the ongoing challenges and propose future prospects to guide the development of high-performance biomedical Mg-RE alloys.

Design of low-alloying and high-performance solid solution-strengthened copper alloys with element substitution for sustainable development

Solid solution-strengthened copper alloys have the advantages of a simple composition and manufacturing process,high mechanical and electrical comprehensive performances,and low cost;thus,they are widely used in high-speed rail contact wires,electronic component connectors,and other devices.Overcoming the contradiction between low alloying and high performance is an important challenge in the development of solid solution-strengthened copper alloys.Taking the typical solid solution-strengthened alloy Cu-4Zn-1Sn as the research object,we proposed using the element In to replace Zn and Sn to achieve low alloying in this work.Two new alloys,Cu-1.5Zn-1Sn-0.4In and Cu-1.5Zn-0.9Sn-0.6In,were designed and prepared.The total weight percentage content of alloying elements decreased by 43%and 41%,respectively,while the product of ultimate tensile strength(UTS)and electrical conductivity(EC)of the annealed state increased by 14%and 15%.After cold rolling with a 90%reduction,the UTS of the two new alloys reached 576 and 627MPa,respectively,the EC was 44.9%IACS and 42.0%IACS,and the product of UTS and EC(UTS×EC)was 97%and 99%higher than that of the annealed state alloy.The dislocations proliferated greatly in cold-rolled alloys,and the strengthening effects of dislocations reached 332 and 356 MPa,respectively,which is the main reason for the considerable improvement in mechanical properties.

Improving creep strength of the fine-grained heat-affected zone of novel 9Cr martensitic heat-resistant steel via modified thermo-mechanical treatment

The infamous type Ⅳ failure within the fine-grained heat-affected zone (FGHAZ) in G115 steel weldments seriously threatens the safe operation of ultra-supercritical (USC) power plants.In this work,the traditional thermo-mechanical treatment was modified via the replacement of hot-rolling with cold rolling,i.e.,normalizing,cold rolling,and tempering (NCT),which was developed to improve the creep strength of the FGHAZ in G115 steel weldments.The NCT treatment effectively promoted the dissolution of preformed M_(23)C_(6)particles and relieved the boundary segregation of C and Cr during welding thermal cycling,which accelerated the dispersed reprecipitation of M_(23)C_(6) particles within the fresh reaustenitized grains during post-weld heat treatment.In addition,the precipitation of Cu-rich phases and MX particles was promoted evidently due to the deformation-induced dislocations.As a result,the interacting actions between precipitates,dislocations,and boundaries during creep were reinforced considerably.Following this strategy,the creep rupture life of the FGHAZ in G115 steel weldments can be prolonged by 18.6%,which can further push the application of G115 steel in USC power plants.

Method to measure tree-ring width,density,elemental composition,and stable carbon and oxygen isotopes using one sample

Tree-ring width(RW),density,elemental com-position,and stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)are widely used as proxies to assess climate change,ecology,and environmental pollution;however,a specific pretreat-ment has been needed for each proxy.Here,we developed a method by which each proxy can be measured in the same sample.First,the sample is polished for ring width meas-urement.After obtaining the ring width data,the sample is cut to form a 1-mm-thick wood plate.The sample is then mounted in a vertical sample holder,and gradually scanned by an X-ray beam.Simultaneously,the count rates of the fluorescent photons of elements(for chemical characteriza-tion)and a radiographic grayscale image(for wood density)are obtained,i.e.the density and the element content are obtained.Then,cellulose is isolated from the 1-mm wood plate by removal of lignin,and hemicellulose.After producing this cellulose plate,cellulose subsamples are separated by knife under the microscope for inter-annual and intra-annual stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)analysis.Based on this method,RW,density,elemental composition,δ^(13)C,and δ^(18)O can be measured from the same sample,which reduces sample amount and treatment time,and is helpful for multi-proxy comparison and combination research.

Numerical analysis of moving train induced vibrations on tunnel,surrounding ground and structure

This study is focused on the effect of vibration induced by moving trains in tunnels on the surrounding ground and structures.A three-dimensional finite element model is established for a one-track railway tunnel and an adjacent twelve-storey building frame by using commercial software Midas GTS-NX(2019)and Midas Gen.This study considered the moving load effect of a complete train,which varies with space as well as with time.The effect of factors such as train speed,overburden pressure on the tunnel and variation in soil properties are studied in the time domain.As a result,the variations in horizontal and vertical acceleration for two different sites,i.e.,the free ground surface(without structure)and the area containing the structure,are compared.Also,the displacement pattern of the raft foundation is plotted for different train velocities.At lower speeds,the heaving phenomenon is negligible,but as the speed increases,both the heaving and differential settlement increase in the foundation.This study demonstrates that the effect of moving train vibrations should be considered in the design of new nearby structures and proper ground improvement should be considered for existing structures.

Effect of Ti and Ta content on the oxidation resistance of Co-Ni-based superalloys

Co-Ni-based superalloys are known for their capability to function at elevated temperatures and superior hot corrosion and thermal fatigue resistance.Therefore,these alloys show potential as crucial high-temperature structural materials for aeroengine and gas turbine hot-end components.Our previous work elucidated the influence of Ti and Ta on the high-temperature mechanical properties of alloys.However,the intricate interaction among elements considerably affects the oxidation resistance of alloys.In this paper,Co-35Ni-10Al-2W-5Cr-2Mo-1Nb-xTi-(5−x)Ta alloys(x=1,2,3,4)with varying Ti and Ta contents were designed and compounded,and their oxidation resistance was investigated at the temperature range from 800 to 1000℃.After oxidation at three test conditions,namely,800℃for 200 h,900℃for 200 h,and 1000℃for 50 h,the main structure of the oxide layer of the alloy consisted of spinel,Cr_(2)O_(3),and Al_(2)O_(3)from outside to inside.Oxides consisting of Ta,W,and Mo formed below the Cr_(2)O_(3)layer.The interaction of Ti and Ta imparted the highest oxidation resistance to 3Ti2Ta alloy.Conversely,an excessive amount of Ti or Ta resulted in an adverse effect on the oxidation resistance of the alloys.This study reports the volatilization of W and Mo oxides during the oxidation process of Co-Ni-based cast superalloys with a high Al content for the first time and explains the formation mechanism of holes in the oxide layer.The results provide a basis for gaining insights into the effects of the interaction of alloying elements on the oxidation resistance of the alloys they form.

Numerical Study of the Biomechanical Behavior of a 3D Printed Polymer Esophageal Stent in the Esophagus by BP Neural Network Algorithm

Esophageal disease is a common disorder of the digestive system that can severely affect the quality of life andprognosis of patients. Esophageal stenting is an effective treatment that has been widely used in clinical practice.However, esophageal stents of different types and parameters have varying adaptability and effectiveness forpatients, and they need to be individually selected according to the patient’s specific situation. The purposeof this study was to provide a reference for clinical doctors to choose suitable esophageal stents. We used 3Dprinting technology to fabricate esophageal stents with different ratios of thermoplastic polyurethane (TPU)/(Poly-ε-caprolactone) PCL polymer, and established an artificial neural network model that could predict the radial forceof esophageal stents based on the content of TPU, PCL and print parameter. We selected three optimal ratios formechanical performance tests and evaluated the biomechanical effects of different ratios of stents on esophagealimplantation, swallowing, and stent migration processes through finite element numerical simulation and in vitrosimulation tests. The results showed that different ratios of polymer stents had different mechanical properties,affecting the effectiveness of stent expansion treatment and the possibility of postoperative complications of stentimplantation.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Analysis of mesoscopic mechanical dynamic characteristics of ballast bed with under sleeper pads

The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.

Computer vision-aided DEM study on the compaction characteristics of graded subgrade filler considering realistic coarse particle shapes

The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.

A particle-resolved heat-particle-fluid coupling model by DEM-IMB-LBM

Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.