Numerical Simulation of Stochastic Di erential Equations

Stochastic differential equation (SDE) is a differential equation in which some of the terms and its solution are stochastic processes. SDEs play a central role in modeling physical systems like finance, Biology, Engineering, to mention some. In modeling process, the computation of the trajectories (sample paths) of solutions to SDEs is very important. However, the exact solution to a SDE is generally difficult to obtain due to non-differentiability character of realizations of the Brownian motion. There exist approximation methods of solutions of SDE. The solutions will be continuous stochastic processes that represent diffusive dynamics, a common modeling assumption for financial, Biology, physical, environmental systems. This Masters' thesis is an introduction and survey of numerical solution methods for stochastic differential equations. Standard numerical methods, local linearization methods and filtering methods are well described. We compute the root mean square errors for each method from which we propose a better numerical scheme. Stochastic differential equations can be formulated from a given ordinary differential equations. In this thesis, we describe two kind of formulations: parametric and non-parametric techniques. The formulation is based on epidemiological SEIR model. This methods have a tendency of increasing parameters in the constructed SDEs, hence, it requires more data. We compare the two techniques numerically.

Conventional and nonconventional Kramers-Kronig analysis in optical spectroscopy

This Master’s Thesis is dedicated to the investigation and testing conventional and nonconventional Kramers-Kronig relations on simulated and experimentally measured spectra. It is done for both linear and nonlinear optical spectral data. Big part of attention is paid to the new method of obtaining complex refractive index from a transmittance spectrum without direct information of the sample thickness. The latter method is coupled with terahertz tome-domain spectroscopy and Kramers-Kronig analysis applied for testing the validity of complex refractive index. In this research precision of data inversion is evaluated by root-mean square error. Testing of methods is made over different spectral range and implementation of this methods in future is considered.

Electrical status epilepticus during sleep. Continuous spikes and waves during sleep

