We simulate ultrashort pulse propagation through water by numerical methods, which is a kind of optical communication research. Ultrashort pulses have been known to have non Beer-Lambert behavior, whereas continuous waves (CW) obey the Beer-Lambert law. People have expected that the ultrashort pulse loses less intensity for a given distance in water than CW which implies that the pulse can travel over longer distances. In order to understand this characteristic of the pulse, we model numerically its spectral and temporal evolution as a function of traveling distance through water. We achieve the pulse intensity attenuation with traveling distance, obtain the temporal envelope of the pulse and compare them with experimental data. This research proves that the spectral and temporal profile of a pulse can be predicted knowing only the intensity spectrum of the input pulse and the refractive index spectrum of water in the linear regime. The real feasibility and the advantage of using an ultrashort pulse as a communication carrier will also be discussed.

It is estimated that mind-boggling 14.1 million new cases of cancer occurred worldwide in 2012 alone. This number is alarming. Although healthy lifestyle may reduce a risk of developing cancer, there is always some probability that cancer would develop even in an absolutely fit individual. There are two main conditions for successful treatment of cancer. Firstly, early diagnostic is absolutely crucial. Secondly, there is a need for suitable surgical methods for affected tissue removal. Ultrasound has a great potential to be used for both purposes as a non-invasive method. Photoacoustic spectroscopy is imaging method for tumor detection of great properties making the use of ultrasound while High-Intensity Focused Ultrasound (HIFU) is non-invasive surgical method. These methods would be impossible without precise ultrasound propagation simulations. The k-Wave is an open source MATLAB toolbox implementing such simulations. So, why are not these methods already deployed in treatment? Unfortunately, the simulation of ultrasound propagation is a very time consuming task, which makes it ineffective for medical purposes. However, there are a few options how to accelerate these simulations. The use of GPU is a very promising way to accelerate simulation.   The main topic of this thesis is the acceleration of the simulation of soundwaves propagation in bones and hard tissue. The implementation developed as a part of this thesis was benchmarked on various supercomputers including Anselm in Ostrava and Piz Daint in Lugano. The implemented solution provides remarkable acceleration compared to the original MATLAB prototype. It was able to accelerate the simulation around 160 times in the best case. It means that the simulation, which would otherwise last for 6.5 days, can be now computed in one hour. This acceleration was achieved using an NVIDIA Tesla P100 to run the simulation with the domain size of 416x416x416 grid points. The thesis includes performance benchmarks on different GPUs to provide complex image acceleration capabilities of developed implementation and provides discussion about memory usage and numerical accuracy. Thanks to the implemented solution harnessing the power of modern GPUs, doctors and researchers all around the world have a powerful tool in hands.

This research was a collaborative effort between the Air Force Research Laboratory (AFRL) and the University of Texas at Austin to study laser induced thermal lensing. Analysis with a high-frame rate thermal camera system in an optical cuvette determined rates of heating and cooling in water exposed to 1200, 1310 and 1318-nm wavelengths. Thermal modeling software (BTEC) developed at AFRL was used to simulate the beam parameters that were used experimentally. The simulation was also used to compute axial temperature for various power levels, beam diameters, and pulse durations. Laser-induced optical pathlength modulation or phase-shift was computed to study the thermal lensing effect. Power and irradiance damage thresholds were calculated for collimated and focused geometries to study the effect in a focusing eye on retinal damage thresholds. / text

