Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗
DOE OSTI2021

Core-collapse contamination in photometric samples of Type Ia Supernovae

This is an exciting time for cosmology with type Ia supernovae (SNe Ia). The recentlyconcluded Dark Energy Survey SN programme (DES-SN) has obtained the largest anddeepest high-redshift cosmological SN Ia sample, and the Vera Rubin Observatory isexpected to observe at least one order of magnitude more SNe Ia in the next decade.In both these experiments, only a limited fraction (.10 per cent) of the SNe can bespectroscopically classified. This leaves us with large ‘photometric’ SN samples, withthe potential for significant contamination by core-collapse SNe that may bias SN Iacosmological measurements. This thesis demonstrates how this contamination can bemodelled and accounted for in current and future cosmological analyses.First, I present state-of-the-art simulations of the SN universe. These are designedto accurately model the population of SNe Ia, peculiar SNe Ia and core-collapse SNe, aswell as their host galaxies. To improve the diversity and quality of the simulated corecollapseSNe, I build a new library of core-collapse SN templates using spectroscopicand photometric (optical and near-ultraviolet) data of 67 core-collapse SNe from theliterature. I account for our incomplete knowledge of core-collapse SN properties bygenerating a set of SN simulations (rather than a single one), each exploring differentmodelling choices and template libraries. I then characterise selection effects in theDES-SN survey and incorporate them in the simulations, thus obtaining a series ofDES-like simulated SN samples that can be compared to the observed DES-SN data.The agreement between the simulations and data is excellent across many observed SNproperties, including Hubble residuals. These simulations are the first to reproduce theobserved photometric SN and host galaxy properties in high-redshift surveys with no fine-tuning of the input parameters.I use my simulation framework to train and test the performance of SuperNNova,a photometric SN classifier based on recurrent neural networks. I explore differenttraining and validation strategies and show that, across all the DES-SN simulationstested, SuperNNova reduces core-collapse SN contamination to 0.8–3.5 per cent. Ithen show that biases due to contamination on the equation-of-state of dark energy,w, are < 0:008 when using our reference SuperNNova model. This compares to anexpected statistical uncertainty on w from the DES-SN sample of 0:039, thus showing that contamination is not a limiting systematic for the cosmological analysis of theDES-SN sample.The results presented in this thesis are the foundation of the DES SN Ia cosmologicalanalysis; they also provide important implications for the future of SN cosmology,as they demonstrate that contamination is not expected to significantly degrade thecosmological figure of merit of the Rubin SN Ia analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electric Dipole Polarizability of Neutron Rich Nuclei

We report that insights into the equation of state of neutron rich matter obtained from the neutron skin thickness of 208 Pb are in sharp conflict with earlier measurements of the electric dipole polarizability. A set of accurately calibrated energy density functionals are used to highlight the tension and to articulate how a highly-intense Gamma Factory may help resolve the present dilemma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Current challenges in hydrate‐based desalination: Kinetic and thermodynamic perspective

Abstract Water scarcity is becoming a severe problem worldwide due to inadequate freshwater resources and swift population growth. Seawater desalination is one of the vital approaches to meet the demand for freshwater. However, energy and associated costs with conventional seawater desalination techniques are incentivizing non‐conventional water desalination processes. Water desalination using gas hydrates formation is one of the emerging non‐conventional processes. In this perspective article, recent advances in hydrate‐based seawater desalination (HBSD) have been critically analyzed to outline a future path towards a clean and efficient hydrate‐based desalination process. It provides a detailed comparison of various processes developed over decades, and measured desalination efficiencies with their process details. Moreover, the current challenges, limitations, and future perspectives of hydrate‐based desalination are also discussed. The study also recapitulates the thermodynamics and kinetics aspects of the hydrate‐based desalination process. In addition, various factors controlling the desalination efficiencies, such as control of the separation of hydrate crystals, salt deposition on hydrate particles, and hydrate morphology, were thoroughly investigated with their proposed process designs. The kinetics of hydrate formation is also assessed, with the possibility of a zero‐induction regime and its consequent impact on hydrate morphology. The current capabilities of the thermodynamics models (Gibbs energy minimization + electrolyte equation of state) were discussed using various commercially available software. Additionally, the role of hydrate promotors is also discussed, which can reduce the higher cost associated with the hydrate‐based desalination process.

