SEARCH · Engineering Papers

Results for “Shape function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Confinement-Controlled Rearrangements in Dioxolane Upgrading on H-ZSM-5 Revealed by Periodic DFT

Zeolitic Brønsted acid sites catalyze carbocation rearrangements central to upgrading biomass-derived oxygenates. Here we elucidate the mechanism of dioxolane conversion to methyl ethyl ketone and isobutanal on H-ZSM-5 using periodic density functional theory on the MFI model, complemented by ab initio molecular dynamics to probe confinement effects. Dioxolane adsorption at the Brønsted site is followed by protonation-assisted ring opening to form an oxocarbenium intermediate stabilized by the deprotonated framework. From this common intermediate, selectivity is governed by two competing rearrangements, namely, the 1,2-hydride shift with a free-energy barrier of 18.05 kcal mol –1 at 498 K leading toward MEK, and the 1,2-methyl shift with a higher barrier of 25.40 kcal mol –1 leading toward isobutanal. The hydride-shift channel is kinetically preferred over the methyl-shift channel, lowering the isobutanal/MEK ratio below the 3:1 limit expected for equal branching. Adsorption thermodynamics further indicate stronger stabilization of MEK than isobutanal within ZSM-5 channels, suggesting that confinement-controlled binding can bias product distributions in addition to intrinsic rearrangement barriers. These results highlight how Brønsted acidity and pore confinement jointly shape the rearrangement landscape in MFI zeolites.

adsorption↗

DOE OSTI2021

Optimal deformation procedure for hybrid adaptive x-ray mirror based on mechanical and piezo-driven bending system

Here, a hybrid deformable x-ray mirror consisting of a mechanical bender and a bimorph deformable mirror has been developed to realize adaptive optical systems, such as zoom condenser optics, for synchrotron-radiation-based x-ray microscopy. In the developed system, both bending mechanisms comprehensively contribute to the formation of the target mirror shape and can narrow the role of piezoelectric actuators, thereby enabling a more stable operation. In this study, the behavior of the bimorph mirror under the clamped condition was investigated, and the sharing of the deformation amount for each bending mechanism was optimized to minimize the amplitude of the voltage distribution of the bimorph mirror.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Microstructure development in a cold worked SRF Nb sheet undergoing 700°C-900°C/3h heat treatments.

Bulk Nb for TESLA shaped SRF cavities is a mature technology, with high gradient (>35 MV/m) and high Q cavities (>1010, ~2K) routinely fabricated. Significant advances are in order to push Q?s to 1011 (~2K), and involve modifications to the sub-surface Nb layers by impurity doping (N, O). In order to achieve the lowest surface resistance any trapped flux needs to be expelled for cavities to reach high Q?s. There is clear evidence that cavities fabricated from polycrystalline sheets meeting current specifications require higher temperatures beyond 800oC leads to better flux expulsion, and hence improves Q. Recently, cavities fabricated with a non-traditional Nb sheet with initial cold work due to cold rolling expelled flux better after 800°C/3h heat treatment than cavities fabricated using fine-grain polycrystalline Nb sheets. Here, we analyze the microstructure development in Nb from the vendor supplied cold work non-annealed sheet that was fabricated into an SRF cavity as a function of heat treatment building upon the methodology development to analyze microstructure being developed by the FSU-MSU-UT-Austin- JLAB collaboration. The results indicate correlation between full recrystallization and better flux expulsion.

Balachandran, S.↗

DOE OSTI2022

Galaxy morphologies revealed with Subaru HSC and super-resolution techniques. I. Major merger fractions of L UV ∼ 3–15 L*UV dropout galaxies at z ∼ 4–7

Abstract We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z ∼ 4–7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ∼1″ to $\lesssim \!\!{0{^{\prime \prime }_{.}}1}$, which is comparable to that of the Hubble Space Telescope, allowing us to identify z ∼ 4–7 bright major mergers at a high completeness value of $\gtrsim \!\!90\%$. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z ∼ 4, 5, 6, and 7, respectively, in a UV luminosity range of LUV ∼ 3–$15\, L_{\rm UV}^*$ corresponding to −24 ≲ MUV ≲ −22. The major merger fractions are estimated to be $f_{\rm merger}\sim 10\%$–$20\%$ at z ∼ 4 and $\sim 50\%$–$70\%$ at z ∼ 5–7, which shows no fmerger difference compared to those of a control faint galaxy sample. Based on the fmerger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z ∼ 4–7 galaxy UV LFs. While these two effects partly explain the DPL shape at LUV ∼ 3–$10\, L_{\rm UV}^*$, the DPL shape cannot be explained at the very bright end of $L_{\rm UV}\gtrsim 10\, L_{\rm UV}^*$, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs.

Shibuya, Takatoshi↗

DOE OSTI2025

Non-Gaussianity from explicit U(1)-breaking interactions

