• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.66 no.2
• /
• pp.261-284
• /
• 2017

This paper describes the results of a low-frequency oscillation analysis using data measured in PMU installed in the KEPCO system, and the comparison with eigenvalues computed from the linear model. The dominant oscillation modes are estimated by applying various algorithms. The algorithms are: the extended Prony method; multiple time interval parameter estimation method; subspace system identification method; and spectral analysis. From the measurement data, modes of frequency 0.68[Hz] and 0.92[Hz] were estimated, and modes of frequency 0.63[Hz] and 0.80[Hz] were computed from the eigenvalue calculation. There was a difference between the mode estimated from measurement data and that from the linear model. This is possibly because of an error in the dynamic data of the KEPCO system used in eigenvalue calculation. Because wide area modes exist in the KEPCO system, these modes should be monitored continuously for the reliable operation of the system. In order to prevent total blackouts caused by wide area oscillation, moreover, contingency analysis should be performed in relation to this mode and appropriate measures should be established.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.7 no.1
• /
• pp.108-116
• /
• 1997

A laboratory experiment was made of critical transition flow modes in Czocllralski convection. Numerical computation was also made to delineate the dynamic transition. The period of temperature oscillation ($t_p$) and the interval of temperature oscillation ($\Delta\theta$) were scrutinized to capture the critical transition regime. The mixed convection parameter was varied in the range of $0.134\le Ra/PrRe^2 \le3.804$. The data from calculation were in good agreement with ones from experiment. The influence of the Prandtl number on the transition was examined for Pr ＝ 910, 4445 and 8889. To understand the transition mechanism, the detailed temperature oscillation modes, the isolines of meridional temperature and the axial velocity profiles were investigated.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.2
• /
• pp.401-409
• /
• 2010

