Search NASA⌕ Search

SEARCH · Search NASA

Results for “Velocity Model”

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 235 records · Page 13

Autoignition of cyclohexane at low-to-intermediate temperatures: Rapid compression machine experiments and improved comprehensive chemical kinetic model

Cycloalkanes are an important class of hydrocarbons found in conventional and alternative transportation fuels. Being one of the simplest practical cycloalkanes, cyclohexane is often used as a surrogate in studies to understand the oxidation, ignition, and sooting characteristics of cycloalkanes more broadly. Although cyclohexane has been extensively studied, very limited datasets for cyclohexane oxidation at low-to-intermediate temperatures are available in the literature. Moreover, the literature autoignition data from rapid compression machines (RCMs) did not provide information on volume histories or experimental heat transfer characteristics which can significantly impact model predictions, thereby hindering the detailed validation of chemical kinetic models. It has also been observed that the cyclohexane kinetic models developed so far are not able to predict the literature data across a wide range of operating and mixture conditions. Hence, further studies are recommended to understand and improve model predictions for cyclohexane oxidation. In this study, RCM experiments were conducted using cyclohexane/air mixtures at varying compressed pressures (P C =15, 20, and 30 bar), equivalence ratios (ϕ=0.5, 1.0, and 2.0), and compressed temperatures (T C =670–960 K). The current RCM results of cyclohexane were found to complement well with the literature data, with the first-stage IDTs at high pressures and volume histories at all test conditions being first reported for better development and validation of chemical kinetic models. Additionally, a detailed chemical kinetic model for the cyclohexane oxidation has been developed and validated using both newly acquired experimental IDTs in this study and literature experimental studies (IDTs in RCMs and STs, species profiles in jet-stirred reactors, and laminar burning velocities). Unlike literature models, the proposed model was observed to provide an overall better prediction for experimental results over a wide range of conditions. Chemical kinetic analysis of the developed model was further conducted to identify the reaction pathways controlling the autoignition of cyclohexane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Modeled Performance of the ORNL Matched Filter Detector for Seismic Signatures of Multi-Axle Tractor Trailers

The MINOS program seeks methods to passively monitor the transfer of special materials, which are reproducible and quantitatively defensible. To date, we have demonstrated that short-period seismometers deployed near vehicle tracks can detect and monitor signatures of multi-axle tractor trailers impacting small road obstacles as means to detect and identify some transfer activities. In this report, we model the seismic velocity signal output from a tractor-trailer impacting a small obstacle (a template waveform) to parameterize a matched filter that identifies signals with similar shapes buried in noise, and that indicate repeating events (target waveforms). We then use this model to perform thousands of synthetic experiments, in which we apply our detector to identify noisy waveforms output from a similarly modeled vehicle data set; we process all data in the same manner as we process real data collected from the Oak Ridge National Laboratory (ORNL) seismic network. Our experiments, while limited, reveal that the matched filter shows a high discrimination capability, and at very low SNR values (at visibility thresholds). This means that template waveforms can match target waveforms that are sourced by models that share the same vehicle speed (target vehicles), and have a low probability of mis-identification. In particular, a matched filter designed to detect a five-axle tractor-trailer will reliably identify that target vehicle down to an SNR of 11.25 dB, and within ±0.1 km/hr, while successfully screening target vehicles with three axles.

42 ENGINEERING↗

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Here, we systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in Ref.Yang et al. (2023) with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and subhalos. We focus on the accuracy of the model in predicting halo density profiles at z = 0 and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within 10%–50% for most of the simulated (sub)halos, 50%–100% for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts.