We investigate primordial non-Gaussianity (NG) arising from the explicit U(1) symmetry-breaking interactions during inflation involving a nearly massless axial component of a complex scalar field P. We analyze the induced NG parameter f NL under scenarios where the axial field functions as either a curvaton or cold dark matter (CDM). In the curvaton framework, there is a conventional contribution to the local NG of f NL ≃ -O(1). Additional positive local NG can result from either the self-interactions of axial field fluctuations, their interactions with a light radial partner, or kinetic mixing with the inflaton via U(1) symmetry-breaking terms. We identify parameter regions where the interactions lead to cancellations, suppressing the overall local NG to |f loc NL | ≲ O(0.1). In the CDM scenario, these interactions enhance the NG in the isocurvature fluctuations. Moreover, interactions between the axial field and another light scalar, such as a curvaton, can generate O(1) curvature NG signals and significant mixed curvature-isocurvature NGs that are within the reach of future experiments with σ(f loc NL ) ∼ 1. We also explore the role of a heavy radial field in generating oscillating correlation signals, noting that such signals can dominate the shape of the mixed adiabatic-isocurvature bispectrum. In certain cases, an oscillatory isocurvature bispectrum signal may be observable in the future, aiding in distinguishing between certain types of the U(1)-breaking self-interactions of the axial field.

axions↗

DOE OSTI2023

Crystallographic variant mapping using precession electron diffraction data

In this work, we developed three methods to map crystallographic variants of samples at the nanoscale by analyzing precession electron diffraction data using a high-temperature shape memory alloy and a VO2 thin film on sapphire as the model systems. The three methods are (I) a user-selecting-reference pattern approach, (II) an algorithm-selecting-reference-pattern approach, and (III) a k-means approach. In the first two approaches, Euclidean distance, Cosine, and Structural Similarity (SSIM) algorithms were assessed for the diffraction pattern similarity quantification. We demonstrated that the Euclidean distance and SSIM methods outperform the Cosine algorithm. We further revealed that the random noise in the diffraction data can dramatically affect similarity quantification. Denoising processes could improve the crystallographic mapping quality. With the three methods mentioned above, we were able to map the crystallographic variants in different materials systems, thus enabling fast variant number quantification and clear variant distribution visualization. The advantages and disadvantages of each approach are also discussed. We expect these methods to benefit researchers who work on martensitic materials, in which the variant information is critical to understand their properties and functionalities.

Crystallographic variant mapping↗

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 OSTI2024

3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources

3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Unconventional gate-induced superconductivity in transition-metal dichalcogenides

Superconductivity in few-layer semiconducting transition-metal dichalcogenides (TMDs) can be induced by field-effect doping through ionic-liquid gating. While several experimental observations have been collected over the years, a fully consistent theoretical picture is still missing. Here we develop a realistic framework that combines the predictive power of first-principles simulations with the versatility and insight of Bardeen-Cooper-Schrieffer gap equations to rationalize such experiments. The multivalley nature of semiconducting TMDs is taken into account, together with the doping- and momentum-dependent electron-phonon and Coulomb interactions. Consistently with experiments, we find that superconductivity occurs when the electron density is large enough that the 𝑄 valleys get occupied, as a result of a large enhancement of electron-phonon interactions. Despite being phonon driven, the superconducting state is predicted to be sensitive to Coulomb interactions, which can lead to the appearance of a relative sign difference between valleys and thus to a 𝑠 +− character. We discuss qualitatively how such scenario may account for many of the observed physical phenomena for which no microscopic explanation has been found so far, including in particular the presence of a large subgap density of states, and the sample-dependent dome-shaped dependence of 𝑇 𝑐 on accumulated electron density. Our results provide a comprehensive analysis of gate-induced superconductivity in semiconducting TMDs, and introduce an approach that will likely be valuable for other multivalley electronic systems, in which superconductivity occurs at relatively low electron density.

2-dimensional systems↗

DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

DOE OSTI2025

The three-point energy correlator in the coplanar limit

Energy correlators are a type of observables that measure how energy is distributed across multiple detectors as a function of the angles between pairs of detectors. In this paper, we study the three-point energy correlator (EEEC) at lepton colliders in the three-particle near-to-plane (coplanar) limit. The leading-power contribution in this limit is governed by the three-jet (trijet) configuration. We introduce a new approach by projecting the EEEC onto the volume of the parallelepiped formed by the unit vectors aligned with three detected final-state particles. Analogous to the back-to-back limit of the two-point energy correlator probing the dijet configuration, the small-volume limit of the EEEC probes the trijet configuration. We derive a transverse momentum dependent (TMD) based factorization theorem that captures the soft and collinear logarithms in the coplanar limit, which enables us to achieve the next-to-next-to-next-to-leading logarithm (N 3 LL) resummation. To our knowledge, this is the first N 3 LL result for a trijet event shape. Additionally, we demonstrate that a similar factorization theorem can be applied to the fully differential EEEC in the three-particle coplanar limit, which provides a clean environment for studying different coplanar trijet shapes.

