Search NASA⌕ Search

SEARCH · Search NASA

Results for “curves”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dispersion curve regression – assisted wideband local wavenumber analysis for characterizing three-dimensional (3D) profile of hidden corrosion damage

Corrosion as common damage in civil, petrochemical, nuclear, and aerospace structures affects the integrity and safety of the structures and may lead to catastrophic failures. This paper presents dispersion curve regression-assisted local wavenumber analysis method, which can analyze the time–space wavefield containing wideband information of wave-damage interaction and further extract the structural information carried by such wavefield for characterizing hidden corrosion damage in an isotropic plate. To acquire the time–space wavefield, a noncontact scanning laser Doppler vibrometer is used. In our analysis method, the acquired time–space wavefield is firstly processed to generate a local wavenumber function. The derived wavenumber function is further analyzed with a dispersion curve regression step, which searches in a set of theoretical frequency-wavenumber dispersion curves for different plate thicknesses and identifies the optimal plate thickness whose theoretical curve best matches the frequency-wavenumber relation contained in the local wavenumber function. By this means, a 3D profile including both in plane and thickness dimensions of the structure can be constructed for corrosion visualization and quantification. The experimental study demonstrates that our method can quantify the profile of hidden uniform corrosion in metal plates in terms of its in-plane shape and size as well as its out-of-plane depth in the subwavelength scale. Moreover, compared to previous Lamb wave-based corrosion inspection methods, our technique allows for noncontact 3D characterization of hidden corrosion from a far distance to the structure. Finally, we believe this work will inspire new 3D damage quantification methods that are based on wavefield analysis as well as enable potential applications for the quantification of hidden corrosion in civil, petrochemical, nuclear, and aerospace structures.

42 ENGINEERING↗
DOE OSTI2022

Phase curves of small bodies from the SLOAN Moving Objects Catalog

Extensive photometric surveys continue to produce enormous stores of data on small bodies. These data are typically sparsely obtained at arbitrary (or unknown) rotational phases. Therefore, new methods for processing such data need to be developed to make the most of these vast catalogs. We aim to produce a method of recreating the phase curves of small bodies by considering the uncertainties introduced by the nominal errors in the magnitudes and the effect introduced by rotational variations. Here, we use the SLOAN Moving Objects Catalog data as a benchmark to construct phase curves of all small bodies in u', g', r', i', and z' filters. From the phase curves, we obtain the absolute magnitudes and we use them to set up the absolute colors, which are the colors of the asteroids that are not affected by changes in the phase angle. We selected objects with ≥3 observations taken in at least one filter and spanning over a minimum of 5 degrees in the phase angle. We developed a method that combines Monte Carlo simulations and Bayesian inference to estimate the absolute magnitudes using the HG 12 * photometric system. We obtained almost 15 000 phase curves, with about 12 000 of these including all five filters. The absolute magnitudes and absolute colors are compatible with previously published data that support our method. The method we developed is fully automatic and well suited for a run based on large amounts of data. Moreover, it includes the nominal uncertainties in the magnitudes and the whole distribution of possible rotational states of the objects producing what are possibly less precise values, that is, larger uncertainties, but more accurate, namely, closer to the actual value. To our knowledge, this work is the first to include the effect of rotational variations in such a manner.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Are light curve classification metrics good proxies for SN Ia cosmological constraining power?

Context. When selecting a light curve classifier for use as part of a photometric supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, such as the contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would eliminate the computational expense of a full cosmology forecast in the analysis pipeline design process. Aims. This study tests the assumption that light curve classification metrics are an appropriate proxy for cosmology metrics. Methods. We emulated photometric SN Ia cosmology light curve samples with controlled contamination rates of individual contaminant classes and evaluated each of them under a set of classification metrics. We then derived cosmological parameter constraints from all samples under two common analysis approaches and quantified the impact of contamination by each contaminant class on the resulting cosmological parameter estimates. Results. We observe that cosmology metrics are sensitive to both the contamination rate and the class of the contaminating population, whereas the classification metrics are shown to be insensitive to the latter. Conclusions. Based on these findings, we discourage any exclusive reliance on light curve classification-based metrics for analysis design decisions, which (counterintuitively) include but are not limited to the classifier choice. Instead, we recommend optimising science analysis pipeline design choices using a metric of the information gained about the physical parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Solitary wave solutions of the 2+1 and 3+1 dimensional nonlinear Dirac equation constrained to planar and space curves

