Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Permitting”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗
DOE OSTI2020

Getting to the Core of PARAFAC2, A Nonnegative Approach

In this paper, the authors present a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Using small building blocks to assemble ultra-complex, multifaceted metal-organic frameworks with zeolitic, mesoporous subnetwork

The assembly of ultra-complex structures from simple building units remains a long-term challenge in chemistry. Using small molecular building blocks (MBBs) in a mixed-ligand approach permitted the assembly of unprecedented metal-organic frameworks (MOFs), M-kum-MOF-1 (M = Y, Tb), exhibiting extra-large mesoporous cavities with small access windows. The ultra-complex cage of M-kum-MOF-1 consists of 240 vertices bridged by 432 edges, leading to a 194 faces-containing tile. This tile exhibits more faces than in any periodic structures (zeolites, MOFs, metal-organic polyhedra [MOPs], etc.) known to date. M-kum-MOF-1 not only possess zeolitic features (anionic framework), but they also contain an underlying wse zeolitic topology, which is observed for the first time.

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

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2023

Experimental and numerical investigation of flame stabilization and pollutant formation in matrix stabilized ammonia-hydrogen combustion

Ammonia (NH 3 ) is a carbon-free fuel that offers an attractive alternative for reducing greenhouse gas emissions. However, the slow flame speed, low heating value, and emissions of nitrogen-containing pollutants present significant issues for practical combustion applications. To address these issues, we investigate the use of matrix stabilized combustion. In this type of burner, combustion is performed within an inert porous ceramic foam, heat is recirculated by solid conduction and radiation, which enhances flame speed and combustion stabilization, thereby permitting combustion over a wide range of equivalence ratio conditions. We present a new porous media burner (PMB) capable of stabilizing NH 3 /air flames at ambient conditions. An extensive experimental characterization of the stability of this burner is conducted with up to 30% by volume of hydrogen (H 2 ) in the fuel stream. A 15:1 turndown ratio is demonstrated, with a high thermal power density of 62 MW m -3 . Concentrations of NO, unburnt NH 3 , and H 2 in the exhaust stream are measured. Two regimes are identified for low NO operation: rich and very lean. For rich conditions, NO emissions decrease with increasing equivalence ratio and decreasing H 2 blending. Unburnt NH 3 emissions follow opposite trends. These measurements are complemented by simulations in which the burner is represented by a coupled solid-gas reactor network. In conclusion, this model captures the burner’s pollutant emissions to good accuracy and is used to analyze the mechanisms of pollutant formation.

42 ENGINEERING↗
DOE OSTI2024

Finding a common ground for RCM experiments. Part B: Benchmark study on ethanol ignition

Rapid compression machines (RCMs) are widely used to investigate gas phase reaction kinetics of various kind of fuels at application relevant conditions. In principle, the operation of an RCM is based on the idea of compressing a homogenous pre-mixed fuel-air mixture by a piston. Usually creviced pistons ensure a homogenous adiabatic core in the center of the reaction chamber which permits the assumption of an isentropic relation between the measured pressure and gas temperature. Despite the ideal core gas compression, non-ideal effects such as heat loss, differences in the compression behavior, and ultimately non-standardized design and operation of rapid compression machines lead to different experimental results in different facilities at nominally the same end of compression conditions. In this study ignition delay times of ethanol are investigated at four different conditions in five independent RCMs. As expected, the raw results of the different facilities indeed show notable differences at the same end of compression conditions. However, according to the adiabatic core hypothesis the agreement between kinetic simulations and experiments should be consistent for all facilities provided that the facility effects are correctly accounted for. To elaborate upon this hypothesis, a kinetic mechanism is optimized to reflect the experimental results of all facilities. In the end, the optimized mechanism predicts all experimental data within the expected uncertainty. This confirms the reliability of RCM experiments for kinetic investigations and the validity of the effective volume approach in simulating RCM data.

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

Two-dimensional burn velocity and analysis of burn products of nitromethane at high pressure

We have studied the laser-initiated deflagration of nitromethane (CH 3 NO 2 ) under high static pressure in the diamond anvil cell. Time-resolved images of the deflagrations have been obtained using intensified CCDs (ICCD). In contrast with previous work, we rely on spontaneous emission from the reaction, rather than changes in the speckle pattern produced by artificial illumination of the sample. Furthermore, as opposed to the 1D records obtained previously with streak cameras, ICCDs permit imaging of the burn in two dimensions providing the ability to directly observe anisotropic deflagration behavior. Here, we report several examples of this behavior and discuss its possible origins. We have also investigated the products of the reaction using Raman spectroscopy. At pressures below 25 GPa, the burn product is observed to be opaque and has a Raman spectrum consistent with a carbonaceous soot. At pressures above 25 GPa, the burn product is observed to be transparent (consistent with earlier reports) and its Raman spectrum reveals the presence of molecular N 2 and a nitrogen-carbonate species. The latter species can be recovered to ambient pressure and has a Raman spectrum consistent with ammonium carbonate or ammonium bicarbonate. To aid interpretation of the Raman spectra, measurements have also been made on the isotopologues CH 3 15 NO 2 and 13 CH 3 NO 2 . This work establishes the importance of advanced time-resolved imaging to reveal the details of deflagration under high pressure and further advances our understanding of the corresponding chemistry.

