SEARCH · Engineering Papers

Results for “variables”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,279 records · Page 71

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Variability in modelled reservoir greenhouse gas emissions: comparison of select US hydropower reservoirs against global estimates

Abstract Greenhouse gas (GHG) emissions from reservoirs have most often been evaluated on a global extent through areal scaling or linear-regression models. These models typically rely on a limited number of characteristics such as age, size, and average temperature to estimate per reservoir or areal flux. Such approaches may not be sufficient for describing conditions at all types of reservoirs. Emissions from hydropower reservoirs have received increasing attention as industry and policy makers seek to better understand the role of hydropower in sustainable energy solutions. In the United States (US), hydropower reservoirs span a wide range of climate regions and have diverse design and operational characteristics compared to those most heavily represented in model literature (i.e., large, tropical reservoirs). It is not clear whether estimates based on measurements and modeling of other subsets of reservoirs describe the diverse types of hydropower reservoirs in the US. We applied the Greenhouse Gas from Reservoirs (G‐res) emissions model to 28 hydropower reservoirs located in a variety of ecological, hydrological, and climate settings that represent the range of sizes and types of facilities within the US hydropower fleet. The dominant pathways for resulting GHG emissions estimates in the case-study reservoirs were diffusion of carbon dioxide, followed by methane ebullition. Among these case-study reservoirs, total post-impoundment areal flux of carbon ranges from 84 to 767 mgCm −2 d −1 , which is less variable than what has been reported through measurements at other US and global reservoirs. The net GHG reservoir footprint was less variable and towards the lower end of the range observed from modeling larger global reservoirs, with a range of 138 to 1,052 g CO 2 eq m −2 y −1 , while the global study reported a range of 115 to 145,472 g CO 2 eq m −2 y −1 . High variation in emissions normalized with respect to area and generation highlights the need to be cautious when using area or generation in predicting or communicating emissions footprints for reservoirs relative to those of other energy sources, especially given that many of the hydropower reservoirs in the US serve multiple purposes beyond power generation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗

DOE OSTI2020

Candidate Periodically Variable Quasars from the Dark Energy Survey and the Sloan Digital Sky Survey

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg 2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-year long multi-color ( griz ) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ~23.5 mag) powered by less massive black holes (with masses ≳ 10 8.5 M ⊙ ) at high redshift for the first time. We find five candidates with significant (at >99.74% single-frequency significance in at least two bands with a global p-value of ~7 × 10 -4 –3× 10 -3 accounting for the look-elsewhere effect) periodicity with observed periods of ~3–5 years (i.e., 1–2 years in rest frame) having ~4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of |${\sim }0.8^{+0.5}_{-0.3}$| % or |${\sim }1.1^{+0.7}_{-0.5}$| quasar per deg 2 . Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. Furthermore, we discuss implications on the future direct detection of low-frequency gravitational waves. Additionally, continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

First detection of VHE gamma-ray emission from TXS 1515–273, study of its X-ray variability and spectral energy distribution

ABSTRACT We report here on the first multiwavelength (MWL) campaign on the blazar TXS 1515–273, undertaken in 2019 and extending from radio to very-high-energy gamma-rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotron-peaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6σ. The average integral VHE flux of the source is 6 ± 1 per cent of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMM–Newton, and NuSTAR. The long continuous X-ray observations were separated by ∼9 h, both showing clear hour scale flares. In the XMM–Newton data, both the rise and decay time-scales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability time-scales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme synchrotron-peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band.

Acciari, V. A.↗

DOE OSTI2022

Radio transients and variables in the tenth Deeper, Wider, Faster observing run

The Deeper, Wider, Faster (DWF) program coordinates observations with telescopes across the electromagnetic spectrum, searching for transients on time-scales of milliseconds to days. The tenth DWF observing run was carried out in near real-time during September 2021, and consisted of six consecutive days of observations of the NGC 6744 galaxy group and a field containing the repeating fast radio burst FRB 190711 with the Australian Square Kilometre Array Pathfinder, the Dark Energy Camera, the Hard X-ray Modulation Telescope, and the Parkes 64-m ‘Murriyang’ radio telescope. Here, in this work, we present the results of an image-domain search for transient, variable, and circularly polarized sources carried out with ASKAP, using data from the observing run along with test observations prior to the run and follow-up observations carried out during and after the run. We identified eight variable radio sources, consisting of one pulsar, six stellar systems (five of which exhibit circularly polarized emission), and one previously uncatalogued source. Of particular interest is the detection of pulses from the ultra-cool dwarf SCR J1845–6357 with a period of 14.2 ± 0.3 h in good agreement with the known optical rotation period, making this the slowest rotating radio-loud ultra-cool dwarf discovered.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A step towards unveiling the nature of three cataclysmic variables: LS Cam, V902 Mon, and SWIFT J0746.3-1608