Here, we study the effect of curvature and torsion on the solitons of the nonlinear Dirac equation considered on planar and space curves. Since the spin connection is zero for the curves considered here, the arc variable provides a natural setting to understand the role of curvature. We obtain for various curves in two and three dimensions the transformation of the Dirac equation to the 1 + 1 dimensional flat space Dirac equation written in terms of the arc length. Depending on the curvature, the soliton profile either narrows or expands. Our results may be applicable to yet-to-be-synthesized curved quasi-one dimensional Bose condensates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The effect of jet–ejecta interaction on the viewing angle dependence of kilonova light curves

ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $1.5\, \mathrm{h}$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

Klion, Hannah↗
DOE OSTI2023

Evidence for a luminosity–decay correlation in GRB GeV light curves

ABSTRACT Correlations between intrinsic properties of gamma-ray burst (GRB) light curves provide clues to the nature of the central engine, the jet, and a possible means to standardize GRBs for cosmological use. Here, we report on the discovery of a correlation between the intrinsic early-time luminosity, LG,10 s, measured at rest frame 10 s, and the average decay rate measured from rest frame 10 s onward, $\alpha _{\mathrm{G,avg\gt 10\, s}}$, in a sample of 13 Fermi Large Area Telescope long GRB light curves. We note that our selection criteria, in particular the requirement for a redshift to construct luminosity light curves, naturally limits our sample to energetic GRBs. A Spearman’s rank correlation gives a coefficient of –0.74, corresponding to a confidence level of 99.6 per cent, indicating that brighter afterglows decay faster than less luminous ones. Assuming a linear relation with log(LG,10s), we find $\alpha _{\mathrm{G,avg\gt 10\, s}}$$= -0.31_{-0.09}^{+0.12}\log ($LG,10s$) + 14.43_{-5.97}^{+4.55}$. The slope of −0.31 is consistent at 1σ with previously identified correlations in the optical/ultraviolet and X-ray light curves. We speculate that differences in the rate at which energy is released by the central engine or differences in observer viewing angle may be responsible for the correlation.

Hinds, K. R. (ORCID:000000020129806X)↗
DOE OSTI2024

Deep drilling in the time domain with DECam – II: characterizing the light curves of candidates in the extragalactic fields

ABSTRACT In this second paper on the DECam deep-drilling field (DDF) program, we release 2020 optical $gri$-band light curves for transients and variables in the extragalactic COSMOS and ELAIS fields based on time series observations with a 3-d cadence from semester 2021A through 2023A. In order to demonstrate the wide variety of time domain events detected by the program and encourage others to use the data set, we characterize the sample by presenting a brief analysis of the light-curve parameters such as time span, amplitude, and peak brightness. We also present preliminary light-curve categorizations, and identify potential stellar variables, active galactic nuclei, tidal disruption events, supernovae (SNe) (such as Type Ia, Type IIP, superluminous, and gravitationally lensed SNe), and fast transients. Where relevant, the number of identified transients is compared to the predictions of the original proposal. We also discuss the challenges of analyzing DDF data in the context of the upcoming Vera C. Rubin Observatory and its Legacy Survey of Space and Time, which will include DDFs. Images from the Dark Energy Camera DDF program are available without proprietary period and the light curves presented in this work are publicly available for analysis.

Graham, Melissa L. (ORCID:0000000291543136)↗
DOE OSTI2023

Designing dynamically corrected gates robust to multiple noise sources using geometric space curves