ABSTRACT Maria Peltola Electrical status epilepticus during sleep – Continuous spikes and waves during sleep Department of Clinical Neurophysiology, University of Turku Department of Clinical Neurophysiology and Department of Pediatric Neurology, Children’s Hospital, Helsinki University Central Hospital Annales Universitatis Turkuensis, Medica-Odontologica, Turku, Finland, 2014 Background: Electrical status epilepticus during sleep (ESES) is an EEG phenomenon of frequent spikes and waves occurring in slow sleep. ESES relates to cognitive deterioration in heterogeneous childhood epilepsies. Validated methods to quantitate ESES are missing. The clinical syndrome, called epileptic encephalopathy with continuous spikes and waves during sleep (CSWS) is pharmacoresistant in half of the patients. Limited data exists on surgical treatment of CSWS. Aims and methods: The effects of surgical treatment were studied by investigating electroclinical outcomes in 13 operated patients (nine callosotomies, four resections) with pharmacoresistant CSWS and cognitive decline. Secondly, an objective paradigm was searched for assessing ESES by the semiautomatic quantification of spike index (SI) and measuring spike strength from EEG. Results: Postoperatively, cognitive deterioration was stopped in 12 (92%) patients. Three out of four patients became seizure-free after resective surgery. Callosotomy resulted in greater than 90% reduction of atypical absences in six out of eight patients. The preoperative propagation of ESES from one hemisphere to the other was associated with a good response. Semiautomatic quantification of SI was a robust method when the maximal interspike interval of three seconds was used to determine the “continuous” discharge in ten EEGs. SI of the first hour of sleep appeared representative of the whole night SI. Furthermore, the spikes’ root mean square was found to be a stable measure of spike strength when spatially integrated over multiple electrodes during steady NREM sleep. Conclusions: Patients with pharmacoresistant CSWS, based on structural etiology, may benefit from resective surgery or corpus callosotomy regarding both seizure outcome and cognitive prognosis. The semiautomated SI quantification, with proper userdefined settings and the new spatially integrated measure of spike strength, are robust and promising tools for quantifying ESES. Keywords: Electrical status epilepticus during sleep, ESES, continuous spikes and waves during sleep, CSWS, epilepsy surgery, spike index, spike strength, RMS TIIVISTELMÄ Maria Peltola Unenaikainen sähköinen status epilepticus Kliininen neurofysiologia, Turun yliopisto Kliininen neurofysiologia ja lastenneurologia, Lasten ja nuorten sairaala, Helsingin yliopistollinen keskussairaala Annales Universitatis Turkuensis, Medica-Odontologica, Turku, Suomi, 2014 Tausta: Sähköinen status epilepticus unessa (ESES) on aivosähkökäyrä (EEG)-ilmiö, jossa hidasaaltounen aikana esiintyy tiheä piikkihidasaaltopurkaus. ESES:n kvantifioimiseen ei ole olemassa validoituja menetelmiä. ESES on liitetty kognitiivisen tason laskuun ja tällöin puhutaan CSWS (continuous spikes and waves during sleep) - oireyhtymästä. CSWS ei vastaa lääkehoitoon puolella potilaista ja sen epilepsiakirurgisesta hoidosta on olemassa vain vähän tietoa. Tavoitteet ja menetelmät: Selvitimme retrospektiivisesti epilepsiakirurgian vaikusta elektrokliinisiin löydöksiin 13:lla lääkeresistenttiä CSWS-oireyhtymää sairastavalla lapsella, joilla oli rakenteellinen aivojen poikkeavuus. Toinen tavoite oli löytää objektiivinen puoliautomaattinen tapa mitata purkauksen määrää ja piikkien voimakkuutta EEG:stä. Tulokset: Kognitiivisen tason jatkuva heikentyminen loppui 12 (92 %) potilaalla leikkauksen jälkeen. Kolme neljästä resektiopotilaasta tuli kohtauksettomaksi. Kallosotomian jälkeen kuudella kahdeksasta potilaasta päivittäiset kohtaukset vähenivät yli 90 %:lla. Purkauksen leviäminen leikkausta edeltävästi vain yhdestä hemisfääristä toiseen liittyi hyvään leikkaushoitovasteeseen. Piikki-indeksi, jossa käytetään jatkuvan purkauksen määritelmänä maksimissaan kolmea sekuntia piikkien välillä, osoittautui luotettavaksi menetelmäksi ESES:n kvantifioimiseen. Useammasta elektrodista integroitu piikkien neliöllinen keskiarvo oli piikin voimakkuuden vakaa mitta häiriintymättömässä NREM-unessa. Päätelmät: Lääkehoidolle vastaamatonta CSWS:ää sairastavat potilaat, joilla on rakenteellinen aivopoikkeavuus ja yhdensuuntainen purkauksen leviämismalli, näyttävät kohtausten vähenemisen lisäksi hyötyvän epilepsiakirurgiasta kognitiivisesti. Puoliautomaattinen piikki-indeksin kvantifiointi sopivilla käyttäjäasetuksilla ja uusi spatiaalisesti integroitu piikin voimakkuuden mittari ovat stabiileja ja lupaavia ESES:n kvantitatiivisia mittareita. Avainsanat: Unenaikainen sähköinen status epilepticus, ESES, CSWS, epilepsiakirurgia, piikki-indeksi, piikin voimakkuus, neliöllinen keskiarvo

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Kohinan poisto spektrikuvista pehmeillä morfologisilla suodattimilla

Tämän diplomityön tavoitteena oli tutkia kohinan poistoa spektrikuvista käyttäen pehmeitä morfologisia suodattimia. Työssä painotettiin impulssimaisen kohinan suodattamista. Suodattimien toimintaa arvioitiin numeerisesti keskimääräisen itseisarvovirheen, neliövirheen sekä signaali-kohinasuhteen avulla ja visuaalisesti tarkastelemalla suodatettuja kuvia sekä niiden yksittäisiä spektritasoja. Käytettyjä suodatusmenetelmiä olivat suodatus kuvapisteittäin spektrin suunnassa, suodatus koko spektrissä sekä kuutiomenetelmä ja komponenteittainen suodatus. Suodatettavat kuvat sisälsivät joko suola ja pippuri- tai bittivirhekohinaa. Parhaimmat suodatustulokset sekä numeeristen virhekriteerien että visuaalisen tarkastelun perusteella saatiin komponenteittaisella sekä kuvapisteittäisellä menetelmällä. Työssä käytetyt menetelmät on esitetty algoritmimuodossa. Suodatinalgoritmien toteutukset ja suodatuskokeet tehtiin Matlab-ohjelmistolla.

