SEARCH · Search NASA

Results for “gradient based optimization”

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 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗

DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

DOE OSTI2023

Quantum Chemistry-Driven Machine Learning Approach for the Prediction of the Surface Tension and Speed of Sound in Ionic Liquids

Ionic liquids (ILs) have unique solvent properties and have thus garnered significant interest. However, exhaustive experimental determination of the physicochemical properties of ILs is unrealistic due to the large structural diversity of anions and cations, their high cost, the requirements of elevated temperature and pressure, and the time required. To circumvent these experimental costs, computational approaches to accurately calculate these properties have emerged. Here in the present study, we present a demonstration of two machine learning (ML) models for the prediction of two critical IL physical properties, the surface tension and the speed of sound, across a wide range of temperatures and pressures. The models make use of molecular descriptors derived from the COSMO-RS, a quantum chemical-based model. The ML models show excellent agreement with experimental observations, with an R2 value of 0.96–0.99 and RMSE of 1.71 mN/m and 16.12 m/s for the surface tension and speed of sound, respectively. This work paves the way for the development of COSMO-RS-informed ML models for the prediction of IL properties which can help to further optimize and accelerate technology development for ILs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

A self-seeded SRS system for the generation of 1.54-micron eye-safe radiation

A light source is described for the efficient generation of 1.54-micron for eye-safe aerosol lidar operation. The system is based upon an Nd:YAG laser at 1.06 micron which is then Raman-shifted in methane to produce light at the first Stokes wavelength of 1.54 micron. First Stokes light generated in the backward direction was retroreflected back into the Raman cell for amplification in the tail of the 10 ns pump beam. The energy conversion efficiency and the spatial beam quality of the amplified first Stokes were found to be adversely affected when operated at higher repetition rates due to a thermal gradient produced in the generation region. A Stokes energy of 25 mJ was obtained for a pumping energy of 140 mJ at a repetition rate of 10 Hz. The beam divergence of the amplified Stokes radiation was measured to be less than 1 mrad. The optimized results demonstrate the applicability of this radiation for eye-safe lidar measurements.

Chu, Z.↗

DOE OSTI2026

Review of SiC material development for nuclear fusion applications: Cross-cutting research and emerging opportunities

The SiC-based materials, particularly SiC-fiber-reinforced SiC matrix (SiC/SiC) composites, show strong potential for structural and functional applications in future fusion power plants because they can operate at high temperatures with a range of coolants and breeders, thereby enabling higher energy conversion efficiency. Here, this paper presents recent advancements in the development of SiC-based materials, focusing on processing techniques and material performance and resistance under fusion-relevant environments. The processing activities have emphasized near-net-shape fabrication and the joining of SiC subcomponents, with processing methods and material compositions informed by previous irradiation experiments on various grades of SiC. Research on irradiation effects has remained focused on degradation mechanisms and the microstructural optimization of SiC/SiC composites irradiated to high neutron damage levels. Analysis of irradiation defects in SiC has advanced via the application of cutting-edge characterization methods, among which Raman spectroscopy is becoming a common tool to assess atomic-scale chemical disorder. Fusion–fission crosscutting irradiation research has explored combined effects in SiC/SiC composites with application-relevant geometries, including bowing of SiC/SiC composite channels under neutron flux gradients, stress evolution in SiC/SiC composite tubes under through-thickness temperature gradients, and irradiation-enhanced corrosion in SiC. Finally, research opportunities for component testing and assessment under fusion-relevant conditions, in support of emerging concepts from the private fusion sector, are discussed.

Advanced manufacturing↗

DOE OSTI2024

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2017

Development of Multi-Sensor Global Cloud and Radiance Composites for DSCOVR EPIC Imager with Subpixel Definition

The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). EPIC delivers adequate spatial resolution imagery but only in shortwave bands (317-780 nm), while NISTAR measures the top-of-atmosphere (TOA) whole-disk radiance in shortwave and longwave broadband windows. Accurate calculation of albedo and outgoing longwave flux requires a high-resolution scene identification such as the radiance observations and cloud properties retrievals from low earth orbit (LEO, including NASA Terra and Aqua MODIS, Suomi-NPP VIIRS, and NOAA AVHRR) and geosynchronous (GEO, including GOES east and west, METEOSAT, INSAT-3D, MTSAT-2, and Himawari-8) satellite imagers. The cloud properties are derived using the Clouds and the Earth's Radiant Energy System (CERES) mission Cloud Subsystem group algorithms. These properties have to be co-located with EPIC pixels to provide the scene identification and to select anisotropic directional models (ADMs), which are then used to adjust the NISTAR-measured radiance and subsequently obtain the global daytime shortwave and longwave fluxes. This work presents an algorithm for optimal merging of selected radiance and cloud property parameters derived from multiple satellite imagers to obtain seamless global hourly composites at 5-km resolution. Selection of satellite data for each 5-km pixel is based on an aggregated rating that incorporates five parameters: nominal satellite resolution, pixel time relative to the EPIC time, viewing zenith angle, distance from day/night terminator, and probability of sun glint. To provide a smoother transition in the merged output, in regions where candidate pixel data from two satellite sources have comparable aggregated rating, the selection decision is defined by the cumulative function of the normal distribution so that abrupt changes in the visual appearance of the composite data are avoided. Higher spatial accuracy in the composite product is achieved by using the inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling.

