Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer Stability”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design and Stability Analysis of Control System in Multiport Autonomous Reconfigurable Solar Power Plants (MARS)

The Multiport Autonomous Reconfigurable Solar Power Plant (MARS) is an integrated photovoltaic (PV) power generation and energy storage system (ESS), that is designed to connect to both alternating current (AC) transmission grids and high-voltage direct current (HVDC) links. It is a three-phase plant consisting of numerous components with a complex hardware and hierarchical control architecture. This paper presents an approach to decouple the multivariable system of MARS using a recursive reduced-order and boundary layer system methodology. This approach enables efficient computation of the control parameters for the Ll, L2, and L3 controllers. To validate the effectiveness of the proposed control strategy, cyclic tests in accordance with pre-defined performance criteria using controller Hardware-in-the-Loop (cHIL) experiments are conducted. The results demonstrate that the MARS system operates consistently under steady-state conditions. Furthermore, the dynamic response of the MARS system to various grid events is analyzed, underlining the resilience of MARS in presence of faults or loss of generation within the connected WECC system.

Xia, Qian↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2021

Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D

Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. Cryogenically frozen deuterium pellets are injected from the low-field side of the DIII-D tokamak at frequencies lower than the natural ELM frequency, and heat flux is measured by infrared cameras. Ideal MHD pedestal stability calculations show that without pellet injection, these low collisionality pedestals were limited by their current density (peeling-limited) rather than their pressure gradient (ballooning-limited). ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. Cameras observing both heat flux and D-alpha emission often find significant toroidally asymmetric striations in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets. Toroidal asymmetries at the inner strike point are similar between natural and pellet-triggered ELMs, indicating that the reduction in peak heat flux and total fluence at that location is robust for the conditions reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗
DOE OSTI2026

Understanding the Performance Gap between Polycrystalline and Single-Crystal Nickel-Rich Layered Oxide Cathodes

Singe-crystal (SC) nickel-rich layered oxide cathodes, composed of boundary-free particles with high tap density, offer significant advantages in volumetric energy density and mechanical strength compared with polycrystalline (PC) cathode materials. However, as the nickel content increases (≥80%), SC Ni-rich cathodes often suffer from faster performance degradation than PC cathodes of the same composition, and the underlying causes of this discrepancy remain poorly understood. Herein, we reveal the distinct Ni redox behaviors that govern the electrochemical performance of SC and PC Ni-rich cathodes using multiscale and operando characterization techniques. Our results indicate that the increasingly heterogeneous Ni oxidation process in SC cathodes leads to the additional irreversible oxygen redox activity that deteriorates both the mechanical and chemical structures. In contrast, PC cathodes, despite with more pronounced surface reconstruction, exhibit greater chemomechanical stability due to homogeneous redox reactions during charging. Consequently, we find that bulk degradation, more than surface reactions, ultimately leads to fast capacity decay of SC Ni-rich cathodes during cycling. In conclusion, this work offers a comprehensive view on the impact of Ni redox evolutions on the chemomechanical stability in Ni-rich layered oxide cathodes, providing new insights into the longstanding performance gap between SC and PC cathodes, and guiding the rational design of Ni-rich cathode architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mass transport in a highly immiscible alloy on extended shear deformation

Forced mixing to a single phase or supersaturated solid solution (SSS) and its prerequisite microstructure evolution in immiscible systems has been a focus of research for both fundamental science and a variety of applications. Controlling the formation of SSS by shear deformation assisted processing could enable a material design beyond conventional equilibrium microstructure in immiscible systems. Here, a highly immiscible (mixing enthalpy of ~ 20 kJ·mol -1 ) Cu-50 at. % Cr binary alloy was employed to investigate the microstructure evolution and localized tendencies of SSS during severe shear deformation. The present results demonstrate distinctive microstructure refinement process in each phase and how they lead to localized SSS as a function of shear strain. Preferential dynamic recrystallization occurs in the softer Cu phase due to strain localization, leading to substantial grain refinement. The refinement of Cr phase in the top-most layer, however, is enabled by the progressive evolution of grain lamination, splitting, and fragmentation as a function of shear strain. The eventual SSS is found to be strongly dependent on the local environments that affect the dislocation activity, including the level of microstructure refinement, the interfacial orientation relationship, hardness difference, and supposed stability of oxidation. Ab initio simulations comfirm that it is more favorable to oxidize Cr than Cu at incoherent Cu/Cr interfaces which then limits the mass transport on an incoherent boundary. Our results shed light on the underpinning mechanism for non-equilibrium mass transport in immiscible systems upon severe deformation that can be applicable to a variety of processing techniques aimed at producing immiscible alloys with superior mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultranarrow TaS 2 Nanoribbons