Deflagration↗
DOE OSTI2023

On the use of a multigrid-reduction-in-time algorithm for multiscale convergence of turbulence simulations

Simulations of turbulent flow present challenges in terms of accuracy and affordability on modern highly-parallel computer architectures. A multigrid-reduction-in-time algorithm is used to provide a framework for separately evolving different scales of turbulence and for parallelizing the temporal domain, thereby increasing the concurrency. It is hypothesized that the space–time locality of the small scales of turbulence can be used to circumvent difficulties in applying temporal multigrid to flows dominated by inertial physics. For algorithms that fall well short of spectral accuracy (fourth-order is used in this work) attention must be paid to the accuracy of features on scales transferred between multigrid levels. Numerical experiments were performed using implicit large-eddy simulation. Results from applying the approach to an infinite-Reynolds number Taylor–Green flow and a double-shear flow at a Reynolds number of 11650 provide strong evidence that the approach has merit. The multigrid-reduction-in-time framework can be used to parallelize the temporal domain of a high-Reynolds-number turbulent flow and permit independent convergence of different scales. Establishing this foundation allows for future research in reducing the wall-clock time to solve turbulent flows while retaining the same accuracy as sequential solvers. In conclusion, current performance results from parallelizing the temporal domain are not competitive with those from sequential-in-time methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Recurrent neural network-based multiaxial plasticity model with regularization for physics-informed constraints

We report a recurrent neural network (RNN) based model is developed as a surrogate to predict nonlinear plastic response under multiaxial loading. The RNN-based model is trained and tested on stress versus strain curves generated using a numerical solution based on the classical radial return method. Besides simply learning the basic constitutive relationship, a novel approach is taken to enforce certain physical conditions. Specifically, regularization is employed to maintain non-negative plastic power density throughout the loading history thereby ensuring monotonically increasing plastic work and thermodynamic consistency. Enforcing physics in this manner permits coupling of the data-driven RNN approach with physics-based knowledge and laws. This has the effect of reducing the necessary amount of data and ensuring known physical laws are not violated. Since, once trained, the model need not perform the expensive task of solving nonlinear equations, its efficiency is orders of magnitude greater than its numerical counterpart. The RNN-based model has been trained on varied sets of data and the accuracy on test datasets validated. The developed model is general and robust and has widespread application such as in the simulation of metal forming, large scale plasticity, and part life prediction.

42 ENGINEERING↗
DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Testing SOAR tools in use

Investigations within Security Operation Centers (SOCs) are tedious as they rely on manual efforts to query diverse data sources, overlay related logs, correlate the data into information, and then document results in a ticketing system. Security Orchestration, Automation, and Response (SOAR) tools are a relatively new technology that promise, with appropriate configuration, to collect, filter, and display needed diverse information; automate many of the common tasks that unnecessarily require SOC analysts’ time; facilitate SOC collaboration; and, in doing so, improve both efficiency and consistency of SOCs. There has been no prior research to test SOAR tools in practice; hence, understanding and evaluation of their effect is nascent and needed. Here, in this paper, we design and administer the first hands-on user study of SOAR tools, involving 24 participants and six commercial SOAR tools. Our contributions include the experimental design, itemizing six characteristics of SOAR tools, and a methodology for testing them. We describe configuration of a cyber range test environment, including network, user, and threat emulation; a full SOC tool suite; and creation of artifacts allowing multiple representative investigation scenarios to permit testing. We present the first research results on SOAR tools. Concisely, our findings are that: per-SOC SOAR configuration is extremely important; SOAR tools increase efficiency and reduce context switching, although with potentially decreased ticketing accuracy/completeness; user preference is slightly negatively correlated with their performance with the tool; internet dependence varies widely among SOAR tools; and balance of automation with assisting decision making is preferred by senior participants. We deliver a public user- and tool-anonymized and -obfuscated version of the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

ToF-SIMS spectral data analysis of Paenibacillus sp. 300A biofilms and planktonic cells