Current propogation models are no longer sufficient for wireless network planning. They are neither accurate (empirical) nor fast enough (deterministic) to be applicable in the applications of Automated Cell Planning. This thesis focuses on the development of a new method, namely Intelligent Ray Launching Algorithm (IRLA), which is based on a fast, accurate and robust algorithm that is especially suitable for wireless network planning. The infrastructure of IRLA is thoroughly analysed in this thesis and the results are presented. Foster's design methodology has been used to parallelise the new model. Various scenarios for outdoor, indoor, indoor-to-outdoor and outdoor-to-indoor settings have been employed to test the effectiveness and efficiency of IRLA. The field strengths (path loss) and multipath information were calculated, which were used to demonstrate the application of IRLA. The accuracy of IRLA is guaranteed via the use of a meta-based heuristics calibration procedure. In order to achieve a simulation within a realistic time scale, acceleration techniques such as avoid double marking, multi-threading and the use of Parallel Object-Oriented Programming C++ have been employed. Since multipath for a large number of receiver locations can be easily obtained via IRLA, the study of delay spread has been presented. The success of the integration with a wireless network planning platform exemplifies that IRLA is suitable for wireless network planning and optimisation, which is beneficial to relevant academics and industries. Testing demonstrated that depending on various scenarios, IRLA obtains industrially-recognised accuracy ranging from 5 to 8 dB Root-Mean-Square-Error. This model is highly-efficient because its required runtime for most simulations is from a few seconds to a few minutes.

L’imagerie ultrasonore harmonique, qui repose sur la non linéarité du milieu de propagation, est une technique d’imagerie clinique qui améliore la résolution des images. La mesure ultrasonore du paramètre local de non linéarité d'un milieu est une voie de recherche qui amènerait de nouvelles perspectives dans le domaine de la caractérisation des tissus. Cependant, l'accès à cette information se heurte à deux écueils : d'une part il n’existe pas actuellement de méthode de mesure de ce paramètre à partir du mode écho classique et d'autre part, les outils de simulation prenant en compte la non-linéarité du milieu sont peu développés. Une méthode de spectre angulaire a donc été proposée afin de calculer le champ de pression dans des milieux de non linéarité inhomogène. Ce champ de pression est ensuite utilisé pour engendrer des images échographiques contenant l’information harmonique. Cette méthode spectrale a été portée sur GPU afin d’accélérer le calcul et a été intégrée dans un logiciel libre : CREANUIS. Dans un deuxième temps, une extension d’une méthode comparative (ECM) a été proposée pour prendre en compte des milieux de non linéarité non homogène, fonctionnant en mode écho. Grâce aux outils de simulation développés, différentes configurations ont été utilisées pour la mise au point de l’ECM qui a ensuite été validée à partir d'objets tests et in vitro sur foies d’animaux. Même si la méthode de mesure présente une résolution relativement faible, les images obtenues démontrent le potentiel de l’imagerie du paramètre de non linéarité des tissus. / Harmonic imaging, based on the propagated medium nonlinearity, is a clinical imaging technique which increases the resolution of ultrasound images. The ultrasound measure of the local nonlinear parameter brings new perspectives in term tissues characterization. However, access to this information suffers from two strong points: from one hand, there is no current measurement method of this parameter in echo mode configuration and on the other hand, the simulation tools taking into account the nonlinearity are not many developed. An angular spectrum method has been proposed to compute the nonlinear pressure field with inhomogeneous nonlinear parameter. This pressure field is then used to generate ultrasound images containing the harmonic component. This spectral approach has been implemented on a GPU in order to accelerate the computation and package in a free software made available to the scientific community under the name CREANUIS. In a second time, a extension of a comparative method (ECM) has been proposed to take into account media with inhomogeneous nonlinearity, working an echo mode configuration. Thanks the developed simulation tools, different configurations have been used to parameterize and to evaluate the ECM which has then be validated on test objects and in vitro animal’s livers. Even if the measure presents a relatively weak resolution, the obtained images demonstrated a high potential in the nonlinear parameter imaging of tissues.

