Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase estimation”

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

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Boundary layer roll circulations during FIRE

The probable mechanism underlying the development of boundary layer roll circulations are studied using wind and temperature profiles measured by the National Center for Atmospheric Research (NCAR) Electra during the stratocumulus phase of the First ISCCP Regional Experiment (FIRE). The expected, or preferred, roll orientations, horizontal wavelengths, and propagation periods are determined by finding the minimum values of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds number (Re) and moist Rayleigh number (Ra sub m). These minimum values depend on the height z sub T of the capping temperature inversion and on the values of the Fourier coefficients of the background height-dependent vector wind profile. As input to our nonlinear spectral model, descent and ascent runs by the Electra provide for initial estimates of the inversion height and the wind profiles. In the first phase of the investigation presented here, a mechanism is said to be a probable contributor to the development of roll circulations within the stratocumulus-topped boundary layer if the modeled roll orientation and wavelengths agree with their observed values. Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a sharp cloud boundary that was within a LANDSAT/SPOT scene. The stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud streets were oriented nearly parallel to the mean wind direction in the layer, which was approximately 340 degrees. The hypothesis that roll circulations occurred in both the relatively clear and the cloudy regions is investigated using as model input a descent profile obtained in the relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air case are that the pure inflection point mode is not possible and the pure thermal mode was oriented 35 degrees to the right of the mean wind direction. The origin of this unacceptably large discrepancy between the observed and modeled results will be investigated further and the conclusions reported at the next FIRE workshop.

Shirer, Hampton N.↗
NASA NTRS2006

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗
DOE OSTI2020

A Test Bed for Evaluating Frequency Estimation Algorithms in Synthetic Inertia Control: User Manual

As penetration of converter interfaced generators (CIGs) increases, the need for CIG frequency control participation increases. Traditionally, research in this area has been performed using positive sequence simulation software, which provides voltage magnitude and phase measurements, but not point-on-wave (POW) measurements. This means that the effect of frequency estimation algorithms cannot be accurately modeled, especially when the voltage waveform is distorted by faults or load connection events. This report serves as a user manual for an electromagnetic transient simulation testbed, which allows for accurate modeling of frequency estimation and control techniques.

42 ENGINEERING↗
DOE OSTI2022

Between Harmonic Crystal and Glass: Solids with Dimpled Potential-Energy Surfaces Having Multiple Local Energy Minima

Solids with dimpled potential-energy surfaces are ubiquitous in nature and, typically, exhibit structural (elastic or phonon) instabilities. Dimpled potentials are not harmonic; thus, the conventional quasiharmonic approximation at finite temperatures fails to describe anharmonic vibrations in such solids. At sufficiently high temperatures, their crystal structure is stabilized by entropy; in this phase, a diffraction pattern of a periodic crystal is combined with vibrational properties of a phonon glass. As temperature is lowered, the solid undergoes a symmetry-breaking transition and transforms into a lower-symmetry phase with lower lattice entropy. Here, we identify specific features in the potential-energy surface that lead to such polymorphic behavior; we establish reliable estimates for the relative energies and temperatures associated with the anharmonic vibrations and the solid–solid symmetry-breaking phase transitions. We show that computational phonon methods can be applied to address anharmonic vibrations in a polymorphic solid at fixed temperature. To illustrate the ubiquity of this class of materials, we present a range of examples (elemental metals, a shape-memory alloy, and a layered charge-density-wave system); we show that our theoretical predictions compare well with known experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparison of methods to estimate aerosol effective radiative forcings in climate models