Khan, Muhammad Naveed↗
DOE OSTI2022

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An integrated EOS, pore‐crush, strength and damage model framework for near‐field ground‐shock

Abstract An integrated Equation of State (EOS) and strength/pore‐crush/damage model framework is provided for modeling near to source (near‐field) ground‐shock response, where large deformations and pressures necessitate coupling EOS with pressure‐dependent plastic yield and damage. Nonlinear pressure‐dependence of strength up to high‐pressures is combined with a Modified Cam‐Clay‐like cap‐plasticity model in a way to allow degradation of strength from pore‐crush damage, what we call the “Yp‐Cap” model. Nonlinear hardening under compaction allows modeling the crush‐out of pores in combination with a fully saturated EOS, that is, for modeling partially saturated ground‐shock response, where air‐filled voids crush. Attention is given to algorithmic clarity and efficiency of the provided model, and the model is employed in example numerical simulations, including finite element simulations of underground explosions to exemplify its robustness and utility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Cookoff of Powdered and Pressed Explosives using a Micromechanics Pressurization Model

Cookoff experiments of powdered and pressed TATB-based plastic bonded explosives (PBXs) have been modeled using a pressure-dependent universal cookoff model (UCM) in combination with a micromechanics pressurization (MMP) model described in a companion paper. The MMP model is based on the accumulation of decomposition gases at nucleation sites that load the surrounding TATB crystals and binder. This is the first cookoff model to use an analytical mechanics solution for compressibility and thermal expansion to describe internal pressurization caused by both temperature and decomposition occurring within closed-pore explosives. This approach produces more accurate predictions of ignition time and pressurization within high-density explosives than simple equation-of-state models. The current paper gives details of the reaction chemistry, model parameters, predicted uncertainty, and validation using experiments from multiple laboratories with errors less than 6 %. The UCM/MMP model framework gives more accurate thermal ignition predictions for high density explosives that are initially impermeable to decomposition gases.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Thermodynamically consistent versions of approximations used in modelling moist air

Abstract Some existing approaches to modelling the thermodynamics of moist air make approximations that break thermodynamic consistency , such that the resulting thermodynamics does not obey the first and second laws or has other inconsistencies. Recently, an approach to avoid such inconsistency has been suggested: the use of thermodynamic potentials in terms of their natural variables , from which all thermodynamic quantities and relationships (equations of state) are derived. In this article, we develop this approach for unapproximated moist‐air thermodynamics and two widely used approximations: the constant‐ approximation and the dry heat capacities approximation. The (consistent) constant‐ approximation is particularly attractive because it leads to, with the appropriate choice of thermodynamic variable, adiabatic dynamics that depend only on total mass and are independent of the breakdown between water forms. Additionally, a wide variety of material from different sources in the literature on thermodynamics in atmospheric modelling is brought together. It is hoped that this article provides a comprehensive reference for the use of thermodynamic potentials in atmospheric modelling, especially for the three systems considered here.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Islands in the fluid: islands are common in cosmology

We discuss the possibility of entanglement islands in cosmological spacetimes with a general perfect fluid with an equation of state w. We find that flat universes with time-symmetric slices where the Hubble parameter vanishes always have islands on that slice. We then move away from such slices, considering still universes with a general perfect fluid. Under the local thermal equilibrium assumption, the comoving entropy density s c is constant. As a result, the conditions for an island become an inequality between the energy density (or Hubble parameter) and the temperature at some time of normalization. The consequences are that islands can exist for practically all fluids that are not radiation, i.e. w ≠ 1/3. We also discuss the ramifications of our results for universes with spatial curvature. Finally, we show that islands occur in the Simple Harmonic Universe model which has no classical singularity at the background level, in contrast to all previous examples where islands occurred only in space-times with singularities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Integrating out beyond tree level and relativistic superfluids

