Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of local structure and stoichiometry on the dynamic behavior of bi-metal interfaces

Molecular dynamics (MD) simulations are performed to investigate the failure behavior of Al/Al, Al/Fe, and Al/AlFe bi-layers under high strain-rate shock loading. Additional simulations are performed to invesigate solely the tensile response of these systems and to isolate the effects of loading history on damage and failure. Specifically, the compression stage of the shock is excluded from the tensile simulations. This study shows that local structure variation, including the introduction of serrations, plays different roles in controlling damage nucleation depending on the specific loading condition. Under shock loading, the stress for void nucleation and the resulting void distribution is insensitive to an interface structure, whereas under pure tensile loading, the opposite is true, even though the strain rate is comparable. The above difference can be explained based on the corresponding change in the total deformation prior to damage nucleation under shock loading, which is missing from the pure tensile loading. Therefore, it is concluded that whether local interface structure matters under high strain-rate loading depends on the specific loading history and the stress state: the local structure variation matters only when it alters the deformation behavior, and subsequently, damage nucleation and evolution.

36 MATERIALS SCIENCE↗

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗

Modeling of the thermo-mechanical response and texture evolution of WE43 Mg alloy in the dynamic recrystallization regime using a viscoplastic self-consistent formulation

This paper presents a microstructure sensitive model for predicting mechanical response and texture evolution of metals in the dynamic recrystallization regime. A recently proposed viscoplastic self-consistent (VPSC) formulation for the prediction of recrystallization driven by strain energy and intragranular misorientation is extended to hexagonal close-packed (hcp) metals. The model is applied to the dynamic recrystallization of magnesium alloy WE43 at different temperatures and strain rates. Model predictions in terms of stress-strain response and texture evolution are compared to the experimental measurements and acceptable agreement is achieved. According to the model predictions, superplastic behavior of nuclei was found to be the dominant softening mechanism at high temperatures and low strain rates. High concentration of precipitates at the grain boundaries and presence of alloying elements are the likely causes of low boundary mobility, resulting in nucleation dominated dynamic recrystallization. Relatively strong basal compression textures indicate dominant activity of basal slip, which can be achieved only through large difference in slip resistance between soft basal and hard prismatic and pyramidal modes.

36 MATERIALS SCIENCE↗

Optimal measurement of ice-sheet deformation from surface-marker arrays

Surface strain rate is best observed by fitting a strain-rate ellipsoid to the measured movement of a stake network or other collection of surface features, using a least squares procedure. Error of the resulting fit varies as 1/(L delta t square root of N), where L is the stake separation, delta is the time period between initial and final stake survey, and n is the number of stakes in the network. This relation suggests that if n is sufficiently high, the traditional practice of revisiting stake-network sites on successive field seasons may be replaced by a less costly single year operation. A demonstration using Ross Ice Shelf data shows that reasonably accurate measurements are obtained from 12 stakes after only 4 days of deformation. It is possible for the least squares procedure to aid airborne photogrammetric surveys because reducing the time interval between survey and re-survey permits better surface feature recognition.

Macayeal, D. R.↗

Optimal measurement of ice-sheet deformation from surface-marker arrays

Surface strain rate is best observed by fitting a strain-rate ellipsoid to the measured movement of a stake network or other collection of surface features, using a least squares procedure. Error of the resulting fit varies as 1/(L delta t square root of N), where L is the stake separation, delta is the time period between initial and final stake survey, and n is the number of stakes in the network. This relation suggests that if n is sufficiently high, the traditional practice of revisiting stake-network sites on successive field seasons may be replaced by a less costly single year operation. A demonstration using Ross Ice Shelf data shows that reasonably accurate measurements are obtained from 12 stakes after only 4 days of deformation. It is possible for the least squares procedure to aid airborne photogrammetric surveys because reducing the time interval between survey and re-survey permits better surface feature recognition.

Macayeal, D. R.↗

Determining damping mechanisms in a composite beam by experimental modal analysis

A method of estimating the distributed damping parameters of a beam based on measured modal parameters (frequency and damping ratios) is derived. Three different mathematical models were used to model the damping mechanism of a quasi-isotropic pultruded cantilevered beam. These three models were (1) viscous (air) damping, (2) strain-rate damping, and (3) both viscous and strain-rate damping. Values were obtained for each of the damping models as well as the modulus of elasticity of the beam. It is found that the two-parameter damping model provides the best fit to measure modal data. However, the two-parameter damping model can only reproduce the measured damping ratios to within 85 percent.

Cudney, H. H.↗

1300 K compressive properties of a reaction milled NiAl-AlN composites

When B2 crystal-structure nickel aluminide is subjected to high-intensity mechanical ball milling in a liquid nitrogen bath, or 'cryomilling', an NiAl composite is obtained which contains about 10 vol pct AlN particles. This composition arises from the incorporation of N during cryomilling; during subsequent thermomechanical processing, the N reacts with Al. While compressive testing of extruded or isostatically pressed specimens at 1300 K indicated that strength at relatively fast strain rates is slightly dependent on consolidation method, slower strain rates indicate no clear dependency on densification technique: four different consolidation methods were found to yield similar creep strengths. The creep properties of NiAl-AlN are similar to those of the single-crystal Ni-base superalloy NASAIR 100.

Whittenberger, J. Daniel↗

Microstructure and Mechanical Properties of Extruded Gamma Met PX