ABSTRACT We have carried out detailed time-resolved timing analyses of three cataclysmic variables (CVs) namely LS Cam, V902 Mon, and SWIFT J0746.3-1608, using the long-baseline, high-cadence optical photometric data from the Transiting Exoplanet Survey Satellite. Our analysis of LS Cam observations hints the presence of a superorbital period of ∼4.025 ± 0.007 d along with negative and positive superhump periods of ∼3.30 and 3.70 h, respectively. These results can be explained as an interaction of nodal and apsidal precession of the accretion disc with orbital motion. For the other two sources, V902 Mon and SWIFT J0746.3-1608, we have found evidence of a beat period of 2387.0 ± 0.6 and 2409.5 ± 0.7 s, respectively, which were not found in earlier studies. Our results presented in this study indicate the change in the accretion mode during the entire observing period for both sources. For V902 Mon, an apparent orbital period derivative of (6.09 ± 0.60) × 10−10 was also found. Moreover, the second harmonic of orbital frequency dominates the power spectrum of SWIFT J0746.3-1608, suggestive of ellipsoidal modulation of the secondary star. Present analyses suggest that LS Cam could be a superhumping CV, whereas V902 Mon and SWIFT J0746.3-1608 are likely to be variable disc-overflow accreting intermediate polars.

Rawat, Nikita (ORCID:0000000246336832)↗

DOE OSTI2023

Pulse profile modelling of thermonuclear burst oscillations – II. Handling variability

ABSTRACT Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts.

Astronomy & Astrophysics↗

DOE OSTI2024

Impact of survey spatial variability on galaxy redshift distributions and the cosmological 3 × 2-point statistics for the Rubin Legacy Survey of Space and Time (LSST)

