Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solidification”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Laser spot welding of additive manufactured 304L stainless steel

Here, the goal of this work is to understand if an additively manufactured 304L stainless steel exhibits similar spot-welding behavior as wrought 304L stainless steel. Due to the many differences between an additively manufactured component and wrought product, it is important to determine how the material interacts with the laser and how it affects the weld bead morphology. In this paper, the laser coupling efficiency, weld size, and solidification of spot welds produced in wrought and additively manufactured 304L stainless steel were investigated. The coupling efficiency of wrought and additively manufactured 304L stainless steel of similar surface condition were approximately the same over a range of applied laser energies. Laser welding of the untreated (rougher) surface of additively manufactured 304L, however, showed improved coupling efficiency ranging between 3.3 and 100%. The rougher surface traps the incoming light and increases the coupling efficiency at lower laser energies, while at higher energy, the absorption efficiency is dominated by intrinsic absorption from the keyhole formation rather than surface roughness. The resulting spot weld microstructures differed from welds made in wrought 304L and additively manufactured 304L. Welds made in wrought 304L were fully austenitic containing what is suspected to be massive austenite, which suggests that these welds solidified as primary ferrite. Welds made in additively manufactured 304L were also fully austenitic and contained both cellular austenite and what is suspected to be massive austenite. These observations mean that welds made in additively manufactured 304L solidified as primary ferrite and primary austenite. The differences in weld microstructures made in wrought and AM 304L can be attributed to differences in the composition and solidification rate.

304L stainless steel↗
DOE OSTI2018

Strengthening mechanisms in directed energy deposited austenitic stainless steel

Microstructures and mechanical properties are evaluated in austenitic stainless steel structures fabricated by directed energy deposition (DED) considering the effects of applied loading orientation, build geometry, and distance from the deposition baseplate. Locations within an as-deposited build with different thermomechanical history display different yield strength, while those locations with similar history have approximately the same yield strength, regardless of test specimen orientation. Thermal expansion of deposited material near the baseplate is inhibited by the mechanical constraint imposed by the baseplate, promoting plastic deformation and producing a high density of dislocations. Concurrently, high initial cooling rates decrease away from the baseplate as the build is heated, causing an increased spacing of cellular solidification features. An analysis of strengthening mechanisms quantitatively established for the first time the important strengthening contribution of high dislocation densities in the materials (166–191 MPa) to yield strength that ranged from 438 to 553 MPa in the present DED fabricated structures. A newly adopted mechanistic relationship for microsegregation strengthening from the literature indicated an additional important contribution to strengthening (123–135 MPa) due to the cellular solidification features. These findings are corroborated by the measured evolution of microstructure and hardness caused by annealing the DED material. Finally, these results suggest that the mechanical properties of deposited austenitic stainless steels can be influenced by controlling thermomechanical history during the manufacturing process to alter the character of compositional microsegregation and the amount of induced plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An additively manufactured near-eutectic Al-Ce-Ni-Mn-Zr alloy with high creep resistance

A new additively manufactured (AM) Al-7.5Ce-4.5Ni-0.4Mn-0.7Zr (wt.%) near-eutectic alloy is reported, which shows unprecedented creep resistance up to 400 °C (a homologous temperature of 0.72). The eutectic solidification microstructure comprises ~ 27 vol% of coarsening-resistant second phase network with an ultrafine (<100 nm) inter-phase spacing. Both Mn and Zr contribute to creep resistance of the alloy. Small amount of Mn addition promotes selection of coarsening resistant phases without compromising the alloy processability. Zr not only improves hot-tearing resistance, but further enhances the second phase coarsening resistance resulting in improved creep resistance. Neutron diffraction performed during creep deformation reveals that the underlying mechanism for creep resistance in this alloy is impedance to dislocation motion stemming from the ultrafine eutectic solidification microstructure, whereas load transfer strengthening becomes less effective as the creep temperature increases. The second phase forms a continuous network in the as-fabricated condition, which is maintained during long-term creep at 300 °C. However, this network is fragmented into fine dispersoids at higher temperatures. It is proposed that the rate-limiting deformation mechanism at 300–400 °C is (i) dislocation climb for the alloy with fragmented second phase dispersoids and (ii) Orowan looping for the alloy with a continuous second phase network. In conclusion, the present design of an AM-processable multicomponent eutectic alloy with high creep resistance can be applied to other metallic systems exhibiting eutectic reactions, with expected extreme creep resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory-guided design of duplex-phase multi-principal-element alloys