Imposing additional confinement in two-dimensional (2D) materials can yield further control over the associated electronic, optical, and topological properties. However, synthesis of ultra-narrow nanoribbons (NRs) remains a challenge, particularly for the transition metal dichalcogenides (TMDs), and synthesizing TMD NRs narrower than 50 nm has remained elusive. In this paper, we report the vapor-phase synthesis of ultra-narrow TaS 2 NRs. The NRs are grown within the hollow cavity of carbon nanotubes, thereby limiting their lateral dimensions and layer number, while simultaneously stabilizing them against the environment. The NRs reach the monolayer (ML) limit and exhibit widths as low as 2.5 nm. Atomic-resolution scanning transmission electron microscopy (STEM) reveals the detailed atomic structure of the ultra-narrow NRs and we observe a hitherto unseen atomic structure supermodulation phenomenon of ordered defect arrays within the NRs. First-principles calculations based on density functional theory (DFT) show the presence of flat bands, as well as edge- and boundary-localized states, and help identify the atomic configuration of the supermodulation. Nanotube-templated synthesis represents a unique, transferable, and broadly deployable route toward ultra-narrow TMD NR growth.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

X-ray Diffraction Studies of Single-Crystal Materials for Broad Battery Applications

Single-crystal materials have attracted growing interest in battery research due to their well-defined crystallographic orientation, absence of grain boundaries, and enhanced mechanical and electrochemical stability. This Review provides a comprehensive overview of recent advances in the synthesis, structural evolution, and performance optimization of single-crystal electrodes and solid electrolytes. Particular focus is placed on the application of advanced X-ray diffraction (XRD) techniques, including operando synchrotron diffraction, reciprocal space mapping, and Bragg coherent diffraction imaging, which have enabled in-depth investigations of lattice strain, cation disorder, phase transitions, and defect formation. Representative case studies across Ni-rich layered oxides, spinel-type cathodes, and garnet-based electrolytes are examined to highlight the structural features unique to single crystals. Additionally, the synergistic integration of XRD with machine learning, tomography, and spectroscopy is discussed as a powerful direction for real-time analysis and predictive modeling. Furthermore, these insights provide critical guidance for the rational design of high-performance single-crystal materials in lithium, sodium, and solid-state battery systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Insights into the Growth Orientation and Phase Stability of Chemical-Vapor-Deposited Two-Dimensional Hybrid Halide Perovskite Films

Chemical vapor deposition (CVD) offers a large-area, scalable, and conformal growth of perovskite thin films without the use of solvents. Low-dimensional organic–inorganic halide perovskites, with alternating layers of organic spacer groups and inorganic perovskite layers, are promising for enhancing the stability of optoelectronic devices. Moreover, their multiple quantum-well structures provide a powerful platform for tuning excitonic physics. Here, in this work, we show that the CVD process is conducive to the growth of 2D hybrid halide perovskite films. Using butylammonium (BA) and phenylethylammonium (PEA) cations, the growth parameters of BA 2 PbI 4 and PEA 2 PbI 4 and mixed halide perovskite films were first optimized. These films are characterized by well-defined grain boundaries and display characteristic absorption and emission features of the 2D quantum wells. X-ray diffraction (XRD) and a noninteger dimensionality model of the absorption spectrum provide insights into the orientation of the crystalline planes. Unlike BA 2 PbI 4 , temperature-dependent photoluminescence measurements from PEA 2 PbI 4 show a single excitonic peak throughout the temperature range from 20 to 350 K, highlighting the lack of defect states. These results further corroborate the temperature-dependent synchrotron-based XRD results. Furthermore, the nonlinear optical properties of the CVD-grown perovskite films are investigated, and a high third harmonic generation efficiency is observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Scrape-off layer stability, turbulence and transport in MAST-U (Final Report)

U.S. Department of Energy (DOE) supported work at Lodestar Research Corporation under grant DE-SC0019270 is summarized in this final report. The research conducted under this grant explored the properties of the edge plasma in spherical torus magnetic fusion devices. Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The edge and scrape-off layer (SOL) region of fusion plasmas is the region where the hot core plasma transitions to a cooler surrounding boundary plasma. The properties of this boundary plasma determine how interactions with material surfaces will impact both the overall device performance and the lifetime of crucial components such as the walls of the device chamber and other hardware. The research carried out studied filamentary plasma instabilities and turbulent transport in the edge and SOL of the MAST U spherical torus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

New insight into tuning magnetic phases of RMn6Sn6 kagome metals

Abstract Predicting magnetic ordering in kagome compounds offers the possibility of harnessing topological or flat-band physical properties through tuning of the magnetism. Here, we examine the magnetic interactions and phases of ErMn 6 Sn 6 which belongs to a family of R Mn 6 Sn 6 , R = Sc, Y, Gd–Lu, compounds with magnetic kagome Mn layers, triangular R layers, and signatures of topological properties. Using results from single-crystal neutron diffraction and mean-field analysis, we find that ErMn 6 Sn 6 sits close to the critical boundary separating the spiral-magnetic and ferrimagnetic ordered states typical for non-magnetic versus magnetic R layers, respectively. Finding interlayer magnetic interactions and easy-plane Mn magnetic anisotropy consistent with other members of the family, we predict the existence of a number of temperature and field dependent collinear, noncollinear, and noncoplanar magnetic phases. We show that thermal fluctuations of the Er magnetic moment, which act to weaken the Mn-Er interlayer magnetic interaction and quench the Er magnetic anisotropy, dictate magnetic phase stability. Our results provide a starting point and outline a multitude of possibilities for studying the behavior of Dirac fermions in R Mn 6 Sn 6 compounds with control of the Mn spin orientation and real-space spin chirality.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-Entropy Ceramic Coatings: Transformative New Materials for Environmentally-Compatible Thin-Film Insulators against High-T Molten Salts