We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry data set degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-z estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i\lt 25.3$ for 10 yr LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd i-band depth for 1-yr (Y1), 3-yr (Y3), and 5-yr (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10\!-\!40~{{\rm per\ cent}}$ at extreme depth values. Using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times 2$ pt analysis, and assess the impact on cosmological parameters via a Fisher forecast. We find that galaxy clustering is most susceptible to variable depth, and non-uniformity needs to be mitigated below 3 per cent to recover unbiased cosmological constraints. There is little impact on galaxy–shear and shear–shear power spectra, given the expected LSST Y3 noise.

cosmology↗

DOE OSTI2024

Probing the origin of the two-component structure of broad-line region by reverberation mapping of an extremely variable quasar

ABSTRACT The physical origins of quasar components, such as the broad-line region (BLR) and dust torus, remain under debate. To gain insights into them, we focused on changing-state quasars (CSQs) which provide a unique perspective through structural changes associated with accretion disc state transitions. We targeted SDSS J125809.31+351943.0, an extremely variable CSQ, to study its central core structure and kinematics. We conducted reverberation mapping with optical spectroscopy to explore the structure of the BLR and estimate the black hole mass. The results from Hβ reverberation mapping indicated a black hole mass of $10^{9.64^{+0.11}_{-0.20}}\rm {M_\odot }$. Additionally, we analysed variations in the optical to X-ray spectral indices, αox, before and after the state transition, to investigate the accretion disc. These variations in αox and the Eddington ratio (from 0.4 per cent to 2.4 per cent) exhibitied behaviour similar to state transitions observed in X-ray binary systems. Spectral analysis of Hβ revealed a predominantly double-peaked profile during dim periods, transitioning to include a single-peaked component as the quasar brightened, suggesting that Hβ contains a mixture of two components. Each of these components has its distinct characteristics: the first is a double-peaked profile that remains stable despite changes in the accretion rate, while the second is a variable single-peaked profile. Using time-lags from reverberation mapping, we estimated the spatial relationships between these BLR components, the accretion disc, and the dust torus. Our results suggest that the BLR consists of two distinct components, each differing in location and origin.

Nagoshi, Shumpei (ORCID:0000000326504322)↗

DOE OSTI2025

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗

DOE OSTI2024

Preserving isohydricity: vertical environmental variability explains Amazon forest water-use strategies

Abstract Increases in hydrological extremes, including drought, are expected for Amazon forests. A fundamental challenge for predicting forest responses lies in identifying ecological strategies which underlie such responses. Characterization of species-specific hydraulic strategies for regulating water-use, thought to be arrayed along an ‘isohydric–anisohydric’ spectrum, is a widely used approach. However, recent studies have questioned the usefulness of this classification scheme, because its metrics are strongly influenced by environments, and hence can lead to divergent classifications even within the same species. Here, we propose an alternative approach positing that individual hydraulic regulation strategies emerge from the interaction of environments with traits. Specifically, we hypothesize that the vertical forest profile represents a key gradient in drought-related environments (atmospheric vapor pressure deficit, soil water availability) that drives divergent tree water-use strategies for coordinated regulation of stomatal conductance (gs) and leaf water potentials (ΨL) with tree rooting depth, a proxy for water availability. Testing this hypothesis in a seasonal eastern Amazon forest in Brazil, we found that hydraulic strategies indeed depend on height-associated environments. Upper canopy trees, experiencing high vapor pressure deficit (VPD), but stable soil water access through deep rooting, exhibited isohydric strategies, defined by little seasonal change in the diurnal pattern of gs and steady seasonal minimum ΨL. In contrast, understory trees, exposed to less variable VPD but highly variable soil water availability, exhibited anisohydric strategies, with fluctuations in diurnal gs that increased in the dry season along with increasing variation in ΨL. Our finding that canopy height structures the coordination between drought-related environmental stressors and hydraulic traits provides a basis for preserving the applicability of the isohydric-to-anisohydric spectrum, which we show here may consistently emerge from environmental context. Our work highlights the importance of understanding how environmental heterogeneity structures forest responses to climate change, providing a mechanistic basis for improving models of tropical ecosystems.

Forestry↗

DOE OSTI2025

A rigorous framework for an improved Messinger/Myers model of ice accretion under conditions of variable property and unsteady aircraft icing

We analyse the Messinger/Myers model by critically evaluating simplifying assumptions through a rigorous formulation of the rime ice accretion process. We explore the effects of both constant and variable ice density and thermal conductivity, along with the effects of sublimation from the ice surface. The effects of key factors such as droplet impact rate, ambient temperature relative to the freezing temperature and the temperature difference between the ambient air and the airfoil surface are examined. Under these varying conditions, the present rigorous formulation is used to assess the significance of unsteady effects, variable ice properties and sublimation. We observe that the Myers model performs remarkably well in certain icing situations and analyse the reasons for this strong performance. We also show that partially relaxing the model’s assumptions can lead to poorer performance. The Myers model can lead to overprediction of ice surface temperature and correspondingly underprediction of transition time under conditions of relatively weak sublimation and surface cooling. A modified Myers model is presented, which can be used to recover near-perfect results under widely varying icing conditions of relevance. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

DOE OSTI2026

Concatenated dual displacement code for continuous-variable quantum error correction

The continuous-variable (CV) Gaussian no-go theorem fundamentally limits the suppression of Gaussian displacement errors using only Gaussian gates and states. Prior studies have employed Gottesman-Kitaev-Preskill (GKP) states as ancillary qumodes to suppress small Gaussian displacement errors. However, when the displacement magnitude becomes large, inevitable lattice-crossing errors arise beyond the correctable range of the GKP state. To address this issue, we concatenate the Gaussian-noise-suppression circuit with an outer analog Steane code that corrects such occasional lattice-crossing events as well as other abrupt displacement errors. Contrary to conventional concatenation, which primarily aims to reduce logical error rates, the Steane-GKP duality in encoding provides complementary protection against displacement errors at different scales: The inner GKP layer employs non-Gaussian resources to suppress continuous Gaussian noise and reduce residual variance, while the outer analog Steane code corrects discrete lattice-crossing events that exceed the GKP correctable range. It is precisely this separation of error-mitigation roles that enables CV error correction. In contrast to prior work on concatenating GKP and repetition codes to establish error correction for discrete qubit/qudit encoding, we provide correction in the continuous encoding space. Analytical studies show that, under infinite squeezing, the concatenated code suppresses the variance of Gaussian displacement errors acting on all qumodes by up to 50%, while enabling unbiased correction of lattice-crossing errors with a success probability determined by the ratio between the residual Gaussian error standard deviation and the lattice-crossing magnitude. Even with finite squeezing, the proposed architecture still provides Gaussian-error suppression and lattice-crossing correction. Moreover, the presence of the outer analog Steane code relaxes the squeezing requirement of the inner GKP states, indicating near-term experimental feasibility. This work establishes a viable route toward fault-tolerant continuous-variable quantum computation and provides insight into the design of concatenated CV error-correcting architectures.

quantum error correction↗

DOE OSTI2024

Exploration of the polarization angle variability of the Crab Nebula with POLARBEAR and its application to the search for axionlike particles

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau A, and for understanding the validity of Tau A as a calibrator. One intriguing application of such observation is to use it for the search of axionlike particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz P instrument. We place a median 95% upper bound of polarization oscillation amplitude A < 0.06 5 ° over the oscillation frequencies from 0.75 year − 1 to 0.66 hour − 1 . Assuming that no sources other than ALP are causing Tau A’s polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling g a γ γ < 2.16 × 10 − 12 GeV − 1 × ( m a / 10 − 21 eV ) in the mass range from 9.9 × 10 − 23 eV to 7.7 × 10 − 19 eV . This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs. Published by the American Physical Society 2024

Adachi, Shunsuke (ORCID:0000000204007555)↗

DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

A Dual-Active-Bridge Converter Employing a Variable Inductor Without an Auxiliary Circuit

This paper proposes a dual active bridge (DAB) converter employing a variable inductor (VI) without an auxiliary circuit. Unlike conventional VI-based designs that require an external DC bias circuit, the proposed DAB utilizes the input DC current itself as the bias source, enabling automatic inductance variation with load conditions. The inductance naturally increases at low power and decreases at high power, effectively extending the zero voltage switching (ZVS) range and reducing circulating current, respectively. The VI was experimentally implemented and characterized to obtain its inductance-current profile, which was then integrated into a PLECS model of the DAB converter for circuit and thermal simulations. Simulation results confirm that the proposed auxiliary-free VI-DAB converter achieves a wider ZVS range and lower circulating current compared with a conventional fixed-inductor DAB converter. By realizing a variable inductor without any auxiliary bias circuitry, the proposed approach maintains soft switching and reduces reactive current losses across a wide load range, leading to improved efficiency and simplified implementation.

Jo, Cheolhui [ORNL] (ORCID:0000000322692434)↗

DOE OSTI2020

Capacitor Voltage Ripple Estimation and Optimal Sizing of Modular Multi-Level Converters for Variable-Speed Drives

This article deals with the sizing of submodule capacitors for modular multilevel converters (MMCs) used for variable speed drives for medium voltage and high voltage electric machines. Variable speed drives are required to run over a wide frequency range from zero to a few hundred Hz or even over 1 kHz. Under low frequency operation (<30 Hz), a large current ripple flowing through the MMC submodule capacitors will cause significant capacitor voltage fluctuation. This article proposes an optimal submodule capacitor selection method based on the common-mode voltage injection. The proposed sizing model includes two functions: 1) accurately estimating the capacitor voltage ripple based on the drive system parameters; and 2) identifying the minimum capacitance of the submodule capacitors based on the allowed maximum capacitor voltage fluctuations. Simulation and experimental results show that the error between the actual and estimated values of the voltage fluctuation is within 3% using the proposed model.

42 ENGINEERING↗

DOE OSTI2021

Multi-Time Scale Coordinated Control and Scheduling of Inverter-Based TCLs With Variable Wind Generation

To address microgrid tie flow errors caused by wind generation variability, here we propose and develop a multi-time scale coordinated control and scheduling strategy for inverter-based thermostatically controlled loads (TCLs). First, in hour-time scale, inverter-based TCLs with adjusting temperature set-point are modeled as virtual generators to compensate tie flow deviations in the day-ahead plan. Next, in minute-time scale, virtual batteries representing operating behaviors of inverter-based TCLs with frequency control are scheduled determined by the control of virtual generators in hour-time scale. The virtual batteries are scheduled to smooth out tie flow errors corresponding to day-ahead plan and hour-time scale schedules. The multi-time scale control methods are coordinated to employ the response potential of inverter-based TCLs and response curve-based methods are proposed to control inverter-based TCLs considering the customer privacy. The multi-time scale stochastic schedules which are based on response curves of inverter-based TCLs are coordinated to accommodate wind generation variability. Simulation results demonstrate that the microgrid tie flow errors are effectively mitigated by the proposed multi-time scale coordinated control and scheduling of inverter-based TCLs.

24 POWER TRANSMISSION AND DISTRIBUTION↗