Bland, PA

Nazarova, K

Spiral - Imperial College Digital Repository

NUMERICAL MODELLING FOR STRENGTH ESTIMATION OF FRAGMENTING METEOROIDS.  K. Nazarova and P. A. Bland. Impacts and Astromaterials Research Centre (IARC), Department of Earth Science and Engineering, Imperial College London, South Kensington Campus, SW7 2AZ, London, UK, e-mail: katsi-aryna.nazarava@imperial.ac.uk.   Introduction:  It has been estimated that the total number of meteoroids entering the Earth’s atmosphere with masses above 1 kg is about 37000 per year [1]. And approximately 4400 events with mass more than 1 kg strike the Earth every year [1]. Among these im-pactors, those in the mass range ~104-3.106 kg frag-ment and form crater fields, with the largest crater resulting from impact of the largest fragment [2]. Fragmented impactors in the mass range over 106 – 108 kg strike together to form single craters [2]. Ac-cording to previous estimates [2], eleven 105 kg events every 100 years will impact and form a crater field. It is possible to numerically model the formation of many terrestrial crater fields. The solutions, however, are not necessarily unique. A large uncertainty in nu-merical modelling is in the altitude of breakup, which is related to the yield strength of the meteoroid. Since the cross-range width of a given crater field depends primarily upon the breakup altitude and the down-range length is related to the angle of entry, it is pos-sible to approximately derive these parameters from the terrestrial crater fields.  Modelling: Fragmentation of meteoroids by aero-dynamic forces is closely connected with their strength – improving models of meteoroid strength is funda-mental to understanding of meteoroid fragmentation and to improving our understanding of impact rates at planetary surfaces [2, 3, 4, 5].  Meteorite strength. The comparison of the meas-ured and estimated values of material strength of me-teoroids was performed.  We estimated the strength of the main body (1 kg sample) for meteorites Lost City, Pribram, Innisfree, Sikhote-Alin, Benesov, Moravka, and Sumava according to Weibull distribution [6] and using the criterion of the meteoroid breakup [7]. Our estimates are in the same range as the measured val-ues of strength for stony and iron impactors. Maxi-mum strength of the fragmenting body was 4.4·107 H/m2 for iron and 1.5·107 H/m2 for stony bolides. However, some of the bodies have strength as low as 0.13·106 H/m2. There is a great uncertainty and scatter in the measured strength: even for different samples of a common meteorite, values of ultimate strength may differ by a factor of three [8]. Moreover, although irons usually have higher ultimate strength than stones, it is not the case for every meteorite. However, by ignoring the extreme values, the ultimate tensile strength of meteorites can be narrowed to a range from 2·107 to 5·107 H/m2 (the average value is about 3.5·107 H/m2). It should be emphasized that the ulti-mate strength of a large specimen can be much lower than that of its constituents (Weibull law). The ulti-mate tensile strength for polycrystallic specimens of Sikhote-Alin meteorite was about 4.4·107 H/m2, whereas for their constituent monocrystals it was 10 times higher [9]. Modelling terrestrial crater fields. Here we used the model of separate fragments (SF [4, 5, 10, 11] ) to model the size distribution of craters and the size of the Morasko crater field. We performed 30 sets of SF model simulations. The initial flight parameters for the modelling were pre-entry mass, 100 t; initial ve-locity, 13 km/s; trajectory angle, 40˚; drag coefficient, 0.8 [12] for spherical shape of the impactor; heat transfer coefficient, 0.1; strength coefficient, 0.09, and strength of the impactor material, 4.4•107 H/m2. In our runs the mass of the main fragment varied from 50 to 200 t. The best results are obtained for the mass range 100 – 200 t, the velocity range 13 – 14 km/s, and the trajectory angle range 30 - 35˚. The density of the body is assumed to be 7700 kg/m3.  Figure 1 demonstrates the schematic crater field of the Morasko meteoroid, where the radii of the circles correspond to the size of the craters, and an example of the modelled field. The size of the crater field agrees with the observations; however the positions of the individual fragments may vary as they were de-fined by the random procedure. The final velocities of the fragments are in the range from 3 to 7 km/s.   Lunar and Planetary Science XXXVII (2006) 1825.pdf Figure 1. Modelled (above) and real (below) schematic map of the Morasko crater field (all distances are in me-ters).  The results obtained in different runs with the same initial data do not coincide because of the use of random procedures in the process of calculations. The impactors having the same density, velocity, and shape would fragment in a different way, because of the difference in their structure (defect distribution). Moreover, the real strength of the impactor and its fragments is a unique value and can deviate from the strength defined by statistical strength theory [5, 6]. These deviations can result from the features of frag-ment shape (drag coefficient, aerodynamical loading), flow conditions, etc. To model this idea, the strength of the fragments can be described with the normal statistical distribution, where the mean value is de-fined from Weibull law and dispersion δ=1-2 [5]. Unlike the Sikhote-Alin iron meteorite shower, most of the craters at Morasko probably have been destroyed by agriculture, and the meteorites buried. That is why on the real map (Figure 1) there are no small craters. It has become impossible to recognize the ploughed craters at Morasko. Therefore, it is im-possible to collect the meteorites in a systematic man-ner such as at Sikhote-Alin. An important outcome of the SF model simula-tions is the constraint they place on the crater field formation and impact rate. The simulated values of mass, velocity and angle for small terrestrial craters allow us to understand the processes of fragmentation and formation of strewn fields. Also, extending the model to other solid planets with atmospheres, it is expected that crater fields on Mars are at least an or-der of magnitude smaller than terrestrial crater fields (assuming that the present Martian atmospheric con-ditions have persisted throughout the period of mete-orite bombardment [4, 5]). For Venus, one would ex-pect to find crater fields that are approximately an order of magnitude larger than crater fields on Earth [4, 10].  References: [1] Halliday I., Griffin A. A., Black-well A.T. (1996) M&PS, 31, 185-217. [2] Bland P.A. and Artemieva N.A. (2003) 3d Int’l Conference on LMI, abstract 4047. [3] Passey Q.R. and Melosh H.J. (1980) Icarus, 42, 211–233. [4] Bland P.A. and Ar-temieva N.A. (2003) Nature, 424, 288-291. [5] Ar-temieva N. A. and Shuvalov V. V. (2001) J. Geophys. Research, 106, 3297-3309. [6] Weibull W. (1951) J. Apple. Mech., 18, 140–147. [7] Svetsov V. V. and Nemchinov I. V. (1995) Icarus, 116, 131-153. [8] Medvedev R.V. (1974) Meteoritica, 33, 100-104. [9] Yavnel’ A.A. (1963) The Sikhote-Alin Iron Mete-orite shower (V.G. Fesnkov, ed.)  2, 362-371. [10] Artemieva N.A. and Bland P.A. (2003) 34th LPSC. CD-ROM. [11] Artemieva N. A. and Shuvalov V. V. (1996) Shock Waves, 5, 359 [12] Chyba C.F., Thomas P J. and Zahnle K.J. (1993) Nature, 361, 40-44.  Lunar and Planetary Science XXXVII (2006) 1825.pdf

Numerical modeling for strength estimation of fragmenting meteoroids

Abstract

Similar works

Full text

Available Versions

Spiral - Imperial College Digital Repository