Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies ...Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies between 1 keV and 150 keV.The results indicate that a smaller grain size leads to more defects forming in grain boundary regions during cascade processes.The impact of high-energy PKA may cause a certain degree of distortion of the grain boundaries,which has a higher probability in systems with smaller grain sizes and becomes more pronounced as the PKA energy increases.The direction of PKA can affect the formation and diffusion pathways of defects.When the PKA direction is perpendicular to the grain boundary,defects preferentially form near the grain boundary regions;by contrast,defects are more inclined to form in the interior of the grains.These results are of great significance for comprehending the changes in the performance of polycrystalline W under the high-energy fusion environments and can provide theoretical guidance for further optimization and application of W-based plasma materials.展开更多
Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for...Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for simulating solvation structures.However,mainstream force fields often lack accuracy in describing strong ion-solvent interactions,causing disparities between CMD simulations and experimental observations.Although some empirical methods have been employed in some of the studies to address this issue,their effectiveness has been limited.Our CMD research,supported by quantum chemical calculations and experimental data,reveals that the solvation structure is influenced not only by the charge model but also by the polarization description.Previous empirical approaches that focused solely on adjusting ion-solvent interaction strengths overlooked the importance of polarization effects.Building on this insight,we propose integrating the Drude polarization model into mainstream force fields and verify its feasibility in carbonate,ether,and nitrile electrolytes.Our experimental results demonstrate that this approach significantly enhances the accuracy of CMD-simulated solvation structures.This work is expected to provide a more reliable CMD method for electrolyte design,shielding researchers from the pitfalls of erroneous simulation outcomes.展开更多
As coal mining depth increases,the combined effects of high stress,mining stress,and fault structures make dynamic impact hazards more frequent.The reproduction of dynamic impact phenomena is basis for studying their ...As coal mining depth increases,the combined effects of high stress,mining stress,and fault structures make dynamic impact hazards more frequent.The reproduction of dynamic impact phenomena is basis for studying their occurrence patterns and control mechanisms.Physical simulation test represents an efficacious methodology.However,there is currently a lack of simulation devices that can effectively simulate two types of dynamic impact phenomena,including high stress and fault slip dynamic impact.To solve aforementioned issues,the physical simulation test system for dynamic impact in deep roadways developed by authors is employed to carry out comparative tests of high stress and fault slip dynamic impact.The phenomena of high stress and fault slip dynamic impact are reproduced successfully.A comparative analysis is conducted on dynamic phenomena,stress evolution,roadway deformation,and support force.The high stress dynamic impact roadway instability mode,which is characterized by the release of high energy accompanied by symmetric damage,and the fault slip dynamic impact roadway instability mode,which is characterized by the propagation of unilateral stress waves accompanied by asymmetric damage,are clarified.On the basis,the differentiated control concepts for different types of dynamic impact in deep roadways are proposed.展开更多
The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular ...The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular dynamics method through the model of γ-γ′phases containing hole defect.The addition of Re makes the dislocation distribution tend towards the γ phase.The higher the Re content,the earlier theγphase yields,while the γ′phase yields later.Dislocation bends under the combined action of the applied force and the resistance of the Re atoms to form a bend point.The Re atoms are located at the bend points and strengthen the alloy by fixing the dislocation and preventing it from cutting the γ′phase.Dislocations nucleate first in the γ phase,causing theγphase to deform plastically before the γ′phase.As the strain increases,the dislocation length first remains unchanged,then increases rapidly,and finally fluctuates and changes.The dislocation lengths in the γ phase are larger than those in the γ′phase at different temperatures.The dislocation length shows a decreasing tendency with the increase of the temperature.Temperature can affect movement of the dislocation,and superalloys have different plastic deformation mechanisms at low,medium and high temperatures.展开更多
Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the p...Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.展开更多
In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the re...In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.展开更多
To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitro...To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.展开更多
According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the p...According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.展开更多
Bit patterned recording(BPR)has attracted much attention due to its promising potential in achieving high densities in magnetic storage devices.The materials with strong perpendicular magnetic anisotropy(PMA)are alway...Bit patterned recording(BPR)has attracted much attention due to its promising potential in achieving high densities in magnetic storage devices.The materials with strong perpendicular magnetic anisotropy(PMA)are always preferred in designing the BPR.Here,the patterned Co/Ni multilayers showing d-d hybridization induced PMA was studied.In particular,we record the ultrafast spin dynamics by means of time-resolved scanning magneto-optical Kerr effect(TRMOKE)microscopy.We are able to acquire the“snapshot”magnetic maps of the sample surface because of both the femtosecond temporal and submicrometer spatial resolution in our TRMOKE microscopy.Furthermore,the spatially inhomogeneous ultrafast demagnetization was observed in experiment,and this has been evidenced by simulations.展开更多
Rotational dynamics simulations of neutral O_(2)molecules driven by linearly,elliptically and circularly polarized femtosecond pulsed lasers are carried out using a full quantum time-dependent wave packet evolution me...Rotational dynamics simulations of neutral O_(2)molecules driven by linearly,elliptically and circularly polarized femtosecond pulsed lasers are carried out using a full quantum time-dependent wave packet evolution method.Here,the direction of laser propagation is set along the z axis,and the polarization plane is restricted to the xy plane.The results indicate that the alignment of O_(2)molecules in the z direction is weakly affected by varying the ellipticity when the total laser intensity is held constant.For rotation within the xy plane,the linearly polarized laser significantly excites rotational motion,with the degree of excitation increasing as the ellipticity increases.In contrast,under the influence of a circularly polarized laser,the angular distribution of O_(2)molecules in the xy plane remains isotropic.Additionally,the effects of the initial rotational quantum number,the temperature of the O_(2)molecules and the nuclear spin on laser-induced alignment are discussed.展开更多
We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molec...We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molecular dynamics(EMD)simulations.D_(G)and D_(E)are computed for a broad domain of screening length(κ)and coupling parameters(Г)along with different system sizes.It is observed that both D_(G)and D_(E)decrease linearly with increasing G in warm liquid states and increase with increasingκ.In cold liquid states,the Einstein relation accurately predicts D_(E)in 2D SC-DPs because diffusion motion is close to normal diffusion,but the G-K relation provides overestimations of D_(G),because VACF indicates anomalous diffusion;thus,D_(G)is not accurate.Our new simulation outcomes reveal that D_(G)and D_(E)remain independent of system sizes.Furthermore,our investigations demonstrate that at higher temperatures,D_(G)and D_(E)converge,suggesting diffusion motion close to normal diffusion,while at lower temperatures,these two values diverge.We find reasonable agreement by comparing current and existing numerical,theoretical and experimental data.Moreover,when normalizing diffusion coefficients by the Einstein frequency and testing against the universal temperature scaling law,D_(G)deviates from theoretical curves at low temperatures and k,whereas D_(E)only disagrees with theory at very smallκ(■0.10).These findings provide valuable insight into diagnosing dust component parameters within 2D DP systems and contribute to the broader understanding of diffusion processes in DP environments.展开更多
In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforatio...In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.展开更多
Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility ...Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.展开更多
Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear...Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.展开更多
Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced...Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.展开更多
This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing...This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.展开更多
In this study,we focused on a novel parallel mechanism for utilizing the motion simulator of a high-speed boat(HSB).First,we expressed the real behavior of the HSB based on a seakeeping trial.For this purpose,we recor...In this study,we focused on a novel parallel mechanism for utilizing the motion simulator of a high-speed boat(HSB).First,we expressed the real behavior of the HSB based on a seakeeping trial.For this purpose,we recorded the motion parameters of the HSB by gyroscope and accelerometer sensors,while using a special data acquisition technique.Additionally,a Chebychev highpass filter was applied as a noise filter to the accelerometer sensor.Then,a novel 3 degrees of freedom(DoF)parallel mechanism(1T2R)with prismatic actuators is proposed and analyses were performed on its inverse kinematics,velocity,and acceleration.Finally,the inverse dynamic analysis is presented by the principle of virtual work,and the validation of the analytical equations was compared by the ADAMS simulation software package.Additionally,according to the recorded experimental data of the HSB,the feasibility of the proposed novel parallel mechanism motion simulator of the HSB,as well as the necessity of using of the washout filters,was explored.展开更多
Based on the dynamical characteristic parameters of the real vehicle, the modehng approach ancl procedure of dynamics of vehicles are expatiated. The layout of vehicle dynamics is proposed, and the sub-models of the d...Based on the dynamical characteristic parameters of the real vehicle, the modehng approach ancl procedure of dynamics of vehicles are expatiated. The layout of vehicle dynamics is proposed, and the sub-models of the diesel engine, drivetrain system and vehicle multi-body dynamics are introduced. Finally, the running characteristic data of the virtual and real vehicles are compared, which shows that the dynamics model is similar closely to the real vehicle system.展开更多
基金Project supported by the National MCF Energy Research and Development Program of China(Grant No.2018YFE0308101)the National Key Research and Development Program of China(Grant No.2018YFB0704000)+1 种基金the Suqian Science and Technology Program(Grant No.K202337)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.23KJD490001).
文摘Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies between 1 keV and 150 keV.The results indicate that a smaller grain size leads to more defects forming in grain boundary regions during cascade processes.The impact of high-energy PKA may cause a certain degree of distortion of the grain boundaries,which has a higher probability in systems with smaller grain sizes and becomes more pronounced as the PKA energy increases.The direction of PKA can affect the formation and diffusion pathways of defects.When the PKA direction is perpendicular to the grain boundary,defects preferentially form near the grain boundary regions;by contrast,defects are more inclined to form in the interior of the grains.These results are of great significance for comprehending the changes in the performance of polycrystalline W under the high-energy fusion environments and can provide theoretical guidance for further optimization and application of W-based plasma materials.
基金supported by the Science and Technology Project of State Grid Corporation of China(5419-202199552A-0-5-ZN).
文摘Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for simulating solvation structures.However,mainstream force fields often lack accuracy in describing strong ion-solvent interactions,causing disparities between CMD simulations and experimental observations.Although some empirical methods have been employed in some of the studies to address this issue,their effectiveness has been limited.Our CMD research,supported by quantum chemical calculations and experimental data,reveals that the solvation structure is influenced not only by the charge model but also by the polarization description.Previous empirical approaches that focused solely on adjusting ion-solvent interaction strengths overlooked the importance of polarization effects.Building on this insight,we propose integrating the Drude polarization model into mainstream force fields and verify its feasibility in carbonate,ether,and nitrile electrolytes.Our experimental results demonstrate that this approach significantly enhances the accuracy of CMD-simulated solvation structures.This work is expected to provide a more reliable CMD method for electrolyte design,shielding researchers from the pitfalls of erroneous simulation outcomes.
基金supported by the National Natural Science Foundation of China(Nos.U24A2088,42177130,42277174,and 42477166).
文摘As coal mining depth increases,the combined effects of high stress,mining stress,and fault structures make dynamic impact hazards more frequent.The reproduction of dynamic impact phenomena is basis for studying their occurrence patterns and control mechanisms.Physical simulation test represents an efficacious methodology.However,there is currently a lack of simulation devices that can effectively simulate two types of dynamic impact phenomena,including high stress and fault slip dynamic impact.To solve aforementioned issues,the physical simulation test system for dynamic impact in deep roadways developed by authors is employed to carry out comparative tests of high stress and fault slip dynamic impact.The phenomena of high stress and fault slip dynamic impact are reproduced successfully.A comparative analysis is conducted on dynamic phenomena,stress evolution,roadway deformation,and support force.The high stress dynamic impact roadway instability mode,which is characterized by the release of high energy accompanied by symmetric damage,and the fault slip dynamic impact roadway instability mode,which is characterized by the propagation of unilateral stress waves accompanied by asymmetric damage,are clarified.On the basis,the differentiated control concepts for different types of dynamic impact in deep roadways are proposed.
基金Project supported by the Xi’an Science and Technology Plan Project of Shaanxi Province of China(Grant No.23GXFW0086).
文摘The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular dynamics method through the model of γ-γ′phases containing hole defect.The addition of Re makes the dislocation distribution tend towards the γ phase.The higher the Re content,the earlier theγphase yields,while the γ′phase yields later.Dislocation bends under the combined action of the applied force and the resistance of the Re atoms to form a bend point.The Re atoms are located at the bend points and strengthen the alloy by fixing the dislocation and preventing it from cutting the γ′phase.Dislocations nucleate first in the γ phase,causing theγphase to deform plastically before the γ′phase.As the strain increases,the dislocation length first remains unchanged,then increases rapidly,and finally fluctuates and changes.The dislocation lengths in the γ phase are larger than those in the γ′phase at different temperatures.The dislocation length shows a decreasing tendency with the increase of the temperature.Temperature can affect movement of the dislocation,and superalloys have different plastic deformation mechanisms at low,medium and high temperatures.
基金Project supported by the National Natural Science Foundation of China (Grant No. 52271105)。
文摘Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.
文摘In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.
基金supported by the National Natural Science Foundation of China (22275018)the Project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)(Grant No.QNKT20-04)。
文摘To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.
基金Supported by the Sichuan Province Regional Innovation Cooperation Project(21QYCX0048)Sinopec Science and Technology Department Project(P21048-3)。
文摘According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.
基金Project supported by the National Key Research and Development Program of China(Grant No.2022YFA1403302)the National Natural Science Foundation of China(Grant Nos.52031015,U22A20115,and 12104030)+1 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.LZ25A040007)the Natural Science Foundation of Beijing(Grant No.1252026).
文摘Bit patterned recording(BPR)has attracted much attention due to its promising potential in achieving high densities in magnetic storage devices.The materials with strong perpendicular magnetic anisotropy(PMA)are always preferred in designing the BPR.Here,the patterned Co/Ni multilayers showing d-d hybridization induced PMA was studied.In particular,we record the ultrafast spin dynamics by means of time-resolved scanning magneto-optical Kerr effect(TRMOKE)microscopy.We are able to acquire the“snapshot”magnetic maps of the sample surface because of both the femtosecond temporal and submicrometer spatial resolution in our TRMOKE microscopy.Furthermore,the spatially inhomogeneous ultrafast demagnetization was observed in experiment,and this has been evidenced by simulations.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFA1602502)the National Natural Science Foundation of China(Grant No.12450404).
文摘Rotational dynamics simulations of neutral O_(2)molecules driven by linearly,elliptically and circularly polarized femtosecond pulsed lasers are carried out using a full quantum time-dependent wave packet evolution method.Here,the direction of laser propagation is set along the z axis,and the polarization plane is restricted to the xy plane.The results indicate that the alignment of O_(2)molecules in the z direction is weakly affected by varying the ellipticity when the total laser intensity is held constant.For rotation within the xy plane,the linearly polarized laser significantly excites rotational motion,with the degree of excitation increasing as the ellipticity increases.In contrast,under the influence of a circularly polarized laser,the angular distribution of O_(2)molecules in the xy plane remains isotropic.Additionally,the effects of the initial rotational quantum number,the temperature of the O_(2)molecules and the nuclear spin on laser-induced alignment are discussed.
基金support of the Fundamental Research Funds for the Central Universities of China(Grant No.2019ZDPY16).
文摘We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molecular dynamics(EMD)simulations.D_(G)and D_(E)are computed for a broad domain of screening length(κ)and coupling parameters(Г)along with different system sizes.It is observed that both D_(G)and D_(E)decrease linearly with increasing G in warm liquid states and increase with increasingκ.In cold liquid states,the Einstein relation accurately predicts D_(E)in 2D SC-DPs because diffusion motion is close to normal diffusion,but the G-K relation provides overestimations of D_(G),because VACF indicates anomalous diffusion;thus,D_(G)is not accurate.Our new simulation outcomes reveal that D_(G)and D_(E)remain independent of system sizes.Furthermore,our investigations demonstrate that at higher temperatures,D_(G)and D_(E)converge,suggesting diffusion motion close to normal diffusion,while at lower temperatures,these two values diverge.We find reasonable agreement by comparing current and existing numerical,theoretical and experimental data.Moreover,when normalizing diffusion coefficients by the Einstein frequency and testing against the universal temperature scaling law,D_(G)deviates from theoretical curves at low temperatures and k,whereas D_(E)only disagrees with theory at very smallκ(■0.10).These findings provide valuable insight into diagnosing dust component parameters within 2D DP systems and contribute to the broader understanding of diffusion processes in DP environments.
基金the support of the Foundation of Natural Science Foundation of Shaanxi Province, Grant/ Award nos. 2023-JC-YB-361National Natural Science Foundation (Number 52104033)。
文摘In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.
基金Project supported by the National Natural Science Foundation of China(Grant No.52072132).
文摘Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.51775077 and 51909023)。
文摘Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.
基金supported by the National Natural Science Foundation of China (Grant No.11572090)。
文摘Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.
基金supported by the Major Scientific and Technological Projects of CNPC under grant ZD2019-183-006the National Science and Technology Major Project of China(2016ZX05014002-006)the National Natural Science Foundation of China(42072234)。
文摘This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.
文摘In this study,we focused on a novel parallel mechanism for utilizing the motion simulator of a high-speed boat(HSB).First,we expressed the real behavior of the HSB based on a seakeeping trial.For this purpose,we recorded the motion parameters of the HSB by gyroscope and accelerometer sensors,while using a special data acquisition technique.Additionally,a Chebychev highpass filter was applied as a noise filter to the accelerometer sensor.Then,a novel 3 degrees of freedom(DoF)parallel mechanism(1T2R)with prismatic actuators is proposed and analyses were performed on its inverse kinematics,velocity,and acceleration.Finally,the inverse dynamic analysis is presented by the principle of virtual work,and the validation of the analytical equations was compared by the ADAMS simulation software package.Additionally,according to the recorded experimental data of the HSB,the feasibility of the proposed novel parallel mechanism motion simulator of the HSB,as well as the necessity of using of the washout filters,was explored.
基金Sponsored by the Ministerial Level Foundation (070312)
文摘Based on the dynamical characteristic parameters of the real vehicle, the modehng approach ancl procedure of dynamics of vehicles are expatiated. The layout of vehicle dynamics is proposed, and the sub-models of the diesel engine, drivetrain system and vehicle multi-body dynamics are introduced. Finally, the running characteristic data of the virtual and real vehicles are compared, which shows that the dynamics model is similar closely to the real vehicle system.
