SEARCH · Engineering Papers

Results for “heat transfer coefficient”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Measurement of the Natural Convection Heat Transfer in a Magnesium Oxide Electrolytic Cell Concept

The rate of heat transfer by natural convection between the wall and electrolyte of an electrolytic cell that produces magnesium (Mg) from magnesium oxide (MgO) at temperatures near 1000 °C in a molten fluoride salt electrolyte is presented. An experimental model of the cell was developed that enabled measurements of the heat transfer in the absence of electrolysis and at temperatures less than 100 °C over ranges of Rayleigh numbers from 1 × 10 -7 to 7 × 10 -8 and Prandtl numbers from 2 to 6200, ranges that include those anticipated in the operation of the MgO electrolytic cell. The model avoids the substantial experimental challenges associated with the high-temperature, corrosive molten salt to enable a conservative estimate of the heat transfer at a lower cost and greater accuracy than would otherwise be possible. The results are correlated by the expression Nu = 0.412Ra 0.23 Pr 0.15 with Nusselt numbers spanning 30–80. The application of the correlation shows that the heat transfer between the cell wall and the molten fluoride electrolyte at ≈1000 °C is characterized by convection coefficients between 100 and 600 W/m 2 -K and is fast enough to enable heat fluxes up to 10 W/cm 2 without compromising the structural integrity of the steel cell wall.

Engineering↗

DOE OSTI2020

Evaporation due to infrared heating and natural convection

The results obtained from a drying experiment under infrared heating are compared to those obtained from correlations for the mass transfer coefficients, which are commonly used to obtain evaporation rates for drying applications. It was shown that the mass transfer coefficients evaluated by natural convection correlations are approximately ten times lower than those measured under infrared (IR) exposure. In this work, an attempt to explain this discrepancy, which has been also observed in other studies and explained on the basis of thermal radiation effects, optical properties of NMP are presented and the implications of the IR and near infrared (NIR) spectrum on the thermal radiation effects and ensuing evaporation rates are discussed. The results show that expected raise in the film temperature, which would be due to the thermal radiation effects, cannot be used to describe the measured evaporation rates for this solvent. Therefore, this study provides important evaporation mass flux data, which cannot be predicted using traditional correlations, for the design of IR assisted drying systems used in membrane fabrication.

42 ENGINEERING↗

DOE OSTI2024

Compatibility of LaFe13−x−yMnxSiyH1.6 and Eutectic Liquid GaInSn Alloy

The heat transfer rate of magnetocaloric regenerators is a topic of extensive research and the cyclability of these regenerators is critical to the operation of systems with a high coefficient of performance (e.g., potentially >22, significantly higher than typical vapor compression cooling technologies). To enable a high operating frequency that will result in a high specific cooling power, the heat transfer fluid should have high thermal conductivity and lower specific heat, i.e., higher thermal diffusivity. Eutectic metal alloys possess these qualities, such as gallium–indium–tin (Galinstan), whose thermal diffusivity has been found to be approximately an order of magnitude higher than water. For this study, the effects of eutectic liquid Galinstan exposure on the phase stability of LaFe13−x−yMnxSiyH1.6 magnetocaloric powders in an active magnetic regenerator device were investigated. The powders were characterized before and after exposure to Galinstan using X-ray diffraction, in which the phases were determined using the Rietveld refinement technique and X-ray fluorescence. It was found that after Galinstan exposure, hydrogen containing phases were present in the powder, suggesting that the hydrogen was lost from the magnetocaloric phase. The magnetocaloric phase degradation indicates that the powder was incompatible with the Galinstan metal in an environment with moisture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Multi-objective optimization of low-GWP mixture composition and heat exchanger circuitry configuration for improved system performance and reduced refrigerant flammability