Density-functional theory (DFT) is used to identify phase-equilibria in multi-principal-element and high-entropy alloys (MPEAs/HEAs), including duplex-phase and eutectic microstructures. Here, a combination of composition-dependent formation energy and electronic-structure-based ordering parameters were used to identify a transition from FCC to BCC favoring mixtures, and these predictions experimentally validated in the Al-Co-Cr-Cu-Fe-Ni system. A sharp crossover in lattice structure and dual-phase stability as a function of composition were predicted via DFT and validated experimentally. The impact of solidification kinetics and thermodynamic stability was explored experimentally using a range of techniques, from slow (castings) to rapid (laser remelting), which showed a decoupling of phase fraction from thermal history, i.e., phase fraction was found to be solidification rate-independent, enabling tuning of a multi-modal cell and grain size ranging from nanoscale through macroscale. Strength and ductility tradeoffs for select processing parameters were investigated via uniaxial tension and small-punch testing on specimens manufactured via powder-based additive manufacturing (directed-energy deposition). This work establishes a pathway for design and optimization of next-generation multiphase superalloys via tailoring of structural and chemical ordering in concentrated solid solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Additively-manufactured Al-0.3Zr-0.2Ce-0.2Cu alloy with high creep resistance and electrical conductivity

Here, a new, solute-lean Al-0.3Zr-0.2Ce-0.2Cu (wt.%) alloy is developed for additive manufacturing that overcomes the classical tradeoff between conductivity and creep resistance. The rapid-cooling-enabled supersaturation of Zr, and its uniform distribution in α-Al matrix, along with formation of submicron (Ce,Cu)-rich intermetallic particles on solidification lead to unusually high creep resistance at 200 °C. Near-zero secondary creep rates are achieved up to the alloy yield stress (YS) of 65 MPa at 200 °C in as-fabricated state. The Zr-solute-induced dislocation-climb suppression mechanism underlying this improvement also restricts dynamic recovery above YS, as noted from appreciable primary creep and its transitioning to near-zero secondary creep rates. A combination of relatively coarse, epitaxially-grown α-Al grains, low Zr concentration in α-Al, and the impurity-scavenging effect of Ce to purify α-Al matrix produces high electrical conductivity of ∼48 %IACS. Aging precipitation of L1 2 -Al 3 Zr nanoprecipitates doubles the YS (to ∼150 MPa) at room temperature and increases alloy conductivity to ∼58 %IACS, but loss of solid-solution Zr out of α-Al matrix leads to activation of dislocation climb, degrading the creep properties as compared to the supersaturated Al-Zr solid solution in the as-fabricated state. Compared to L1 2 -Al 3 Zr nanoprecipitates, submicron (Ce,Cu)-rich particles formed on solidification are more effective at impeding dislocation climb, producing a threshold stress for dislocation creep of ∼ 50 MPa at 200 °C. The new alloy design concepts, especially solute-induced dislocation-climb suppression for creep resistance, explored here may pave way for the design of new metallic alloys for thermal/electrical conductors and other high-temperature applications.

Additive Manufacturing↗
DOE OSTI2025

Nanoscale interfacial melting enables bonding during high velocity microparticle impacts