We revisit certain subtleties of renormalization that arise when one derives a low-energy effective action by integrating out the heavy fields of a more complete theory. Usually these subtleties are circumvented by matching some physical observables, such as scattering amplitudes, but a more involved procedure is required if one is interested in deriving the effective theory to all orders in the light fields (but still to fixed order in the derivative expansion). As a concrete example, we study the U(1) Goldstone low-energy effective theory that describes the spontaneously broken phase of a $\phi$ 4 theory for a complex scalar. Working to lowest order in the derivative expansion, but to all orders in the Goldstones, we integrate out the radial mode at one loop and express the low-energy effective action in terms of the renormalized couplings of the UV completion. This yields the one-loop equation of state for the superfluid phase of (complex) $\phi$ 4 . We perform the same analysis for a renormalizable scalar SO(N) theory at finite chemical potential, integrating out the gapped Goldstones as well, and confirm that the effective theory for the gapless Goldstone exhibits no obvious sign of the original SO(N) symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Horizon radiation reaction forces

Using Effective Field Theory (EFT) methods, we compute the effects of horizon dissipation on the gravitational interactions of relativistic binary black hole systems. We assume that the dynamics is perturbative, i.e it admits an expansion in powers of Newton’s constant (post-Minkowskian, or PM, approximation). As applications, we compute corrections to the scattering angle in a black hole collision due to dissipative effects to leading PM order, as well as the post-Newtonian (PN) corrections to the equations of motion of binary black holes in non-relativistic orbits, which represents the leading order finite size effect in the equations of motion. The methods developed here are also applicable to the case of more general compact objects, eg. neutron stars, where the magnitude of the dissipative effects depends on non-gravitational physics (e.g, the equation of state for nuclear matter).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Cumulants of multiple conserved charges and global conservation laws

We analyze the behavior of cumulants of conserved charges in a subvolume of a thermal system with exact global conservation laws by extending a recently developed subensemble acceptance method (SAM) to multiple conserved charges. Explicit expressions for all diagonal and off-diagonal cumulants up to sixth order that relate them to the grand canonical susceptibilities are obtained. The derivation is presented for an arbitrary equation of state with an arbitrary number of different conserved charges. The global conservation effects cancel out in any ratio of two second order cumulants, in any ratio of two third order cumulants, as well as in a ratio of strongly intensive measures Σ and Δ involving any two conserved charges, making all these quantities particularly suitable for theory-to-experiment comparisons in heavy-ion collisions. We also show that the same cancellation occurs in correlators of a conserved charge, like the electric charge, with any non-conserved quantity such as net proton or net kaon number. The main results of the SAM are illustrated in the framework of the hadron resonance gas model. We also elucidate how net-proton and net-Λ fluctuations are affected by conservation of electric charge and strangeness in addition to baryon number.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

High pressure raman spectroscopy and X-ray diffraction of K2Ca(CO3)2 bütschliite: multiple pressure-induced phase transitions in a double carbonate

Abstract The crystal structure and bonding environment of K 2 Ca(CO 3 ) 2 bütschliite were probed under isothermal compression via Raman spectroscopy to 95 GPa and single crystal and powder X-ray diffraction to 12 and 68 GPa, respectively. A second order Birch-Murnaghan equation of state fit to the X-ray data yields a bulk modulus, $${K}_{0}=46.9$$ K 0 = 46.9 GPa with an imposed value of $${K}_{0}^{\prime}= 4$$ K 0 ′ = 4 for the ambient pressure phase. Compression of bütschliite is highly anisotropic, with contraction along the c -axis accounting for most of the volume change. Bütschliite undergoes a phase transition to a monoclinic C 2/ m structure at around 6 GPa, mirroring polymorphism within isostructural borates. A fit to the compression data of the monoclinic phase yields $${V}_{0}=322.2$$ V 0 = 322.2 Å 3 $$,$$ , $${K}_{0}=24.8$$ K 0 = 24.8 GPa and $${K}_{0}^{\prime}=4.0$$ K 0 ′ = 4.0 using a third order fit; the ability to access different compression mechanisms gives rise to a more compressible material than the low-pressure phase. In particular, compression of the C 2/ m phase involves interlayer displacement and twisting of the [CO 3 ] units, and an increase in coordination number of the K + ion. Three more phase transitions, at ~ 28, 34, and 37 GPa occur based on the Raman spectra and powder diffraction data: these give rise to new [CO 3 ] bonding environments within the structure.

36 MATERIALS SCIENCE↗