Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global 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.

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

At least 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Cutting edge research and new technologies in heat and mass transfer processes of refrigeration and air conditioning systems

This special edition of the Science and Technology for the Built Environment (STBE) presents eight selected articles originally presented at the 17th International Refrigeration and Air Conditioning Conference held at Purdue University, West Lafayette, IN, USA, during July 9–12, 2018. All eight articles were significantly expanded with additional technical content and received full peer review before being finally accepted for publication in this special issue of the STBE. The eight articles in this special issue cover a wide range of topics from condensation and evaporation of low global warming potential (LGWP) refrigerants to phase change material storage for personal cooling systems, from frosting on hydrophilic and hydrophobic coated surfaces to geometry optimization of heat exchanger fins. Furthermore, these articles describe cutting-edge research and feasibility studies of new technologies that were recently investigated for refrigeration and air conditioning applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Integrated operation of steady-state long-pulse H-mode in Experimental Advanced Superconducting Tokamak

Recent Experimental Advanced Superconducting Tokamak (EAST) experiments have successfully demonstrated a long-pulse steady-state scenario with improved plasma performance through integrated operation since the last IAEA FEC in 2016. A discharge with a duration over 100 s using pure radio frequency (RF) power heating and current drive has been obtained with the required characteristics for future long-pulse tokamak reactors such as good energy confinement quality (H 98y2 ~ 1.1) with electron internal transport barrier inside ρ < 0.4, small ELMs (frequency ~100–200 Hz), and good control of impurity and heat exhaust with the tungsten divertor. The optimization of X-point, plasma shape, the outer gap and local gas puffing near the low hybrid wave (LHW) antenna were integrated with global parameters of BT and line-averaged electron density for higher current drive efficiency of LHW and on-axis deposition of electron cyclotron heating in the long-pulse operation. More recently, a high βP RF-only discharge (β P ~ 1.9 and β N ~ 1.5, $\langle$n e $\rangle$/n GW ~ 0.80, f bs ~ 45% at q 95 ~ 6.8) was successfully maintained over 24 s with improved hardware capabilities, demonstrating performance levels needed for the China Fusion Engineering Test Reactor steady-state operation. In this work, a higher energy confinement is observed at higher β P and with favorable toroidal field direction. Towards the next goal (≥400 s long-pulse H-mode operations with ~50% bootstrap current fraction) on EAST, an integrated control of the current density profile, pressure profile and radiated divertor will be addressed in the near future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

cMonkey2: Automated, systematic, integrated detection of co-regulated gene modules for any organism

The cMonkey integrated biclustering algorithm identifies conditionally co-regulated modules of genes (biclusters). cMonkey integrates various orthogonal pieces of information which support evidence of gene co-regulation, and optimizes biclusters to be supported simultaneously by one or more of these prior constraints. The algorithm served as the cornerstone for constructing the first global, predictive Environmental Gene Regulatory Influence Network (EGRIN) model for a free-living cell, and has now been applied to many more organisms. However, due to its computational inefficiencies, long run-time and complexity of various input data types, cMonkey was not readily usable by the wider community. To address these primary concerns, we have significantly updated the cMonkey algorithm and refactored its implementation, improving its usability and extendibility. These improvements provide a fully functioning and user-friendly platform for building co-regulated gene modules and the tools necessary for their exploration and interpretation. We show, via three separate analyses of data for E. coli, M. tuberculosis and H. sapiens, that the updated algorithm and inclusion of novel scoring functions for new data types (e.g. ChIP-seq and transcription factor over-expression [TFOE]) improve discovery of biologically informative co-regulated modules. The complete cMonkey2 software package, including source code, is available at https://github.com/baliga-lab/cmonkey2.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

WE-G-207-04: Non-Local Total-Variation (NLTV) Combined with Reweighted L1-Norm for Compressed Sensing Based CT Reconstruction

Purpose: Compressed sensing (CS) has been used for CT (4DCT/CBCT) reconstruction with few projections to reduce dose of radiation. Total-variation (TV) in L1-minimization (min.) with local information is the prevalent technique in CS, while it can be prone to noise. To address the problem, this work proposes to apply a new image processing technique, called non-local TV (NLTV), to CS based CT reconstruction, and incorporate reweighted L1-norm into it for more precise reconstruction. Methods: TV minimizes intensity variations by considering two local neighboring voxels, which can be prone to noise, possibly damaging the reconstructed CT image. NLTV, contrarily, utilizes more global information by computing a weight function of current voxel relative to surrounding search area. In fact, it might be challenging to obtain an optimal solution due to difficulty in defining the weight function with appropriate parameters. Introducing reweighted L1-min., designed for approximation to ideal L0-min., can reduce the dependence on defining the weight function, therefore improving accuracy of the solution. This work implemented the NLTV combined with reweighted L1-min. by Split Bregman Iterative method. For evaluation, a noisy digital phantom and a pelvic CT images are employed to compare the quality of images reconstructed by TV, NLTV and reweighted NLTV. Results: In both cases, conventional and reweighted NLTV outperform TV min. in signal-to-noise ratio (SNR) and root-mean squared errors of the reconstructed images. Relative to conventional NLTV, NLTV with reweighted L1-norm was able to slightly improve SNR, while greatly increasing the contrast between tissues due to additional iterative reweighting process. Conclusion: NLTV min. can provide more precise compressed sensing based CT image reconstruction by incorporating the reweighted L1-norm, while maintaining greater robustness to the noise effect than TV min.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Enhancing Least-Squares Finite Element Methods Through a Quantity-of-Interest

