In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary unde...In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.展开更多
YAG (Y_(3)Al_(5)O_(12)) transparent ceramics have attractive application prospects for transparent armor protection modules because of their excellent light transmittance and anti-ballistic capability. Understanding t...YAG (Y_(3)Al_(5)O_(12)) transparent ceramics have attractive application prospects for transparent armor protection modules because of their excellent light transmittance and anti-ballistic capability. Understanding the fracture behavior and damage mechanism of YAG is necessary for armor design. To explore the damage characteristics of YAG under compression and tension, shock compression and shockless spalling experiments with soft recovery technique are conducted. The spall strength of YAG is obtained and the recovered samples are observed by CT and SEM. It is shown that the macroscopic damage characteristic of YAG under compression is vertical split cracks with oblique fine cracks distributed in the entire sample, while that under tension is horizontal transgranular cracks concentrated near the main spall surface. The cracks generated by macroscopic compression, tension and shear stress extend in similar tensile form at the microscale. The proportion of transgranular fractures on spall surfaces is higher than that of cracks induced by macroscopic compression. Meanwhile, higher loading rate and longer loading duration increase the transgranular fracture percentage.展开更多
In this study,the dynamic compressive response behavior of a body-centered cubic(BCC)single-phase TiZrNbV refractory high-entropy alloy(RHEA)was investigated under impact at speeds of 313-1584 m s^(-1)using two-stage,...In this study,the dynamic compressive response behavior of a body-centered cubic(BCC)single-phase TiZrNbV refractory high-entropy alloy(RHEA)was investigated under impact at speeds of 313-1584 m s^(-1)using two-stage,gas-gun-driven,high-speed plate-impact experiments;recovery sample analysis;and theoretical calculations.The strain rate and pressure were approximately 10^(7) s^(−1) and 5.07-29.37 GPa,respectively.The results showed that the TiZrNbV RHEA had a Hugoniot elastic limit of 4.12-5.86 GPa and a spall strength of 1.84-2.03 GPa.The initial yield strength of the alloy showed a strong strain-rate dependence and could be described by the modified Zerilli-Armstrong model,while the phonon-damping effect was the main reason for its high strain-rate sensitivity.Microstructural analysis showed that the dynamic deformation of the TiZrNbV RHEA was controlled by the dislocation slip,dislocation proliferation,intersection of the deformation bands,and grain refinement.The analysis also showed that the intergranular,transgranular,and mixed-type cracks dominated the spall failure of the material.The dynamic Hall-Petch effect and pinning from the lattice distortion led to high dynamic yield strength.The critical strain rate for the phonon drag effect was positively related to the relative atomic mass and local strain field of the metals.Within the experimental loading range,the RHEA showed good structural stability,and simultaneously,the theoretical calculation method for the equation of state based on a cold-energy mixture could accurately predict its shock-response behavior.The valence-electron concentration(VEC)had a direct effect on the shock-compression properties of the HEAs;higher VEC implied more difficulty in compressing the HEAs.The findings of this study provide insights into understanding the mechanical response characteristics of RHEAs under extreme conditions such as high-speed impact and ultrahigh strain-rate loading.展开更多
Single crystals subjected to shock compression exhibit responses with distinct two-wave structures for certain crystal orientations.However,little is known to date regarding how the shock response depends on crystal o...Single crystals subjected to shock compression exhibit responses with distinct two-wave structures for certain crystal orientations.However,little is known to date regarding how the shock response depends on crystal orientation,and especially why the two-wave structure depends on the crystal orientation.In this work,molecular dynamics simulations of shock compressions in copper single crystals are performed to investigate the orientation dependence of shock responses and the corresponding deformation mechanisms.Four copper single crystals with[001],[011],[012],and[123]crystal orientations along the depth direction are investigated.The[011],[012],and[123]crystal orientations of copper single crystals show distinct two-wave structures in their shock responses,while such a two-wave structure in the shock response is not seen for those orientations having a[001]crystal orientation.The potential causes are analyzed by considering the propagation velocities of both elastic and plastic waves.We develop a technique for identifying twin structures in face-centered cubic crystals and this technique can effectively identify the twin structure.The morphology of shock-induced defects(e.g.,dislocations and twins)shows the significant dependence of crystal orientation and the mechanisms behind these are discussed in detail.Finally,the Johnson-Cook constitutive model describing dynamic deformations at high temperatures and high strain rates is used to analyze the relationships between the shock responses and microscopic defects.The predictions of the Johnson-Cook constitutive model are consistent with the results of the molecular dynamics simulations.展开更多
Shock responses of Mg-Al-Zn alloy are investigated by the molecular dynamics(MD)method.The wave propagation,plastic deformation behavior and failure mechanism along the[0001]and[1010]orientations are analyzed.For both...Shock responses of Mg-Al-Zn alloy are investigated by the molecular dynamics(MD)method.The wave propagation,plastic deformation behavior and failure mechanism along the[0001]and[1010]orientations are analyzed.For both orientations,simulation results show that the shock wave has an obvious double-wave structure(plastic-elastic)under a piston velocity of 1200 m/s.A higher Hugoniot elastic limit(HEL)is observed for[0001]-oriented shock.When the shock pressure is along the[1010]direction,the distance between plastic and elastic waves is closer,and higher dislocation density and more twins are observed.Moreover,the spall strength for[1010]-oriented shock is predicted to be higher.In addition,the wave interactions,HEL and spall strength predicted for Mg-Al-Zn alloy are compared with the experimental results and MD simulation results of Mg single crystal in the literature.It is concluded that the shock performance of Mg-Al-Zn is better than that of Mg single crystal.展开更多
Tungsten has promising applications in high-radiation,high-erosion and high-impact environments.Laser peening is an effective method to enhance the surface mechanical properties of tungsten materials.However,the ultra...Tungsten has promising applications in high-radiation,high-erosion and high-impact environments.Laser peening is an effective method to enhance the surface mechanical properties of tungsten materials.However,the ultrafast dynamic mechanism of defect evolutions induced by laser shockwave in tungsten lattice is unclear.Here,we investigated the evolutions and interactions of various defects under ultrafast compressive process in tungsten lattice using molecular dynamic method.The results confirm the brittleness of tungsten and reveal that void can reduce the yield strain and strength of the tungsten lattice by accelerating defect mesh extension and promoting the dislocation nucleation around itself.Dislocation density is increased with compressive strain rate.Meanwhile,dislocation multiplication and motion reduce the elastic stage and play a dominant role during the plastic deformation of tungsten lattice.Additionally,void can disrupt the dislocation displacement and promote the pinning effect on dislocations by defect mesh extension.展开更多
基金Project supported by the National Natural Science Foundation of China (Nos. 11901555, 11901213,11871448, and 11732016)the National Numerical Windtunnel Project (No. NNW2019ZT4-B10)。
文摘In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.
基金This work is funded by the National Natural Science Foundation of China(No.11772159)the NSAF Joint Fund(No.U1730101)the Fundamental Research Funds for the Central Universities(No.30917011104).
文摘YAG (Y_(3)Al_(5)O_(12)) transparent ceramics have attractive application prospects for transparent armor protection modules because of their excellent light transmittance and anti-ballistic capability. Understanding the fracture behavior and damage mechanism of YAG is necessary for armor design. To explore the damage characteristics of YAG under compression and tension, shock compression and shockless spalling experiments with soft recovery technique are conducted. The spall strength of YAG is obtained and the recovered samples are observed by CT and SEM. It is shown that the macroscopic damage characteristic of YAG under compression is vertical split cracks with oblique fine cracks distributed in the entire sample, while that under tension is horizontal transgranular cracks concentrated near the main spall surface. The cracks generated by macroscopic compression, tension and shear stress extend in similar tensile form at the microscale. The proportion of transgranular fractures on spall surfaces is higher than that of cracks induced by macroscopic compression. Meanwhile, higher loading rate and longer loading duration increase the transgranular fracture percentage.
基金This study was financially supported by the Hunan Provin-cial Natural Science Foundation of China(Grant No.2022JJ10058)the National Natural Science Foundation of China(Grant Nos.12072369 and 52171166).The authors thank Dr.Xuehao Zheng from the ZKKF(Beijing)Science and Technology Company of China for supporting the SEM and TEM analyses.The authors would like to acknowledge Dr.Xiang Wu from KAIPLE Centre for Microscopy,Characterisation&Analysis(CMCA)of China for supporting the EBSD analyzes.
文摘In this study,the dynamic compressive response behavior of a body-centered cubic(BCC)single-phase TiZrNbV refractory high-entropy alloy(RHEA)was investigated under impact at speeds of 313-1584 m s^(-1)using two-stage,gas-gun-driven,high-speed plate-impact experiments;recovery sample analysis;and theoretical calculations.The strain rate and pressure were approximately 10^(7) s^(−1) and 5.07-29.37 GPa,respectively.The results showed that the TiZrNbV RHEA had a Hugoniot elastic limit of 4.12-5.86 GPa and a spall strength of 1.84-2.03 GPa.The initial yield strength of the alloy showed a strong strain-rate dependence and could be described by the modified Zerilli-Armstrong model,while the phonon-damping effect was the main reason for its high strain-rate sensitivity.Microstructural analysis showed that the dynamic deformation of the TiZrNbV RHEA was controlled by the dislocation slip,dislocation proliferation,intersection of the deformation bands,and grain refinement.The analysis also showed that the intergranular,transgranular,and mixed-type cracks dominated the spall failure of the material.The dynamic Hall-Petch effect and pinning from the lattice distortion led to high dynamic yield strength.The critical strain rate for the phonon drag effect was positively related to the relative atomic mass and local strain field of the metals.Within the experimental loading range,the RHEA showed good structural stability,and simultaneously,the theoretical calculation method for the equation of state based on a cold-energy mixture could accurately predict its shock-response behavior.The valence-electron concentration(VEC)had a direct effect on the shock-compression properties of the HEAs;higher VEC implied more difficulty in compressing the HEAs.The findings of this study provide insights into understanding the mechanical response characteristics of RHEAs under extreme conditions such as high-speed impact and ultrahigh strain-rate loading.
基金supported by the National Natural Science Foundation of China(Grant Nos.11972165,and 11502085)the Japan Society for the Promotion of Science(Grant No.P18067)the Fundamental Research Funds for the Central Universities of China(Grant No.2016YXMS097)。
文摘Single crystals subjected to shock compression exhibit responses with distinct two-wave structures for certain crystal orientations.However,little is known to date regarding how the shock response depends on crystal orientation,and especially why the two-wave structure depends on the crystal orientation.In this work,molecular dynamics simulations of shock compressions in copper single crystals are performed to investigate the orientation dependence of shock responses and the corresponding deformation mechanisms.Four copper single crystals with[001],[011],[012],and[123]crystal orientations along the depth direction are investigated.The[011],[012],and[123]crystal orientations of copper single crystals show distinct two-wave structures in their shock responses,while such a two-wave structure in the shock response is not seen for those orientations having a[001]crystal orientation.The potential causes are analyzed by considering the propagation velocities of both elastic and plastic waves.We develop a technique for identifying twin structures in face-centered cubic crystals and this technique can effectively identify the twin structure.The morphology of shock-induced defects(e.g.,dislocations and twins)shows the significant dependence of crystal orientation and the mechanisms behind these are discussed in detail.Finally,the Johnson-Cook constitutive model describing dynamic deformations at high temperatures and high strain rates is used to analyze the relationships between the shock responses and microscopic defects.The predictions of the Johnson-Cook constitutive model are consistent with the results of the molecular dynamics simulations.
基金This research are funded by the National Natural Science Foundation of China(11402183,51604206 and 51974217)the Fundamental Research Funds for the Central Universities of China(WUT:2017IA002)National Defense Science and technology foundation strengthening program.
文摘Shock responses of Mg-Al-Zn alloy are investigated by the molecular dynamics(MD)method.The wave propagation,plastic deformation behavior and failure mechanism along the[0001]and[1010]orientations are analyzed.For both orientations,simulation results show that the shock wave has an obvious double-wave structure(plastic-elastic)under a piston velocity of 1200 m/s.A higher Hugoniot elastic limit(HEL)is observed for[0001]-oriented shock.When the shock pressure is along the[1010]direction,the distance between plastic and elastic waves is closer,and higher dislocation density and more twins are observed.Moreover,the spall strength for[1010]-oriented shock is predicted to be higher.In addition,the wave interactions,HEL and spall strength predicted for Mg-Al-Zn alloy are compared with the experimental results and MD simulation results of Mg single crystal in the literature.It is concluded that the shock performance of Mg-Al-Zn is better than that of Mg single crystal.
基金financially support from the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA25040201)the National Natural Science Foundation of China(Grant No.51727901)support provided by the Deanship of Scientific Research(DSR)at King Fahd University of Petroleum&Minerals(KFUPM)(Grant No.DF201020)
文摘Tungsten has promising applications in high-radiation,high-erosion and high-impact environments.Laser peening is an effective method to enhance the surface mechanical properties of tungsten materials.However,the ultrafast dynamic mechanism of defect evolutions induced by laser shockwave in tungsten lattice is unclear.Here,we investigated the evolutions and interactions of various defects under ultrafast compressive process in tungsten lattice using molecular dynamic method.The results confirm the brittleness of tungsten and reveal that void can reduce the yield strain and strength of the tungsten lattice by accelerating defect mesh extension and promoting the dislocation nucleation around itself.Dislocation density is increased with compressive strain rate.Meanwhile,dislocation multiplication and motion reduce the elastic stage and play a dominant role during the plastic deformation of tungsten lattice.Additionally,void can disrupt the dislocation displacement and promote the pinning effect on dislocations by defect mesh extension.
