OTM/HOTM极限力学仿真在小行星防御中的应用

  • 廖祜明 ,
  • 龚自正 ,
  • 宋光明 ,
  • 樊宗岳 ,
  • 杨宏涛 ,
  • 黎波
展开
  • 1.北京航空航天大学能源与动力工程学院,北京 102206;
    2.北京卫星环境工程研究所,北京 100094;
    3.北京大学工学院,北京 100871;
    4.云翼超算(北京)软件科技有限公司,北京 100027
廖祜明(1987—),男,博士后,主要研究方向:计算力学,无网格方法,流固与热流固耦合,超高速撞击与增材制造仿真,以及航空发动机热端部件结构强度与可靠性设计。电子邮箱:liaohuming@buaa.edu.cn
黎波(1979—),男,长聘副教授,主要研究方向:材料在高温、高压、高应变率等极端热力学条件下动态行为的理论与仿真计算研究,多物理场强耦合数值模拟、先进计算方法、基于物理的多尺度材料建模、断裂力学、增材制造以及超大型并行计算。 电子邮箱:bo.li@pku.edu.cn

网络出版日期: 2023-12-11

OTM/HOTM Extreme Thermomechanical Simulations in Asteroid Defense

  • LIAO Hu-ming ,
  • GONG Zi-zheng ,
  • SONG Guang-ming ,
  • FAN Zong-yue ,
  • YANG Hong-tao ,
  • LI Bo
Expand
  • 1. School of Energy and Power Engineering,Beihang University, Beijng 102206,China;
    2. Beijing Institute of Spacecraft Environment Engineering, Beijing 100094, China;
    3. College of Engineering, Peking University, Beijing 100871, China;
    4. ESCAAS Co., Ltd., Beijing 100027, China

Online published: 2023-12-11

摘要

小行星撞击地球是人类生存面临的重大威胁之一,减缓或避免小行星撞击是学术界和国际社会关注的热点问题。动能撞击、侵彻爆炸或激光烧蚀是当前小行星防御的几个主要手段,其涉及超高速撞击、极限条件下的热-力-化学强耦合等现象。高效高精度预测行星材料在极限热力学条件下的动态响应机制是开展计算仿真的难点,其对计算方法、材料模型与高性能计算都提出了极大的挑战。本文介绍了描述极限条件下材料超大变形、熔化气化、流固与热流固耦合、自由表面、多体接触、碎裂、层裂与碎片云等复杂物理现象的多物理场、多尺度强耦合的OTM/HOTM最优运输无网格系列方法,并初步探讨了基于此方法发展的ESCAAS极限力学仿真平台在小行星防御中的应用,为动能撞击等小行星防御手段提供相关理论与技术支撑。

本文引用格式

廖祜明 , 龚自正 , 宋光明 , 樊宗岳 , 杨宏涛 , 黎波 . OTM/HOTM极限力学仿真在小行星防御中的应用[J]. 空间科学与试验学报, 2022 , 22(3) : 5 -18 . DOI: 10.19963/j.cnki.2096-4099.2022.03.002

Abstract

Asteroid impact on Earth is one of the major threats to human being.Asteroid defense is a hot issue in the academic sector and the international community.Kinetic impact deflecting asteroid orbits or destroying asteroids is one of the main means of asteroid defense.However, these measures involve strong thermal-mechanical-chemical interactive phenomena under extreme conditions such as hypervelocity impact, penetrating explosion and laser ablation.High-efficiency and high-precision prediction of the dynamic response mechanism of materials under extreme thermodynamic conditions is a difficult issue in computational simulation, which poses great challenges to computational methods, material models, and high-performance computing.This paper introduces a series of OTM/HOTM methods for describing complex physical phenomena such as ultra-large deformation, melting, gasification, fluid-structure interaction, thermal-fluid-structure interaction, free surface, multi-body contact, fragmentation, spallation and debris cloud of materials under extreme conditions.Finally, the application of ESCAAS, an extreme mechanics simulation platform developed based on OTM/HOTM method, in asteroid defense measures is preliminarily discussed, which provides relevant theoretical and technical support for asteroid defense methods such as kinetic impact.

参考文献