Noise-induced gate errors remain one of the main obstacles to realizing a broad range of quantum information technologies. Dynamical error suppression using carefully designed control schemes is critical for overcoming this challenge. Such schemes must be able to correct against multiple noise sources simultaneously afflicting a qubit to reach error-correction thresholds. Here we present a general framework for designing control fields that simultaneously suppress both noise in the fields themselves as well as transverse dephasing noise. Using the recently developed space curve quantum control formalism, in which robust quantum evolution is mapped to closed geometric curves in a multidimensional Euclidean space, we derive the minimal conditions necessary to guarantee the simultaneous cancellation of both types of noise to leading order. In particular, we find that the cancellation of control field noise requires the derivative of the space curve to have zero-area projections, which is a much more subtle property compared to the closed-curve condition needed to suppress transverse dephasing. We present several techniques for solving both these conditions simultaneously and provide explicit examples of error-resistant control fields. Finally, our work also sheds light on the relation between holonomic evolution and the suppression of control field errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Proposed Material Properties, Allowable Stresses, and Design Curves of Diffusion Bonded Alloy 800H for the ASME Code Section III Division 5

Abstract Increased interest in compact heat exchangers (CHXs) to serve as intermediate heat exchangers of very high temperature reactors resulted in significant research and development on their design, analysis, and construction. Printed circuit heat exchangers are a type of CHXs with high thermal efficiency and compactness achieved through diffusion bonding a stack of etched plates with millimeter scaled channels. The diffusion bonding process changes the microstructural and mechanical properties of the wrought metal plates. The current nonnuclear design code ASME section VIII, division 1 captures the material property change through a “joint efficiency factor.” However, the current nuclear design code ASME section III, division 5 does not address or support the diffusion bonded material properties. Hence, there is a need to develop allowable stresses, isochronous curves, and fatigue life curves for various diffusion bonded alloys. In this study, Alloy 800H material was selected to establish the diffusion bonded material properties under tension, creep, fatigue, and creep-fatigue loads at elevated temperatures in the range 550–760 °C. A set of tests on diffusion bonded Alloy 800H (DB 800H) were performed and the acquired data are used in developing allowable stresses Sy, Su, Sr, Sm, St, Smt, So, isochronous curves and fatigue life curves according to the ASME section III, division 5 requirements. This paper also presents detailed procedures used in developing the ASME code section III division 5 design provisions for diffusion bonded Alloy 800H.

Engineering↗
DOE OSTI2022

Methodology to Calibrate Fragility Curves Using Limited Real-World Data

With increasing human dependence on electricity and increasing energy demand, electrical infrastructure has emerged to be one of the most critical services. This is one of the main reasons for making the power grid more resilient to extreme events during which uninterrupted power supply is crucial in keeping the consequences of the extreme event limited. The first step towards making the overall power grid more resilient is to evaluate the probabilities of failure for all assets at risk from the excitation caused by extreme events. Thus, this paper presents the methodology to calibrate fragility curves. The strength of the proposed approach is its ability to calibrate the fragility curves utilizing limited data, which is the most common constraint in carrying out such analyses. This paper describes the calibration of transmission tower fragility curves for Puerto Rico utilizing only damage reports for hurricane Maria. That, combined with Puerto Rico's wind modeling & geographic information system, this study also calibrates fragility curves for transmission towers made with four different construction materials.

Bereta dos Reis, Fernando↗
DOE OSTI2020

Smaller than Expected Bright-spot Offsets in Spitzer Phase Curves of the Hot Jupiter Qatar-1b

We present Spitzer full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature—and intermediate surface gravity and orbital period—as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse depths of 0.21 ± 0.02% at 3.6 μm and 0.30 ± 0.02% at 4.5 μm, corresponding to dayside brightness temperatures of 1542{sub −31}{sup +32} K and 1557{sub −36}{sup +35} K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are 1117{sub −71}{sup +76} and 1167{sub −74}{sup +69} K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of 0.12{sub −0.16}{sup +0.14} and a moderate day–night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with Spitzer full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: 11° ± 7° at 3.6 μm and −4° ± 7° at 4.5 μm. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Astraea: Predicting Long Rotation Periods with 27 Day Light Curves

The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27 day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4 yr light curves with 13% uncertainty overall (and a 9% uncertainty for periods >30 days). By training on 27 day Kepler light-curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods >30 days). After training this tool on these 27 day Kepler light-curve segments, we applied Astraea to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Mapping Stellar Surfaces. II. An Interpretable Gaussian Process Model for Light Curves