A gamma TiAl alloy with a high Nb content is being assessed as a compressor blade material. The microstructure and mechanical properties of extruded Ti-45Al-X(Nb,B,C) (at %) were evaluated in both an as-extruded condition and after a lamellar heat treatment. Tensile behavior of both as-extruded and lamellar heat treated specimens was studied in the temperature range of RT to 926 C. In general, the yield stress and ultimate tensile strength reached relatively high values at room temperature and decreased with increasing deformation temperature. The fatigue strength of both microstructures was characterized at 650 C and compared to a baseline TiAl alloy and to a Ni-base superalloy. Tensile and fatigue specimens were also exposed to 800 C for 200 h in air to evaluate the alloy's environmental resistance. A decrease in ductility was observed at room temperature due to the 800 C exposure but the 650 C fatigue properties were unaffected. Compressive and tensile creep testing between 727 and 1027 C revealed that the creep deformation was reproducible and predictable. Creep strengths reached superalloy-like levels at fast strain rates and lower temperatures but deformation at slower strain rates and/or higher temperature indicated significant weakening for the as-extruded condition. At high temperatures and low stresses, the lamellar microstructure had improved creep properties when compared to the as-extruded material. Microstructural evolution during heat treatment, identification of various phases, and the effect of microstructure on the tensile, fatigue, and creep behaviors is discussed.

Draper, S. L.↗

Microstructure and Mechanical Properties of Extruded Gamma Microstructure Met PX

A gamma TiAl alloy with a high Nb content is being assessed as a compressor blade material. The microstructure and mechanical properties of extruded Ti-45Al-X(Nb,B,C) (at.%) were evaluated in both an as-extruded condition and after a lamellar heat treatment. Tensile behavior of both as-extruded and lamellar heat treated specimens was studied in the temperature range of RT to 926 C. In general, the yield stress and ultimate tensile strength reached relatively high values at room temperature and decreased with increasing deformation temperature. The fatigue strength of both microstructures was characterized at 650 C and compared to a baseline TiAl alloy and to a Ni-base superalloy. Tensile and fatigue specimens were also exposed to 800 C for 200 h in air to evaluate the alloy's environmental resistance. A decrease in ductility was observed at room temperature due to the 800 C. exposure but the 650 C fatigue properties were unaffected. Compressive and tensile creep testing between 727 and 1027 C revealed that the creep deformation was reproducible and predictable. Creep strengths reached superalloy-like levels at fast strain rates and lower temperatures but deformation at slower strain rates and/or higher temperature indicated significant weakening for the as-extruded condition. At high temperatures and low stresses, the lamellar microstructure had improved creep properties when compared to the as-extruded material. Microstructural evolution during heat treatment, identification of various phases, and the effect of microstructure on the tensile, fatigue, and creep behaviors is discussed.

Draper, S. L.↗

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 6061-T651 Aluminum Alloy from Plate

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

36 MATERIALS SCIENCE↗

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗

Rate and time dependent failure of structural adhesives

Studies on two adhesives (Metlbond 1113 and 1113-2) identified as having important applications in the bonding of composite materials are presented. A testing program to ascertain stress-strain, strain-rate, time, yield, and/or failure behavior of these materials in bulk form using uniaxial tensile constant strain-rate, creep, and relaxation tests is described. The stress-strain behavior of each material is shown to be significantly rate dependent. A rate dependent stress whitening (crazing) phenomenon occurs prior to either yield or fracture. A region of linear elasticity, a region of viscoelasticity, and the onset of yielding are identified in the stress-strain behavior. The linear elastic limit and the yield point are shown to be rate dependent and agree well with an empirical equation proposed by Ludwik. A creep to failure phenomenon is shown to exist and is correlated with a delayed yield equation proposed by Crochet. Analytical predictions based on a modified Bingham model are shown to agree well with experimental stress-strain strain-rate data. Analytical predictions based on a modified Ramberg-Osgood equation are also shown for comparison purposes.

Brinson, H. F.↗

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

The hypothesis is explored that the crust-mantle boundary of Venus is not in phase equilibrium but rather is rate-limited by the temperature-dependent volume diffusion of the slowest ionic species. The 1D thermal evolution problem is solved assuming that the mountains formed by uniform horizontal shortening of the crust and the lithospheric mantle at a constant rate. The time-dependent density structure and surface elevation are calculated by assuming a temperature-dependent reaction rate and local Airy isostatic compensation. For a horizontal strain rate of 10 exp -15/s or greater, the temperature increase at the base of the crust during mountain formation is modest to negligible, the deepening lower crust is metastable, and the surface elevation increases as the crust thickens. For strain rates less than 10 exp -16/s, crustal temperature increases with time because of internal heat production and the lower crust is more readily transformed to the dense eclogite assemblage. For such models, a maximum elevation is reached during crustal shortening.

Namiki, Noriyuki↗

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the viscoplastic behavior of a semicrystalline polymer

In this study, a complex constitutive relation is identified using inverse modeling with the nominal mechanical response as sole experimental input. The methodology is illustrated for a semicrystalline thermoplastic in the presence of strain localization at finite deformations. The experimental database includes cylindrical tensile bars, compression pins and round notched bars loaded at strain rates spanning up to five decades and temperatures below and above T g . The data is organized into a calibration set and a validation set. The response of tensile specimens is determined using finite element analyses and a two-phase constitutive relation for semicrystalline polymers that accounts for temperature- and rate-sensitive plastic flow, pressure-sensitivity, small-strain softening and large-strain orientational hardening of the amorphous phase, along with the evolution of crystallinity. The large number of constitutive parameters is identified using an optimization tool coupled with the finite element solver and the calibration set from experiments. The methodology is shown to be successful in predicting the response of round notched bars and replicating the effects of temperature and strain rate on the severity of necking in tensile bars. The proposed model identification strategy is both simple and effective in comparison with other elaborate methods that attempt to access intrinsic behavior directly from high-fidelity experimental measurements.

36 MATERIALS SCIENCE↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