Towards a unified understanding of event-related changes in the EEG:the Firefly model of synchronization through cross-frequency phase modulation

Although event-related potentials (ERPs) are widely used to study sensory, perceptual and cognitive processes, it remains unknown whether they are phase-locked signals superimposed upon the ongoing electroencephalogram (EEG) or result from phase-alignment of the EEG. Previous attempts to discriminate between these hypotheses have been unsuccessful but here a new test is presented based on the prediction that ERPs generated by phase-alignment will be associated with event-related changes in frequency whereas evoked-ERPs will not. Using empirical mode decomposition (EMD), which allows measurement of narrow-band changes in the EEG without predefining frequency bands, evidence was found for transient frequency slowing in recognition memory ERPs but not in simulated data derived from the evoked model. Furthermore, the timing of phase-alignment was frequency dependent with the earliest alignment occurring at high frequencies. Based on these findings, the Firefly model was developed, which proposes that both evoked and induced power changes derive from frequency-dependent phase-alignment of the ongoing EEG. Simulated data derived from the Firefly model provided a close match with empirical data and the model was able to account for i) the shape and timing of ERPs at different scalp sites, ii) the event-related desynchronization in alpha and synchronization in theta, and iii) changes in the power density spectrum from the pre-stimulus baseline to the post-stimulus period. The Firefly Model, therefore, provides not only a unifying account of event-related changes in the EEG but also a possible mechanism for cross-frequency information processing.

Atenuação de ruídos coerentes utilizando decomposição em modos empíricos

The seismic processing technique has the main objective to provide adequate picture of geological structures from subsurface of sedimentary basins. Among the key steps of this process is the enhancement of seismic reflections by filtering unwanted signals, called seismic noise, the improvement of signals of interest and the application of imaging procedures. The seismic noise may appear random or coherent. This dissertation will present a technique to attenuate coherent noise, such as ground roll and multiple reflections, based on Empirical Mode Decomposition method. This method will be applied to decompose the seismic trace into Intrinsic Mode Functions. These functions have the properties of being symmetric, with local mean equals zero and the same number of zero-crossing and extremes. The developed technique was tested on synthetic and real data, and the results were considered encouraging

Advancements in jet turbulence and noise modeling: accurate one-way solutions and empirical evaluation of the nonlinear forcing of wavepackets

Jet noise reduction is an important goal within both commercial and military aviation. Although large-scale numerical simulations are now able to simultaneously compute turbulent jets and their radiated sound, lost-cost, physically-motivated models are needed to guide noise-reduction efforts. A particularly promising modeling approach centers around certain large-scale coherent structures, called wavepackets, that are observed in jets and their radiated sound. The typical approach to modeling wavepackets is to approximate them as linear modal solutions of the Euler or Navier-Stokes equations linearized about the long-time mean of the turbulent flow field. The near-field wavepackets obtained from these models show compelling agreement with those educed from experimental and simulation data for both subsonic and supersonic jets, but the acoustic radiation is severely under-predicted in the subsonic case. This thesis contributes to two aspects of these models. First, two new solution methods are developed that can be used to efficiently compute wavepackets and their acoustic radiation, reducing the computational cost of the model by more than an order of magnitude. The new techniques are spatial integration methods and constitute a well-posed, convergent alternative to the frequently used parabolized stability equations. Using concepts related to well-posed boundary conditions, the methods are formulated for general hyperbolic equations and thus have potential applications in many fields of physics and engineering. Second, the nonlinear and stochastic forcing of wavepackets is investigated with the goal of identifying and characterizing the missing dynamics responsible for the under-prediction of acoustic radiation by linear wavepacket models for subsonic jets. Specifically, we use ensembles of large-eddy-simulation flow and force data along with two data decomposition techniques to educe the actual nonlinear forcing experienced by wavepackets in a Mach 0.9 turbulent jet. Modes with high energy are extracted using proper orthogonal decomposition, while high gain modes are identified using a novel technique called empirical resolvent-mode decomposition. In contrast to the flow and acoustic fields, the forcing field is characterized by a lack of energetic coherent structures. Furthermore, the structures that do exist are largely uncorrelated with the acoustic field. Instead, the forces that most efficiently excite an acoustic response appear to take the form of random turbulent fluctuations, implying that direct feedback from nonlinear interactions amongst wavepackets is not an essential noise source mechanism. This suggests that the essential ingredients of sound generation in high Reynolds number jets are contained within the linearized Navier-Stokes operator rather than in the nonlinear forcing terms, a conclusion that has important implications for jet noise modeling.

Fadiga em compósitos 5HS carbono/epóxi processados via RTM: relação entre cargas axiais no plano e modos de delaminação