Effective Field Theories of QCD↗

DOE OSTI2024

Self-calibration strategies for reducing systematic slope measurement errors of autocollimators in deflectometric profilometry

Deflectometric profilometers are used to precisely measure the form of beam shaping optics of synchrotrons and X-ray free-electron lasers. They often utilize autocollimators which measure slope by evaluating the displacement of a reticle image on a detector. Based on our privileged access to the raw image data of an autocollimator, novel strategies to reduce the systematic measurement errors by using a set of overlapping images of the reticle obtained at different positions on the detector are discussed. It is demonstrated that imaging properties such as, for example, geometrical distortions and vignetting, can be extracted from this redundant set of images without recourse to external calibration facilities. This approach is based on the fact that the properties of the reticle itself do not change – all changes in the reticle image are due to the imaging process. Firstly, by combining interpolation and correlation, it is possible to determine the shift of a reticle image relative to a reference image with minimal error propagation. Secondly, the intensity of the reticle image is analysed as a function of its position on the CCD and a vignetting correction is calculated. Thirdly, the size of the reticle image is analysed as a function of its position and an imaging distortion correction is derived. It is demonstrated that, for different measurement ranges and aperture diameters of the autocollimator, reductions in the systematic errors of up to a factor of four to five can be achieved without recourse to external measurements.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Metabolic Redox Coupling Controls Methane Production in Permafrost‐Affected Peatlands Through Organic Matter Quality‐Dependent Energy Allocation

ABSTRACT Permafrost thaw represents one of Earth's largest climate feedback risks, potentially releasing vast carbon (C) stores as greenhouse gases (GHG). However, our ability to predict emissions remains limited by poor understanding of how changing organic matter (OM) composition affects microbial carbon processing. We test a metabolism‐centered redox framework, which views microbial processes as coupled oxidative‐reductive reactions, to mechanistically explain how organic matter metabolite quality controls greenhouse gas production in permafrost‐affected peatland ecosystems. Rather than relying solely on geochemical redox measurements, our approach examines how microbes balance electron flow through metabolic pathways. Using active layer peat (9–19 cm) from contrasting environments (bog and fen), we employed multi‐omics approaches, including metabolomics, metagenomics, and metatranscriptomics, to link OM chemistry to microbial function. Our results reveal distinct dissolved organic matter metabolite composition, with fen systems enriched in compounds with higher substrate quality (low molecular weight (MW) sugars with high H:C ratios and low aromaticity) and bog systems dominated by compounds with lower substrate quality (high MW phenols with lower H:C ratios and higher aromaticity). In fen samples, these sugar‐like compounds correlated with higher oxidative metabolism and methanogenesis, supported by increased glycolysis gene expression. Initially, electrons from increased oxidative metabolism were balanced through nitrate and sulfate reduction, but as these electron acceptors were depleted, methanogenesis increased to maintain redox balance. Fen samples showed rapid degradation of both high‐ and low‐substrate‐quality compounds, suggesting sufficient energy for efficient C cycling. Conversely, bog samples exhibited more polyphenolic compounds, lower glycolysis activity, and higher stress‐related gene expression, suggesting energy was diverted towards cell maintenance under acidic conditions rather than C processing. This approach suggests that predicting greenhouse gas emissions requires an understanding of how organic matter quality shapes microbial energy allocation strategies, providing a mechanistic framework for improving emission predictions from permafrost‐affected peatlands and similar ecosystems.

Biodiversity & Conservation↗

DOE OSTI2020

Sensitivity Studies of Premixed Flame Response to Transverse, High Frequency Disturbances

This study analyzes the response of premixed flames to high frequency transverse acoustic disturbances. A significant body of experiment, modeling and simulation efforts has been performed to understand the response of the spatially integrated heat release of premixed flames to low-frequency disturbances. These studies quantified this sensitivity through the Flame Transfer Function (FTF), as the global heat release is the key quantity influencing thermoacoustic stability for acoustically compact flames. Recent efforts have expanded these models to consider high-frequency acoustic disturbances, and the sensitivity of premixed flames through the Rayleigh integral (RI). The full Rayleigh Integral, as opposed to simply the spatially integrated heat release, must be considered in these cases where the flame is not acoustically compact. Prior modeling studies considered several effects - stability of axisymmetric flames centered in a cylindrical combustion chamber and the effect of flame offset from the combustor axis and the effect of mean flame asymmetries in the flame centered in the combustor. These studies showed that the flame's stability was sensitive to all these parameters, as well as to the transverse mode in question. In this paper, we quantify these different parametric sensitivities of premixed flame stability in the presence of high-frequency transverse acoustic modes - specifically sensitivity to flame length and nozzle location within the acoustic mode shape for 2 different transverse modes. The sensitivity of the RI to these control parameters provides an improved understanding of which flame parameters are critical to improve the overall thermoacoustic stability of the system.