Analysis of bacterial biofilms is particularly challenging and important with diverse applications from systems biology to biotechnology. Among the variety of techniques that have been applied, time-of-flight secondary ion mass spectrometry (ToF-SIMS) has many promising features in studying the surface characteristics of biofilms. ToF-SIMS offers high spatial resolution and high mass accuracy, which permit surface sensitive analysis of biofilm components. Thus, ToF-SIMS provides a powerful solution to addressing the challenge of bacterial biofilm analysis. This dataset covers ToF-SIMS analysis of Paenibacillus sp. 300A (300A) isolated from the Hanford site in Richland, WA. The strain is known to have metal and sulfur reducing properties and can be used for bioremediation, wastewater treatment, bioengineering and technology development. There is a current need to identify small molecules and fragments produced from bacterial biofilms. Static ToF-SIMS spectra of 300A were obtained using an IONTOF TOF-SIMS V instrument equipped with a 25 keV Bi 3 + metal ion gun. Identified molecules and molecular fragments are compared against known biological databases and the reported peaks have at least 65 ppm mass accuracy. These molecules range from lipids and fatty acids to flavonoids, quinolones, and other naturally occurring organic compounds. It is anticipated that the spectral identification of key peaks will assist detection of metabolites, extracellular polymeric substance molecules like polysaccharides, and biologically relevant small molecules using ToF-SIMS in future surface and interface research of bacterial biofilms.

Biofilms↗
DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗
DOE OSTI2021

Insights for Canadian electricity generation planning from an integrated assessment model: Should we be more cautious about hydropower cost overruns?

Hydropower accounts for approximately 60% of electricity generation in Canada, with growth expected in the coming decades as part of renewable energy transitions; however, frequent cost overruns threaten the viability of this growth. Here, using the integrated assessment model GCAM, we develop an endogenous representation of hydropower for Canada that accounts for market dynamics, thus permitting analysis of hydropower competition with other electricity generation technologies, both with and without cost overruns. Results show that modelling hydropower resources endogenously increases Canadian hydropower deployment relative to an assumption of fixed hydropower production, from 417 to 495 TWh annually by 2050. In scenarios that apply cost overruns at historical levels, hydropower loses market share to more easily scalable technologies like wind power. When including high cost overrun assumptions, the model determines that hydropower falls from about 73% to 65% of Canadian electricity generation by 2050, while wind power increases from about 8% to 11%. Countries may be better able to achieve electrification and renewable energy targets at lower cost by avoiding large-scale, overrun-prone hydropower and nuclear generation projects. Model results support that cost overruns are important considerations for policy decisions related to electricity sector development in Canada and elsewhere.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Transport of coexisting Ni-Cu sulfide liquid and silicate melt in partially molten peridotite

Transport of coexisting sulfide and silicate melts in partially molten peridotite contributes to the redistribution of chalcophile elements within the upper mantle as well as the genesis of magmatic Ni-Cu sulfide deposits, but has not been investigated systematically. Using laboratory experiments, theoretical calculations, and X-ray synchrotron microtomography, this study documents the topology and considers controls on the extraction of two immiscible liquids during partial melting of mantle peridotite. Under hydrostatic conditions, the measured dihedral angle at silicate melt-mineral-mineral junctions is 13.7–21.3°. Silicate melt is distributed along grain edges forming incompletely interconnected melt channels that disconnected by some dead ends at melt fractions ~7–9 vol%. Application of theoretically predicted permeability (k ~10 -14 -10 -16 m 2 ) permits estimation of the extraction velocity of silicate melt of 0.7-11.1 μm/day within a single interconnected melt channel. In the absence of silicate melt, isolated sulfide droplets (3.77 vol%) show a sulfide-olivine-olivine dihedral angle of 91.5–101.3°. However, in the presence of silicate melt, sulfide droplets (average size ~2.53±2.14μm, 1σ) are partially surrounded by silicate melt and stranded in triple junctions or melt pockets due to the limitation of the smallest dimension (0.3μm) of melt channels. Thus, the extraction of sulfide liquid is highly restricted by these dead ends and the smallest dimension of melt channels during porous flow of silicate melt. In contrast, during large-strain shear deformation (strain ~1.6–2.5), initially stranded sulfide droplets were elongated and extracted with silicate melt into liquid-rich sheets with a length of several hundred microns, constantly oriented at 14.3±4.5°to the shear plane and antithetic to the shear direction. The angle is lower than that (18-30°) of those sheets containing sulfide liquid only, indicating that silicate melt dominates liquid-rich sheets. Driven by stress, silicate melt-dominated liquid-rich sheets open the appropriately oriented grain boundaries between silicate minerals, thus providing an efficient pathway for the extraction of sulfide liquid during deformation. When such sheet-like channels remain open, sulfide droplets (> millimeter-scale) can be potentially mobile through high strain domains of the upper mantle, contributing to the addition of chalcophile elements and the fertilization of the lithospheric mantle.

high P-T experiment↗