Results for “Composite models” — Page 9

DOE OSTI2024

Chemical composition based machine learning model to predict defect formation in additive manufacturing

With a goal of exploiting additive manufacturing to improve the manufacturing of existing reactor materials, we developed a chemical composition-based machine learning model to predict the printability of any given alloy in laser powder bed fusion (L-PBF) using experimental data from peer-reviewed literature. We defined printability as the ability to avoid defects like cracking, balling, porosity, and lack of fusion, that are caused by thermal stresses (during solidification or liquation), molten pool disintegration into disconnected small beads or lack of heat input respectively. Our models predict the tendency of balling defect formation and porosity percentage for a given composition, under a given set of processing conditions. To predict the likelihood of balling defect, three models: a random forest classifier, a gradient boost regressor and a neural network were trained on a dataset containing both traditional alloys and high entropy alloys. The neural network model showed the highest accuracy of 92.3 % in predicting the balling defect formation. A random forest regressor, gradient boost regressor and neural network were trained and tested on a dataset of various alloys to predict porosity. The random forest regressor showed the best predictions with an R 2 score of 0.97. The models also revealed the relative importance of the input descriptors on defect-formation tendency. Of particular significance was the identification of carbon as an important element in determining the occurrence of balling and percent porosity in alloys like steel, as well as being moderately important to the percentage porosity in other alloys as well as steel. Manganese was also identified as a key descriptor for the percentage of porosity in steel and other alloys. Manganese’s low thermal conductivity and consistent presence in the dataset is the likely cause for its contribution. Carbon’s role is attributable to its relatively high specific heat and high melting temperature. In conclusion, our model serves as a swift, chemistry-based tool to design experiments and find modified compositions better suited for additive manufacturing.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Neutrino Dipole Moments and Radiative Signatures from Partial Compositeness

We investigate composite neutrino models where heavy neutrinos emerge as bound states from a near-conformal strongly coupled sector. Standard Model neutrinos mix with these composite singlets via an inverse seesaw mechanism, where the anomalous scaling dimensions of the composite-sector operators naturally suppress light neutrino masses to sub-eV scales. Matching the conformal dynamics onto low-energy theory yields enhanced electromagnetic transition dipole operators with couplings $d_{μN} \sim 10^{-6}$-$10^{-8}\,\mathrm{GeV}^{-1}$, parametrically larger than the loop-level predictions of minimal Dirac or Majorana models. We carry out a dedicated event-level simulation of the production-and-decay chain $νX \to \mathcal{U} X \to νγX$ and compute the resulting event rates at MiniBooNE and MINERvA within the model, accounting for the composite production cross section and decay kinematics in detail. We further present predictions for the photon energy, angular, and multiplicity distributions. For the benchmark scenarios accessible at these experiments the radiative signal is predominantly single-photon; the composite structure additionally permits fragmentation of the up-scattered state into multiple heavy neutrinos, each decaying as $N\toνγ$, with multi-photon final states emerging for lighter compositeness scales or higher beam energies as a qualitatively new probe of the composite dynamics.

