The construction of models for the internal constitution and the temporal evolution of large planetesimals, the parent bodies of chondrites, requires information on the heat conductivity of the complex mixture of minerals and iron metal found in chondrites. It is attempted to evaluate the heat conductivity of a multi-component mineral mixture and granular medium from the heat conductivities of its mixture components. Random mixtures of solids with chondritic composition and packings of spheres are numerically generated. The heat conduction equation is solved in high spatial resolution for a test cube filled with such matter. From the heat flux through the cube the heat conductivity of the mixture is derived. The model results for porous material are consistent with data for compacted sandstone, but are at odds with measurements for H and L chondrites. The discrepancy is traced back to shock modification of the currently available meteoritic material by impacts on the parent body over the last 4.5 Ga. This causes numerous micro-cracks that act as additional barriers for heat transfer. The void structure in meteorites is different from that which probably existed in the pristine material of the parent bodies. The results obtained for the heat conductivity of the pristine material are used for calculating models for the evolution of the H chondrite parent body which are fitted to the cooling data of a number of H chondrites. The fit to the data good.

Reconstruction of the thermal history of individual meteorites which can be assigned to the same parent body allows to derive general characteristics of the parent body, which hold important clues on the planetary formation process. This requires to construct models for the heating of such bodies by short lived radioactives, in particular by 26Al, and its cooling by heat conduction, which then are compared with the cooling histories of the meteorites. The heat conductivity of the planetesimal material depends critically on the porosity of the chondritic material and changes by sintering of the material at elevated temperatures and pressures. Compaction of a granular material is a key process for the thermal history of the parent bodies of meteorites. The modelling of the compaction process is improved by applying concepts originally developed for hot isostatic pressing in metallurgical processes. It is extended to a binary mixture of matrix and chondrules, as observed in chondrites. By comparison with published data on sintering experiments it is shown that the algorithm allows a sufficiently accurate modelling of the compaction of silicate material. It is shown that the characteristic temperature at which sintering occurs is different for matrix or chondrule dominated precursor material. We apply the new method to model the thermal evolution of the parent body of the H chondrites and determine an improved optimized set of model parameters for this body.

The cooling histories of individual meteorites can be empirically reconstructed by using ages from different radioisotopic chronometers with distinct closure temperatures. For a group of meteorites derived from a single parent body such data permit the reconstruction of the cooling history and properties of that body. Particularly suited are H chondrites because precise radiometric ages over a wide range of closure temperatures are available. A thermal evolution model for the H chondrite parent body is constructed by using all H chondrites for which at least three different radiometric ages are available. Several key parameters determining the thermal evolution of the H chondrite parent body and the unknown burial depths of the H chondrites are varied until an optimal fit is obtained. The fit is performed by an 'evolution algorithm'. Empirical data for eight samples are used for which radiometric ages are available for at least three different closure temperatures. A set of parameters for the H chondrite parent body is found that yields excellent agreement (within error bounds) between the thermal evolution model and empirical data of six of the examined eight chondrites. The new thermal model constrains the radius and formation time of the H chondrite parent body (possibly (6) Hebe), the initial burial depths of the individual H chondrites, the average surface temperature of the body, the average initial porosity of the material the body accreted from, and the initial 60Fe content of the H chondrite parent body.

The major aim of this study is to assess the effects of sintering of initially porous material on the thermal evolution of planetesimals, and to constrain the values of basic parameters that determined the structure and evolution of the H chondrite parent body. A new code is presented for modeling the thermal evolution of ordinary chondrite parent bodies that initially are highly porous and undergo sintering by hot pressing as they are heated by decay of radioactive nuclei. The pressure and temperature stratification in the interior of the bodies is calculated by solving the equations of hydrostatic equilibrium and energy transport. The decrease of porosity of the granular material by hot pressing due to self-gravity is followed by solving a set of equations for the sintering of powder materials. For the heat conductivity of granular material we combine recently measured data for highly porous powder materials, relevant for the surface layers of planetesimals, with data for heat conductivity of chondrite material, relevant for the strongly sintered material in deeper layers. To demonstrate the capability of our new model, the thermal evolution of the H chondrite parent body was reconstructed. The model starts with a porous body that is later compacted first by 'cold pressing' at low temperatures and then by 'hot pressing' for temperatures above \approx 700 K, i.e., the threshold temperature for sintering of silicates. The thermal model was fitted to the well constrained cooling histories of the two H chondrites Kernouve (H6) and Richardton (H5).