L’intégration du monde physique et des systèmes d’information est l’un des faits majeurs de la récente décennie. Cette intégration est rendue possible par les systèmes de communication, en particulier les systèmes sans fil, et par les recherches à l’interface physique dans le domaine des capteurs.L’usage de cette interface peut se concevoir de manière locale, par exemple le contrôle d’un objet, ou une aide dirigée vers une personne. Cette interface peut aussi exister au niveau distribué, en permettant la synthèse d’informations éparses et des prises de décisions concernant par exemple des économies de ressources : éteindre l’éclairage d’une rue, guider vers une place de parking libre.Les réseaux de capteurs sans fil contribuent à cette seconde catégorie et cette thèse s’est consacrée à la simulation de la propagation d’informations sonores en milieu urbain, avec une caractérisation distribuée.Le support décisif pour cette simulation est la représentation cellulaire de la ville et de ses zones : rues, jardins, rocades, bâtiments, rivières. Cette représentation cellulaire est produite par l’analyse d’imagesgéo-localisées et peut être complétée par la consultation de bases de données, pour l’élévation, par exemple. La représentation cellulaire est ensuite transformée en systèmes de processus interconnectés, à même de reproduire quantité de comportements collectifs, physiques ou numériques. / Physical world and information systems unification is a major achievement during this last decade. Communication systems, especially wireless technology, and research in physical sensors field have taken a major part in this integration.The interface between the physical word and a computer can be used in two ways. First, from a local perspective, where such a system can control for instance a single object or provide a personal assistance.Moreover, in a distributed system, this interface collects and analyses sparse data to make decisions, for instance about resource savings: to turn off streetlights, to guide to a parking spot.Wireless sensor networks are part of this second category. Thus, this thesis is focused on sound propagation simulation in urban environment, in a distributed manner.The key point of this work is a cellular representation of the city in order to model streets, gardens, ring roads, buildings, and rivers. Georeferenced image analysis, complemented by database consultation, for example to retrieve elevation data, generates this cellular model. Then, this model is converted into a system made of interconnected processes, that can reproduce many collective behaviors, whether physical or digital.We have developed a cellular automata that models sound propagation, including reflection and refraction, working on a graphics accelerator. Finally, we have produced a coupling method between observation systems by sensor networks and physical systems.

Simulácia šírenia ultrazvuku prostredníctvom mäkkých biologických tkanív má širokú škálu praktických aplikácií. Patria sem dizajn prevodníkov pre diagnostický a terapeutický ultrazvuk, vývoj nových metód spracovania signálov a zobrazovacích techník, štúdium anomálií ultrazvukových lúčov v heterogénnych médiách, ultrazvuková klasifikácia tkanív, učenie rádiológov používať ultrazvukové zariadenia a interpretáciu ultrazvukových obrazov, modelové vrstvenie medicínskeho obrazu a plánovanie liečby pre ultrazvuk s vysokou intenzitou. Ultrazvuková simulácia však predstavuje výpočtovo zložitý problém, pretože simulačné domény sú veľmi veľké v porovnaní s akustickými vlnovými dĺžkami, ktoré sú predmetom záujmu. Ale ak je problém osovo symetrický, problém môže byť riešený v 2D.To umožňuje spúšťanie simulácií na mriežke s väčším počtom bodov, s menším využitím výpoč- tových zdrojov za kratšiu dobu. Táto práca modeluje a implementuje zrýchlenie vlnovej nelineárnej ultrazvukovej simulácie v axisymetrickom súradnicovom systéme realizovanom v Matlabe pomocou Mex súborov pre diskrétne sínové a kosínové transformácie. Axisymetrická simulácia bola implementovaná v C++ ako open source rozšírenie K-WAVE toolboxu. Kód je optimalizovaný na beh na jednom uzle superpočítaču Salomon (IT4Innovations, Ostrava, Česká republika) s dvoma dvanásť-jadrovými procesormi Intel Xeon E5-2680v3. Na maximalizáciu výpočtovej efektívnosti boli vykonané viaceré optimalizácie kódu. Po prvé, fourierové tramsformácie boli vypočítané pomocou real-to-complex FFT z knižnice FFTW. V porovnaní s complex-to-complex FFT to znížilo čas výpočtu a pamäť spojenú s výpočtom FFT o takmer 50%. Taktiež diskrétne sínové a kosínové transformácie sa počítali pomocou knižnice FFTW, ktoré v Matlab verzii museli byť vyvolané z dynamicky načítaných MEX súborov. Po druhé, aby sa znížilo zaťaženie priepustnosti pamäte, boli všetky operácie počítané jednoduchej presnosti pohyblivej rádovej čiarky. Po tretie, elementárne operá- cie boli paralelizované pomocou OpenMP a potom vektorizované pomocou rozšírení SIMD (SSE). Celkový výpočet C++ verzie je až do 34-násobne rýchlejší a využíva menej ako tretinu pamäte ako Matlab verzia simulácie. Simulácia ktorá by trvala takmer dva dni tak môže byť vypočítaná za jeden a pol hodinu. Toto všetko umožňuje počítať simuláciu na výpočetnej mriežke s veľkosťou 16384 × 8192 bodov v primeranom čase.