Assi, Benoît [Cincinnati U.; Fermilab] (ORCID:0000↗

DOE OSTI2023

Glass formulation and composition optimization with property models: A review

Abstract Glass is a versatile material with a remarkable history and many practical applications. It plays a critical role in our everyday lives, the advancement of science, and the development of many technologies. The Edisonian type trial‐and‐error method was commonly used for conventional design of glass compositions, which was time‐consuming and costly. With the urgent need to develop new glass compositions for technology applications rapidly, it has become necessary to develop precise property models with predictive powers using large databases and efficient formulation approaches. This paper reviews the design of glass compositions using these analytical and numerical models of composition–structure–property relations of glasses, some based on large databases and machine learning approaches. Aspects of data collection, model fitting, feature extraction, model evaluation, and uncertainty quantification will be covered. Furthermore, advances in the glass optimization framework and available tools are summarized with examples. The outlook and perspective for further glass property model development and formulation approaches are discussed.

Lu, Xiaonan↗

DOE OSTI2020

Vadose Zone Model for REDOX Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Vadose Zone Model for PFP Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Vadose Zone Model for PUREX Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Vadose Zone Model for B Complex for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3. Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B Complex vadose zone model. The B Complex model encompasses the region around single-shell tank farms 241-B, 241-BX, and 241-BY, and adjacent cribs, trenches, and burial grounds in the northwestern part of the 200 East Area. Sources of waste to this area were primarily from B Plant (221-B) and U Plant (221-U).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Weakly supervised anomaly detection with event-level variables

We introduce a new topology for weakly supervised anomaly detection searches, diobject plus X. In this topology, one looks for a resonance decaying to two standard model particles produced in association with other anomalous event activity (X). This additional activity is used for classification. We demonstrate how anomaly detection techniques which have been developed for dijet searches focusing on jet substructure anomalies can be applied to event-level anomaly detection in this topology. To robustly capture event-level features of multiparticle kinematics, we employ new physically motivated variables derived from the geometric structure of a collision’s phase space manifold. As a proof of concept, we explore the application of this approach to several benchmark signals in the di-𝜏 and di-𝜇 plus X final states. We demonstrate that our anomaly detection approach can reach discovery-level significances for signals that would be missed in a conventional bump-hunt approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Model Validation for the FY2021 SRS Composite Analysis Monitoring Plan

Using a projected end-state date of 2065 (SRNS 2015b), the Savannah River Site (SRS) Composite Analysis (CA) modeling for each facility and waste site began on the inventory year assigned to it so that source depletion and radionuclide transport out of the system could be appropriately captured. Some SRS waste sites that have already achieved their end states (i.e., end-state inventories and end-state configuration) are currently contributing to the potential off-site public dose through source release, groundwater transport, discharge to on-site surface streams, and stream transport to the CA point of assessments (POAs). The inventory year assigned to these waste sites is 2002 or before. This means that SRS CA results from 2002 and beyond are a reasonable representation for these waste sites that have already achieved their end states and are currently contributing to the potential off-site public dose. The SRS Annual Environmental Report (AER) monitoring can differentiate and separate liquid pathway data allowing the data representing only waste sites at their end state to be produced. Because the SRS CA has projected reasonable end-state impacts from 2002 and beyond, and the AER monitoring can differentiate and separate operating and end-state contributions to annual liquid pathway release, an opportunity exists to use the AER monitoring data to validate the SRS CA model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Multiscale modeling-enabled design of multifunctional composites

This study aims to create a comprehensive model that considers multiple scales and physics for predicting the electromechanical behavior of fiber-reinforced composites enhanced with barium titanate (BaTiO3). In our earlier work, we have demonstrated that depositing BaTiO3 microparticles of 200-nm-diameter, on fiber surfaces during fiber-reinforced composite fabrication enhances mechanical strength, passive self-sensing, and energy harvesting properties. The key is to carefully control the microparticle concentration to prevent agglomeration. Since the particles are micron-sized, understanding how agglomeration affects the composites' electromechanical properties is crucial for guiding such multifunctional materials’ design. This study introduces a micromechanics-based approach to explore the impact of microparticle dispersion on the bulk composites' electromechanical properties. Insights gained from this investigation are applied in experiments, enabling accurate predictions of mechanical and self-sensing responses in BaTiO3-enhanced fiber-reinforced composites. Micro-level findings from this computational approach can be integrated into larger continuum models to comprehensively capture the electromechanical behavior of the composite structures at bulk scale. The proposed model is validated by comparing predictions with experimental results, accounting for the nonlinear mechanical and electromechanical behaviors of constituent materials. Consequently, this computational model serves as a digital platform for efficiently designing multifunctional composites.

Gupta, Sumit↗

DOE OSTI2021

Composite fermion nonlinear sigma models

We study particle-hole symmetry at the integer quantum Hall plateau transition using composite fermion mean-field theory. Because this theory implicitly includes some electron-electron interactions, it also has applications to certain fractional quantum Hall plateau transitions. Previous work [P. Kumar et al., Phys. Rev. B 100, 235124 (2019)] using this approach showed that the diffusive quantum criticality of this transition is described by a nonlinear sigma model with topological θ = π term. This result, which holds for both the Dirac and Halperin, Lee, and Read composite fermion theories, signifies an emergent particle-hole (reflection) symmetry of the integer (fractional) quantum Hall transition. Here we consider the stability of this result to various particle-hole symmetry-violating perturbations. In the presence of quenched disorder that preserves particle-hole symmetry, we find that finite longitudinal conductivity at this transition requires the vanishing of a symmetry-violating composite fermion effective mass, which if present would generally lead to θ ≠ π and a corresponding violation of particle-hole symmetric electrical transport σ xy ≠ $\frac1{2}$$\frac{e^2}{h}$ . When the disorder does not preserve particle-hole symmetry, we find that θ can vary continuously within the diffusive regime. Our results call for further study of the universality of the quantum Hall plateau transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

A composite Higgs with a heavy composite axion

We consider the strong dynamics associated with a composite Higgs model that simultaneously produces dynamical axions and solves the strong CP problem. The strong dynamics arises from a new Sp or SU(4) hypercolor gauge group containing QCD colored hyperfermions that confines at a high scale. The hypercolor global symmetry is weakly gauged by the Standard Model electroweak gauge group and an enlarged color group, SU(N + 3) × SU(N)'. When hyperfermion condensates form, they not only lead to an SU(5)/SO(5) composite Higgs model but also spontaneously break the enlarged color group to SU(3) c × SU(N) D . At lower energies, the SU(N) D group confines, producing two dynamical axions that eliminates all CP violation. Furthermore, small instantons from the SU(N)' group can enhance the axion mass, giving rise to TeV scale axion masses that can be detected at collider experiments. Our model provides a way to unify the composite Higgs with dynamical axions, without introducing new elementary scalar fields, while also extending the range of axion masses that addresses the strong CP problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

An ICME Modeling Framework for Titanium/Tungsten-Carbide Metal Matrix Composites

This report describes a collaborative project to develop a validated, predictive model for the high temperature mechanical properties of a titanium-matrix, tungsten-carbide/cobalt-reinforced metal matrix composite. The modeling approach was to first develop a detailed, microstructural model linking the material structure and the interfacial debonding properties to the effective properties of the material. The project then completed a throughput simulation campaign to generate a large number of simulations for discrete microstructures and different debonding parameters. Finally, the project trained a fast, Gaussian process surrogate model against this simulation database to provide a quick model linking the material compositions, structure, and processing parameters to the resulting material properties. This model was validated against high temperature tensile test data on a few particular composite compositions. The tests validate the model predictions for ultimate tensile strength and uniform elongation/ductility, meaning the final surrogate model can now be used to tune the material composition and processing parameters to identify optimal composite compositions for particular applications.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Analytical Modeling of Metal Foam Composite Phase Change Materials (PCM) in Thermal Energy Storage Using Asymptotic Analysis

The use of phase change materials (PCMs) for thermal energy storage can release or absorb a significant amount of latent heat during the freezing or melting process, offering a higher energy storage density. One of the main drawbacks of PCMs is their low thermal conductivity, resulting in poor thermal performance. Recent research has attempted to enhance heat transfer and increase the thermal conductivity of PCMs, including the use of metal foams. However, modeling the metal foam composite PCM using conventional methods is computationally expensive. This paper proposes an asymptotic solution for a Stefan-like problem subject to a convective boundary for outward solidification in a hollow cylinder, capable of predicting the freeze-melt cycle of the metal foam composite PCM. Specifically, three temporal regimes and four spatial layers are considered in the asymptotic analysis for each phase change process. The thermal conductivity is calculated by a theoretical three-dimensional tetrakaidecahedron model, while other thermophysical properties are obtained using the method of volume averaging. The results are verified with numerical data and validated against experimental data in the literature. The presented analytical modeling framework could have the potential to be applied to other types of composite PCMs with considerably lower computational costs compared with conventional methods.

analytical model↗

DOE OSTI2021

Methods and devices for preventing computationally explosive calculations in a computer for model parameters distributed on a hierarchy of geometric simplices

A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.

Source record↗

DOE OSTI2025

Uncertainty-Based Design: Finite Element and Explainable Machine Learning Modeling of Carbon–Carbon Composites for Ultra-High Temperature Solar Receivers

Design under uncertainty has significantly grown in research developments during the past decade. Additionally, machine learning (ML) and explainable ML (XML) have offered various opportunities to provide reliable predictable models. The current article investigates the use of finite element modeling (FEM), ML and XML predictions, and uncertain-based design of carbon-carbon (C-C) composites for use in ultra-high temperatures. A C-C composite concentrating solar power (CSP) as a microvascular receiver is considered as a case study. These C-C composites are fiber composites with directly integrated carbonized microchannels to form a lightweight, high-absorptivity material that includes an embedded microvascular network of channels. The topology of these microchannels is engineered to optimize heat transfer to a supercritical carbon dioxide (sCO2) heat transfer fluid. The mechanical characterization of C-C composites is highly challenging. Thus, designing every component made of C-C composites for ultra-high temperature applications needs an uncertainty-based analysis. As a part of a comprehensive project on the development of a novel carbonized microvascular C-C composite, this paper explores C-C composite sensitivity analysis, FEM, ML prediction, and XML analysis. The resulting composite can then be carbonized and coated with an oxidation-resistant coating to form a thermally efficient and mechanically robust C-C composite. An ANSYS 3-D-FE model was used to analyze the CSP’s stress/strain. To consider the variability in the mechanical and thermal properties of C-C composites, various mechanical properties are considered as the ANSYS FEM’s input. A synthetic dataset from 730 ANSYS runs was produced to feed into the ML and XML algorithms for uncertainty analysis and prediction. The ML and XML algorithms could accurately predict the CSP stresses/strains.

Daghigh, Vahid (ORCID:0000000298941620)↗

DOE OSTI2020

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗

DOE OSTI2022

Mixing of Condensable Constituents with H–He during the Formation and Evolution of Jupiter

Simulations of Jupiter's formation are presented that incorporate mixing of H–He with denser material entering the planet as solids. Heavy compounds and gas mix substantially when the planet becomes roughly as massive as Earth, because incoming planetesimals can fully vaporize. Supersaturation of vaporized silicates causes the excess to sink as droplets, but water remains at higher altitudes. Because the mean molecular weight decreases rapidly outward, some of the compositional inhomogeneities produced during formation can survive for billions of years. After 4.57 Gyr, our Jupiter model retains compositional gradients; proceeding outward, one finds (i) an inner heavy-element core, the outer part derived from hot supersaturated rain-out; (ii) a composition-gradient region, containing most of the heavy elements, where H–He abundance increases outward, reaching about 0.9 mass fraction at 0.3 of the radius, with silicates enhanced relative to water in the lower parts and depleted in the upper parts; (iii) a uniform-composition region (neglecting He immiscibility) that is enriched over protosolar and contains most of the planet’s mass; and (iv) an outer region where cloud formation (condensation) of heavy constituents occurs. This radial compositional profile has heavy elements more broadly distributed than predicted by classical formation models but less diluted than suggested by Juno-constrained gravity models. The compositional gradients in the region containing the bulk of the heavy elements prevent convection, in both our models and those fitting current gravity, resulting in a hot interior where much of the accretion energy remains trapped.

79 ASTRONOMY AND ASTROPHYSICS↗