An emerging new class of promising materials, High-Entropy Ceramics (HECs) or Compositionally-Complex Ceramics (CCCs), have been successfully fabricated in the last three years. This project evaluates HECs/CCCs in the concentrated solar power (CSP) context. In this project, we have identified a couple of simple, and key descriptors governing thermal conductivity in HECs. Size disorder among the cations was found to have a strong influence on the thermal conductivity across multiple systems and throughout a vast compositional space. Since identifying this, we have optimized previous compositions to achieve low thermal conductivities of $k$~$1-1.4Wm{-1}K{-1}$. We have further reduced our selection of promising compositions based on project guidelines and applied a theory that can potentially predict reactivity with molten salts to identify a small set of HECs that are suitable for further molten salt heat transfer media (HTM) stability testing. Many compositions reveal poor resistance to molten chloride salt after being exposed for 24 h compared to the reference structural alloy, Haynes 230. However, after 100 h, a few HECs/CCCs show strong parabolic rate constant effects and provide greater resistance to chloride salt compared to Haynes 230 alloy, which shows more linear kinetics. Postmortem characterization of the microstructure reveals the parabolic rate constant originates from a dense reaction layer form at the interface. This feature is absent in Haynes 230, which only shows grain boundary infiltration. The HECs/CCCs were also put in contact with the metal alloy under Ar for 800 h at 800°C and show negligible interaction. Lastly, parabolic rate constant HECs/CCCs also have low thermal conductivity $(<2 Wm^{-1} K^{-1})$ at 800°C. These specimens satisfy all project requirements and demonstrate their promising application in next-generation CSP.

14 SOLAR ENERGY↗
DOE OSTI2023

Stabilizing Ni-rich layered cathode for high-voltage operation through hierarchically heterogeneous doping with concentration gradient

High-nickel LiNi x Mn y Co 1-x-y O 2 (NMC) cathodes have demonstrated superior energy density, yet their stability is compromised under high voltage conditions. To address this, here we propose a strategy of heterogeneous doping with a concentration gradient, specifically through Sr–Zr co-modification. We synthesized Ni-rich NMC particles featuring several micron-sized secondary particles composed of micron-sized primary grains. This design aims to harness the structural robustness of single-crystalline grains and the favorable diffusion kinetics of polycrystalline secondary particles. Systematic characterization using a combination of electrochemical measurements and synchrotron analytics reveals an intriguing pattern of hierarchically heterogeneous Sr–Zr co-doping. It demonstrates a depth-dependent concentration gradient at the secondary particle level and competing dopant segregation over the buried grain boundaries. This unique characteristic creates opportunities for enhancing battery performance, particularly by optimizing precursors and implementing advanced modulation techniques. We also investigate the dissolution and precipitation of the cathode's transition metal cations upon high-voltage cycling. These insights suggest that a tailored compositional variation can be a viable approach to effectively design the next-generation high-Ni NMC cathode materials for high-voltage lithium batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interface velocity dependent solute trapping and phase selection during rapid solidification of laser melted hypo-eutectic Al-11at.%Cu alloy

Microstructure formation of a hypo-eutectic Al-11Cu (atom percent) alloy during solidification after laser melting has been studied by transmission electron microscopy (TEM). The evolution of the solid-liquid-interface velocity, V SL , during the solidification process has been determined from direct observation by in-situ TEM. This enabled correlating V SL with development of four distinct microstructure zones. Crystal growth mode transitions from planar to cellular, cellular to dendritic, dendritic to cellular, and cellular to planar, have been observed for the accelerating solid-liquid-interface. The transition from coupled two-phase growth to the single-phase growth occurred for V SL = V a = (0.80±0.05) m/s, where V a is the velocity of absolute stability, at the onset of banded morphology grain formation. The in-situ and post-mortem TEM uniquely permitted determination of the non-equilibrium solidus for the rapidly solidifying Al-11Cu alloy. Experimental evidence for solute clustering and chemical ordering tendencies at length scales on the order of ≤ 5nm has been detected in the single-phase regions of the banded grains. The structural features of the single-phase bands have been interpreted as signatures of ‘frozen in’ configurations present in the liquid boundary layer adjacent to the growing crystal, which has a width of about 3nm. Finally, the nano-scale spatiotemporal resolution experimental TEM studies performed here provided quantitative metrics, e.g. the solidification interface velocity dependent solute concentration of the α-Al phase and the near-atomic scale structure in the single-phase bands, that are uniquely suitable for comparison with theory and model predictions for solidification microstructure development in multicomponent alloys after laser melting.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure : Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.

diopside↗