Here, we introduce an approach that augments least-squares finite element formulations with user-specified quantities-of-interest. The method incorporates the quantity-of-interest into the least-squares functional and inherits the global approximation properties of the standard formulation as well as increased resolution of the quantity-of-interest. We establish theoretical properties such as optimality and enhanced convergence under a set of general assumptions. Central to the approach is that it offers an element-level estimate of the error in the quantity-of-interest. As a result, we introduce an adaptive approach that yields efficient, adaptively refined approximations. Several numerical experiments for a range of situations are presented to support the theory and highlight the effectiveness of our methodology. Notably, the results show that the new approach is effective at improving the accuracy per total computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Lattice Element Error Solver for Modeling of Accelerators, Storage Ring, Transport Systems, and Insertion Devices

The goals of CRADA 350 were achieved in this Phase I SBIR work. Specifically, a software environment was developed which included a nonlinear solver, OASIS. During the Phase I work, we applied the algorithm to a number of nonlinear-solver benchmark problems, including SLAC accelerator problems. By searching over only two optimization-control parameters, this Phase I work has enabled the development of a systematic procedure to facilitate the finding of a “best-guessed” global minimum solution for various types of non-linear problems with different size and complexity. This set of rules was codified and led to a versatile nonlinear program, IMIGO.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
NASA NTRS1977

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

The described investigation is concerned with he development of a finite element scheme which can be used in a study of the acoustics of aircraft-engine ducts. In the absence of suitable variational principles for acoustic fluctuations within an aircraft fan engine, an acoustic analysis must proceed directly from the differential equations which describe compressible flow. The derived equations cannot be solved algebraically. The numerical technique used for solving them makes use of a linear rectangular element of a type considered by Zienkiewicz (1971). Attention is given to aspects of element derivation, the global matrix assembly, the solution of the matrix equation, questions of acoustic attenuation, and illustrations of the potential of the current model in duct optimization.

Abrahamson, A. L.↗
NASA NTRS1998

GCM Hindcasts for SST Forced Climate Variability over Agriculturally Intensive Regions

The ability to forecast seasonal climate is of great practical interest. One of the most obvious benefits would be agriculture, for which various preparations (planting, machinery, irrigation, manpower) would be enabled. The expectation of being able to make such forecasts far enough in advance (on the order of 9 months) hinges on components of the system with the longest persistence or predictability. The mixed results of El Nino forecasts has raised the hope that tropical Pacific sea surface temperatures (SST) fall into this category. For agriculturally-relevant forecasts to be made, and utilized, requires several conditions. The SST in the regions that affect agricultural areas must be forecast successfully, many months in advance. The climate response to such sea surface temperatures must then be ascertained, either through the use of historical empirical studies or models (e.g., GCMS). For practical applications, the agricultural production must be strongly influenced by climate, and farmers on either the local level or through commercial concerns must be able to adjust to using such forecasts. In a continuing series of papers, we will explore each of these components. This article concerns the question of utilizing SST to forecast the climate in several regions of agricultural production. We optimize the possibility of doing so successfully by using observed SST in a hindcast mode (i.e., a perfect forecast), and we also use the globally observed values (rather than just those from the tropical Pacific, for which predictability has been shown). This then is the ideal situation; in subsequent papers we will explore degrading the results by using only tropical Pacific SSTs, and then using only

Druyan, Leonard M.↗
NASA NTRSDate not supplied

Space Weather Investigation Frontier (SWIFT) Mission Concept: Continuous Distributed Observations of Geo-Effective, Heliospheric Structures from the Vantage Points of Sun-Earth L1 and Sub- L1

