Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable predictive control”

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 235 records · Page 13

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah)

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

60 APPLIED LIFE SCIENCES↗

Engine identification for adaptive control

An attempt to obtain a dynamic model for a turbofan gas turbine engine for the purpose of adaptive control is described. The requirements for adaptive control indicate that a dynamic model should be identified from data sampled during engine operation. The dynamic model identified was of the form of linear differential equations with time varying coefficients. A turbine engine is, however, a highly nonlinear system, so the identified model would be valid only over a small area near the operating point, thus requiring frequent updating of the coefficients in the model. Therefore it is necessary that the identifier use only recent information to perform its function. The identifier selected minimized the square of the equation errors. Known linear systems were used to test the characteristics of the identifier. It was found that the performance was dependent on the number of data points used in the computations and upon the time interval over which the data points were obtained. Preliminary results using an engine deck for the quiet, clean, shorthaul experimental engine indicate that the identified model predicts the engine motion well when there is sufficient dynamic information, that is when the engine is in transient operation.

Leonard, R. G.↗

Metric Learning to Accelerate Convergence of Operator Splitting Methods

Recent developments in machine learning have led to promising advances in accelerating the solution of constrained optimization problems. Increasing demand for real-time decision-making capabilities in applications such as artificial intelligence and optimal control has led to a variety of proposed strategies for learning to produce fast solutions to optimization problems. For example, recent works have shown that it is possible to accelerate the convergence of optimization algorithms by learning to select their parameters, such as gradient descent stepsizes. This work proposes a new approach, in which the underlying metric spaces of proximal operator splitting algorithms are learned to maximize convergence rate. While prior works in optimization theory have derived optimal metrics in simple cases, no such result exists for many practical problem forms including general Quadratic Programming (QP). This paper shows how differentiable optimization can enable the end-to-end learning of proximal metrics, enhancing the convergence of proximal algorithms for QP problems beyond what is possible based on known theory. Additionally, the results illustrate a strong connection between the learned proximal metrics and active constraints at the optima, leading to an interpretation in which the predicted proximal metrics can be viewed as a form of active set prediction.

King, Ethan [BATTELLE (PACIFIC NW LAB)]↗

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗

Detached Bridgman Growth of Germanium and Germanium-Silicon Alloy Crystals

Earth based experiments on the science of detached crystal growth are being conducted on germanium and germanium-silicon alloys (2at% Si average composition) in preparation for a series of experiments aboard the International Space Station (ISS) to differentiate among proposed mechanisms contributing to detachment. Sessile drop measurements were first carried out for a large number of substrates made of potential ampoule materials to determine the contact angles and the surface tension as a function of temperature and composition. The process atmosphere and duration of the experiment (for some cases) were also found to have significant influence on the wetting angle. Growth experiments have used pyrolytic boron nitride (pBN) and fused silica ampoules with the majority of the detached results occurring predictably in the pBN. The contact angles were 173 deg (Ge) and 165 deg (GeSi) for pBN. For fused silica, the contact angle decreases to an equilibrium value with duration of measurement ranging from 150 to 117 deg (Ge), 129 to 100 deg (GeSi). Forming gas (Ar + 2% H2) and vacuum have been used in the growth ampoules. With gas in the ampoule, a variation of the temperature profile during growth has been used to control the pressure difference between the top of the melt and the volume below the melt caused by detachment of the growing crystal. The stability of detachment has been modeled and substantial insight has been gained into the reasons that detachment has most often been observed in reduced gravity but nonetheless has occurred randomly even there. An empirical model for the conditions necessary to achieve sufficient stability to maintain detached growth for extended periods has been developed and will be presented. Methods for determining the nature and extent of detachment include profilometry and optical and electron microscopy. This surface study is the subject of another presentation at this Congress. Results in this presentation will show that we have established the effects of different ampoule materials, temperature profiles, pressure differences, and silicon concentrations and that samples that are nearly completely detached can be grown repeatedly.

Szofran, F. R.↗

Structural characterization and regulatory element analysis of the heart isoform of cytochrome c oxidase VIa