The recently introduced hydrofluoroolefin (HFO) refrigerants, including R1234yf and R1234ze(E), have significantly lower global warming potentials (GWPs) than traditional hydrofluorocarbon (HFC) refrigerants like R410A. However, prior tests show that direct drop-in of pure R1234yf or R1234ze(E) into equipment designed for R410A results in a decrease in heat exchanger capacity and the system coefficient of performance. The primary reason is the lower in-tube heat transfer performance of R1234yf and R1234ze(E) compared with that of R410A. To address this issue, previous studies have mixed the mildly flammable HFC R32 with HFOs to improve system performance, with HFC R125 also added to suppress flammability. Previous studies selected compositions based on simple cycle analyses and did not consider modifications of the heat exchanger circuitry configuration to adapt to the new refrigerants. This study presents a novel multi-objective optimization approach to design a refrigerant composition that maximizes energy efficiency within flammability and GWP limits. The approach in this work simultaneously optimizes mixture composition and heat exchanger circuitry configuration. A case study on a rooftop unit indicates that, compared with mixture-only optimization, simultaneous optimization of mixture and heat exchanger circuitry yields a 5.9% improvement in cycle efficiency and a 48.6% reduction in refrigerant flammability with a GWP of 268. Circuitry optimization using refrigerants with different temperature glides shows that the larger the temperature glide is, the larger EER improvement is obtained. The results show that zeotropic blends with a large temperature glide are more sensitive to the refrigerant circuitry than pure refrigerants and may suffer significant performance degradation with subpar heat exchanger circuitry design. The proposed optimization approach is generally applicable to mixtures with any number of components. Using this approach to design a HVAC system can yield higher system efficiency within flammability and GWP constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Optimization of Refrigerant Compositions for Low-GWP Refrigerant Mixtures Using Segment-by-segment Heat Exchanger and Detailed System Models

Li, Zhenning↗

DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗

DOE OSTI2023

A thermodynamic model of integrated liquid-to-liquid thermoelectric heat pump systems

Thermoelectric (TE) heat pumps (TEHPs) are advantageous for heating and cooling in various applications because of their modularity and simple design. A TEHP system includes the TE modules with p- and n-type materials bonded to substrates, plus heat exchangers, thermal interfaces to the heat exchangers, and heat transfer fluids. Although modeling an individual TE module has been extensively studied, limited studies have reported performance at the larger system-level. Furthermore, no prior study has addressed the impact of temperature-dependent TE material properties (e.g., electric resistivity, thermal conductivity, and Seebeck coefficient) on overall heat-pump-system-level performance. Here, this work presents a mathematical model for TEHP system performance based on Goldsmid's approach for TE material performance, “effective” TE material properties, Gnielinski's correlation for convective heat transfer, and thermal balance theory for a heat exchange network. This combined approach provides an accurate model of the liquid-to-liquid TEHP system. Three different approaches—one empirical, one based on the manufacturer's specifications, and one drawn from the literature—were then used to determine values for TE material properties. The first two methods treated properties as constants, while the last approach treated properties as surface-temperature-based functions. Finally, experimental TEHP data was used to validate the models, all with relative absolute deviations of approximately 10% when predicting heating capacity and 10%–25% when forecasting cooling capacity up to a 30 K surface temperature lift. The results demonstrated that, at the TEHP system level, the TE material properties could be treated as constants, avoiding solver iterations and reducing the performance uncertainty by up to 95%.

42 ENGINEERING↗

DOE OSTI2020

Developing Equivalent Surface Heat Transfer Input Parameters for Updated Standard 140/BESTEST Thermal Fabric Test Cases

ANSI/ASHRAE Standard 140, Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs, is widely cited for evaluating software for use with performance-path energy efficiency analysis in conjunction with well-known energy-efficiency standards, including: ASHRAE Standard 90.1, the International Energy Conservation Code, and other international standards. This paper summarizes recent work to develop improved equivalent surface heat transfer input parameters for the building thermal fabric comparative test cases that comprise a portion of Standard 140. The software-to-software comparative test method allows comparison of simulation results from a wide variety of building energy simulation programs to each other. In this type of test method, it is of great importance to specify inputs that are as well-equivalenced as possible to accommodate a wide variety of input schemes for the tested programs, and thereby reduce the potential for disagreements among them in the tests. Results include new alternative constant surface coefficients that are appropriate to the test cases and demonstrate their ability to approximate the more detailed surface heat transfer models from which they are derived.