Στην παρούσα εργασία παρουσιάζεται μια μέθοδος για την αριθμητική προσομοίωση δισδιάστατης, μη συνεκτικής ροής με ελεύθερη επιφάνεια, που προκύπτει από τη διάδοση κυμάτων βαρύτητας πάνω από πυθμένα με τυχαία μορφολογία. Η μέθοδος βασίζεται στην αριθμητική επίλυση των εξισώσεων Euler, που υπόκεινται σε πλήρως μη γραμμικές οριακές συνθήκες ελεύθερης επιφάνειας και κατάλληλες οριακές συνθήκες πυθμένα, εισόδου και εξόδου, χρησιμοποιώντας ένα υβριδικό σχήμα πεπερασμένων διαφορών και ψευδό-φασματικής μεθόδου. Οι εξισώσεις ροής μετασχηματίζονται έτσι ώστε τα όρια του υπολογιστικού πεδίου να είναι ανεξάρτητα του χρόνου. Η επαλήθευση της μεθόδου επίλυσης γίνεται με την εφαρμογή της στο πρόβλημα της κατανομής θερμοκρασίας σε λεπτή ορθογωνική πλάκα, υπό σταθερές συνθήκες. Για την ελαχιστοποίηση της ανάκλασης χρησιμοποιείται ζώνη απορρόφησης στην περιοχή απορροής. Η αποτελεσματικότητα της ζώνης απορρόφησης τεκμηριώνεται με την παρουσίαση αποτελεσμάτων προσομοίωσης διάδοσης γραμμικών κυματισμών σε πυθμένα σταθερού βάθους. Προκύπτει ότι η ζώνη απορρόφησης που βασίζεται στην επιβολή εξωτερικής δυναμικής πίεσης στην ελεύθερη επιφάνεια εμφανίζει την ελάχιστη ανάκλαση. Αποτελέσματα παρουσιάζονται για την προσομοίωση ροής με ελεύθερη επιφάνεια πάνω από πυθμένα σταθερής κλίσης 1:10 και 1:50, για διαφορετικά μήκη και ύψη κυμάτων εισόδου. Ο μετασχηματισμός των γραμμικών κυμάτων πάνω από την περιοχή σταθερής κλίσης συμφωνεί με τη θεωρία γραμμικής διασποράς για ροή με μικρά ύψη κύματος. Για μη γραμμικούς κυματισμούς, το μήκος κύματος μειώνεται πάνω από την περιοχή σταθερής κλίσης, ενώ η ανύψωση της ελεύθερης επιφάνειας αποκλίνει από την αρχική ημιτονοειδή μορφή και το ύψος κύματος αυξάνει λόγω της ρήχωσης. / A method for the numerical simulation of two-dimensional, inviscid, free-surface flow resulting from the propagation of regular gravity water waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid scheme of finite-differences and pseudo-spectral method. The formulation includes a boundary-fitted transformation. The validation of the pressure solver is accomplished by applying it to the two-dimensional, temperature, steady problem, under differing boundary conditions. For the free-surface flow, a wave absorption zone is attached at the outflow domain in order to minimize reflection effects. The absorption zone effectiveness is validated by the simulation of linear waves propagation over constant-depth bottom. Minimal reflection occurs when an appropriate external dynamic pressure is imposed on the free surface of the absorption zone. Results are presented for cases of wave propagation over constant slope bottom, with slopes 1:10 and 1:50, for various incoming wavelengths and wave heights. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape and the wave height increases due to shoaling.

