Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thresholds”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Study of hydride reorientation and ductility of SRA ZIRLO and PRXA ZIRLO using ring compression testing (RCT)

Before dry storage, spent nuclear fuel is dried at high temperatures, which induces hydride dissolution. Radial hydride particles precipitate when cooled under a stress above a threshold stress. This process can severely degrade the ductility of cladding tubes. Zirconium alloys with different microstructures, such as grain boundary density, were quantified and studied to examine their effect on hydride reorientation. The ring compression loading mode creates a range of stress states depending on the location within the ring with respect to the applied load, producing various hydride morphologies that can be used to study the effects of microstructures on hydride reorientation. Finite element modeling was performed to calculate the hoop stress, which is the dominant stress, during hydride reorientation. These results were coupled with hydride morphology to determine the threshold stress for hydride reorientation, depending on the microstructure. In conclusion, the samples were then mechanically tested using ring compression to evaluate the ductility of zirconium alloys as a function of different hydride morphologies.

Hydride embrittlement↗
DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Beta-delayed charged-particle spectroscopy using TexAT

β-delayed charged-particle emission is a sensitive probe of three-body decays in light nuclei. Time Projection Chambers (TPCs) offer a significant advantage over traditional charged-particle spectroscopy techniques due to a low-energy threshold and a high-geometric efficiency (≈ 4π) which are essential for use with radioactive ion beams where the beam intensities are limited. The technique for high-sensitivity spectroscopy of β-delayed charged-particle emission is shown to be possible using the Texas Active Target (TexAT) TPC in conjunction with the General Electronics for TPCs (GET) system. The benchmark case studied was that of 12 N β-decay to the first α-unbound state in 12 C, the Hoyle state. Here, half-life and branching ratio measurements are presented and are in good agreement with previous studies. The efficacy of using TPCs to study such a near-threshold state and disentangle the three-body dynamics of the decay products is demonstrated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Parameter dependencies of the experimental nitrogen concentration required for detachment on ASDEX Upgrade and JET

While current tokamak experiments are beginning to use real-time feedback control systems to manage the plasma exhaust, future tokamaks still require validation of theoretical models used to predict the threshold impurity concentration required to sufficiently reduce the power and particle fluxes to the divertor. This work exploits new spectroscopic measurements of the divertor nitrogen concentration, c n , in partially detached N 2 -seeded H-mode plasmas on ASDEX Upgrade (AUG) and JET with the ITER-Like Wall (JET-ILW) to test the parameter dependencies of the power flowing to the outer divertor, P div,outer , and the separatrix electron density, n e.sep . A least-squares regression of the AUG measurements demonstrates that the threshold c n required for detachment scales as $c_n∝P^{1.19±0.32}_{div,outer}n^{-2.77±0.36}_{e,sep}$. This scaling of n e,sep is also consistent with the measurements from JET which, at constant P div,outer , show $c_n∝n^{-2.43±0.27}_{e,sep}$. The dependencies of P div,outer and n e.sep is demonstrated over at least a factor of two change in both parameters and indicates a stronger dependence on n e.sep in comparison to the Lengyel model, which could be due to the assumption in this model that the heat flux channel width is independent of density. This first assessment of detachment scaling with impurity seeding highlights the need for further analysis of the systematic uncertainties of the measurement and more consistent scenarios from more tokamaks to investigate the machine size scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Integral LOCA fragmentation test on high-burnup fuel

Increasing the fuel burnup limit in light water reactors is sought to enhance fuel cycle economics and requires establishing a technical basis. Experimental observations of severe fuel fragmentation under loss-of-coolant-accident conditions at Halden and Studsvik had raised the need for additional considerations during the development of this technical basis. These test data suggested that the burnup threshold for high burnup fuel fragmentation may be influenced by pre-transient power. Additional loss-of-coolant-accident test data is therefore valuable to complement these tests and enhance the current state of understating. Oak Ridge National Lab has developed the Severe Accident Test Station capable of examining the oxidation kinetics and accident response of irradiated fuel and cladding materials for design basis accident and beyond design basis accident scenarios. Severe Accident Test Station provides various temperature profiles, steam, and the thermal shock conditions necessary for integral loss of coolant accident testing, defueled oxidation quench testing, and high-temperature beyond design basis accident testing. Severe Accident Test Station has been successfully installed and demonstrated in the Irradiated Fuels Examination Laboratory at Oak Ridge National Lab. Furthermore, descriptions of the in-cell re-fabrication capabilities and assembly of the loss-of-coolant-accident test train are provided. Installation of the Severe Accident Test Station system and in-cell re-fabrication restores United States capability to examine postulated and extended loss-of-coolant-accident conditions on spent fuel and cladding and provides a platform for evaluating advanced fuel and accident-tolerant fuel cladding concepts. Lastly, three in-cell integral loss-of-coolant-accident test were performed in the Severe Accident Test Station and compared to the Nuclear Regulatory Commission sponsored loss-of-coolant-accident test as well as the original loss-of-coolant-accident test performed at Argonne National Lab. Finally, the results of these tests as well as all publicly available integral loss-of-coolant-accident test were used to validate a threshold for high burnup fuel fragmentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Note on the definitions of branching ratios of overlapping resonances