BESTEST↗

DOE OSTI2022

Photothermal Atomic Force Microscopy Coupled with Infrared Spectroscopy (AFM-IR) Analysis of High Extinction Coefficient Materials: A Case Study with Silica and Silicate Glasses

Photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) brings significant value as a spatially resolved surface analysis technique for disordered oxide materials such as glasses, but additional development and fundamental understanding of governing principles is needed to interpret AFM-IR spectra, since the existing theory described for organic materials does not work for materials with high extinction coefficients for infrared (IR) absorption. This paper describes theoretical calculation of a transient temperature profile inside the IR-absorbing material considering IR refraction at the interface as well as IR adsorption and heat transfer inside the sample. Furthermore, this calculation explains the differences in peak positions and amplitudes of AFM-IR spectra from those of specular reflectance and extinction coefficient spectra. It also addresses the information depth of the AFM-IR characterization of bulk materials. AFM-IR applied to silica and silicate glass surfaces has demonstrated novel capability of characterizing subsurface structural changes and surface heterogeneity due to mechanical stresses from physical contacts, as well as chemical alterations manifested in surface layers through aqueous corrosion.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Dynamic modeling of near isothermal compressor for Transcritical carbon dioxide cycle to support efficiency improvement

Compressors are the primary components of energy consumption in vapor compression systems, drawing considerable research effort to reduce carbon emissions and improve energy efficiency. This work presents an isothermal compressor that integrates the compression chamber with the gas cooler to achieve near-isothermal compression, resulting in up to 30% energy reduction compared to traditional isentropic compression, according to the experiment results. A detailed dynamic model for the isothermal compression cycle of carbon dioxide (CO 2 ) is first presented, coupled with a liquid-piston model for chamber design. The model is experimentally validated on a 3.5 kW prototype, yielding a 6.2% relative error in chamber pressure and a 0.5 K deviation in chamber temperature. The liquid piston exhibits reduced frictional heating and geometric adaptability relative to the mechanical piston. Leveraging these attributes, the study further evaluates approaches that enhance the heat transfer to achieve isothermal compression by increasing the effective heat transfer area during compression. As a result, simulations demonstrated up to a 5% improvement in the coefficient of performance over the baseline case through chamber design, providing guidelines for design of isothermal compression cycle.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Alpert multi-wavelets for functional inverse problems: direct optimization and deep learning

Computational engineering models often contain unknown entities (e.g. parameters, initial and boundary conditions) that require estimation from other measured observable data. Estimating such unknown entities is challenging when they involve spatio-temporal fields because such functional variables often require an infinite-dimensional representation. Here, we address this problem by transforming an unknown functional field using Alpert wavelet bases and truncating the resulting spectrum. Hence the problem reduces to the estimation of few coefficients that can be performed using common optimization methods. We apply this method on a one-dimensional heat transfer problem where we estimate the heat source field varying in both time and space. The observable data is comprised of temperature measured at several thermocouples in the domain. This latter is composed of either copper or stainless steel. The optimization using our method based on wavelets is able to estimate the heat source with an error between 5% and 7%. We analyze the effect of the domain material and number of thermocouples as well as the sensitivity to the initial guess of the heat source. Finally, we estimate the unknown heat source using a different approach based on deep learning techniques where we consider the input and output of a multi-layer perceptron in wavelet form. We find that this deep learning approach is more accurate than the optimization approach with errors below 4%.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Impact of uranium oxide (UO 2 ) fuel with molybdenum (Mo) inserts on pressurized water reactor performance and safety