Continuous, in-situ, multi-point observations along the Sun-Earth line at and inside the Lagrange point L1 (sub-L1) will enable a better understanding of the three-dimensional structure and temporal evolution of heliospheric structures that drive terrestrial space weather. The proposed SWIFT (Space Weather Investigation Frontier) mission will use a new solar sail propulsion system developed by NASA to enable a suite of science instruments onboard a smallsat to maintain observations along the Sun-Earth line, sub-L1, for extended periods. Three identically instrumented small spacecraft at L1 will fly in concert with the sailcraft at sub-L1 in an optimized tetrahedron constellation, covering distances between 10 to 100’s of Earth radii. This viewing geometry will enable scientists to distinguish between local and global processes driving space weather by revealing the spatial characteristics, temporal evolution, and geo-effectiveness of small- to meso-scale solar wind structures and substructures of macro-scale structures, such as interplanetary coronal mass ejections (ICMEs) and stream interaction regions (SIRs). In addition, real time measurements of earth-bound heliospheric structures from sub-L1 will improve our current forecasting lead-times by up to 35 percent. This paper will provide an overview of the proposed SWIFT science and mission concept.

Spacecraft Propulsion↗
NASA NTRSDate not supplied

Space Weather Investigation Frontier (SWIFT)

Continuous, in-situ, multi-point observations along the Sun-Earth line at and inside the Lagrange point L1 (sub-L1) will enable a better understanding of the three-dimensional structure and temporal evolution of heliospheric structures that drive terrestrial space weather. The proposed SWIFT (Space Weather Investigation Frontier) mission will use a new solar sail propulsion system developed by NASA to enable a suite of science instruments onboard a smallsat to maintain observations along the Sun-Earth line, sub-L1, for extended periods. Three identically instrumented small spacecraft at L1 will fly in concert with the sailcraft at sub-L1in an optimized tetrahedron constellation, covering distances between 10 to 100’s of Earth radii. This viewing geometry will enable scientists to distinguish between local and global processes driving space weather by revealing the spatial characteristics, temporal evolution, and geo-effectiveness of small-to meso-scale solar wind structures and substructures of macro-scale structures, such as interplanetary coronal mass ejections (ICMEs) and stream interaction regions (SIRs). In addition, real time measurements of earth-bound heliospheric structures from sub-L1 will improve our current forecasting lead-times by up to 35 percent. This paper will provide an overview of the proposed SWIFT science and mission concept.

Heliophysics↗
DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗
DOE OSTI2023

Learning Provably Stable Local Volt/Var Controllers for Efficient Network Operation

Here this paper develops a data-driven framework to synthesize local Volt/Var control strategies for distributed energy resources (DERs) in power distribution grids (DGs). Aiming to improve DG operational efficiency, as quantified by a generic optimal reactive power flow (ORPF) problem, we propose a two-stage approach. The first stage involves learning the manifold of optimal operating points determined by an ORPF instance. To synthesize local Volt/Var controllers, the learning task is partitioned into learning local surrogates (one per DER) of the optimal manifold with voltage input and reactive power output. Since these surrogates characterize efficient DG operating points, in the second stage, we develop local control schemes that steer the DG to these operating points. We identify the conditions on the surrogates and control parameters to ensure that the locally acting controllers collectively converge, in a global asymptotic sense, to a DG operating point agreeing with the local surrogates. We use neural networks to model the surrogates and enforce the identified conditions in the training phase. AC power flow simulations on the IEEE 37-bus network empirically bolster the theoretical stability guarantees obtained under linearized power flow assumptions. The tests further highlight the optimality improvement compared to prevalent benchmark methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Harnessing the power of gradient-based simulations for multi-objective optimization in particle accelerators

Abstract Particle accelerator operation requires simultaneous optimization of multiple objectives. Multi-objective optimization (MOO) is particularly challenging due to trade-offs between the objectives. Evolutionary algorithms, such as genetic algorithms (GAs), have been leveraged for many optimization problems, however, they do not apply to complex control problems by design. This paper demonstrates the power of differentiability for solving MOO problems in particle accelerators using a deep differentiable reinforcement learning (DDRL) algorithm. We compare the DDRL algorithm with model-free reinforcement learning (MFRL), GA, and Bayesian optimization (BO) for simultaneous optimization of heat load and trip rates in the continuous electron beam accelerator facility. The underlying problem enforces strict constraints on both individual states and actions as well as cumulative (global) constraints on energy requirements of the beam. Using historical accelerator data, we develop a physics-based surrogate model which is differentiable and allows for back-propagation of gradients. The results are evaluated in the form of a Pareto-front with two objectives. We show that the DDRL outperforms MFRL, BO, and GA on high dimensional problems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Global Value Chain and Manufacturing Analysis on Geothermal Power Plant Turbines