Melting during high-velocity particle impact has been understood to be typically detrimental to bonding by lowering the strength at the interface and promoting rebound before solidification can occur. Here we establish a possible remedy to this challenge: by dramatically restricting the volume of molten material, its resolidification is accelerated, effectively forming a nanoscale, braze-type joint during impact. In-situ single particle impact imaging is combined with post-mortem structural and chemical analyses to reveal a regime where adhesion is governed not by extensive plastic deformation, but by the kinetics of melt layer resolidification. Furthermore, these findings redefine the role of melting in impact-based processes, establishing transient melting and rapid solidification as a viable strategy for engineering successful adhesion events.

Additive manufacturing↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Synthesis of precipitation-strengthened Al-Sc, Al-Zr and Al-Sc-Zr alloys via selective laser melting of elemental powder blends

Selective laser melting is used to create Al-1.5Sc, Al-1.5Zr and Al-0.75Sc-0.75Zr (at.%) alloys from blends of elemental Al, Sc, and Zr powders. This study investigates elemental alloying elements (Sc and Zr) which are high-melting and highly reactive, unlike previous work which focused on more concentrated elemental additions of lower-melting, lower-reactivity Cu and Si to aluminum. High-speed in situ synchrotron x-ray imaging and diffraction show that the 20–30 μm Al, Sc, and Zr powders fully melt and sufficiently mix in the molten state to create, on solidification, a homogeneous distribution of primary, micron-size L12 precipitates (Al 3 Sc, Al 3 Zr, and Al 3 (Sc,Zr), respectively) and nucleate micron-size Al matrix grains, as confirmed by SEM imaging of cross-sections. Here, a second laser pass, simulating a realistic additive-manufacturing build condition, fully remelts the initial volume which shows, after solidification, the same Al 3 (Sc,Zr) L1 2 primary micro-precipitates and very fine Al grains. After aging at 300–400°C, the alloys show large increases in hardness, consistent with an exceptionally high number density (1.4 × 10 24 m –3 ) and volume fraction (2.5%) of secondary Al 3 (Sc,Zr) nano-precipitates with a Sc-rich core and Zr-rich shell, as measured via atom-probe tomography.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Laser powder bed fusion of Inconel 718 on 316 stainless steel

The joining of Inconel 718 and 316 stainless steel has notable industrial importance. However, traditional arc-welding has many issues. In particular, the formation of the detrimental Nb-rich Laves phase in the heat affected zone can promote crack initiation and propagation, which leads to poor mechanical properties. The small heat input and the rapid solidification of laser powder bed fusion (PBF) has the potential to suppress this problem. Despite the infancy of the multi-material PBF, it is of interest to investigate the relationship among the process parameter, the microstructure at the materials joining point, and the mechanical performance. In this study, we used high-speed x-ray imaging to in-situ monitor the printing of Inconel 718 tracks on a 316H SS substrate via laser powder bed fusion, followed by post-build characterization and testing with SEM, XRD and nanoindentation. We observed that a higher linear energy density resulted in a more diffused interface and a smaller compositional gradient, but as a trade-off, generated more keyhole pores and thermal cracks. Furthermore, we also showed that the built layer exhibited chemical inhomogeneity at two length scales: a tens-of-mu m-scale compositional vortex due to the melt flow, and a sub-mu m-scale Nb-rich dendritic structure due to the non-equilibrium solidification that led to the formation of the Laves phase. It was inferred that the chemical inhomogeneity may benefit the mechanical properties at the interface by providing more interlocking between the two materials and also by limiting the effect of the brittle Laves phase to potentially short range.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of processing parameters and strut dimensions on the microstructures and hardness of stainless steel 316L lattice-emulating structures made by powder bed fusion