This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

A molecular dynamics study of collisional heat transfer to nanoclusters in the gas phase

Nanoclusters in the gas phase grow by condensation and coagulation, which heat up nanoclusters, creating thermal non-equilibrium with the surrounding gas. The extent of non-equilibrium for nanoclusters is dependent upon the rate of heating relative to the rate of cooling brought about by nanocluster-gas molecule collisions, hence the latter process plays an important role in nanocluster energy evolution during growth. We apply molecular dynamics simulations to investigate heat transfer between metal nanoclusters and gas molecules relevant to high temperature aerosol systems. In analysis, we first define a thermal reemission coefficient (TRC), which is a correction factor for the portion of energy possessed by reemitted gas molecules after collision. Here, the TRC differs from the thermal accommodation coefficient (TAC), which is a correction factor for the maximum possible energy transfer rate from nanoclusters to gas molecules based on the diffusive collision model. We find the TRC is insensitive to nanocluster size but is affected by the nanocluster to gas temperature ratio and the nanocluster to gas atomic mass ratio. Comparison to the traditionally-defined TAC shows that the TAC is insensitive to temperature, a phenomenon that has been widely indicated in literature but not quantitatively explained. We show how the temperature insensitivity of the TAC arises because of the manner in which TRCs vary with both nanocluster and gas temperature, i.e. the manner in which the reemitted gas molecule energy from a surface varies with the temperature of that surface relative to that of the gas. In addition, while calculations ultimately reveal temperature insensitive TACs and justify their continued use in modeling gas-phase heat transfer, we argue that the TRC is more appropriate to calculate a priori in analyzing simulations; it is physically more realistic to first introduce a correction factor solely for the energy of reemitted gas molecules after a collision (as TRCs) instead of a correction factor applying to both the reemission energy and the initial energy (as TACs), as the latter is exactly calculable and independent of nanocluster properties. In line with prior work, for monoatomic gases we find the TAC decreases as the nanocluster to gas atomic mass increases, approaching 0.08 for gold nanoclusters in helium.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2020

Modeling of gap conductance for LWR fuel rods applied in the BISON code

A new gap conductance model is proposed in this study as a combination of Toptan’s model and the Ross-Stoute model. A variance-based sensitivity analysis is performed to understand how simulation results depend on all input parameters of the proposed model. Additionally, new modeling options (e.g. fill gas thermal conductivity, temperature jump distance, thermal accommodation coefficient, etc.) are added into the nuclear fuel performance code, BISON. The need for further investigation of the gap heat transfer between fuel and cladding in BISON motivated this study to evaluate its impact on the code’s predictions. New gap conductance modeling is proposed. A series of integral-effects validation tests is performed: (1) to demonstrate the impact of the proposed model on the code’s fuel temperature predictions at the beginning of life and through the reactor’s life; (2) to ensure that the proposed model is capable of accurately modeling gap heat transfer characteristics in real-world problems; and (3) to investigate the impact of the estimation of fission gas release on the fuel temperature predictions with the proposed model. Overall, the results indicate that the proposed gap conductance model improves BISON’s predictions.

42 ENGINEERING↗

DOE OSTI2024

Thermal diode and thermal switch for bi-directional heat transfer in building envelopes

The present disclosure relates to a directional heat transfer using thermal control devices, including a dual phase change thermal diode and an active contact-based thermal switch. The thermal diode includes a positive temperature coefficient switching material and a negative temperature coefficient switching material arranged in series. The thermal switch includes two thermally conducting surfaces which may be moved to contact (i.e., having a distance between them of substantially zero) creating minimal thermal contact resistance. Both thermal control devices may be used to control heat flow into and/or out of a building.

Kishore, Ravi Anant↗

DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance: Preprint

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer, rather than also comparing control systems responses. HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study shows that HEATNETS, while not intended to replace other modeling tools, can be a unique modeling tool for the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