Pós-graduação em Engenharia Mecânica - FEG

Twist phase in Gaussian-beam optics

The recently discovered twist phase is studied in the context of the full ten-parameter family of partially coherent general anisotropic Gaussian Schell-model beams. It is shown that the nonnegativity requirement on the cross-spectral density of the beam demands that the strength of the twist phase be bounded from above by the inverse of the transverse coherence area of the beam. The twist phase as a two-point function is shown to have the structure of the generalized Huygens kernel or Green's function of a first-order system. The ray-transfer matrix of this system is exhibited. Wolf-type coherent-mode decomposition of the twist phase is carried out. Imposition of the twist phase on an otherwise untwisted beam is shown to result in a linear transformation in the ray phase space of the Wigner distribution. Though this transformation preserves the four-dimensional phase-space volume, it is not symplectic and hence it can, when impressed on a Wigner distribution, push it out of the convex set of all bona fide Wigner distributions unless the original Wigner distribution was sufficiently deep into the interior of the set.

Interdependence between Greece and other European stock markets: A comparison of wavelet and VMD copula, and the portfolio implications

The interdependence of Greece and other European stock markets and the subsequent portfolio implications are examined in wavelet and variational mode decomposition domain. In applying the decomposition techniques, we analyze the structural properties of data and distinguish between short and long term dynamics of stock market returns. First, the GARCH-type models are fitted to obtain the standardized residuals. Next, different copula functions are evaluated, and based on the conventional information criteria and time varying parameter, Joe-Clayton copula is chosen to model the tail dependence between the stock markets. The short-run lower tail dependence time paths show a sudden increase in comovement during the global financial crises. The results of the long-run dependence suggest that European stock markets have higher interdependence with Greece stock market. Individual country’s Value at Risk (VaR) separates the countries into two distinct groups. Finally, the two-asset portfolio VaR measures provide potential markets for Greece stock market investment diversification.

Dynamic Stall on Vertical Axis Wind Turbines

In this study the dynamics of flow over the blades of vertical axis wind turbines was investigated using a simplified periodic motion to uncover the fundamental flow physics and provide insight into the design of more efficient turbines. Time-resolved, two-dimensional velocity measurements were made with particle image velocimetry on a wing undergoing pitching and surging motion to mimic the flow on a turbine blade in a non-rotating frame. Dynamic stall prior to maximum angle of attack and a leading edge vortex development were identified in the phase-averaged flow field and captured by a simple model with five modes, including the first two harmonics of the pitch/surge frequency identified using the dynamic mode decomposition. Analysis of these modes identified vortical structures corresponding to both frequencies that led the separation and reattachment processes, while their phase relationship determined the evolution of the flow.

Detailed analysis of the leading edge vortex found multiple regimes of vortex development coupled to the time-varying flow field on the airfoil. The vortex was shown to grow on the airfoil for four convection times, before shedding and causing dynamic stall in agreement with 'optimal' vortex formation theory. Vortex shedding from the trailing edge was identified from instantaneous velocity fields prior to separation. This shedding was found to be in agreement with classical Strouhal frequency scaling and was removed by phase averaging, which indicates that it is not exactly coupled to the phase of the airfoil motion.

The flow field over an airfoil undergoing solely pitch motion was shown to develop similarly to the pitch/surge motion; however, flow separation took place earlier, corresponding to the earlier formation of the leading edge vortex. A similar reduced-order model to the pitch/surge case was developed, with similar vortical structures leading separation and reattachment; however, the relative phase lead of the separation mode, corresponding to earlier separation, necessitated that a third frequency to be incorporated into the reattachment mode to provide a relative lag in reattachment.

Finally, the results are returned to the rotating frame and the effects of each flow phenomena on the turbine are estimated, suggesting kinematic criteria for the design of improved turbines.

Effect of low-frequency Rossby wave on thermal structure of the upper southwestern tropical Indian Ocean

We investigate the influence of low-frequency Rossby waves on the thermal structure of the upper southwestern tropical Indian Ocean (SWTIO) using Argo profiles, satellite altimetric data, sea surface temperature, wind field data and the theory of linear vertical normal mode decomposition. Our results show that the SWTIO is generally dominated by the first baroclinic mode motion. As strong downwelling Rossby waves reach the SWTIO, the contribution of the second baroclinic mode motion in this region can be increased mainly because of the reduction in the vertical stratification of the upper layer above thermocline, and the enhancement in the vertical stratification of the lower layer under thermocline also contributes to it. The vertical displacement of each isothermal is enlarged and the thermal structure of the upper level is modulated, which is indicative of strong vertical mixing. However, the cold Rossby waves increase the vertical stratification of the upper level, restricting the variability related to the second baroclinic mode. On the other hand, during decaying phase of warm Rossby waves, Ekman upwelling and advection processes associated with the surface cyclonic wind circulation can restrain the downwelling processes, carrying the relatively colder water to the near-surface, which results in an out-of-phase phenomenon between sea surface temperature anomaly (SSTA) and sea surface height anomaly (SSHA) in the SWTIO.