Forecasting of wind speeds and directions with artificial neural networks

In this master’s thesis, wind speeds and directions were modeled with the aim of developing suitable models for hourly, daily, weekly and monthly forecasting. Artificial Neural Networks implemented in MATLAB software were used to perform the forecasts. Three main types of artificial neural network were built, namely: Feed forward neural networks, Jordan Elman neural networks and Cascade forward neural networks. Four sub models of each of these neural networks were also built, corresponding to the four forecast horizons, for both wind speeds and directions. A single neural network topology was used for each of the forecast horizons, regardless of the model type. All the models were then trained with real data of wind speeds and directions collected over a period of two years in the municipal region of Puumala in Finland. Only 70% of the data was used for training, validation and testing of the models, while the second last 15% of the data was presented to the trained models for verification. The model outputs were then compared to the last 15% of the original data, by measuring the mean square errors and sum square errors between them. Based on the results, the feed forward networks returned the lowest generalization errors for hourly, weekly and monthly forecasts of wind speeds; Jordan Elman networks returned the lowest errors when used for forecasting of daily wind speeds. Cascade forward networks gave the lowest errors when used for forecasting daily, weekly and monthly wind directions; Jordan Elman networks returned the lowest errors when used for hourly forecasting. The errors were relatively low during training of the models, but shot up upon simulation with new inputs. In addition, a combination of hyperbolic tangent transfer functions for both hidden and output layers returned better results compared to other combinations of transfer functions. In general, wind speeds were more predictable as compared to wind directions, opening up opportunities for further research into building better models for wind direction forecasting.

Korteweg–de Vries solitoniyhtälön diskretointi- ja aikaintegrointimenetelmien vertailua

Solitoni on tunnettu ilmiönä jo 1800-luvun alkupuolelta lähtien. Se on eräänlainen muotonsa säilyttävä ja vakionopeudella etenevä aalto. 1800-luvun loppupuolella esitettiin osittaisdifferentiaaliyhtälön kuvaamaan tällaista matalassa ja kapeassa kanavassa esiintynyttä solitoniaaltoa. Tätä Kortewegin ja de Vries’n mukaan nimettyä osittaisdifferentiaaliyhtälöä tutkivat numeerisesti ensimmäistä kertaa Zabusky ja Kruskal vuonna 1965. Osittaisdifferentiaaliyhtälöiden ratkaisuun tarvitsee usein käyttää numeerisia menetelmiä. Tämän työn alkupuoli käsittelee yleisesti tarvittavia matemaattisia menetelmiä sekä KdV-yhtälön analyyttistä tarkastelua. Loppupuolella tutkitaan KdV-yhtälön mallintamista tietokoneen avulla. Zabuskyn ja Kruskalin käyttämien menetelmien lisäksi kokeillaan montaa muutakin tapaa KdV-yhtälön mallintamiseen. Näistä menetelmistä vertaillaan laskentatehokkuutta sekä menetelmän tarkkuutta. Zabuskyn ja Kruskalin käyttämä paikkadiskretointi todettiin mallinnuksissa tarkimmaksi, mutta ei kuitenkaan mallinnusaikaa tarkastellen tehokkaimmaksi. Aikaintegroinneista Runge-Kutta-menetelmät todettiin parhaiksi. Menetelmien vertailun lisäksi niistä parhaiksi havaittuja sovellettiin muutaman erikoistapauksen, kuten kolmen aallon törmäyksen, mallintamiseen.