Branching ratios for the decay of hadrons with large width or near thresholds depend on their definition. We test different definitions and show that rather different branching ratios can be obtained. For wide resonances and for sequential decays with wide intermediate resonances, integration over the spectral functions is mandatory. The tests are performed exploiting the latest solution of the Bonn-Gatchina multi-channel analysis and published values for residues of light scalar mesons. For a resonance overlapping with a threshold, in case its pole lies in a non-adjacent sheet, we show how the total width, needed for the branching ratios, does not correspond to the imaginary part of the pole position. We use the Madrid-Krakow dispersive parameterizations to illustrate this situation with the f 0 (980).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Response of PETN detonators to elevated temperatures

In the present study, commercially available detonators with pentaerythritol tetranitrate (PETN) were subjected to elevated temperatures. The detonators were thermally ignited over a range of heating rates to measure ignition delay time and assess detonator violence. The violence of the detonator was quantified by measuring the velocity of the detonator closure disc (or “flyer”). The maximum flyer velocity of a thermally ignited detonator was comparable in magnitude to that obtained by initiating a room temperature pristine detonator with an exploding bridge wire (under the same confinement); however, the high flyer velocity was not an indication of deflagration to detonation transition (DDT) in the thermally ignited detonator. The detonator responded more violently than a thermally ignited detonator when initiated at 95% of the ignition delay time. Furthermore, inoperability thresholds were also measured by varying the detonator temperature and the threshold was found to be sensitive at detonator temperatures below the melting point of PETN.

42 ENGINEERING↗
DOE OSTI2022

Hydrodynamics and Load Shedding Behavior of a Variable-Geometry Oscillating Surge Wave Energy Converter (OSWEC)

In order to improve their long-term viability, wave energy converters (WECs) need to be able to shed loads when a threshold wave condition is exceeded. As shown by Tom et al. (2016) [1], provision of adjustable flaps within the body of an oscillating surge wave energy converter (OSWEC) allows wave energy to pass through the device. A control system may then be able to open and close the flaps when waves approaching the device exceed preset thresholds. The variable-geometry OSWEC (VG-OSWEC) concept studied in this paper is a bottom-hinged, rectangular wave paddle with five flaps of elliptical cross-section embedded into the face of the paddle. System ID tests were conducted on this VG-OSWEC device at a 1:14 scale in a wave basin. Free decay tests showed that the damping was distinctly nonlinear when the flaps were fully open, and the natural frequency increased by 40% when compared with the flaps in a fully closed configuration. Tests with regular wave conditions were used to develop the response amplitude operator for the rotational motion about the hinge. These response amplitude operator results when compared with numerical simulations run using WEC-Sim/WAMIT and ANSYS AQWA, show strong agreement with the flap open and closed conditions. The regular-wave condition measurements also show that the wave excitation moment about the hinge was reduced by up to 60% when the flaps were fully open. The experiments serve to demonstrate the potential of the variable geometry design to shed loads and survive harsh ocean environments.

load shedding↗
DOE OSTI2025

Velocity- and pressure-based metrics for estimating strike injuries during fish passage through hydro turbines