Currently, geothermal project developers customize the size of the power plant to fit the resource being developed. The turbine is designed and sized to optimize efficiency and resource utilization for electricity production; most often, other power plant components are then chosen to complement the turbine design. These custom turbine designs demand one-off manufacturing processes, which result in higher manufacturing setup costs, longer lead-times, and higher capital costs overall in comparison to larger-volume line manufacturing processes. In contrast, turbines produced in standard increments, manufactured in larger volumes, could result in lower costs per turbine. In this study, we developed a manufacturing cost model to identify requirements for equipment, facilities, raw materials, and labor. Then we used discounted cash flow (DCF) analysis to calculate the minimum sustainable price (MSP). We analyzed three different cases 1) 1 MWe geothermal ORC turboexpander 2) 5 MWe ORC turboexpander and 3) 20 MWe geothermal steam turbine, and calculated the cost of manufacturing the major components, such as the impellers/blades, shaft/rotor, nozzles, inlet guide lanes, disks, and casings. The results showed that MSP could highly vary between 893 /kW and 30 /kW based on turbine size, standardization and volume of manufacturing. The analysis also showed that the economy of scale applies both to the size of the turbine and the number manufactured in a single run. As an example, the unit price of a 5 MW standard design turbine could be 150 /W cheaper than the custom design. Sensitivity analysis indicated that these savings come largely from reduced labor costs for design and engineering and manufacturing setup. In addition to manufacturing cost savings, there is a delivery time saving up to 10 months, which could have a positive effect to construction financing operation time. Another advantage of these standard turbines is their adaptability to different geothermal systems by operating at off-design conditions. We also modeled off-design turbine efficiencies within the commercially-favorable operating range of a standard ORC compared to custom-designed ORC equipment. Than we used these performance calculations and power output capacities in a discounted cash flow (DCF) analysis of plant operations and financing, thereby creating representative techno economic models of a total geothermal power plant. I this scope we used NREL's System Advisor model (SAM) and performed DCF analysis of standard and custom design turbines using results from IPSEpro over range of temperatures of interest. In this scope 63 different off-design cases were analyzed. This helped us to answer the question; 'Can today's capital cost savings compensate the future revenue losses due to lower electricity generation?' The results showed that, the net capital cost savings for the modeled 5 MW case study may reach up to plus 1,298,672. Also, the difference between NPV of standard design and custom design turbine could reach up to 733,940. To sum up, when all the things mentioned above are taken into consideration, it can be said that, the standard turbines are competitive over a wide range of temperatures and flow rates and the give positive NPV differences for certain cases near their design point.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

2018 Research Highlight in Manufacturing: Global Value Chain and Manufacturing Analysis on Geothermal Power Plant Turbines

Currently, geothermal project developers customize the size of the power plant to fit the resource being developed. The turbine is designed and sized to optimize efficiency and resource utilization for electricity production; most often, other power plant components are then chosen to complement the turbine design. These custom turbine designs demand one-off manufacturing processes, which result in higher manufacturing setup costs, longer lead-times, and higher capital costs overall in comparison to larger-volume line manufacturing processes. In contrast, turbines produced in standard increments, manufactured in larger volumes, could result in lower costs per turbine. In this study, we developed a manufacturing cost model to identify requirements for equipment, facilities, raw materials, and labor. Then we used discounted cash flow (DCF) analysis to calculate the minimum sustainable price (MSP). We analyzed three different cases 1) 1 MWe geothermal ORC turboexpander 2) 5 MWe ORC turboexpander and 3) 20 MWe geothermal steam turbine, and calculated the cost of manufacturing the major components, such as the impellers/blades, shaft/rotor, nozzles, inlet guide lanes, disks, and casings. The results showed that MSP could highly vary between 893 $/kW and 30 $/kW based on turbine size, standardization and volume of manufacturing. The analysis also showed that the economy of scale applies both to the size of the turbine and the number manufactured in a single run. As an example, the unit price of a 5 MW standard design turbine could be 150 $/W cheaper than the custom design. Sensitivity analysis indicated that these savings come largely from reduced labor costs for design and engineering and manufacturing setup. In addition to manufacturing cost savings, there is a delivery time saving up to 10 months, which could have a positive effect to construction financing operation time. Another advantage of these standard turbines is their adaptability to different geothermal systems by operating at off-design conditions. We also modeled off-design turbine efficiencies within the commercially-favorable operating range of a standard ORC compared to custom-designed ORC equipment. Than we used these performance calculations and power output capacities in a discounted cash flow (DCF) analysis of plant operations and financing, thereby creating representative techno economic models of a total geothermal power plant. I this scope we used NREL's System Advisor model (SAM) and performed DCF analysis of standard and custom design turbines using results from IPSEpro over range of temperatures of interest. In this scope 63 different off-design cases were analyzed. This helped us to answer the question; 'Can today`s capital cost savings compensate the future revenue losses due to lower electricity generation?' The results showed that, the net capital cost savings for the modeled 5 MW case study may reach up to +$1,298,672. Also, the difference between NPV of standard design and custom design turbine could reach up to +$733,940. To sum up, when all the things mentioned above are taken into consideration, it can be said that, the standard turbines are competitive over a wide range of temperatures and flow rates and the give positive NPV differences for certain cases near their design point.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