In this study, we present the effects of input processing parameters and strut thickness (in square struts) on microstructure and properties in laser powder bed fusion additively manufactured stainless steel 316L lattice-emulating structures. Lattice-emulating X-structures with square cross-sections of 1.5, 1.0, and 0.5 mm were fabricated using three different parameter sets with varying power, speed, and therefore, linear energy density. Grain size and morphology were shown to be dictated by epitaxial growth, which was dependent on weld pool morphology. Additionally, grain size and morphology were shown to change across the thickness direction of the struts (from the bottom inclined surface to the top inclined surface). The spatial variation in grain size was reflected by changes in hardness through the thickness of each strut. The 0.5 mm struts exhibited more significant grain elongation in the strut direction and larger sub-grain solidification cell diameters than their thicker counterparts. The larger sub-grain solidification cell diameters in the 0.5 mm samples resulted in correspondingly lower hardness values when compared to samples of higher thicknesses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Grain structure and texture selection regimes in metal powder bed fusion

Additive manufacturing (AM) offers opportunities to produce complex part geometries not possible with conventional processing and in some cases even improve part performance. However, adoption has been slowed by difficulties assessing microstructure variability and there is no straightforward approach to relate processing to grain structure characteristics. In this study, datasets from AdditiveFOAM heat transport simulations of laser powder bed fusion (LPBF) are used to drive ExaCA simulations of grain structure. The GPU utilization of ExaCA and an algorithmic update for modeling melt pool overlap region solidification enabled rapid and parallel simulation across a wider range of process conditions than previously explored with cellular automata-based solidification models. A texture selection angle $θ_s$ is defined based on melt pool overlap geometry, and the range of $θ_s$ over which a commonly observed texture transition occurs in characterized AM builds was well-reproduced by ExaCA simulations over a wide range of melt pool shape, hatch spacing, and layer height. ExaCA simulations with 90 degree rotation of the scan direction on every other layer reproduced a number of trends from the AM literature including grain refinement, the dominance of layers with larger melt pools on the final grain structure, and the weakening or strengthening of texture depending on odd and even layer melt pool overlap geometry. EBSD data from a benchmark AM part is used to validate the simulated mechanism of a layer rotation-induced texture strengthening effect. Importantly, these results expand the understanding of the mechanisms for texture selection in alloys with cubic crystal symmetry and offer an approach to easily evaluate processing conditions. With this new understanding, these modeling tools will enable anticipation of previously unexpected variations in grain structure and target specific microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Numerical modeling study of a man-made low-permeability barrier for the compressed air energy storage in high-permeability aquifers

Compressed air energy storage (CAES) is a grid-scale energy storage technology for intermittent energy, as proven by the decades-long successful operation of two existing compressed air energy storage in cavern (CAESC) power plants. Because of the limited availability of salt domes appropriate for CAESC, the more widely available aquifers (compressed air energy storage in aquifers, CAESA) have recently attracted considerable attention as candidates for CAES. An ideal aquifer for CAESA is highly permeable around the well to facilitate easy injection and withdrawal of air, but the high-permeability region is surrounded by low-permeability zones to minimize the loss of injected air and decrease in energy efficiency. However, such ideal geological structures are not always available in nature. Therefore, the potential of creating man-made low-permeability barrier in high-permeability aquifers is very interesting. In this paper, we investigate the feasibility of man-made low-permeability barriers in high-permeability aquifers using the numerical simulator TOUGH2/Gel to calculate the three-component flow (including a miscible gelling liquid). The simulation results show that an expected low-permeability barrier can be created by injecting grout with certain properties, and the altered aquifer performs well for CAESA. Additional sensitivity studies are also performed to reveal the effects of the various factors on the success of the low-permeability barrier creation, including the critical solidification concentration, the scale factor of the time dependence of the grout viscosity, the relative density of the grout, and the volume of the follow-up water injection. The results indicate that, in a horizontal aquifer, low critical solidification concentrations, and small scale factors are generally preferred and the density of grout should be close to that of the in situ water. For the given volume of the injected grout, there is an optimal follow-up water injection that will create the largest storage space without damaging the barrier. Furthermore, these results may help to extend the candidate sites for CAESA and the prospect of large scale energy storage.

25 ENERGY STORAGE↗
DOE OSTI2020

The microstructure of a selective laser melting (SLM)-fabricated NiTi shape memory alloy with superior tensile property and shape memory recoverability