In order to investigate the mechanism(s) governing the striated muscle-specific expression of cytochrome c oxidase VIaH we have characterized the murine gene and analyzed its transcriptional regulatory elements in skeletal myogenic cell lines. The gene is single copy, spans 689 base pairs (bp), and is comprised of three exons. The 5'-ends of transcripts from the gene are heterogeneous, but the most abundant transcript includes a 5'-untranslated region of 30 nucleotides. When fused to the luciferase reporter gene, the 3.5-kilobase 5'-flanking region of the gene directed the expression of the heterologous protein selectively in differentiated Sol8 cells and transgenic mice, recapitulating the pattern of expression of the endogenous gene. Deletion analysis identified a 300-bp fragment sufficient to direct the myotube-specific expression of luciferase in Sol8 cells. The region lacks an apparent TATA element, and sequence motifs predicted to bind NRF-1, NRF-2, ox-box, or PPAR factors known to regulate other nuclear genes encoding mitochondrial proteins are not evident. Mutational analysis, however, identified two cis-elements necessary for the high level expression of the reporter protein: a MEF2 consensus element at -90 to -81 bp and an E-box element at -147 to -142 bp. Additional E-box motifs at closely located positions were mutated without loss of transcriptional activity. The dependence of transcriptional activation of cytochrome c oxidase VIaH on cis-elements similar to those found in contractile protein genes suggests that the striated muscle-specific expression is coregulated by mechanisms that control the lineage-specific expression of several contractile and cytosolic proteins.

Non-NASA Center↗

NO(x) Concentrations in the Upper Troposphere as a Result of Lightning

Upper tropospheric NO(x) controls, in part, the distribution of ozone in this greenhouse-sensitive region of the atmosphere. Many factors control NO(x) in this region. As a result it is difficult to assess uncertainties in anthropogenic perturbations to NO from aircraft, for example, without understanding the role of the other major NO(x) sources in the upper troposphere. These include in situ sources (lightning, aircraft), convection from the surface (biomass burning, fossil fuels, soils), stratospheric intrusions, and photochemical recycling from HNO3. This work examines the separate contribution to upper tropospheric "primary" NO(x) from each source category and uses two different chemical transport models (CTMS) to represent a range of possible atmospheric transport. Because aircraft emissions are tied to particular pressure altitudes, it is important to understand whether those emissions are placed in the model stratosphere or troposphere and to assess whether the models can adequately differentiate stratospheric air from tropospheric air. We examine these issues by defining a point-by-point "tracer tropopause" in order to differentiate stratosphere from troposphere in terms of NO(x) perturbations. Both models predict similar zonal average peak enhancements of primary NO(x) due to aircraft (approx. = 10-20 parts per trillion by volume (pptv) in both January and July); however, the placement of this peak is primarily in a region of large stratospheric influence in one model and centered near the level evaluated as the tracer tropopause in the second. Below the tracer tropopause, both models show negligible NO(x) derived directly from the stratospheric source. Also, they predict a typically low background of 1 - 20 pptv NO(x) when tropospheric HNO3 is constrained to be 100 pptv of HNO3. The two models calculate large differences in the total background NO(x) (defined as the source of NO(x) from lightning + stratosphere + surface + HNO3) when using identical loss frequencies for NO(x). This difference is primarily due to differing treatments of vertical transport. An improved diagnosis of this transport that is relevant to NO(x) requires either measurements of a surface-based tracer with a substantially shorter lifetime than Rn-222 or diagnosis and mapping of tracer correlations with different source signatures. Because of differences in transport by the two models we cannot constrain the source of NO(x) from lightning through comparison of average model concentrations with observations of NO(x).

Penner, Joyce E.↗

Sex ratio of Western Bluebirds ( Sialia mexicana ) is mediated by phenology and clutch size

Mothers may produce more of one sex to maximize their fitness if there are differences in the cost of producing each sex or there are differences in their relative reproductive value. Breeding date and clutch size are known to influence offspring sex ratios in birds through sex differences in dispersal, social behaviours, differential mortality and available food resources. We tested whether breeding date, clutch size and drought conditions influenced offspring sex ratios in a sexually size-monomorphic species, the Western Bluebird Sialia mexicana, by interrogating a 21-year dataset. After controlling for differential mortality, we found that hatch dates late in the breeding season were associated with the production of more females, suggesting that the value of producing males declines as the breeding season progresses. When clutch size was taken into account, small clutches yielded significantly more females late in the breeding season than in the early and middle parts of the breeding season, which produced significantly more males. Here, large clutches early in the season tended to produce more females, although this was not significant. Drought severity was not correlated with sex ratio adjustment. We propose and discuss several explanations for these patterns, including male offspring, but not female offspring, acting as helpers, increased female nestling provisioning late in the breeding season, differences in food abundance and egg-laying order. Future work will help to uncover the mechanisms leading to these patterns. Identifying patterns and mechanisms of sex ratio skew from long-term datasets is important for informing predictions regarding life-history trade-offs in wildlife populations.

59 BASIC BIOLOGICAL SCIENCES↗

Integrated Modeling Tools for Thermal Analysis and Applications