Acharya, Vishal S.↗

DOE OSTI2020

Machine learning for nonintrusive model order reduction of the parametric inviscid transonic flow past an airfoil

Fluid flow in the transonic regime finds relevance in aerospace engineering, particularly in the design of commercial air transportation vehicles. Computational fluid dynamics models of transonic flow for aerospace applications are computationally expensive to solve because of the high degrees of freedom as well as the coupled nature of the conservation laws. While these issues pose a bottleneck for the use of such models in aerospace design, computational costs can be significantly minimized by constructing special, structure-preserving surrogate models called reduced-order models. In this work, we propose a machine learning method to construct reduced-order models via deep neural networks and we demonstrate its ability to preserve accuracy with a significantly lower computational cost. In addition, our machine learning methodology is physics-informed and constrained through the utilization of an interpretable encoding by way of proper orthogonal decomposition. Application to the inviscid transonic flow past the RAE2822 airfoil under varying freestream Mach numbers and angles of attack, as well as airfoil shape parameters with a deforming mesh, shows that the proposed approach adapts to high-dimensional parameter variation well. Notably, the proposed framework precludes the knowledge of numerical operators utilized in the data generation phase, thereby demonstrating its potential utility in the fast exploration of design space for diverse engineering applications. Comparison against a projection-based nonintrusive model order reduction method demonstrates that the proposed approach produces comparable accuracy and yet is orders of magnitude computationally cheap to evaluate, despite being agnostic to the physics of the problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Fast Butterfly-Compressed Hadamard–Babich Integrator for High-Frequency Helmholtz Equations in Inhomogeneous Media with Arbitrary Sources

Here we present a butterfly-compressed representation of the Hadamard-Babich (HB) ansatz for the Green's function of the high-frequency Helmholtz equation in smooth inhomogeneous media. For a computational domain discretized with Nv discretization cells, the proposed algorithm first solves and tabulates the phase and HB coefficients via eikonal and transport equations with observation points and point sources located at the Chebyshev nodes using a set of much coarser computation grids, and then butterfly compresses the resulting HB interactions from all Nv cell centers to each other. The overall CPU time and memory requirement scale as O(Nv log2 Nv) for any bounded two-dimensional (2D) domains with arbitrary excitation sources. A direct extension of this scheme to bounded 3D domains yields an O(Nv4/3) CPU complexity, which can be further reduced to quasi-linear complexities with proposed remedies. The scheme can also efficiently handle scattering problems involving inclusions in inhomogeneous media. Although the current construction of our HB integrator does not accommodate caustics, the resulting HB integrator itself can be applied to certain sources, such as concave-shaped sources, to produce caustic effects. Compared to finite-difference frequency domain methods, the proposed HB integrator is free of numerical dispersion and requires fewer discretization points per wavelength. As a result, it can solve wave propagation problems well beyond the capability of existing solvers. Remarkably, the proposed scheme can accurately model wave propagation in 2D domains with 640 wavelengths per direction and in 3D domains with 54 wavelengths per direction on a state-of-the-art supercomputer at Lawrence Berkeley National Laboratory.

Hadamard--Babich ansatz↗

DOE OSTI2020

Time-resolved in situ visualization of the structural response of zeolites during catalysis

Zeolites form an important class of materials that are widely used in various fields. Their unique properties come primarily from their size, shape, and connectivity of their sub-nm crystalline voids, the Si/Al ratio of the anionic framework, and the number and the nature of charge-balancing cations.The inhomogeneous distribution of the Si/Al ratio within crystals, and hence the cation density, the presence of organic residues, and their inhomogeneous distribution within crystals sensitively affect their performances because these factors influencethe intra-crystalline diffusion rates of the reactants and products. However, the information regarding these inhomogeneities in zeolite crystals cannot be obtained by conventional analytical tools because they only provide information regardingthe average properties of zeolite crystals. Here we employin situX-ray free electron laser-based time-resolved coherent X-ray diffraction imaging to investigate the internal deformations originating from the inhomogeneous Cu ion distributions in Cu-exchanged ZSM-5 zeolite crystals during the deoxygenation of nitrogen oxides with propene. We show that the interactions between the reactants and the active sites lead to an unusual strain distribution,which is confirmed by density functional theory simulations. Moreover, this elucidates the factors responsible for distinct active regions,which are essential to maximizing some applications with enhanced performances. These observations provide insights into the role of structural inhomogeneity in zeolites during a catalytic processand will assist the future design of zeolites for their applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy's Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

closed-loop geothermal systems↗