Hydropower is a crucial source of clean and reliable energy worldwide, and its importance will continue to grow. To enhance the environmental sustainability of hydropower development and operations, it is essential to predict the strike injury rates of fish passing through turbines accurately and cost-effectively. However, conventional experiments involving a large number of live fish are still commonly conducted in practice, and previous attempts mainly focused on the dose-response relationships at an individual level or relative comparisons of biological characterization between different conditions. Thus, this study proposes two novel strike metrics based on velocity and pressure (M V and M P ) measured by the cutting-edge Sensor Fish (SF) technology, designed to quantify the biological effects of strikes and collisions between fish and rigid hydraulic structures. SF measurements and 48-hour survival rates of live fish from a study conducted at the Ice Harbor Dam on the Snake River were used to establish thresholds for the two proposed metrics. These metrics and thresholds were then applied to three additional SF studies conducted at hydropower facilities featuring Kaplan turbines within the Columbia River basin. Results from the application of the strike metrics indicate that the estimated survival rates derived from both metrics agree well with previously reported live fish survival rates. Furthermore, the study identified specific passage regions with elevated risks of fish mortality due to strikes and collisions. Overall, the proposed strike metrics present a reliable and cost-effective approach to reducing reliance on live fish and lowering costs in fish passage studies.

Biological characterization↗
DOE OSTI2021

Hydrological projections over the Upper Indus Basin at 1.5 °C and 2.0 °C temperature increase

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Breakdown of the $v=1$ integer quantum Hall effect in the single particle and collective localization regimes

We report large signal current–voltage characteristics measurements near the $v=1$ integer quantum Hall plateau in a high mobility GaAs/AlGaAs sample that exhibits both single particle and collective localization. According to prior results, the central part of the $v=1$ plateau corresponds to single particle localization and hence an Anderson insulator. In contrast, the two regions in the flanks of this Anderson insulator were associated with collective localization and were identified with the integer quantum Hall Wigner solid. We find that, at low temperatures and under large bias currents, both the Anderson insulator and integer quantum Hall Wigner solid have a sharp onset of conduction above a threshold. Furthermore, this threshold shows an approximately monotonic dependence on the filling factor across the boundary between the two localized phases. Finally, to our surprise, we find that these large signal characteristics do not clearly distinguish between the Anderson insulator and the integer quantum Hall Wigner solid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Quantitative Method for In Situ Pump-Beam Metrology in 4D Ultrafast Electron Microscopy

In this work, we report a method for measuring spot size and focusing conditions of the femtosecond (fs) excitation laser in situ at the specimen location in 4D ultrafast electron microscopy (UEM). The method makes use of threshold laser ablation behaviors of thin amorphous carbon membranes. For Gaussian beam profiles and for ablation threshold fluence values, we analytically derive expressions describing the relationship between ablated hole size and the actual laser spot size. Using these expressions, we developed experimental procedures for characterizing the shape and spot size of the pump beam at the specimen. We demonstrate the viability of the approach for incident excitation wavelengths of 343 nm and 515 nm, thus illustrating the applicability of the method to a range of optical wavelengths without modification. Further, we show that by measuring ablated hole size as a function of focusing condition, a full metrological characterization of the Gaussian beam propagation properties can be performed. Finally, we find good agreement for spot sizes determined with this method and with those determined by extrapolation from measurements taken outside the microscope column. Overall, this method is a simple, cost-effective means for accurate and precise determination of key pump-beam parameters in situ at the specimen location in UEM experiments.

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

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗
DOE OSTI2024

Exploring the Influence of Linker Substitution and Ratios on Cooperative Framework Flexibility through the Mixed-Linker Approach

Cooperatively flexible metal–organic frameworks that exhibit step-shaped, or Type V-like, adsorption–desorption profiles can lower requisite pressure–temperature swings and thus energy input, necessary for an array of gas storage, delivery, and separations applications. However, such benefits are lost if the pressure threshold of the adsorption and desorption steps at a given temperature does not match the conditions dictated by the application, such as H 2 storage and delivery or olefin–paraffin separations. Unfortunately, the discovery of cooperatively flexible frameworks remains wholly serendipitous and de novo design remains impossible. Accordingly, there is a great need to further our understanding of flexibility such that we can intuitively derivatize known frameworks and, ultimately, design entirely new ones to meet the requisite conditions of energy-consumptive processes. In this work, we demonstrate that the mixed-linker, or multivariate, approach is a powerful tool for the derivation of a known flexible framework, with variances in linker substitution and ratio giving rise to a family wherein significant changes to the step-shaped adsorption–desorption profiles for multiple adsorbates are observed. Specifically, we report 12 isostructural mixed-linker derivatives of CdIF-13 (sod–Cd(benzimidazolate) 2 ) with six point-modified benzimidazole linkers spontaneously synthesized through prototypical solvothermal conditions. Furthermore, each is shown by PXRD to exhibit similar reversible flexibility to CdIF-13 and by TGA and DSC to be similarly thermally stable. Isothermal gas adsorption measurements with N 2 at 77 K, CO 2 at 195 K, and propane at 298 K demonstrate the wide-ranging, adsorbate-dependent effects of linker substitution and ratio, including dramatic reductions in the adsorption threshold pressure, evolution of multiple steps, and complete absence of observable adsorption. Aided by prior crystallographic characterization of CdIF-13, the observed trends were analyzed in an attempt to resolve the influence of linker substitution and ratios on structural behavior. In total, this work illustrates how the mixed-linker approach enables the synthesis and study of a large catalog of functionally modified frameworks, which can help identify the variables that influence flexibility and sorption behavior.

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