This thesis focuses on the development of tools and methods dedicated for ultra wide band (UWB) localization systems in indoor environment. The thesis work was conducted within the European FP7 project Where2, about the cooperative localization in cellular networks. Data from a measurement campaign conduct during the project are used to validate the proposed algorithms. This thesis is divided in four parts : The first part is focused on the description of an original raytraing tool based on a graph description. In order to be compliant with the requirement of a mobile simulation, a new concept of rays signature enabling incremental computation, and a vectorized formalism for processing rays are described and implemented. The second part is focused on the indoor localization techniques, where a novel technique based on interval analysis approaches is presented and compared to alternative techniques. Advantageously using this approach, a specific processing based on an hypothesis testing method using received power observations to resolve ambiguities appearing in under determined localization problems is described. A third part describes different aspects of the dynamic platform. In particular a realistic mobility model based on ''steering behaviors'', a graph description of the network scene and an inter agents communication protocol are detailed. The fourth section uses measured data obtained from an heterogeneous measurement campaign to validate both the developed software platform and the proposed localization algorithms.

Die vorliegende Dissertation befasst sich mit Entstehungsmechanismen kardialer Arrhythmien auf der Grundlage pathologisch veränderten Myokards. Es wurde eine systematische Analyse der elektrophysiologischen Veränderungen, die als Folge von Gap-Junction- Lateralisierung und zellulärer Hypertrophie auftreten, durchgeführt. Die Analyse beruht auf einem mathematischen Computermodell, das zur Simulation der Aktionspotentialausbreitung innerhalb einer Einzelzellschicht humaner ventrikulärer Kardiomyozyten entwickelt wurde. Ausgehend von bestehenden Einzelzellmodellen wurde ein räumlich und zeitlich hoch aufgelöstes Multizellmodell generiert und in der Programmiersprache Object Pascal implementiert. Nach Validierung des Modells wurde es zur gezielten, an experimentellen Daten orientierten Manipulation geometrischer Eigenschaften der Zellen (Länge, Durchmesser) und des Zellverbandes (Anordnung der Zellen untereinander) sowie der Gap-Junction-Verteilung genutzt. Die Analyse der elektrophysiologischen Effekte im Vergleich zur Kontrolle fand sowohl unter Normalbedingungen als auch unter Bedingungen, die pathologischen Veränderungen entsprechen (Entkopplung der Gap-Junctions, verringerte Aktivität des schnellen Natriumkanals, erhöhte Inhomogenität), statt. Es zeigte sich, dass ein größerer Zelldurchmesser bzw. erhöhte laterale Gap-Junction-Leitfähigkeit (Simulation von kardialer Hypertrophie bzw. Connexin- Lateralisierung) die Entstehungswahrscheinlichkeit eines unidirektionalen Leitungsblocks erhöhte. Die Erregungsausbreitungsgeschwindigkeit in hypertrophierten Zellen war zudem weniger stabil als in normalen Zellen. Beide Effekte gehören zu den Hauptursachen der Entstehung und Aufrechterhaltung ventrikulärer Arrhythmien. Die Ergebnisse der Arbeit erklären somit Ursachen des erhöhten Arrhythmierisikos in pathologisch veränderten und hypertrophierten Herzen und liefern eine theoretische Grundlage für zukünftige Studien.