The use of Gaussian processes (GPs) as models for astronomical time series data sets has recently become almost ubiquitous, given their ease of use and flexibility. In particular, GPs excel at marginalization over the stellar signal when the variability due to starspots is treated as a nuisance, as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest, since it is not obvious how the parameters of the GP relate to physical parameters like the spot size, contrast, and latitudinal distribution. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is ill-posed and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many stars. Here we derive a closed-form expression for a GP that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more light curves, including the typical spot radii and latitudes. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from light curves and radial velocity measurements, as well as for modeling correlated noise in exoplanet searches. Our implementation is efficient, user-friendly, and open-source, available in the package starry-process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Early-time Light Curves of Type Ia Supernovae Observed with TESS

We present the early-time light curves of Type Ia supernovae (SNe Ia) observed in the first six sectors of Transiting Exoplanet Survey Satellite (TESS) data. Ten of these SNe were discovered by ASAS-SN, seven by ATLAS, six by ZTF, and one by Gaia. For nine of these objects with sufficient dynamic range (>3.0 mag from detection to peak), we fit power-law models and searched for signatures of companion stars. We found a diversity of early-time light-curve shapes, although most of our sources are consistent with fireball models where the flux increases as ∝t 2 . Three SNe displayed a flatter rise with flux ∝t. We did not find any obvious evidence for additional structures, such as multiple power-law components, in the early rising light curves. For assumptions about the SN properties and the observer viewing angle (ejecta mass of 1.4 M⊙, expansion velocity of 10 4 km s –1 , opacity of 0.2 cm 2 g –1 , and viewing angle of 45°) and a further assumption that any companion stars would be in Roche lobe overflow, it is possible to place upper limits on the radii of any companion stars. Six of the nine SNe had complete coverage of the early-time light curves, and we placed upper limits on the radii of companion stars of ≲32 R ⊙ for these SNe, ≲20 R ⊙ for five of the six, and ≲4 R ⊙ for two of the six. The small sample size did not allow us to put limits on the occurrence rate of companion stars in the progenitors of SNe Ia. However, we expect that TESS observed enough SNe in its two-year primary mission (26 sectors) to either detect the signature of a large companion (R > 20 R ⊙ ) or constrain the occurrence rate of such systems, at least for the fiducial SN properties adopted here. Here, we also show that TESS is capable of detecting emission from a 1 R ⊙ companion for an SN Ia within 50 Mpc and has a reasonable chance of doing so after about six years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The South Pole Telescope AGN Monitoring Campaign: First Release of SPTpol Bright AGN Light Curves

The South Pole Telescope (SPT) collaboration has recently embarked upon a campaign to monitor the brightness of a sample of active galactic nuclei (AGN), both in real time and in archival SPT data. The original design of the SPT was optimized for observations of the cosmic microwave background (CMB) at arc-minute and larger angular scales, and it has been used for this purpose for nearly twenty years, using three generations of CMB cameras. Recently it has been recognized that data from CMB experiments have the potential to be used for AGN monitoring. In this paper, we present the first public release of data from a full sample of SPT-monitored AGN, comprising 158 AGN light curves and associated data from the SPTpol camera, which was operational from 2012-2016. These light curves were created using observations from the SPTpol 500 deg$^{2}$ survey, in which the instrument was used to scan a 500 deg$^2$ patch of the sky several times per day with detectors sensitive to radiation in bands centered at 90 and 150 GHz. We provide a comprehensive description of the observations, the data processing methods, and the resulting light curve catalog. As an example of analyses that these data enable, we searched for a correlation between variability and spectral index, and we looked for ``bluer-when-brighter'' trends in the sample. Our analysis finds $> 10 σ$ correlation between fractional intrinsic variance and mean spectral index in the sample, but no significant evidence for bluer-when-brighter trends. The datasets from this study can be accessed through the SPT Treasury Record of AGN With Historical Activity and Time-Series or STRAWHAT catalog. This initial data release includes SPTpol light curves at 90 and 150 GHz, focusing on total intensity. In later updates, SPTpol polarization data and new observations from the SPT-3G instrument at 90, 150, and 220 GHz will be included.