Strain Relaxation and Relative Defect Density with Thickness in MBE-Grown Ge 0.85 Sn 0.15 on Ge(001)

Germanium–tin (GeSn) alloys are emerging as promising materials for mid-infrared optoelectronics and silicon-compatible photonic devices, owing to their tunable direct bandgap. However, the growth of high-quality GeSn films with high Sn content remains challenging due to strain-induced defect formation. In this study, we investigate the role of film thickness on strain-induced relaxation, defect density, and Sn segregation. A series of five samples with varying thicknesses and ∼15% Sn-containing GeSn layers were grown, ranging from the critical thickness for strain relaxation to the onset of Sn segregation. All GeSn samples were analyzed using X-ray diffraction reciprocal space mapping (XRD-RSM) to explore the evolution of strain-induced relaxation as a function of thickness. Photoluminescence measurements reveal that increasing the GeSn thickness enhances strain relaxation while reducing defect-related emission, indicating a decrease in effective defect density prior to reaching the threshold thickness of GeSn layer. At a thickness of ∼150 nm, the GeSn layer shows the onset of Sn segregation, evident in the XRD-RSM spectrum, marking the threshold thickness for Sn segregation. This work defines an effective growth window in terms of thickness (35 to 150 nm) for fabricating relaxed, defect-suppressed GeSn layers with 15% Sn content. These findings emphasize the crucial role of thickness control in balancing strain relaxation and defect suppression, advancing the fabrication of high-quality, high Sn-content relaxed GeSn using molecular beam epitaxy.

Defects↗
DOE OSTI2020

Localized High Concentration Electrolytes for High Voltage Lithium-Metal Batteries: Correlation between the Electrolyte Composition and Its Reductive/Oxidative Stability

We demonstrate a first-principles screening methodology as an effective tool to explore electrolyte formulations for the new generation of high energy density rechargeable batteries. We study the liquid structure and electronic properties in dilute electrolytes, high concentration electrolytes (HCE), and localized high concentration electrolytes (LHCE), with focus on electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and bis(2,2,2-trifluoroethyl) ether (BTFE) as a diluent. We describe the solvation complexes in the dilute electrolyte and explore structural changes triggered by the increase in lithium salt concentration for HCEs and the diluent effects in LHCEs. In HCE formulations, there is a 4-fold coordination environment of lithium-ions as in the dilute electrolyte, but the number of lithium-ion interactions with O atoms from FSI- anions dominates. In these solutions, the ability of the FSI- anions to interact with multiple lithium-ions allows complex 3D network formation and influences the reductive/oxidative behavior of the electrolyte. Interestingly, in LHCEs, the BTFE diluent molecules do not change the 3D solution structure when diluting the HCE formulation from 5.49 to 3.83 M. However, there is a composition threshold where the structural and electronic behavior may change. We show that diluting the HCE electrolyte with BTFE down to 1.77 M breaks the three-dimensional solution structure into an island-like solvation complex. We relate these structural changes to the electronic properties of the electrolytes finding a causal relationship between the reductive/oxidative behavior and the lithium-oxygen interaction mechanisms in the solvated complexes. The coordination with lithium-ions lowers the electrolyte LUMO and HOMO levels: the higher is the number of interactions with lithium-ions, the more likely the solvent molecule, FSI- anion, or diluent molecule is to be reduced and the less likely it is to become oxidized. The evolution of the solvated ion structure in HCE and LHCE suggests a close connection to a corresponding change in the lithium-ion transport mechanisms for these electrolytes.

anions, molecules, oxygen, solvents, electrolytes↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