A selective laser melting (SLM)-fabricated NiTi with superior tensile property and shape memory recoverability was obtained by using a unique stripe rotation scanning strategy. Here, the alloy characteristics, formation mechanisms and evolution in terms of twins, dislocations and precipitations of the alloy were systematically studied. Compared to the conventional smelting-followed-by-machining, the SLM fabrication process involves rapid solidification and repeated heating, which confer distinctive characteristics to the microstructures of SLM-fabricated NiTi alloys. Rapid solidification promotes the formation of a supersaturation solid solution matrix containing a high concentration vacancies, which in turn aggregate to generate a high density of dislocations. During subsequent repeated heating stages, these dislocations occur thermal motion along three directions of <001 > , <111> and <110 >, leading to the formation of thermal kinks, helical dislocations and wave morphology. Simultaneously, precipitated particles Ti 3 Ni 4 repeatedly nucleate and heterogeneous grow with the movement of dislocation. Such precipitation behavior, termed repeated precipitation, has not been previously reported in the conventional NiTi alloys, suggesting that it could be a unique characteristic of such alloy. After martensitic transformation, only two twins, {1 1¯ 1} type I twin and compound twin, are detected. The twinning lamellae of these two twins, where precipitation and dislocation pile-ups exist, often have uneven thickness and chaotic arrangement. Besides, the unique self-accommodated microstructures, such as secondary {1 1¯ 1} type I twin and compound twin with “herring-bone’’ lamellae, which often appear in the deformed or nanocrystalline NiTi, can also be observed. These unique microstructures may confer the distinctive properties to the SLM-fabricated NiTi.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A sensitivity analysis of the columnar-to-equiaxed transition for Ni-based superalloys in electron beam additive manufacturing

Solidification conditions experienced during Selective-Electron Beam Melting (S-EBM) of Nickel-Base superalloys generally lead to epitaxial grain growth from previous layers, resulting in highly textured columnar grain structures, and consequently, highly anisotropic mechanical properties. However, in recent work, researchers have been able to produce equiaxed grain structures by changing the scan strategy to influence the solidification conditions, thereby obtaining isotropic properties. The competition between columnar and equiaxed grain formation is governed by both the local thermal conditions and alloy thermodynamics. The above examples demonstrate successful manipulation of the thermal conditions through changes in process parameters. However, there is no guarantee that process conditions that encourage the formation of equiaxed grains for a given alloy will coincide with the production of defect free material. The goal of this work is therefore to understand the influence of alloy composition on the columnar-to-equiaxed transition (CET) so that alloys may be designed such that grain structure control is accessible to a wider range of S-EBM process conditions. A CET model originally developed by Hunt and extended by Gäumann was used to relate processing conditions to resulting microstructure. A sensitivity analysis of the CET model was performed using the Element Effects Methods. Furthermore, the results lead to a comparison of the relative effects of various alloying additions on the CET, as well as an assessment of other model parameters, including, most significantly, the nucleation density of equiaxed grains.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Isotopic evidence for a young lunar magma ocean