Abstract. Uncertainty in the effective radiative forcing (ERF) of climate primarily arises from the unknown contribution of aerosols, which impact radiative fluxes directly and through modifying cloud properties. Climate model simulations with fixed sea surface temperatures but perturbed atmospheric aerosol loadings allow for an estimate of how strongly the planet's radiative energy budget has been perturbed by the increase in aerosols since pre-industrial times. The approximate partial radiative perturbation (APRP) technique further decomposes the contributions to the direct forcing due to aerosol scattering and absorption and to the indirect forcing due to aerosol-induced changes in cloud scattering, amount, and absorption, as well as the effects of aerosols on surface albedo. Here we evaluate previously published APRP-derived estimates of aerosol effective radiative forcings from these simulations conducted in the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and find that they are biased as a result of two large coding errors that – in most cases – fortuitously compensate. The most notable exception is the direct radiative forcing from absorbing aerosols, which is more than 40 % larger averaged across CMIP6 models in the present study. Correcting these biases eliminates the residuals and leads to better agreement with benchmark estimates derived from double calls to the radiation code. The APRP method – when properly implemented – remains a highly accurate and efficient technique for diagnosing aerosol ERF in cases where double radiation calls are not available, and in all cases it provides quantification of the individual contributors to the ERF that are highly useful but not otherwise available.

58 GEOSCIENCES↗
NASA NTRS1994

A comparison of GPS broadcast and DMA precise ephemerides

We compare the broadcast ephemerides from Global Positioning Satellites (GPS) to the postprocessed ephemerides from the Defense Mapping Agency (DMA). We find significant energy in the spectrum of the residuals at 1 cycle/day and higher multiples. We estimate the time variance of the residuals and show that the short term residuals, from 15 min, exhibit power law processes with greater low frequency perturbations than white phase modulation. We discuss the significance of these results for the performance of the GPS Kalman filter which estimates the broadcast orbits.

Weiss, Marc A.↗
NASA NTRS2011

Stardust-NExT, Deep Impact, and the Accelerating Spin of 9P/Tempel 1