In this paper, we introduce an iterative method to study oscillatory properties of delay difference equations of the following form ${\nabla}_{\alpha}\;[x(t)\;-\;r(t)x(t\;-\;k)]\;+\;p(t)x(t\;-\;{\tau})\;-\;q(t)x(t\;-\;{\sigma})\;=\;0$, $t\;{\geq}\;t_0$, where $t_0\;{\in}\;\mathbb{R}$, t varies in the real interval ($t_0,\;{\infty}$), $\alpha$ > 0, $\kappa$, $\tau$, ${\sigma}\;{\geq}\;0$, $r\;{\in}\;C\;([t_0-{\alpha},\;{\infty}),\;\mathbb{R}^+$, p, $q\;{\in}\;C\;([t_0,\;{\infty}),\;\mathbb{R}^+)$ and ${\nabla}_{\alpha}x(t)\;=\;x(t)\;-\;x(t\;-\;{\alpha})$ for $t\;{\geq}\;t_0$.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.54 no.6
• /
• pp.293-305
• /
• 2005

This paper describes an improved Prony method in its speed, accuracy and reliability by efficiently determining the optimal sampling interval with use of DFT (discrete Fourier transformation). In the Prony method the computation time is dominated by the size of the linear prediction matrix, which is given by the number of data times the modeling order The size of the matrix in a general Prony method becomes large because of large number of data and so does the computation time. It is found that the Prony method produces satisfactory results when SNR is greater than three. The maximum sampling interval resulting minimum computation time is determined using the fact that the spectrum in DFT is inversely proportional to sampling interval. Also the process of computing the modes is made efficient by applying Hessenberg method to the companion matrix with complex shift and computing selectively only the dominant modes of interest. The proposed method is tested against the 2003 KEPCO system and found to be efficient and reliable. The proposed method may play a key role in monitoring in real time low frequency oscillations of power systems .

• Journal of Electrical Engineering and Technology
• /
• v.9 no.1
• /
• pp.27-36
• /
• 2014

This paper presents an adaptive wide-area damping controller (WADC) based on generalized predictive control (GPC) and model identification for damping the inter-area low frequency oscillations in large-scale inter-connected power system. A recursive least-squares algorithm (RLSA) with a varying forgetting factor is applied to identify online the reduced-order linearlized model which contains dominant inter-area low frequency oscillations. Based on this linearlized model, the generalized predictive control scheme considering control output constraints is employed to obtain the optimal control signal in each sampling interval. Case studies are undertaken on a two-area four-machine power system and the New England 10-machine 39-bus power system, respectively. Simulation results show that the proposed adaptive WADC not only can damp the inter-area oscillations effectively under a wide range of operation conditions and different disturbances, but also has better robustness against to the time delay existing in the remote signals. The comparison studies with the conventional lead-lag WADC are also provided.

• Journal of Power Electronics
• /
• v.19 no.5
• /
• pp.1270-1277
• /
• 2019

A problem with virtual synchronous generator (VSG) systems is that they are difficult to operate stably with photovoltaic (PV) power as the DC side. With this problem in mind, a PV-VSG control strategy considering the dynamic characteristics of the DC side is proposed after an in-depth analysis of the dynamic characteristics of photovoltaic power with a parallel energy-storage capacitor. The proposed PV-VSG automatically introduces DC side voltage control for the VSG when the PV enters into an unstable working interval, which avoids the phenomenon where an inverter fails to work due to a DC voltage sag. The stability of the original VSG and the proposed PV-VSG were compared by a root locus analysis. It is found that the stability of the PV-VSG is more sensitive to the inertia coefficient J than the VSG, and that a serious power oscillation may occur. According to this, a new rotor model is designed to make the inertial coefficient automatically change to adapt to the operating state. Experimental results show that the PV-VSG control strategy can achieve stable operation and maximum power output when the PV output power is insufficient.

• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.245-253
• /
• 1995

We have analyzed the time series of Ca II H,K and ${\lambda}8498$ line profiles taken for a sunspot (SPO 5007) with the Echelle spectrograph attached to Vacuum Tower Telescope at Sacramento Peak Solar Observatory. Each set of spectra was taken simultaneously for 20 minutes at a time interval of 30 seconds. A total of 40 photographic films for each line was scanned by a PDS at Korea Astronomy Observatory. The central peak intensity of Ca II H ($I_{max}$), the intensity measured at ${\Delta}{\lambda}=-0.1{\AA}$ from the line center of ${\lambda}8498(I_{{\lambda}8489})$, the radial velocity ($V_r$) and the Doppler width (${\Delta}{\lambda}_D$) estimated from Ca II H have been measured to study the dynamical behaviors of the sunspot chromosphere. Fourier analysis has been carried out for these measured quantities. Our main results are as follows: (1) We have confirmed the 3-minute oscillation being dominant throughout the umbra. The period of oscillations jumps from 180 sec in the umbra to 500 to 1000 sec in the penumbra. (2) The nonlinear character of the umbral oscillation is noted from the observed sawtooth shaped radial velocity fluctuations with amplitudes reaching up to $5{\sim}6\;km/sec$. (3) The spatial distribution of the maximum powers shows that the power of oscillations is stronger in the umbra than in the penumbra. (4) The spatial distributions of the time averaged < $I_{max}$ > and < $V_r$ > across the spot are found to be nearly axially symmetric, implying that the physical quantities derived from the line profiles of Ca II H and ${\lambda}8498$ are inherently associated with the geometry of the magnetic field distribution of the spot. (5) The central peaks of the CaII H emission core lead the upward motions of the umbral atmosphere by $90^{\circ}$, while no phase delay is found in intensities between $I_{max}$ and $I_{{\lambda}8498}$, suggesting that the umbral oscillation is of standing waves.

• Korean Journal of Food Science and Technology
• /
• v.28 no.4
• /
• pp.624-631
• /
• 1996

Effect of fluctuation range and intervals of defrosting temperature on quality of frozen foods stored in a domestic refrigerator equipped with an automatic defrost system was evaluated. As defrost system was operated, temperatures of domestic refrigerators were elevated from $-18^{\circ}C\;to\;-5^{\circ}C\;and\;-15^{\circ}C$, and fluctuation intervals were l6 hrs and 30 hrs, respectively. Quality deterioration such as protein denaturation, vitamin loss, exudate production and changes in appearance of frozen foods was minimized by reducing temperature oscillation during storage. Considerable effects of thermal oscillation on ice crystal sizes were observed for frozen beef tissue and ice cream. TTI (time temperature indicator) system also proved that the temperature control of defrost system in domestic refrigerator can improve the quality of foods during storage.

• Journal of Biomedical Engineering Research
• /
• v.28 no.4
• /
• pp.503-511
• /
• 2007

Heart rate variability(HRV) is the clinical consequence of various influences of the autonomic nervous system(ANS) on heart beat. HRV can estimate the potential physiologic rhythm from the interval between consecutive beats(RR interval or HRV data), but cardiovascular system governed by ANS is in relation to respiration and autonomic regulation. It is known as RSA representing respiration-related HR rhythmic oscillation. Because the mechanism linking the variability of HR to respiration is complex, it has so far been unknown well. In this paper, we tried to evaluate 5-min RR interval segments under control of respiration in order to find out a proper respiration rate that can estimate the ANS function. 10 healthy volunteers were included to evaluate 5-min HRV data under 4 different respiration-controlled environments; 0.03Hz, 0.1Hz, 0.2Hz, and 0.4Hz respiration. HRV data were analyzed both in the frequency and the time domain, with cross-correlation coefficient(cross-coeff.) for HRV and respiration signal. The results showed maximum cross-coeff. of 0.84 at 0.1 Hz and minimum that of 0.16 at 0.4Hz respiration. Cross-coeff was decreased at a faster rate from 0.1Hz respiration. All mean SDNN, RMSSD, and pNN50 of time domain measures were 108.7ms, 71.85ms, and 28.47%, respectively, and LF, HF, and TP of frequency domain measures were $12,722ms^2,\;658.8ms^2$, and $7,836.64ms^2$ at 0.1Hz respiration, respectively. In conclusion, 0.1Hz respiration was observed to be very meaningful from time domain and frequency domain analysis in relation to respiration and autonomic regulation of the heart.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2000.06a
• /
• pp.166-172
• /
• 2000

Friction characteristics of phenolic resin-based friction composites containing threedifferent relative amounts of graphite and zirconium silicate were investigated by using a pad-on-disk type friction tester. Constant temperature test and constant interval test at three different initial temperatures(100. 200, 300$^{\circ}C$) were performed to examine the effects of friction heat on friction characteristics at elevated temperature. The friction composite(FMO.7) with higher content of ZrSiO$_4$showed unstable friction force at higher temperature and resulted in larger fluctuations of vibration during friction test. The abrasive action of ZrSiO$_4$in friction composite impeded stable transfer film and induced higher friction heat at friction interface. Friction oscillations according to the temperature were associated with the formation of transfer film(i'd body layer) on the friction composite and the counter part.