Equivalent-barotropic definition of tropospheric mean temperature

An equivalent-barotropic (EB) description of the tropospheric temperature field is derived from the geostrophic empirical mode (GEM) in the form of a scalar function Gamma(p, phi), where p is pressure and phi is 300-850-mb thickness. Baroclinic parameter phi plays the role of latitude at each longitudinal section. Compared with traditional Eulerian-mean methods, GEM defines a mean field in baroclinic streamfunction space with a time scale much longer than synoptic variability. It prompts an EB concept that is only based on a baroclinic field. Monthly GEM fields are diagnosed from NCEP-NCAR reanalysis data and account for more than 90% of the tropospheric thermal variance. The circumglobal composite of GEM fields exhibits seasonal, zonal, and hemispheric asymmetries, with larger rms errors occurring in winter and in the Northern Hemisphere (NH). Zonally asymmetric features and planetary deviation from EB are seen in the NH winter GEM. Reconstruction of synoptic sections and correlation analysis reveal that the tropospheric temperature field is EB at the leading order and has a 1-day phase lag behind barotropic variations in extratropical regions.

PN断面东海黑潮的地转经验模态分析

基于东海PN断面近五十年的历史水文数据，本文应用一种斜压流函数空间投影的方法分析冲绳海槽北部黑潮水团的水文结构及其季节性和年代际变异。通过将历史数据投影到一个斜压流函数空间中，GEM方法能很好地消除黑潮中尺度变异引起的噪音。诊断生成的东海黑潮温度和盐度GEM场以压力和斜压流函数为参数，能反映黑潮基本的温盐结构。通过分析1987年之后黑潮季节平均的温盐GEM场，结果显示黑潮中层水的盐度在夏季相对较低, 在春季达到最高。通过比较1987年之后温盐深仪（CTD）数据生成的GEM场与1987年之前南森瓶（Bottle）数据生成的GEM场，我们证实黑潮的年代际变化特征，即黑潮表层水近二十年来温度升高而盐度减小，同时发现黑潮中层水具有盐度减小而黑潮热带水具有盐度增加的趋势。进一步通过对PN断面上黑潮水团来源的分析讨论，得出中层水盐度最小值靠近黑潮左侧的现象是由于冲绳岛-宫古岛之间入侵的北太平洋中层水受黑潮主轴的平流作用所致。

反假频重建理论及其在储层预测中的应用

To extend the cross-hole seismic 2D data to outside 3D seismic data, reconstructing the low frequency data to high frequency data is necessary. Blind deconvolution method is a key technology. In this paper, an implementation of Blind deconvolution is introduced. And optimized precondition conjugate gradient method is used to improve the stability of the algorithm and reduce the computation. Then high-frequency retrieved Seismic data and the cross-hole seismic data is combined for constraint inversion. Real data processing proved the method is effective. To solve the problem that the seismic data resolution can’t meet the request of reservoir prediction in the river face thin-layers in Chinese eastern oil fields, a high frequency data reconstruction method is proposed. The extrema of the seismic data are used to get the modulation function which operated with the original seismic data to get the high frequency part of the reconstruction data to rebuild the wide band data. This method greatly saves the computation, and easy to adjust the parameters. In the output profile, the original features of the seismic events are kept, the common feint that breaking the events and adding new zeros to produce alias is avoided. And the interbeded details are enhanced compared to the original profiles. The effective band of seismic data is expended and the method is approved by the processing of the field data. Aim to the problem in the exploration and development of Chinese eastern oil field that the high frequency log data and the relative low frequency seismic data can’t be merged, a workflow of log data extrapolation constrained by time-phase model based on local wave decomposition is raised. The seismic instantaneous phase is resolved by local wave decomposition to build time-phase model, the layers beside the well is matched to build the relation of log and seismic data, multiple log info is extrapolated constrained by seismic equiphase map, high precision attributes inverse sections are produced. In the course of resolve the instantaneous phase, a new method of local wave decomposition --Hilbert transform mean mode decomposition(HMMD) is raised to improve the computation speed and noise immunity. The method is applied in the high resolution reservoir prediction in Mao2 survey of Daqing oil field, Multiple attributes profiles of wave impedance, gamma-ray, electrical resistivity, sand membership degree are produced, of which the resolution is high and the horizontal continuous is good. It’s proved to be a effective method for reservoir prediction and estimation.