Integrated modeling of spacecraft systems is a rapidly evolving area in which multidisciplinary models are developed to design and analyze spacecraft configurations. These models are especially important in the early design stages where rapid trades between subsystems can substantially impact design decisions. Integrated modeling is one of the cornerstones of two of NASA's planned missions in the Origins Program -- the Next Generation Space Telescope (NGST) and the Space Interferometry Mission (SIM). Common modeling tools for control design and opto-mechanical analysis have recently emerged and are becoming increasingly widely used. A discipline that has been somewhat less integrated, but is nevertheless of critical concern for high precision optical instruments, is thermal analysis and design. A major factor contributing to this mild estrangement is that the modeling philosophies and objectives for structural and thermal systems typically do not coincide. Consequently the tools that are used in these discplines suffer a degree of incompatibility, each having developed along their own evolutionary path. Although standard thermal tools have worked relatively well in the past. integration with other disciplines requires revisiting modeling assumptions and solution methods. Over the past several years we have been developing a MATLAB based integrated modeling tool called IMOS (Integrated Modeling of Optical Systems) which integrates many aspects of structural, optical, control and dynamical analysis disciplines. Recent efforts have included developing a thermal modeling and analysis capability, which is the subject of this article. Currently, the IMOS thermal suite contains steady state and transient heat equation solvers, and the ability to set up the linear conduction network from an IMOS finite element model. The IMOS code generates linear conduction elements associated with plates and beams/rods of the thermal network directly from the finite element structural model. Conductances for temperature varying materials are accommodated. This capability both streamlines the process of developing the thermal model from the finite element model, and also makes the structural and thermal models compatible in the sense that each structural node is associated with a thermal node. This is particularly useful when the purpose of the analysis is to predict structural deformations due to thermal loads. The steady state solver uses a restricted step size Newton method, and the transient solver is an adaptive step size implicit method applicable to general differential algebraic systems. Temperature dependent conductances and capacitances are accommodated by the solvers. In addition to discussing the modeling and solution methods. applications where the thermal modeling is "in the loop" with sensitivity analysis, optimization and optical performance drawn from our experiences with the Space Interferometry Mission (SIM), and the Next Generation Space Telescope (NGST) are presented.

Milman, Mark H.↗

Maps of land surface phenology derived from PlanetScope data, 2018-2022, Teller, Kougarok, and Council, Seward Peninsula

Remote sensing maps of land surface phenology derived from PlanetScope (Planet Team, 2017) normalized differential greenness index (NDGI; Yang et al., 2019) time series data. These maps include four phenological timing metric - start of spring (SOS), end of spring (EOS), start of fall (SOF), and end of fall (EOF), corresponding NDGI value at the four phenological timings, and annual maximum and minimum NDGI. This package includes maps for Next-Generation Ecosystem Experiment Arctic (NGEE Arctic)’s Teller Mile Marker (MM) 27, Kougarok MM64, and Council MM 71 watersheds. Maps of 5 years from 2018 to 2022 are included in this dataset. The map data and metadata are provided as image (ENVI) and text (*.txt, *hdr) formats. Additional supporting map quicklooks are provided as GIS *.kml files. These datasets are provided in support of Yang et al., (In revision) “Fine-scale Landscape Characteristics, Vegetation Composition, and Snowmelt Timing Control Phenological Heterogeneity across Arctic Tundra Landscapes”The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Reliable Real-Time Solution of Parametrized Partial Differential Equations: Reduced-Basis Output Bound Methods

We present a technique for the rapid and reliable prediction of linear-functional outputs of elliptic (and parabolic) partial differential equations with affine parameter dependence. The essential components are (i) (provably) rapidly convergent global reduced basis approximations, Galerkin projection onto a space W(sub N) spanned by solutions of the governing partial differential equation at N selected points in parameter space; (ii) a posteriori error estimation, relaxations of the error-residual equation that provide inexpensive yet sharp and rigorous bounds for the error in the outputs of interest; and (iii) off-line/on-line computational procedures, methods which decouple the generation and projection stages of the approximation process. The operation count for the on-line stage, in which, given a new parameter value, we calculate the output of interest and associated error bound, depends only on N (typically very small) and the parametric complexity of the problem; the method is thus ideally suited for the repeated and rapid evaluations required in the context of parameter estimation, design, optimization, and real-time control.

Prudhomme, C.↗

FutureTense