Herein, mare basalt sources and ferroan anorthosite suite cumulates define a linear array on a 146 Sm/ 144 Nd versus 142 Nd/ 144 Nd isochron plot demonstrating these materials were derived from a common reservoir at 4336 +31 / –32 Ma. The minimum proportion of the Moon that was in isotopic equilibrium at this time is estimated to be 1-3% of its entire volume based on the geographic extent from which the analyzed samples were collected and the calculated depths from which the samples were derived. Scenarios in which large portions of the Moon were molten to depths of many hundreds of kilometers are required to produce the observed Sm-Nd isotopic equilibrium between the mantle and crustal rocks at 4.34 Ga. This is a consequence of the fact that limited heating of a solid Moon above the blocking temperature of the Sm-Nd isotopic system is insufficient to diffusively homogenize radiogenic Nd throughout the mantle and crust. There are three scenarios that might account for global-scale isotopic equilibrium on the Moon relatively late in Solar System history including: (1) Sm-Nd re-equilibration of a solid Moon resulting from widespread melting in response to mantle overturn or a very large impact, (2) early accretion of the Moon followed by delayed cooling due to the presence of an additional heat source that kept a large portion of the Moon molten until 4.34 Ga, or (3) late accretion of the Moon followed by rapid cooling of the magma ocean late in Solar System history. Neither density-driven overturn of the mantle, nor a large impact, are likely to homogenize the mantle and crust to the extent required by the Sm-Nd isochron. Likewise, secondary heating mechanisms, such as tidal heating or radioactive decay, are not efficient enough to keep the Moon molten to the depth of the mare basalt source regions for many tens to hundreds of millions of years. Instead, the age of equilibrium between such a compositionally diverse set of rocks, produced on a global scale, likely records the time of primordial solidification of the Moon from a magma ocean. This scenario accounts for both the petrogenetic characteristics of lunar rock suites, as well as their Sm-Nd isotopic systematics. It is supported by the preponderance of ~4.35 Ga ages obtained for other hypothetical magma ocean crystallization products, such as ferroan anorthosite suite rocks and K, REE, and P enriched cumulates that are thought to represent flotation cumulates of the magma ocean and the last vestiges of magma ocean solidification, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Accretion timescale and impact history of Mars deduced from the isotopic systematics of martian meteorites

High precision Sm-Nd isotopic analyses have been completed on a suite of 11 martian basaltic meteorites in order to better constrain the age of silicate differentiation on Mars associated with the formation of their mantle sources. Our data is used to evaluate the merits and disadvantages of various mathematical approaches that have been employed in previous work on this topic. Ages determined from the Sm-Nd isotopic systematics of individual samples are strongly dependent on the assumed Nd isotopic composition of the bulk planet. This assumption is problematic given differences observed between the Nd isotopic composition of Earth and chondritic meteorites and the fact that these materials are both commonly used to represent bulk planetary Nd isotopic compositions. Ages determined from the slope of Sm-146-Nd-142 whole rock isochrons are not dependent on the assumed Nd-142/Nd-144 ratio of the planet, but require the sample suite to be derived from complementary, contemporaneously-formed reservoirs. In this work, we present a mathematical expression that defines the age of formation of the source regions of such a suite of samples that is based solely on the slope of a Nd-143-Nd-142 whole rock isochron and is also independent of any a priori assumptions regarding the bulk isotopic composition of the planet. This expression is also applicable to mineral isochrons and has been used to successfully calculate Nd-143-Nd-142 model crystallization ages of early refractory solids as well as lunar samples. This permits ages to be obtained using only Nd isotopic measurements without the need for Sm-147/Nd-144 isotope dilution determinations. When used in conjunction with high-precision Nd isotopic measurements completed on martian meteorites this expression yields an age of formation of the martian basaltic meteorite source regions of 4504 +/- 6 Ma. Because the Sm-Nd model ages for the formation of martian source regions are commonly interpreted to record the age at which large scale mantle reservoirs formed during planetary differentiation associated with magma ocean solidification, the age determined here implies that magma ocean solidification occurred several tens of millions of years after the beginning of the Solar System. Recent thermal models, however, suggest that Mars-sized bodies cool rapidly in less than similar to 5 Ma after accretion ceases, even in the presence of a thick atmosphere. In assuming these models are correct, an extended period of accretion is necessary to provide a mechanism to keep portions of the martian mantle partially molten until 4504 Ma. Late accretional heating of Mars could either be associated with protracted accretion occurring at a quasi-steady state or alternatively be associated with a late giant impact. If this scenario is correct, then accretion of Mars-sized bodies takes up to 60 Ma and is likely to be contemporaneous with the core formation and possibly the onset of silicate differentiation. This further challenges the concept that isotopic equilibrium is attained during primordial evolution of planets, and may help to account for geochemical evidence implying addition of material into planetary interiors after core formation was completed.

58 GEOSCIENCES↗