Power system parameters forecasting using Hilbert-Huang transform and machine learning

A novel hybrid data-driven approach is developed for forecasting power system parameters with the goal of increasing the efficiency of short-term forecasting studies for non-stationary time-series. The proposed approach is based on mode decomposition and a feature analysis of initial retrospective data using the Hilbert-Huang transform and machine learning algorithms. The random forests and gradient boosting trees learning techniques were examined. The decision tree techniques were used to rank the importance of variables employed in the forecasting models. The Mean Decrease Gini index is employed as an impurity function. The resulting hybrid forecasting models employ the radial basis function neural network and support vector regression. A part from introduction and references the paper is organized as follows. The second section presents the background and the review of several approaches for short-term forecasting of power system parameters. In the third section a hybrid machine learningbased algorithm using Hilbert-Huang transform is developed for short-term forecasting of power system parameters. Fourth section describes the decision tree learning algorithms used for the issue of variables importance. Finally in section six the experimental results in the following electric power problems are presented: active power flow forecasting, electricity price forecasting and for the wind speed and direction forecasting.

Global atmospheric energetics under present and future climate conditions

An evaluation of the global atmospheric energetics is presented in the framework of the basic decomposition into the zonal mean and eddy components, the zonal wavenumber decomposition, and the three-dimensional normal mode decomposition. An extension to the normal mode energetics formulation is also presented in the study, which enables the explicit evaluation of the conversion rate between available potential energy and kinetic energy along with their generation and dissipation rates, in both the zonal wavenumber and vertical mode domains. In addition, it has been proposed an extended energy cycle diagram describing the flow of energy among the zonal mean and eddy components, and also among the barotropic and baroclinic components. The energetics is first assessed for three reanalysis datasets and five state-ofthe- art climate models simulations representing the present climate conditions. It is performed a comparative analysis between the observationally based energetics and that based on the climate models' simulations. In order to appraise possible changes in the atmospheric energetics of a future climate scenario relative to that of the present climate conditions, the analysis is extended using the datasets simulated by the same five climate models for a future climate scenario experiment, as defined in the Special Report on Emissions Scenarios (SRES) of the Intergovernmental Panel on Climate Change (IPCC).

What causes violent crime?

This study uses a new data set of crime ratesfor a large sample of countriesfor the period 1970- 1994, based on information from the United Nations World Crime Surveys, to ana/yze the determinants ofnational homicide and robbery rates. A simple model of the incentives to commit crimes is proposed, which explicit/y considers possible causes of the persistence of crime over time (criminal inertia). Several econometric mode/s are estimated, attempting to capture the . determinonts of crime rates across countries and over time. The empirical mode/s are first run for cross-sections and then applie'd to panel data. The former focus on erplanatory variables that do not change markedly over time, while the panel data techniques consider both the eflect of the business cyc1e (i.e., GDP growth rate) on the crime rate and criminal inertia (accountedfor by the inclusion of the /agged crime rate as an explanatory variable). The panel data techniques a/so consider country-specific eflects, the joint endogeneity of some of the erplanatory variables, and lhe existence of some types of measurement e"ors aJjlicting the crime data. The results showthat increases in income inequality raise crime rates, dete"ence eflects are significant, crime tends to be counter-cyclical, and criminal inertia is significant even after controlling for other potential determinants of homicide and robbery rates.

Four Decades of Studying Global Linear Instability: Progress and Challenges

Global linear instability theory is concerned with the temporal or spatial development of small-amplitude perturbations superposed upon laminar steady or time-periodic three-dimensional flows, which are inhomogeneous in two(and periodic in one)or all three spatial directions.After a brief exposition of the theory,some recent advances are reported. First, results are presented on the implementation of a Jacobian-free Newton–Krylov time-stepping method into a standard finite-volume aerodynamic code to obtain global linear instability results in flows of industrial interest. Second, connections are sought between established and more-modern approaches for structure identification in flows, such as proper orthogonal decomposition and Koopman modes analysis (dynamic mode decomposition), and the possibility to connect solutions of the eigenvalue problem obtained by matrix formation or time-stepping with those delivered by dynamic mode decomposition, residual algorithm, and proper orthogonal decomposition analysis is highlighted in the laminar regime; turbulent and three-dimensional flows are identified as open areas for future research. Finally, a new stable very-high-order finite-difference method is implemented for the spatial discretization of the operators describing the spatial biglobal eigenvalue problem, parabolized stability equation three-dimensional analysis, and the triglobal eigenvalue problem; it is shown that, combined with sparse matrix treatment, all these problems may now be solved on standard desktop computers