Protective vaccines and reliable diagnostics are essential tools for controlling viral diseases. However, the efficacy of these tools can be diminished by mutations in viral genomes. The delay between the emergence of new viral strains and the redesign of vaccines and diagnostics allows for continued viral transmission. Is it possible to address this challenge by computationally predicting viral genome sequence evolution? Can we “future-proof” vaccines and diagnostics by targeting both current and anticipated future sequence variants? While predicting viral evolution is still an unsolved, “grand challenge” problem in biology, the large, and rapidly growing, number of SARS-CoV-2 genome sequences provide an opportunity to quantify the ability of machine learning to predict viral genome sequence evolution. Towards this end, we have developed a simple computational model for predicting viral evolution at the level of individual nucleotides. The key metric for quantifying the per-base, prediction accuracy for viral evolution is the Mann-Whitney U statistic (or, equivalently, the area under the receiver operator curve). Since the Mann-Whitney U statistic is not a differentiable function, existing deep leaning packages (like Pytorch and Keras/TensorFlow) are not useful, as they require that the accuracy metric/objective function be analytically differentiable with respect to the model parameters. To overcome this challenge, we have implemented custom software, “FutureTense”, that can train a machine learning model by maximizing the non-differentiable Mann-Whitney U statistic. This software trains a machine learning model by exploring along the direction of the discrete gradient of the Mann-Whitney U statistic in the model parameter space. Parallel computing and genome sequence-specific optimizations are used to accelerate model training. The resulting machine learning model learns the observed high C->U mutation rates in the SARS-CoV-2 genome (which are potentially induced by host defenses) and provides prediction accuracies that are significantly better than one would expect from random chance. While predicting viral evolution is still quite far from a solved problem, the surprising performance of this simple model gives hope that the accuracy of predicting viral genome evolution can be further increased by more sophisticated approaches.

Gans, Jason↗

New Mathematical Functions for Vacuum System Analysis

A new bivariate function has been found that provides solutions of integrals having the form u (sup minus eta) e (sup u) du which arise when developing predictions for the behavior of pressure within a rigid volume under high vacuum conditions in the presence of venting as well as sources characterized by power law transient decay over the range [0,1] for eta and for u greater than or equal to 0. A few properties of the new function are explored in this work. For instance the eta equals 1/2 case reproduces the Dawson function. In addition, a slight variation of the solution technique reproduces the exponential integral for eta equals 1. The technique used to generate these functions leads to an approach for solving a more general class of nonlinear ordinary differential equations, with the potential for identifying other new functions that solve other integrals.

Dawson Function↗

Supercritical CO2 Recuperators, Presented at The ASME Turbo Expo 2017, June 29, 2017, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

High Effectiveness, Compact, High Pressure and Low-Cost Recuperator, Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗

High Effectiveness, Compact, High Pressure and Low Cost Recuperator for Super-Critical CO2 Power Cycles, Paper Presented at The Fifth International Symposium – Super-Critical CO2 Power Cycles, San Antonio, Texas, March 28-31, 2016, by Dr. John Kelly, President, Altex Technologies Corporation

Closed Brayton super-critical CO2 power cycles are well suited to waste heat bottoming cycles, due to their increased efficiency and compactness, relative to Rankine steam bottoming cycles. Since waste heat applications would be retrofits, the power system compactness is important. To achieve high efficiency, these power cycles require high pressure recuperative type heat exchangers, of substantial heat duty. Current Printed Circuit Heat Exchangers (PCHE), originally developed for high pressure gas and oil applications, can be used as recuperators, but costs are higher than desired. Altex is developing a purpose-built high pressure and effectiveness recuperator to provide the reliability, compactness, performance and pressure capability of current recuperators, but at a reduced cost. The High Effectiveness Low Cost (HELC) recuperative heat exchanger design yields volume and weight metrics of .0024 m3/UA and 10.2 kg/UA, which are 4% and 77.3% below recuperator target metrics, respectively. A 50 kW test article was designed and fabricated. Performance tests on water and oil showed that the HELC design model could predict heat transfer, to within 10% of the measured value. This model was then used to project HELC performance, when operating on supercritical CO2. Besides performance tests, the test article was hydrostatically tested, for integrity at up to 4,000 psi pressure. At these conditions, some distortion of channels was encountered. To mitigate distortion, the inserts were redesigned and these results were used to project the cost of 500 kW and 5,083 kW HELC units. The cost metric for the 5,083 kW unit was determined to be $1,349/UA, which is 10% lower than the recuperator desired cost metric desired target of $1,500/UA. In addition, the 5,083 kW unit HELC cost of $61.09/kW is 33.6% lower than the $92/kW estimated cost for a PCHE. Hydrostatic pressure tests, at up to 4,000psi, showed that the unit did not leak. However, channels were distorted at this pressure differential. Design updates to minimize stress concentrations and distortion were prepared and analyzed, to show that distortion could be controlled, but to date tests have not been run to prove the design. Project results show the potential of the HELC approach, but more work is required to confirm this potential, at the larger scales of interest.

14 SOLAR ENERGY↗