Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 757 records · Page 42

Measurements of the Coexistence Curve near the Liquid-Gas Critical Point

The shape of the liquid-gas coexistence curve of He-3 very near the critical point (-2x10(exp -6) < t < -5x10(exp -3) was measured using the quasi-static thermogram method. The study was performed in Earth s gravitational field using two different height calorimetry cells, both originally designed for simultaneous measurements of the isochoric heat capacity, isothermal compressibility, and PVT. The heights of two cells were 0.5 mm and 4.8 cm. The uncertainty in measuring the phase transition temperature was typically +/-2 micro-K. The measured coexistence curve near the critical point was strongly affected by the gravitational field. Away from the critical point, the coexistence curve obtained using this technique was also consistent with the earlier work using the local density measurements of Pittman et al. The recent crossover parametric model of the equation-of-state are used to analyze the height-dependent measured coexistence curves. Data analyses have indicated that microgravity will permit measurements within two additional decades in reduced temperatures beyond the best gravity-free data obtained in Earth-bound experiments.

Hahn, Inseob↗

Development and Characterization of Embedded Sensory Particles Using Multi-Scale 3D Digital Image Correlation

A method for detecting fatigue cracks has been explored at NASA Langley Research Center. Microscopic NiTi shape memory alloy (sensory) particles were embedded in a 7050 aluminum alloy matrix to detect the presence of fatigue cracks. Cracks exhibit an elevated stress field near their tip inducing a martensitic phase transformation in nearby sensory particles. Detectable levels of acoustic energy are emitted upon particle phase transformation such that the existence and location of fatigue cracks can be detected. To test this concept, a fatigue crack was grown in a mode-I single-edge notch fatigue crack growth specimen containing sensory particles. As the crack approached the sensory particles, measurements of particle strain, matrix-particle debonding, and phase transformation behavior of the sensory particles were performed. Full-field deformation measurements were performed using a novel multi-scale optical 3D digital image correlation (DIC) system. This information will be used in a finite element-based study to determine optimal sensory material behavior and density.

Cornell, Stephen R.↗

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

PAHs and the Diffuse Interstellar Bands. What have we Learned from the New Generation of Laboratory and Observational Studies?

Polycyclic Aromatic Hydrocarbons (PAHs) are an important and ubiquitous component of carbon-bearing materials in space. PAHs are the best-known candidates to account for the IR emission bands (UIR bands) and PAH spectral features are now being used as new probes of the ISM. PAHs are also thought to be among the carriers of the diffuse interstellar absorption bands (DIBs). In the model dealing with the interstellar spectral features, PAHs are present as a mixture of radicals, ions and neutral species. PAH ionization states reflect the ionization balance of the medium while PAH size, composition, and structure reflect the energetic and chemical history of the medium. A major challenge for laboratory astrophysics is to reproduce (in a realistic way) the physical conditions that exist in the emission and/or absorption interstellar zones, An extensive laboratory program has been developed at NASA Ames to characterize the physical and chemical properties of PAHs in astrophysical environments and to describe how they influence the radiation and energy balance in space and the interstellar chemistry. In particular, laboratory experiments provide measurements of the spectral characteristics of interstellar PAH analogs from the ultraviolet and visible range to the infrared range for comparison with astronomical data. This paper will focus on the recent progress made in the laboratory to measure the direct absorption spectra of neutral and ionized PAHs in the gas phase in the near-W and visible range in astrophysically relevant environments. These measurements provide data on PAHs and nanometer-sized particles that can now be directly compared to astronomical observations. The harsh physical conditions of the IS medium - characterized by a low temperature, an absence of collisions and strong V W radiation fields - are simulated in the laboratory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature (of the order of 100 K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase. Intrinsic band profiles and band positions of cold gas-phase PAHs can now be measured with high- sensitivity spectroscopy and directly compared to the astronomical data. The electronic bands measured for ionized PAH are found to be intrinsically broad (about 20/cm) while the bands associated with the neutral precursors are narrower (of the order of 2 - 10/cm).

Salama, Farid↗

Accelerated Aging of PV Cables - The Development of Methods Towards Combined-Accelerated Stress Testing

PV cables allow the collection and distribution of electricity from modules to the energy grid. Durable cabling allows continuous operation of PV installations. Cables with a life less than that of the modules must be replaced - reducing electricity generation and adding to the system operating expense. We have examined the durability of a variety of representative cable samples using artificial UV weathering (IEC TS 62788-7-2) and combined accelerated stress testing (C-AST). Jackets in this study include PV-rated (polyolefin), building wire (PVC/PA), and multi-conductor (rubber, as in PV trackers), with some examples of black and red cables. Characterizations at read points include appearance (camera and optical microscope), surface roughness (linear profilometry), instrumented indentation, surface morphology (SEM), chemical composition (EDS), thermal decomposition (TGA), phase transition (DSC), crystalline structure (WAXS), and polymer chemical structure (FTIR). Complimentary examinations confirm damage to the external and interior layer(s) of less durable samples, including changes in morphology (e.g., cracking), chemistry (e.g., oxidation), structure, and mechanical performance (e.g., through-thickness modulus and hardness). The results may be generally compared to other examples from the field (PV installations) or after accelerated testing (IEC 61215). The methods in this study are presently being applied in a follow-on C-AST study, incorporating both dynamic and kinematic fixtures.

balance of systems↗

Quantum phase transition dynamics in the two-dimensional transverse-field Ising model

The quantum Kibble-Zurek mechanism (QKZM) predicts universal dynamical behavior near the quantum phase transitions (QPTs). It is now well understood for the one-dimensional quantum matter. Higher-dimensional systems, however, remain a challenge, complicated by the fundamentally different character of the associated QPTs and their underlying conformal field theories. In this work, we take the first steps toward theoretical exploration of the QKZM in two dimensions for interacting quantum matter. We study the dynamical crossing of the QPT in the paradigmatic Ising model by a joint effort of modern state-of-the-art numerical methods, including artificial neural networks and tensor networks. As a central result, we quantify universal QKZM behavior close to the QPT. We also note that, upon traversing further into the ferromagnetic regime, deviations from the QKZM prediction appear. We explain the observed behavior by proposing an extended QKZM taking into account spectral information as well as phase ordering. Our work provides a testing platform for higher-dimensional quantum simulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Synchrotron-based techniques for characterizing STCH water-splitting materials

Understanding the role of oxygen vacancy–induced atomic and electronic structural changes to complex metal oxides during water-splitting processes is paramount to advancing the field of solar thermochemical hydrogen production (STCH). The formulation and confirmation of a mechanism for these types of chemical reactions necessitate a multifaceted experimental approach, featuring advanced structural characterization methods. Synchrotron X-ray techniques are essential to the rapidly advancing field of STCH in part due to properties such as high brilliance, high coherence, and variable energy that provide sensitivity, resolution, and rapid data acquisition times required for the characterization of complex metal oxides during water-splitting cycles. X-ray diffraction (XRD) is commonly used for determining the structures and phase purity of new materials synthesized by solid-state techniques and monitoring the structural integrity of oxides during water-splitting processes (e.g., oxygen vacancy–induced lattice expansion). X-ray absorption spectroscopy (XAS) is an element-specific technique and is sensitive to local atomic and electronic changes encountered around metal coordination centers during redox. While in operando measurements are desirable, the experimental conditions required for such measurements (high temperatures, controlled oxygen partial pressures, and H 2 O) practically necessitate in situ measurements that do not meet all operating conditions or ex situ measurements. Here, we highlight the application of synchrotron X-ray scattering and spectroscopic techniques using both in situ and ex situ measurements, emphasizing the advantages and limitations of each method as they relate to water-splitting processes. The best practices are discussed for preparing quenched states of reduction and performing synchrotron measurements, which focus on XRD and XAS at soft (e.g., oxygen K-edge, transition metal L-edges, and lanthanide M-edges) and hard (e.g., transition metal K-edges and lanthanide L-edges) X-ray energies. The X-ray absorption spectra of these complex oxides are a convolution of multiple contributions with accurate interpretation being contingent on computational methods. The state-of-the-art methods are discussed that enable peak positions and intensities to be related to material electronic and structural properties. Through careful experimental design, these studies can elucidate complex structure–property relationships as they pertain to nonstoichiometric water splitting. A survey of modern approaches for the evaluation of water-splitting materials at synchrotron sources under various experimental conditions is provided, and available software for data analysis is discussed.

08 HYDROGEN↗

Spectral induced polarization monitoring of induced calcite precipitation in subsurface sediments

Co-precipitation of contaminants within the crystalline structure of calcite is a promising natural attenuation or remedial technology being considered at contaminated sites. We explore the sensitivity of the spectral induced polarization (SIP) method to induced calcite precipitation in natural sediments as a path forward to non-invasively monitor these sites. We performed time-lapse column experiments using phased (I-IV) injections over 40 days on natural sediments from the Hanford Site (Washington State, USA). In the phased injections, abiotic calcite precipitation was induced and confirmed to have occurred. Previous work in glass beads and homogeneous sand was limited to high frequency detection of calcite, however in this work we observed the development of two polarization mechanisms, one at high frequency (>100 Hz) and one at low frequency (< 100 Hz). Based on the characteristic frequencies from the SIP high and low frequency regimes, characteristic length scales (L) were computed where the adsorption mode of Na+ versus Ca2+ was compared by using diffusion coefficients corresponding to Na+ versus an arithmetically averaged value for Na+ and Ca2+. Using the diffusion coefficient of Na+, the high frequency L was found to correlate well with the size of the calcite crystals. The low frequency L correlated well with the individual natural sediment grain sizes within the columns. During late experimental times (day 36 and 40), the characteristic low frequency in two of the experimental columns shifted to lower frequencies (<0.001 Hz) which may signify SIP sensitivity of the formed calcite with the sediment grains. In field applications, the development of a low frequency polarization length scale to monitor calcite precipitation is promising for field monitoring applications, however further laboratory work needs to be performed to examine the SIP sensitivity of calcite formation in the presence of to natural sediments

calcium carbonate, calcite, precipitation, spectra↗

Leveraging intermediate resonances to probe CP violation at colliders

We explore the phenomenological impact of interference in tree-level contributions to three-body final states in $2\rightarrow 3$ scattering processes. This work introduces a novel search strategy leveraging asymmetries to enable sensitivity to CP-violating effects in less well-explored regions of phase space. Analytically, we demonstrate the effectiveness of this observable in probing interference between Standard Model charged-current decays and effective left-handed vector interactions, illustrated in a toy model featuring a scalar leptoquark, $S_1 \sim (3, 1, -\,1/3)$. Numerically, we apply this framework to studying the process $pp\rightarrow b \tau \nu $; unlike traditional high-$p_T$ searches or “bump hunts”, this approach utilizes an intermediate energy regime – where new physics is neither light enough to be produced on shell or heavy enough to justify an effective field theory treatment. A proof-of-principle analysis at parton level demonstrates a percent-level asymmetry, with sensitivity also to BSM weak-CP phase. While the specific phase sensitivity is diminished at particle level due to showering and detector effects, a machine learning classifier can recover sensitively to the presence of SM-BSM interference, significantly outperforming standard analysis methods. Notably discrimination between BSM signal and SM background could be achieved at the 2$\sigma $ level for the current LHC dataset and 8$\sigma $ at the High-Luminosity LHC. Moreover, this asymmetry observable as defined can also be more broadly applied to other searches for CP-violation in $2\rightarrow 3$ processes in present and future collider environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Topological Spin Textures in an Insulating van der Waals Ferromagnet

Generation and control of topological spin textures constitutes one of the most exciting challenges of modern spintronics given their potential applications in information storage technologies. Of particular interest are magnetic insulators, which due to low damping, absence of Joule heating and reduced dissipation can provide energy-efficient spin-textures platform. Here, it is demonstrated that the interplay between sample thickness, external magnetic fields, and optical excitations can generate a prolific paramount of spin textures, and their coexistence in insulating CrBr 3 van der Waals (vdW) ferromagnets. Using high-resolution magnetic force microscopy and large-scale micromagnetic simulation methods, the existence of a large region in T-B phase diagram is demonstrated where different stripe domains, skyrmion crystals, and magnetic domains exist and can be intrinsically selected or transformed to each-other via a phase-switch mechanism. Lorentz transmission electron microscopy unveils the mixed chirality of the magnetic textures that are of Bloch-type at given conditions but can be further manipulated into Néel-type or hybrid-type via thickness-engineering. The topological phase transformation between the different magnetic objects can be further inspected by standard photoluminescence optical probes resolved by circular polarization indicative of an existence of exciton-skyrmion coupling mechanism. The findings identify vdW magnetic insulators as a promising framework of materials for the manipulation and generation of highly ordered skyrmion lattices relevant for device integration at the atomic level.

36 MATERIALS SCIENCE↗

Improved accuracy in degenerate variational integrators for guiding centre and magnetic field line flow

First-order-accurate degenerate variational integration (DVI) was introduced in Ellison et al. ( Phys. Plasmas , vol. 25, 2018, 052502) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introduce second-order-accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second-order scheme with a preserved two-form by composing a first-order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second-order-accurate DVI schemes. We test these second-order schemes numerically on two systems having non-canonical variables, namely the magnetic field line and guiding centre systems. Variational integration for Hamiltonian systems with non-uniform time steps, in terms of an extended phase space Hamiltonian, is generalized to non-canonical variables. It is shown that preservation of proper degeneracy leads to single-step (one-step) methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second-order-accurate as well as first-order schemes, and can be applied to adapt the time stepping to an error estimate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Structure and Phase Behavior of a Block Copolymer Electrolyte under dc Polarization

There is an important consideration when designing lithium battery electrolytes for advanced applications is how the electrolyte facilitates ion transport at fast charge and discharge rates. Large current densities are accompanied by large salt concentration gradients across the electrolyte. Nanostructured composite electrolytes have been proposed to enable the use of high energy density lithium metal anodes, but many questions about the interplay between the electrolyte morphology and the salt concentration gradient that forms under dc polarization remain unanswered. To address these questions, we use an in situ small-angle X-ray scattering technique to examine the nanostructure of a polystyrene- block -poly(ethylene oxide) copolymer electrolyte under dc polarization with spatial and temporal resolution. In the quiescent state, the electrolyte exhibits a lamellar morphology. The passage of ionic current in a lithium symmetric cell leads to the formation of concurrent phases: a disordered morphology near the negative electrode, lamellae in the center of the cell, and coexisting lamellae and gyroid near the positive electrode. The most surprising result of this study was obtained after the applied electric field was turned off: a current-induced gyroid phase grows in volume for 6 h in spite of the absence of an obvious driving force. We show that this reflects the formation of localized pockets of salt-dense electrolyte, termed concentration hotspots, under dc polarization. Our methods may be applied to understand the dynamic structure of composite electrolytes at appreciable current densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Differentiable multiphase flow model for physics-informed machine learning in reservoir pressure management

Accurate subsurface reservoir pressure control is extremely challenging due to geological heterogeneity and multiphase fluid-flow dynamics. Predicting behavior in this setting relies on high-fidelity physics-based simulations that are computationally expensive. Yet, the uncertain, heterogeneous properties that control these flows make it necessary to perform many of these expensive simulations, which is often prohibitive. To address these challenges, we introduce a physics-informed machine learning workflow that couples a fully differentiable multiphase flow simulator, which is implemented in the DPFEHM framework with a convolutional neural network (CNN). The CNN learns to predict fluid extraction rates from heterogeneous permeability fields to enforce pressure limits at critical reservoir locations. By incorporating transient multiphase flow physics into the training process, our method enables more practical and accurate predictions for realistic injection-extraction scenarios compared to previous works. To speed up training, we pretrain the model on single-phase, steady-state simulations and then finetune it on full multiphase scenarios, which dramatically reduces the computational cost. We demonstrate that high-accuracy training can be achieved with fewer than three thousand full-physics multiphase flow simulations – compared to previous estimates requiring up to ten million. This drastic reduction in the number of simulations is achieved by leveraging transfer learning from much less expensive single phase simulations.

25 ENERGY STORAGE↗

Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in 76 Ge

The observation of neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0⁢𝜈⁢𝛽⁢𝛽 experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0⁢𝜈⁢𝛽⁢𝛽-decay NME, in the key nucleus 76 Ge. Here, our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60$^{+1.28}_{−1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187$^{+205}_{−62}$ meV. The result is important for designing next generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.

Ab initio calculations↗

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗

Analytical and experimental studies of an optimum multisegment phased liner noise suppression concept

Results are presented from detailed analytical studies made to define methods for obtaining improved multisegment lining performance by taking advantage of relative placement of each lining segment. Properly phased liner segments reflect and spatially redistribute the incident acoustic energy and thus provide additional attenuation. A mathematical model was developed for rectangular ducts with uniform mean flow. Segmented acoustic fields were represented by duct eigenfunction expansions, and mode-matching was used to ensure continuity of the total field. Parametric studies were performed to identify attenuation mechanisms and define preliminary liner configurations. An optimization procedure was used to determine optimum liner impedance values for a given total lining length, Mach number, and incident modal distribution. Optimal segmented liners are presented and it is shown that, provided the sound source is well-defined and flow environment is known, conventional infinite duct optimum attenuation rates can be improved. To confirm these results, an experimental program was conducted in a laboratory test facility. The measured data are presented in the form of analytical-experimental correlations. Excellent agreement between theory and experiment verifies and substantiates the analytical prediction techniques. The results indicate that phased liners may be of immediate benefit in the development of improved aircraft exhaust duct noise suppressors.

Sawdy, D. T.↗

Thermochemical Stability of Ca2Yb8(SiO4)6O2 Apatite in Presence of Molten Calcium-Magnesium-AluminoSilicate (CMAS)

Thermochemical stability of ytterbium silicon oxyapatite Ca 2 Yb 8 (SiO 4 ) 6 O 2 (CYbS) in the presence of molten calcium-magnesium aluminosilicate (CMAS) has been investigated at elevated temperatures. CYbS apatite powder was synthesized from the constituent oxides via solid state reaction method. Hot pressed apatite substrates were exposed to molten CMAS at 1200, 1300, and 1400 °C for 1, 10, and 50 h. Development of phases in the interaction region of the heat-treated specimens was monitored using scanning electron microscopy, transmission electron microscopy, high angle annular dark field imaging, selected area electron diffraction and energy dispersive X-ray spectroscopy. Monoclinic cyclosilicate Ca 3 Yb 2 (Si 3 O 9 ) 2 formed from interaction of CYbS apatite with CaO in the CMAS melt at the apatite-CMAS reaction front and continued to nucleate and grow within the residual CMAS in diffusion couples annealed for 1-50 h at 1200 °C and those heat treated at 1300 °C for 1 h. Residual CMAS was Ca-depleted when cyclosilicate was present. Dendritic wollastonite CaSiO 3 was observed within the residual CMAS in couples annealed at 1200 and 1300 °C. Ingress of molten CMAS, because of its exponential decrease in viscosity, occurred through open pores and along the grain boundaries of the apatite substrates without any detectable chemical reaction at 1300 and 1400 °C. Results of this study indicate that Ca 2 Yb 8 (SiO 4 ) 6 O 2 apatite has the potential to mitigate the CMAS corrosion up to about 1200 °C but not at higher temperatures.

X-ray diffraction↗