Khlopenkov, Konstantin V.↗

DOE OSTI2023

Sucrose transport inside the phloem: Bridging hydrodynamics and geometric characteristics

In plants, the delivery of the products of photosynthesis is achieved through a hydraulic system labeled as phloem. This semi-permeable plant tissue consists of living cells that contract and expand in response to fluid pressure and flow velocity fluctuations. The Münch pressure flow theory, which is based on osmosis providing the necessary pressure gradient to drive the mass flow of carbohydrates, is currently the most accepted model for such sucrose transport. When this hypothesis is combined with the conservation of fluid mass and momentum as well as sucrose mass, many simplifications must be invoked to mathematically close the problem and to resolve the flow. Here this study revisits such osmotically driven flows by developing a new two-dimensional numerical model in cylindrical coordinates for an elastic membrane and a concentration-dependent viscosity. It is demonstrated that the interaction between the hydrodynamic and externally supplied geometrical characteristic of the phloem has a significant effect on the front speed of sucrose transport. These results offer a novel perspective about the evolutionary adaptation of plant hydraulic traits to optimize phloem soluble compounds transport efficiency.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1995

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗

NASA NTRS2006

Effects of Vegetation and of Heat and Vapor Fluxes from Soil on Snowpack Evolution and Radiobrightness

The radiobrightness of a snowpack is strongly linked to the snow moisture content profile, to the point that the only operational inversion algorithms require dry snow. Forward dynamic models do not include the effects of freezing and thawing of the soil beneath the snowpack and the effect of vegetation within the snow or above the snow. To get a more realistic description of the evolution of the snowpack, we reported an addition to the Snow-Soil-Vegetation-Atmosphere- Transfer (SSVAT) model, wherein we coupled soil processes of the Land Surface Process (LSP) model with the snow model SNTHERM. In the near future we will be adding a radiobrightness prediction based on the modeled moisture, temperature and snow grain size profiles. The initial investigations with this SSVAT for a late winter and early spring snow pack indicate that soil processes warm the snowpack and the soil. Vapor diffusion needs to be considered whenever the ground is thawed. In the early spring, heat flow from the ground into a snow and a strong temperature gradient across the snow lead to thermal convection. The buried vegetation can be ignored for a late winter snow pack. The warmer surface snow temperature will affect radiobrightness since it is most sensitive to snow surface characteristics. Comparison to data shows that SSVAT provides a more realistic representation of the temperature and moisture profiles in the snowpack and its underlying soil than SNTHERM. The radiobrightness module will be optimized for the prediction of brightness when the snow is moist. The liquid water content of snow causes considerable absorption compared to dry snow, and so longer wavelengths are likely to be most revealing as to the state of a moist snowpack. For volumetric moisture contents below about 7% (the pendular regime), the water forms rings around the contact points between snow grains. Electrostatic modeling of these pendular rings shows that the absorption of these rings is significantly higher than a sphere of the same volume. The first implementation of the radiobrightness module will therefore be a simple radiative transfer model without scattering.

soil↗

DOE OSTI2024

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

On the Convergence of the Inexact Running Krasnosel'skii-Mann Method

This letter leverages a framework based on averaged operators to tackle the problem of tracking fixed points associated with maps that evolve over time. In particular, this letter considers the Krasnosel'skii-Mann (KM) method in a settings where: 1) the underlying map may change at each step of the algorithm, thus leading to a 'running' implementation of the KM method, and 2) an imperfect information of the map may be available. An imperfect knowledge of the maps can capture cases where processors feature a finite precision or quantization errors, or the case where (part of) the map is obtained from measurements. The analytical results are applicable to inexact running algorithms for solving optimization problems, whenever the algorithmic steps can be written in the form of (a composition of) averaged operators; examples are provided for inexact running gradient methods and the forward-backward splitting method. Convergence of the average fixed-point residual is investigated for the non-expansive case; linear convergence to a unique fixed-point trajectory is showen in the case of inexact running algorithms emerging from contractive operators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Enabling On-Demand Aerospace Component Manufacturing: Topology Optimization of GE Engine Bracket and Fabrication Using Metal FFF