The evolution of the spin rate of Comet 9P/Tempel 1 through two perihelion passages (in 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13 year as part of a World-Wide observing campaign and from 2888 observations taken over a period of 50 days from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 +/- 0.025deg/dy (41.335 0.005 h) prior to the 2000 perihelion passage, 210.448 +/- 0.016deg/dy (41.055 +/- 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 +/- 0.030deg/dy (40.783 +/- 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 +/- 0.012deg/dy (40.827 +/- 0.002 h) in the interval 2006-2010 following the 2005 perihelion passage. The period decreased by 16.8 +/- 0.3 min during the 2000 passage and by 13.7 +/- 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 +/- 0.003deg/dy/dy at JD2453530.60510 (i.e., 25.134 dy before impact), which provides independent confirmation of the change seen in the Earth-based observations. The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. (Thomas et al. [2007]. Icarus 187, 4-15) pole and longitude system. The possibility of a 180deg error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 at 20:42. We find that a net torque in the range of 0.3-2.5 x 10(exp 7) kg/sq m/sq s acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the Sun. We find that the pole direction is unlikely to precess by more than approx. 1deg per perihelion passage. The trend of the period with time and the fact that the modeled peak torque occurs before perihelion are in agreement with published accounts of trends in water production rate and suggests that widespread H2O out-gassing from the surface is largely responsible for the observed spin-up.

Belton, Michael J. S.↗
NASA NTRS1962

Application of Statistical Filter Theory to the Optimal Estimation of Position and Velocity on Board a Circumlunar Vehicle

Statistical filter theory is employed to develop a method for determining the best possible estimate of the position and velocity of a space vehicle in the midcourse phase of flight. Results of a computer simulation are given to illustrate the performance attainable. An on-board system is visualized in which the source of information is an arbitrary sequence of observations of space angles, corrupted by measurement errors. The scheme is in effect a dynamical time-varying filter, implemented by a digital computer, which processes the incoming data to compute an up-to-date optimal estimate of position and velocity.

Smith, Gerald L.↗
NASA NTRS2000

Analysis of the NASA White Sands Test Facility (WSTF) Test System for Friction-Ignition of Metallic Materials

Friction is a known ignition source for metals in oxygen-enriched atmospheres. The test system developed by the NASA White Sands Test Facility in response to ASTM G-94 has been used successfully to determine the relative ignition from friction of numerous metallic materials and metallic materials pairs. These results have been ranked in terms of a pressure-velocity product (PV) as measured under the prescribed test conditions. A high value of 4.1(exp 8) watts per square meter for Inconel MA 754 is used to imply resistance to friction ignition, whereas a low value of 1.04(exp 8) watts per square meter for stainless steel 304 is taken as indicating material susceptible to friction ignition. No attempt has been made to relate PV values to other material properties. This work reports the analysis of the WSTF friction-ignition test system for producing fundamental properties of metallic materials relating to ignition through friction. Three materials, aluminum, titanium, and nickel were tested in the WSTF frictional ignition instrument system under atmospheres of oxygen or nitrogen. Test conditions were modified to reach a steady state of operation, that is applied, the force was reduced and the rotational speed was reduced. Additional temperature measurements were made on the stator sample. The aluminum immediately galled on contact (reproducible) and the test was stopped. Titanium immediately ignited as a result of non-uniform contact of the stator and rotor. This was reproducible. A portion of the stator sampled burned, but the test continued. Temperature measurements on the stator were used to validate the mathematical model used for estimating the interface (stator/rotor) temperature. These interface temperature measurements and the associate thermal flux into the stator were used to distinguish material-phase transitions, chemical reaction, and mechanical work. The mechanical work was used to analyze surface asperities in the materials and to estimate a coefficient of fiction. The coefficient of fiction was analyzed in terms of material properties that is, hardness, Young's modulus and elasticity/plasticity of the material.

Shoffstall, Michael S.↗
DOE OSTI2020

Excavation Volume Growth Factors at Environmental Remediation Sites - 20318

The U.S. Army Corps of Engineers (USACE) manages numerous environmental remediation projects in accordance with the Comprehensive Environmental Response, Compensation, and Liability Act, otherwise known as CERCLA. These include projects for the Formerly Utilized Sites Remedial Action Program (FUSRAP). Utilizing the CERCLA framework to address these projects requires detailed cost estimates and a significant portion of costs are associated with the volume of material to be excavated and that material's final disposition. Accordingly, understanding the expected volumes and uncertainty associated with those volumes is vital to project cost estimating. These estimates are completed at various stages of the CERCLA process, thus the project team has varying degrees of information available from which to develop these estimates depending on how far along the project is. Understanding the factors contributing to volume growth (increase in volume of material from that originally estimated) is a crucial element to determining volume and cost uncertainty. This paper presents a case study of a FUSRAP site which involves excavation and offsite disposal of contaminated soils from complex commercial, industrial, and residential properties. Original volume estimates were compared to post excavation volumes on a property by property basis. Reasons for differences are determined along with magnitude of impacts to cost estimating. For both residential and commercial property types, volume estimating factors, uncertainty considerations, and methods to account for each are discussed with emphasis on volume growth considerations. Cost estimating model input considerations are discussed with specific emphasis on accounting for volume uncertainty at the Feasibility Study (FS) phase of CERCLA. Underestimating or overestimating costs at the FS phase has repercussions under CERCLA that may impact project close out or completion schedule. Discussions presented will be helpful to both reviewers and preparers of cost estimates involving excavations. Additionally, project planners may be able to use the information to plan for volume growth contingencies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Laboratory for Engineering Man/Machine Systems (LEMS): System identification, model reduction and deconvolution filtering using Fourier based modulating signals and high order statistics

Several important problems in the fields of signal processing and model identification, such as system structure identification, frequency response determination, high order model reduction, high resolution frequency analysis, deconvolution filtering, and etc. Each of these topics involves a wide range of applications and has received considerable attention. Using the Fourier based sinusoidal modulating signals, it is shown that a discrete autoregressive model can be constructed for the least squares identification of continuous systems. Some identification algorithms are presented for both SISO and MIMO systems frequency response determination using only transient data. Also, several new schemes for model reduction were developed. Based upon the complex sinusoidal modulating signals, a parametric least squares algorithm for high resolution frequency estimation is proposed. Numerical examples show that the proposed algorithm gives better performance than the usual. Also, the problem was studied of deconvolution and parameter identification of a general noncausal nonminimum phase ARMA system driven by non-Gaussian stationary random processes. Algorithms are introduced for inverse cumulant estimation, both in the frequency domain via the FFT algorithms and in the domain via the least squares algorithm.

Pan, Jianqiang↗
NASA NTRS1992

Comparison of computer- and human-derived coronary angiographic end-point measures for controlled therapy trials

The Cholesterol Lowering Atherosclerosis Study, a randomized angiographic clinical trial, demonstrated the beneficial effect of niacin/colestipol plus diet therapy on coronary atherosclerosis. Outcome was determined by panel-based estimates (viewed in both still and cine modes) of percent stenosis severity and change in native artery and bypass graft lesions. Computer-based quantitative coronary angiography (QCA) was also used to measure lesion and bypass graft stenosis severity and change in individual frames closely matched in orientation, opacification, and cardiac phase. Both methods jointly evaluated 350 nonoccluded lesions. The correlation between QCA and panel estimates of lesion size was 0.70 (p less than 0.0001) and for change in lesion size was 0.28 (p = 0.002). Agreement between the two methods in classifying lesion changes (i.e., regression, unchanged, or progression) occurred for 60% (210 of 350) of the lesions kappa +/- SEM = 0.20 +/- 0.05, p less than 0.001). The panel identified 442 nonoccluded lesions for which QCA stenosis measurements could not be obtained. Lesions not measurable by QCA included those with stenosis greater than 85% that could not be reliably edge tracked, segments with diffuse or ecstatic disease that had no reliable reference diameter, and segments for which matched frames could not be located. Seventy-nine lesions, the majority between 21% and 40% stenosis, were identified and measured by QCA but were not identified by the panel. This comparison study demonstrates the need to consider available angiographic measurement methods in relation to the goals of their use.

NASA Center JPL↗
DOE OSTI2023

Structure–Transport Relationships of Water–Organic Solvent Co-transport in Carbon Molecular Sieve (CMS) Membranes

We explore the effects of the carbon molecular sieve (CMS) microstructure on the separation performance and transport mechanism of water–organic mixtures. Specifically, we utilize PIM-1 dense films and integrally skinned asymmetric hollow fiber membranes as polymer precursors for the CMS materials. The PIM-1 membranes were pyrolyzed under several different pyrolysis atmospheres (argon, carbon dioxide, and diluted hydrogen gas) and at multiple pyrolysis temperatures. Detailed gas physisorption measurements reveal that membranes pyrolyzed under 4% H 2 and CO 2 had broadened ultramicropore distributions (pore diameter <7 Å) compared to Ar pyrolysis, and pyrolysis under CO 2 increased ultramicropore volume and broadened micropore distributions at increased pyrolysis temperatures. Gravimetric water and p-xylene sorption and diffusion measurements reveal that the PIM-1-derived CMS materials are more hydrophilic than other CMS materials that have been previously studied, which leads to sorption-diffusion estimations showing water-selective permeation. Water permeation in the vapor phase, pervaporation, and liquid-phase hydraulic permeation reveal that the isobaric permeation modes (vapor permeation and pervaporation) are reasonably well predicted by the sorption-diffusion model, whereas the hydraulic permeation mode is significantly underpredicted (>250×). Conversely, the permeation of p-xylene is well predicted by the sorption-diffusion model in all cases. The collection of pore size analysis, vapor sorption and diffusion, and permeation in different modalities creates a picture of a combined transport mechanism in which water–under high transmembrane pressures–permeates via a Poiseuille-style mechanism, whereas p-xylene solutes in the mixture permeate via sorption-diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
DOE OSTI2022

Ambient Noise Attenuation Tomography Reveals an Asymmetric Damage Zone Across San Jacinto Fault Near Anza, California

Here, we perform seismic attenuation tomography of the shallow structure for the San Jacinto Fault in the Ramona Reservation of southern California. The study uses ambient seismic noise recorded by a linear array of 65 3-C sensors across the fault. We extract amplitude decay information, with uncertainty, from noise interferometry functions. To account for strong heterogeneities in the complex shallow fault zone structure, we apply a frequency-dependent amplitude correction for focusing/defocusing effects using inverted phase velocity maps obtained by solving the transport equation. We then estimate the attenuation structure based on the linear station-triplet method for both Love and Rayleigh waves. The attenuation tomography indicates strong attenuation correlated with known San Jacinto Fault surface traces. The Love wave attenuation tomography reveals an asymmetric damage zone that exists primarily on the fault side with faster seismic velocity, consistent with earthquake ruptures on a fault bimaterial interface with preferred propagation direction to the northwest.

58 GEOSCIENCES↗
DOE OSTI2024

Edgewise Structural Damping of a 2.8-MW Land-Based Wind Turbine Rotor Blade

Modern wind turbines push the predictive capabilities of state-of-the-art aero-servo-elastic tools. The existing limitations hide across the numerical tool chain and can result in serious issues, such as missing aeroelastic instabilities during the design phase. Structural damping is an input that is usually hard to estimate, but also has a major impact on the turbine behavior. In this paper, we discuss an experiment that aims to accurately quantify the structural damping characterizing the edgewise modes of modern wind turbine blades. The experiment is carried out on a 2.8-MW land-based wind turbine and features a fast yaw actuation that induces an edgewise motion on one of the three blades. The Covariant-subspace system identification (Cov-SSI) method is then used to post-process the blade root moment to estimate the short-term edgewise structural damping. Despite limitations of the Cov-SSI method, which consistently under-predicts the absolute values of damping, we observe that structural damping decreases across the first three blade edgewise modes, which is different from the stiffness-proportional damping model that assumes that structural damping increases with the modes. This paper argues that a stiffness-proportional damping model, which is implemented in most aeroelastic tools, is therefore not conservative and might hide aeroelastic instabilities that can instead appear in the field.

17 WIND ENERGY↗
DOE OSTI2021

Model-based Cyber-attack Detection for Voltage Source Converters in Island Microgrids

With the upgrading of communication and control in microgrids, cyber threats on the power converters are increasing. Here, in this paper, a novel model-based cyber-attack detection methodology is proposed for each voltage source converter in microgrids. The Harmonics State Space Matrix (H-Matrix) is used to build the closed-loop transfer function, providing model-based estimation with the grid voltage and control reference. Then, the space phase model (SPM) is used to calculate residual to detect cyber-attacks in voltage source converter (VSC). As controller and converter parameters are included in the H-matrix, the proposed method also could detect cyber-attacks when voltage fluctuation occurs in the grid. To verify the feasibility of the proposed method, several attacks targeting VSC in an islanded microgrid are simulated in MATLAB. Simulation and detection results are introduced to demonstrate the resilience and feasibility of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Bayesian Framework for Multi-Timescale State Estimation in Low-Observable Distribution Systems

To support the smart grid paradigm, there has been a significant increase in sensor deployments and metering infrastructure in distribution systems. However, the measurements provided by these sensors and metering devices are typically sampled at different rates and could suffer from losses during the aggregation process. It is crucial to effectively reconcile the time-series measurements for a reliable state estimation. While weighted least squares has been the traditional approach for state estimation, sparsity-based approaches like matrix completion have become popular due to their superior performance in low-observability conditions. This paper proposes a Bayesian framework for both multi-timescale data aggregation and matrix completion based state estimation. Specifically, the multiscale time-series data aggregated from heterogenous sources are reconciled using a multitask Gaussian process that exploits the spatio-temporal correlations. Here, the resulting consistent timeseries alongwith the confidence bound on the imputations are fed into a Bayesian matrix completion method augmented with linearized power-flow constraints to accurately estimate the states in low-observability conditions. Results on three phase unbalanced IEEE 37 and IEEE 123 bus test systems reveal the superior performance of the proposed Bayesian framework. The computational complexity for the proposed Bayesian framework is also quantified.

42 ENGINEERING↗