Hood, J.C., II [Chicago U., KICP; Chicago U., Astr↗
DOE OSTI2022

Adapting Technology Learning Curves for Prospective Techno-Economic and Life Cycle Assessments of Emerging Carbon Capture and Utilization Pathways

Comparisons of emerging carbon capture and utilization (CCU) technologies with equivalent incumbent technologies are necessary to support technology developers and to help policy-makers design appropriate long-term incentives to mitigate climate change through the deployment of CCU. In particular, early-stage CCU technologies must prove their economic viability and environmental reduction potential compared to already-deployed technologies. These comparisons can be misleading, as emerging technologies typically experience a drastic increase in performance and decrease in cost and greenhouse gas emissions as they develop from research to mass-market deployment due to various forms of learning. These changes complicate the interpretation of early techno-economic assessments (TEAs) and life cycle assessments (LCAs) of emerging CCU technologies. The effects of learning over time or cumulative production themselves can be quantitatively described using technology learning curves (TLCs). While learning curve approaches have been developed for various technologies, a harmonized methodology for using TLCs in TEA and LCA for CCU in particular is required. To address this, we describe a methodology that incorporates TLCs into TEA and LCA to forecast the environmental and economic performance of emerging CCU technologies. This methodology is based on both an evaluation of the state of the art of learning curve assessment and a literature review of TLC approaches developed in various manufacturing and energy generation sectors. Additionally, we demonstrate how to implement this methodology using a case study on a CO 2 mineralization pathway. Finally, commentary is provided on how researchers, technology developers, and LCA and TEA practitioners can advance the use of TLCs to allow for consistent, high-resolution modeling of technological learning for CCU going forward and enable holistic assessments and fairer comparisons with other climate technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Solving the duck curve in a smart grid environment using a non-cooperative game theory and dynamic pricing profiles

With the intermittency that comes with electricity generation from renewables, utilizing dynamic pricing will encourage the demand-side to respond in a smart way that would minimize the electricity costs and flatten the net electricity demand curve. Determining the optimal dynamic pricing profile that would leverage distributed storage to flatten the curve is a novel idea that needs to be studied. Moreover, the economic feasibility of utilizing distributed electrical energy storage is still not given in the literature. Therefore, in this paper, a novel way of solving a citywide dynamic model using a bilevel programming algorithm is introduced. The problem is developed as a novel non-cooperative Stackelberg game that utilizes air-conditioning systems and electrical storage through the end-users to determine the optimal dynamic pricing profile. The results show that the combined effect of utilizing demand-side air-conditioning systems and distributed storage together can flatten the curve while employing the optimal dynamic pricing profile. An economic study is performed to determine the economic feasibility of 20 different cases with different battery designs and the level of solar penetration. Three metrics were used to evaluate the economic performance of each case: the levelized cost of storage, the levelized cost of energy, and the simple payback period. Most cases had levelized cost of storage values lower than 0.457 $/kWh, which is the lower bound available in the literature. Seven out of 16 cases have a simple payback period shorter than the lifetime of the system (25 years). The case with a 100 MW PV power plant and a battery storage of size 597 MWh, was found to be the most promising case with a simple payback period of 12.71 years for the photovoltaic plant and 19.86 years for the demand-side investments.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Modeling State Transition and Head-Dependent Efficiency Curve for Pumped Storage Hydro in Look-Ahead Dispatch

As one of the most widely installed utility-scale storage facilities, pumped storage hydro (PSH) plays an essential role in providing flexibility for power systems worldwide. Thus, to accurately quantify the flexibility of PSH units in operational optimization problems is important. The conventional PSH models in the literature rarely consider detailed state transitions due to their hourly-based settings. However, it becomes imperative for operational optimization on short-term intervals, especially with increasing shares of renewable energy in power systems. To this end, this paper presents a novel deterministic PSH model that considers the transition time and trajectory between three states of PSH units in look-ahead dispatch. Moreover, to better characterize the varying efficiency of PSH units with water head and flow rate, this paper proposes to model detailed headdependent efficiency curves (hereafter called input-output curves) in look-ahead dispatch, which encounters heavy computational burdens when short time intervals are applied. In this work, a zig-zag piece-wise linear approximation method is used for input-output curve modeling. This can enable an accurate quantification for variable efficiency and head dependence in a computationally effective manner. Here, numerical results are presented to show performances of the proposed PSH model in both flexibility quantification and computation time.

13 HYDRO ENERGY↗