Additive Manufacturing (AM) offers advantages over conventional manufacturing processes, particularly by reducing the number of parts produced through multistage combined technologies, but these often result in low manufacturing yields or require post-processing. AM facilitates the production of complex geometries with fine features, overhangs, and lattice structures. For instance, Laser Powder Bed Fusion (LPBF) technology enables the fabrication of intricate parts that can be easily post-processed by removing residual powder. Laser powder bed AM technologies are widely discussed in the literature due to their design freedom in creating complex geometries, with and without the need for support generation. However, rapid solidification due to a thermal gradient in the build direction, which leads to the formation of columnar grains and warpage, is one of the challenges. To address this challenge, we propose layer-by-layer metal FFF technology, followed by the debinding and sintering process, as an alternative to powder- and laser-based approaches. Furthermore, design for additive manufacturing (DfAM) principles are discussed to minimize the need for support generation and enable easy post-processing, thereby improving surface finish to meet high tolerances in fabricating components for aerospace and healthcare applications.

Singh, Abhishek [University of Michigan, Ann Arbor↗

DOE OSTI2021

Solid Phase Gradient Alloying Method via ShAPE

Alloying has been used to confer desirable properties to various metallic systems since the bronze age. For example, Ni, Mn, and Cr are alloyed into an Fe matrix to make steel, which has improved strength and corrosion resistance (as displayed in Fig. 1). In addition to the type of alloying elements, the amount of alloying elements added into a matrix is also very critical for a material’s mechanical property, physical property, machinability, and cost. Conventionally, new alloys are discovered by combining alloy precursors and elemental constituents into a mixture and melting them together. The solidified product or ingot represents a usually non-homogeneous combination of the chemistry with a microstructure dictated by the physics of solidification. It is time and energy intensive to optimize the alloying amounts via the casting process and hence, economically not conducive to rapid alloy discovery. In addition, the cast microstructure seen in the as-cast ingot is often not the ideal, nor even desired structure for optimized performance. Cast materials often need further processing such as homogenization, annealing, and deformation work put in through rolling, forging, extruding, etc. to have favorable microstructures and properties. A faster method to discover new alloy combinations and evaluate them in their worked or processed form is needed. Recently, Laser Engineered Net Shaping was applied to fabricate gradient compositional material for alloy designing. However, the mechanical properties of this gradient alloy were poor due to the existence of impurities, oxidation and cracks that are inherent in the melt-solidification process. Because of these limitations, we propose to use a solid-phase process (ShAPE) to create bulk materials (extrudates) that vary in composition from one end of the extruded solid to the other. Various manufacturing methods for alloy design including conventional casting method, laser based combinatorial method and current solid phase gradient alloying method are displayed and compared in Fig. 2. ShAPE machine and a schematic of ShAPE process are displayed in Fig. 3. Starting material is processed by a rotating and plunging die to form an extrudate. We are proposing through this methodology to invent a new, bulk-scale combinatorial technique. With this novel technique, alloying element can be dissolved into the matrix with a continued gradient without bulk melting. Formation of inter-metallics and defects originating from melt processing can be avoided. In addition, the ShAPE process creates the mixed alloy chemistry, and meanwhile subjects the material to severe plastic strain, which induces the favorable “worked” microstructure. Products from the ShAPE process can be tested in hardness directly providing a path to rapid evaluation of mechanical properties. The ability to create graded structures using friction stir processing (FSP) has been demonstrated, however using the ShAPE process we believe will lead to much faster and higher fidelity results. When combined with high throughput screening methods of physical, mechanical and microstructural property characterization, efficiency and accuracy of alloy design can be significantly improved.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Objective and algorithm considerations when optimizing the number and placement of turbines in a wind power plant

Abstract. Optimizing turbine layout is a challenging problem that has been extensively researched in the literature. However, optimizing the number of turbines within a given boundary has not been studied as extensively and is a difficult problem because it introduces discrete design variables and a discontinuous design space. An essential step in performing wind power plant layout optimization is to define the objective function, or value, that is used to express what is valuable to a wind power plant developer, such as annual energy production, cost of energy, or profit. In this paper, we demonstrate the importance of selecting the appropriate objective function when optimizing a wind power plant in a land-constrained site. We optimized several different wind power plants with different wind resources and boundary sizes. Results show that the optimal number of turbines varies drastically depending on the objective function. For a simple, one-dimensional, land-based scenario, we found that a wind power plant optimized for minimal cost of energy produced just 72 % of the profit compared to the wind power plant optimized for maximum profit, which corresponded to a loss of about USD 2 million each year. This paper also compares the performance of several different optimization algorithms, including a novel repeated-sweep algorithm that we developed. We found that the performance of each algorithm depended on the number of design variables in the problem as well as the objective function.

17 WIND ENERGY↗

DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