Yang, Daneng [Chinese Academy of Sciences, Nanjing↗

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗

SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations

We present SIDM Concerto: 14 cosmological zoom-in simulations in cold dark matter (CDM) and self-interacting dark matter (SIDM) models based on the Symphony and Milky Way-est suites. SIDM Concerto includes one Large Magellanic Cloud– (LMC-) mass system (host mass ∼10 11 M ⊙ ), two Milky Way (MW) analogs (∼10 12 M ⊙ ), two group-mass hosts (∼10 13 M ⊙ ), and one low-mass cluster (∼10 14 M ⊙ ). Each host contains ≈2 × 10 7 particles and is run in CDM and one or more strong, velocity-dependent SIDM models. Our analysis of SIDM (sub)halo populations over seven subhalo mass decades reveals that (1) the fraction of core-collapsed isolated halos and subhalos peaks at a maximum circular velocity corresponding to the transition of the SIDM cross section from a v −4 to v 0 scaling; (2) SIDM subhalo mass functions are suppressed by ≈50% relative to CDM in LMC, MW, and group-mass hosts but are consistent with CDM in the low-mass cluster host; (3) subhalos’ inner density profile slopes, which are more diverse in SIDM than in CDM, are sensitive to both the amplitude and shape of the SIDM cross section. Our simulations provide a benchmark for testing SIDM predictions with astrophysical observations of field and satellite galaxies, strong lensing systems, and stellar streams. Data products are publicly available at doi:10.5281/zenodo.14933624.

dark matter↗

Graph Learning in Physical-informed Mesh-reduced Space for Real-world Dynamic Systems

This Git repository contains codes for the 'Graph Learning in Physical-informed Mesh-reduced Space for Real-world Dynamic Systems' paper that will be published in 2023 SIGKDD. This work uses physical-informed prior (PiP) information to learn and predict fluid dynamics in a reduced mesh space. We propose a two-stage graph-based model for fluid velocity field reconstruction and prediction. In the first stage, we learn a subgraph autoencoder to summarize the information in a mesh-reduced space using physical-informed priors. In the second stage, we learn a dynamics predictor to predict subgraph evolution. We demonstrate the effectiveness of our model on two fluid flow datasets: lid-driven cavity flow data and cylinder flow data. This code is also applicable for any other dynamic systems with corresponding data.

Lei, Bo↗

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal observing strategies for velocity-suppressed dark matter annihilation

Numerous particle models for the cosmological dark matter feature a pair-annihilation rate that scales with powers of the relative velocity between the annihilating particles. As a result, the annihilation rate in the central regions of a dark matter halo can be significantly lower than at the halo’s periphery for particular ambient gravitational potentials. While this might be offset by an increasing dark matter pair number density in the inner halo, it raises the question; what angular region for dark matter models with velocity-suppressed annihilation rates optimizes the signal to noise ratio? We consider simplified background models for galactic and extragalactic targets and demonstrate that the optimal observing strategy varies greatly case by case. Generally, a bright central source warrants an annular region of interest, while a flatter background warrants as large as possible an angular region, possibly including the central regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗

Laser-induced-fluorescence imaging of a spin-polarized ultracold neutral plasma in a magnetic field

In this work, we report Doppler-sensitive laser-induced-fluorescence (LIF) imaging of an ultracold neutral plasma in a magnetic field. Local values of ion density, hydrodynamic fluid velocity, temperature, and spin polarization are obtained using a fluorescence model based on velocity-resolved rate equations (REs) including the transfer of ions between states due to laser coupling and spontaneous emission. The RE approach captures optical pumping of ions into states that are not driven by the LIF excitation laser, and this is validated with experimental data. Combined molecular-dynamics and quantum-trajectories simulations verify that velocity-changing collisions have a negligible impact on the state population evolution for typical experimental conditions. Relative intensities of Zeeman components of the LIF spectra provide clear evidence that the ions are electron-spin-polarized when created by photoionization of magnetically trapped 88 Sr atoms. This probe opens many possibilities for studying thermal transport and the equilibration of neutral plasmas in overlapping regimes of strong coupling and magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Lagrangian perturbation theory in the presence of massive neutrinos

We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using ConvolutionLPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote 1 simulation suite, we find good accuracy down to r = 20 Mpc h -1 at redshift z = 0.5, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to r = 20 Mpc h -1 when comparing to simulated halos.

79 ASTRONOMY AND ASTROPHYSICS↗

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗

In Situ Determination of Thermoelastic Properties of Magnesite at High Pressure and Temperature With Implications to Seismic Detectability of Moderately Carbonated Lithologies in the Earth's Mantle

Abstract Subduction of carbonate‐bearing oceanic plates into Earth's interior recycles carbon from the surface to the deep mantle. The subducted carbonates can significantly affect mantle properties and dynamics. Magnesite is recognized as one of the major potential carbon hosts in the deep mantle because of its stability up to deep lower mantle conditions. However, despite many previous studies on the equation of state and elastic properties of magnesite, large discrepancies still exist for its elastic moduli and their pressure and temperature derivatives. Here we report in situ density and elastic wave velocity measurements on a natural magnesite at simultaneous high pressure‐temperature conditions up to ∼8 GPa–1073 K in a multi‐anvil apparatus using ultrasonic and synchrotron X‐ray techniques. Global fitting of the data set to finite strain equations yieldsK S0 = 114.0 ± 1.2 GPa, = 4.9 ± 0.3, = (−0.019 ± 0.002) GPa/K,G 0 = 68.6 ± 0.4 GPa,G′ = 1.6 ± 0.1, and = (−0.018 ± 0.001) GPa/K. Compared to other major upper mantle and transition zone minerals, magnesite has intermediateV P , the lowestV S , and the lowest density. Thus, magnesite possesses higherV P /V S ratio than other mantle minerals in normal mantle regions, whereas this feature is less pronounced in subduction zone environments. Modeling of the velocity profiles of carbonated lithologies along different geotherms suggests that moderately‐enriched magnesite domains are unlikely to be detected seismically in the Earth's mantle.

Geochemistry & Geophysics↗