[1] 柳森,党雷宁,赵君尧,等.小行星撞击地球的超高速问题[J].力学学报,2018,50(6):1311-1327.
[2] NASA Jet Propulsion Laboratory, Center for NEO Studies (CNEOS), Discovery Statistics [EB/OL]. https://cneos.jpl.nasa.gov/stats/totals.html.2022-09-26.
[3] 龚自正,高珂佳,宋光明,等.超高速撞击研究在航天领域的应用现状与展望[J].现代应用物理,2018,9(4):前插1,1-12.
[4] 龚自正,李明,陈川,等.小行星监测预警、安全防御和资源利用的前沿科学问题及关键技术[J].科学通报,2020,65(5):346-372.
[5] 马鹏斌,宝音贺西.近地小行星威胁与防御研究现状[J].深空探测学报,2016,3(1):10-17.
[6] 李飞,孟林智,王彤,等.国外近地小行星撞击地球防御技术研究[J].航天器工程,2015,24(2):87-95.
[7] 刘雪奇,孙海彬,孙胜利.近地小行星防御策略分析[J].深空探测学报,2017,4(6):557-563.
[8] 文德阳.如何防御小行星的来袭[J].科学大观园,2016(12):39-41.
[9] Cheng AF, Rivkin AS, Michel P,et al.AIDA DART asteroid deflection test: Planetary defense and science objectives[J].Planetary and Space Science, 2018, 157: 104-115.
[10] 张韵,刘岩,李俊峰.小行星防御动能撞击效果评估[J].深空探测学报,2017,4(1):51-57.
[11] Jutzi, M.SPH calculations of asteroid disruptions: The role of pressure dependent failure models[J].Planetary and space science, 2015, 107: 3-9.
[12] 姜宇,程彬,宝音贺西,等.潜在威胁小行星碰撞防御的计算与分析[J].深空探测学报,2017,4(2):190-195.
[13] Raducan S, Davison T, Luthe R,et al.The role of asteroid strength, porosity and internal friction in impact momentum transfer[J].Icarus, 2019,329:282-295.
[14] Raducan S, Jutzi M, Davison T, et al.Influence of the projectile geometry on the momentum transfer from a kinetic impactor and implications for the DART mission[J].International Journal of Impact Engineering,2022, 162:104147.
[15] Raducan S, Davison T,Collins G.Ejecta distribution and momentum transfer from oblique impacts on asteroid surfaces[J].2022,Icarus,374:114793.
[16] 马鑫,周琪,李明涛,等.深空撞击载荷总体技术分析与效能仿真[J].中国空间科学技术,2021, 41 (1): 120-130.
[17] 李毅, 陈鸿, 兰胜威, 等.一种提升近地小行星防御中拦截效率的方法[J].航天器环境工程, 2017, 34(6): 585-592.
[18] 汤文辉,张昆,冉宪文.基于核爆炸的小行星偏转方案分析[J].空间碎片研究, 2021, 21 (4): 45-48.
[19] King PK, Syal MB, Dearborn DS, et al.Late-time small body disruptions for planetary defense[J].Acta Astronautica, 2021, 188:367-386.
[20] Michael Owen, J.A compatibly differenced total energy conserving form of SPH[J].International Journal for Numerical Methods in Fluids, 2014,75(11):749-774.
[21] Donea J, Giuliani S, Halleux J P.An arbitrary lagrangian-eulerian finite element method for transient dynamic fluid-structure interactions[J].Computer Methods in Applied Mechanics and Engineering, 1982, 33(1): 689-723.
[22] Donea J, Huerta A, Ponthot JP, et al.Arbitrary lagrangian-rulerian methods[J].Encyclopedia of computational mechanics, 2004, 1: 413-437.
[23] Suetsugu R, Tanaka H, Kobayashi H, et al. Collisional disruption of planetesimals in the gravity regime with iSALE code: Comparison with SPH code for purely hydrodynamic bodies[J].Icarus, 2018, 314: 121-132.
[24] Quan X, Birnbaum N K, Cpwler MS, et al. Numerical simulation of structural deformation under shock and impact loads using a coupled multi-solver approach[C].5th Asia-Pacific Conference on Shock and Impact Loads on Structures. Hunan, China, 2003.
[25] Gingold RA, Minghan JJ.Smoothed particle hydrodynamics: theory and application to non-spherical stars[J].Monthly Notices of the Royal Astronomical Society, 1997, 181(3): 375-389.
[26] Jutzi M, Benz W, Michel P.Numerical simulations of impacts involving porous bodies: I.Implementing sub-resolution porosity in a 3D SPH hydrocode[J].Icarus,2008,198(1):242-255.
[27] Liu WK, Jun S, Zhang YF.Reproducing kernel particle methods[J].International Journal for Numerical Methods in Fluids, 1995, 20(8-9): 1081-1106.
[28] Schwartz SR, Richardson DC,Michel P.An implementation of the soft-sphere discrete element method in a high-performance parallel gravity tree-code[J]. Granular Matter, 2012,14(3):363-380.
[29] Zhang X, Chen Z, Liu Y.The material point method[M].Oxford: Academic Press, 2017.
[30] 林健宇, 罗斌强, 徐名扬, 等.铝弹丸超高速撞击防护结构的研究进展[J] , 高压物理学报, 2019, 33 (3): 030112-1.
[31] Li Bo, Feras Habbal, Michael Ortiz.Optimal transportation meshfree approximation schemes for fluid and plastic flows[J].International Journal for Numerical Methods in Engineering,2010,83 (12): 1541-1579.
[32] Hao Wang, Huming Liao, Zongyue Fan,et al.The Hot Optimal Transportation Meshfree (HOTM) method for materials under extreme dynamic thermomechanical conditions[J].Computer Methods in Applied Mechanics and Engineering, 2020, 364: 112958.
[33] Li Bo, Stalzer M,Ortiz M.A massively parallel implementation of the Optimal Transportation Meshfree method for explicit solid dynamics[J].International Journal for Numerical Methods in Engineering, 2014,100(1):40-61.
[34] Li Bo, Pandolfi ANNA,Ortiz M.Material-point erosion simulation of dynamic fragmentation of metals[J].Mechanics of Materials, 2015,80: 288-297.
[35] Yang Q, Stainier L, Ortiz M.A variational formulation of the coupled thermos-mechanical boundary-value problem for general dissipative solids[J].J Mech Phys Solids,2006, 54(2): 401-424.
[36] 廖祜明.整体拉格日无网格流固耦合计算方法[D]. 北京: 北京航空航天大学,2018.
[37] Arroyo M, Ortiz M.Local maximum-entropy approximation schemes: a seamless bridge between finite elements and meshfree methods[J].Int J Numer Methods Eng,2006, 65(13): 2167-2202.
[38] Fan J, Liao H, Wang H, et al.Local maximum-entropy based surrogate model and its application to structural reliability analysis[J].Struct Multidiscipl Optim,2018, 57(1): 373-392.
[39] Schmidt B, Fraternali F, Ortiz M.Eigenfracture: an eigendeformation approach to tariational fracture[J]. Multiscale Modeling and Simulation, 2009, 7(3): 1237-1266.
[40] DJ Steinberg.Equation of State and Strength Properties of Selected Materials.Technical report[J].Lawrence Livermore National Laboratory, 1996.
[41] 樊江, 袁圆, 廖祜明, 等.基于最优运输无网格法的Whipple 屏超高速撞击数值模拟[J].爆炸与冲击, 2019, 40(7): 94-104.
[42] 廖祜明, 黎波, 樊江, 等.超高速撞击下碎片云的OTM分析[J].爆炸与冲击,2022,42(10):50-60.
[43] 张天龙, 马天宝, 郝莉.基于pOTM方法的超高速碰撞并行数值模拟研究[J].兵器装备工程学报, 2021, 42 (1): 144-149.
[44] Li B, et al.Large scale optimal transportation meshfree (OTM) simulations of hypervelocity impact [J]. Procedia Engineering, 2013, 58: 320-327.
[45] Jiang, Hao, et al."Numerical modeling of compressive failure mechanisms in ceramic materials at high strain rates[J].Computer Methods in Applied Mechanics and Engineering,2019, 347 : 806-826.
[46] Jiang, Hao, Valerie Scott, Bo Li.Modeling and simulations of high and hypervelocity impact of small ice particles[J].International Journal of Impact Engineering,2021, 155: 103906.
[47] 龚自正, 徐坤博, 牟永强, 等.空间碎片环境现状与主动移除技术[J].航天器环境工程, 2014, 31: 129-135.
[48] 陈川, 龚自正, 杨武霖, 等.空间碎片几何形状对激光烧蚀冲量的影响规律[J].高压物理学报, 2018, 3: 040101.
[49] Chen C, Yang W L, Gong Z Z, et al.Research progress in laser active debris removal of CAST[J].Aerospace China, 2018, 19: 3-19.
[50] Phipps C R, Bonnal C.A spaceborne, pulsed UV laser system for re-entering or nudging LEO debris, and re-orbiting GEO debris[J].Acta Astronaut, 2016, 118: 224-236.
[51] Kubiak M, Piekarska W, Stano S.Modelling of laser beam heat source based on experimental research of yb: Yag laser power distribution[J].Int J Heat Mass Transf.2015, 83: 679-689.
[52] Anisimov SI, Khokhlov VA.Instabilities in Laser-Matter Interaction [M]. Boca Raton, FL: CRC Press, 1995.
[53] Fan Z, Wang H, Huang Z, et al.A Lagrangian meshfree mesoscale simulation of powder bed fusion additive manufacturing of metals[J].Int J numer Methods Eng.2021, 122(2): 483-514.
[54] Yang, Gang, et al.Feasibility analysis of SPH method in the simulation of condensed explosives detonation with ignition and growth model[J].Computers & Fluids,2013, 88: 51-59.
[55] Liu MB, et al.Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology[J].Computers & Fluids,2003, 32(3): 305-322.
文章导航

/