Computational imaging has been playing an important role in understanding and analysing the captured images. Both image segmentation and restoration has been in-tegral parts of computational imaging. The studies performed in this thesis has been focussed toward developing novel algorithms for image segmentation and restoration. Study related to usage of Morozov Discrepancy Principle in Di use Optical Imaging was also presented here to show that hyper parameter selection could be performed with ease. The Laplacian of Gaussian (LoG) and Canny operators use Gaussian smoothing be-fore applying the derivative operator for edge detection in real images. The LoG kernel was based on second derivative and is highly sensitive to noise when compared to the Canny edge detector. A new edge detection kernel, called as Helmholtz of Gaussian (HoG), which provides higher di suavity is developed in this thesis and it was shown that it is more robust to noise. The formulation of the developed HoG kernel is similar to LoG. It was also shown both theoretically and experimentally that LoG is a special case of HoG. This kernel when used as an edge detector exhibited superior performance compared to LoG, Canny and wavelet based edge detector for the standard test cases both in one- and two-dimensions. The linear inverse problem encountered in restoration of blurred noisy images is typically solved via Tikhonov minimization. The outcome (restored image) of such min-imitation is highly dependent on the choice of regularization parameter. In the absence of prior information about the noise levels in the blurred image, ending this regular-inaction/hyper parameter in an automated way becomes extremely challenging. The available methods like Generalized Cross Validation (GCV) may not yield optimal re-salts in all cases. A novel method that relies on minimal residual method for ending the regularization parameter automatically was proposed here and was systematically compared with the GCV method. It was shown that the proposed method performance was superior to the GCV method in providing high quality restored images in cases where the noise levels are high Di use optical tomography uses near infrared (NIR) light as the probing media to recover the distributions of tissue optical properties with an ability to provide functional information of the tissue under investigation. As NIR light propagation in the tissue is dominated by scattering, the image reconstruction problem (inverse problem) is non-linear and ill-posed, requiring usage of advanced computational methods to compensate this. An automated method for selection of regularization/hyper parameter that incorporates Morozov discrepancy principle(MDP) into the Tikhonov method was proposed and shown to be a promising method for the dynamic Di use Optical Tomography.

The primary goal of this work is to arrive at direction tomography (DT) algorithms freed from the severe linearization in the formulation, and as-assumptions on variation of the refractive index distribution (RID), involved in the earlier approaches based on Born and Royton approximations and the Fourier di reaction theorem (FDT). To start with, a direct single-step re-covery of RID from intensity measurements is demonstrated, replacing the common two-step procedure involving, rest the recovery of phase from in-density followed by the inversion of scattered led for the RID. The information loss, unavoidable in a two-step procedure is thus successfully addressed. Secondly, an iterative method which works with a forward model obtained directly from the Helmholtz equation is developed. This forward model, though has simplifying assumptions, is more general and can accommodate larger variations in RID than that allowed in the previous linear models. The iterative procedure has an update step which uses a linearization of the forward model and a re-linearization step at the updated RID. The procedure which directly employs the measured intensities is used as part of a deterministic Gauss-Newton algorithm and a stochastic optimization algorithm which uses the ensemble Kalman lter to arrive at the recursive update. The stochastic method is found to be more noise-tolerant and efficient to take care of process model inaccuracies. The proof is seen in better reconstructions from experimental data for two example objects, namely, a graded-index optical bre and a photonic-crystal bre. It is further ob-served that the reconstructions from photonic crystal bre are blurred, noisy and less accurate. Identifying the inaccurate implementation of the modemed Helmholtz equation for large k values employing the current sampling rate as the shortcoming, a new procedure, which splits the bandwidth into smaller components using short-time Fourier Transform is developed. The set of equations arrived at, each t for a narrow frequency band, is solved and the solutions are reassembled to obtain the scattered led for the original problem. The simulated di rated intensities so obtained are better matched to their measured experimental counterparts. However, the impel-mentation of the mode end procedure is computation-intensive, for which a parallel-processing machine can be a good solution. The recovery of RID with this mode cation is not attempted in this work and is left for future implementation.

Dans cette thèse, nous présentons notre travail de recherche dans le domaine du calcul haute performance en mécanique des fluides. Avec la demande croissante de simulations à haute résolution, il est devenu important de développer des solveurs numériques pouvant tirer parti des architectures récentes comprenant des processeurs multi-cœurs et des accélérateurs. Nous nous proposons dans cette thèse de développer un solveur efficace pour la résolution sur architectures hétérogènes CPU/GPU des équations de Navier-Stokes (NS) relatives aux écoulements 3D de fluides incompressibles.Tout d'abord nous présentons un aperçu de la mécanique des fluides avec les équations de NS pour fluides incompressibles et nous présentons les méthodes numériques existantes. Nous décrivons ensuite le modèle mathématique, et la méthode numérique choisie qui repose sur une technique de prédiction-projection incrémentale.Nous obtenons une distribution équilibrée de la charge de calcul en utilisant une méthode de décomposition de domaines. Une parallélisation à deux niveaux combinée avec de la vectorisation SIMD est utilisée dans notre implémentation pour exploiter au mieux les capacités des machines multi-cœurs. Des expérimentations numériques sur différentes architectures parallèles montrent que notre solveur NS obtient des performances satisfaisantes et un bon passage à l'échelle.Pour améliorer encore la performance de notre solveur NS, nous intégrons le calcul sur GPU pour accélérer les tâches les plus coûteuses en temps de calcul. Le solveur qui en résulte peut être configuré et exécuté sur diverses architectures hétérogènes en spécifiant le nombre de processus MPI, de threads, et de GPUs.Nous incluons également dans ce manuscrit des résultats de simulations numériques pour des benchmarks conçus à partir de cas tests physiques réels. Les résultats obtenus par notre solveur sont comparés avec des résultats de référence. Notre solveur a vocation à être intégré dans une future bibliothèque de mécanique des fluides pour le calcul sur architectures parallèles CPU/GPU. / In this PhD thesis, we present our research in the domain of high performance software for computational fluid dynamics (CFD). With the increasing demand of high-resolution simulations, there is a need of numerical solvers that can fully take advantage of current manycore accelerated parallel architectures. In this thesis we focus more specifically on developing an efficient parallel solver for 3D incompressible Navier-Stokes (NS) equations on heterogeneous CPU/GPU architectures. We first present an overview of the CFD domain along with the NS equations for incompressible fluid flows and existing numerical methods. We describe the mathematical model and the numerical method that we chose, based on an incremental prediction-projection method.A balanced distribution of the computational workload is obtained by using a domain decomposition method. A two-level parallelization combined with SIMD vectorization is used in our implementation to take advantage of the current distributed multicore machines. Numerical experiments on various parallel architectures show that this solver provides satisfying performance and good scalability.In order to further improve the performance of the NS solver, we integrate GPU computing to accelerate the most time-consuming tasks. The resulting solver can be configured for running on various heterogeneous architectures by specifying explicitly the numbers of MPI processes, threads and GPUs. This thesis manuscript also includes simulation results for two benchmarks designed from real physical cases. The computed solutions are compared with existing reference results. The code developed in this work will be the base for a future CFD library for parallel CPU/GPU computations.

Cieľom predloženej diplomovej práce je numerické riešenie šírenia elektromagnetických vĺn v tuneli. Za týmto účelom bola sformulovaná integrálna rovnica a numericky riešená pomocou metódy hraničných prvkov (BEM). Implementácia v prostredí MATLAB sľubne poukazuje na nízke výpočtové nároky oproti štandardným diferenciálnych diskretizačným metódam. Súčasťou projektu je vykreslenie rozloženia elektromagnetického poľa pre rôzne profily tunelov. Overenie výsledkov je vykonané pomocou zjednodušeného analytického modelu. V rámci práce je pozorované štúdium vplyvov zmien profilu tunela a rôznych impedančných podmienok na stenách tunela na výsledné rozloženie elektromagnetického poľa vo vnútri tunela.

<div>Accurate, fast, large-scale, and circuit-intuitive electromagnetic analysis is of critical importance to the design of integrated circuits (IC) and systems. Existing methods for the analysis of integrated circuits and systems have not satisfactorily achieved these performance goals. In this work, rapid modeling and simulation methods are developed for large-scale and circuit-intuitive electromagnetic analysis of integrated circuits and systems. The derived model is correct from zero to high frequencies where Maxwell's equations are valid. In addition, in the proposed model, we are able to analytically decompose the layout response into static and full-wave components with neither numerical computation nor approximation. This decomposed yet rigorous model greatly helps circuit diagnoses since now designers are able to analyze each component one by one, and identify which component is the root cause for the design failure. Such a decomposition also facilitates efficient layout modeling and simulation, since if an IC is dominated by RC effects, then we do not have to compute the full-wave component; and vice versa. Meanwhile, it makes parallelization straightforward. In addition, we develop fast algorithms to obtain each component of the inverse rapidly. These algorithms are also applicable for solving general partial differential equations for fast electromagnetic analysis.</div><div><br></div><div>The fast algorithms developed in this work are as follows. First, an analytical method is developed for finding the nullspace of the curl-curl operator in an arbitrary mesh for an arbitrary order of curl-conforming vector basis function. This method has been applied successfully to both a finite-difference and a finite-element based analysis of general 3-D structures. It can be used to obtain the static component of the inverse efficiently. An analytical method for finding the complementary space of the nullspace is also developed. Second, using the analytically found nullspace and its complementary space, a rigorous method is developed to overcome the low-frequency breakdown problem in the full-wave analysis of general lossy problems, where both dielectrics and conductors can be lossy and arbitrarily inhomogeneous. The method is equally valid at high frequencies without any need for changing the formulation. Third, with the static component part solved, the full-wave component is also ready to obtain. There are two ways. In the first way, the full-wave component is efficiently represented by a small number of high-frequency modes, and a fast method is created to find these modes. These modes constitute a significantly reduced order model of the complementary space of the nullspace. The second way is to utilize the relationship between the curl-curl matrix and the Laplacian matrix. An analytical method to decompose the curl-curl operator to a gradient-divergence operator and a Laplacian operator is developed. The derived Laplacian matrix is nothing but the curl-curl matrix's Laplacian counterpart. They share the same set of non-zero eigenvalues and eigenvectors. Therefore, this Laplacian matrix can be used to replace the original curl-curl matrix when operating on the full-wave component without any computational cost, and an iterative solution can converge this modified problem much faster irrespective of the matrix size. The proposed work has been applied to large-scale layout extraction and analysis. Its performance in accuracy, efficiency, and capacity has been demonstrated.</div>

Vessels constructed in electrically conductive materials give rise to frequency-dependent, induced magnetic fields when waves of water cause them to roll in the Earth's magnetic field. These fields, typically referred to as roll-induced magnetic vortex fields, are a component of the ship's overall signature, where signature refers to measurable quantities which can reveal or identify objects. It is crucial for military platforms to keep the signature low and thereby increase the possibilities of operation. For magnetic signatures, this is done through strategic design and construction of the platform or by using magnetic silencing systems. The signature is then decreased to minimize the risk of detection from naval mines and marine detection systems.  This report covers the initial research on the subject of an inverse source model for roll induced magnetic fields. By limiting the analysis to two basic objects and applying a time variant magnetic field to them, we induce a magnetic field which we then model. The inverse modeling is done using magnetic dipoles as sources which are placed around the area of the object. The parameters of the model are then found by applying a least squares algorithm coupled with Tikhonov regularization. The focus of this report is the configuration of this setup in terms of measurements and sources, as well as finding a proper regularization parameter. Since the applied magnetic field is dependent on the roll frequency, also the inverse model depends on a frequency parameter in addition to the geometry and material of the object. The objects here studied are of two simple geometries, a rectangular block and a hollow cylinder. Both objects are constructed in an aluminum alloy with well known material parameters. Measurement data is gathered using a numerical solver utilizing the finite element method for solving the partial differential equations. The numerical measurement data is compared to physical measurements as well. The physical measurement data is gathered by placing the objects in a Helmholtz-cage which is used to apply a homogeneous time variant magnetic field upon them. The project was carried out at the Swedish Defence Research Agency (FOI) at the department of underwater research.

No description available.

Мы исследуем термодинамические свойства магнитных жидкостей с небольшими и умеренными значениями параметра диполь-дипольного взаимодействия в отсутствии внешнего магнитного поля. Магнитные жидкости моделируются монодисперсной системой дипольных твердых сфер. Построена теория, основанная на обратном кумулянтном преобразовании вириального разложения свободной энергии Гельмгольца в ряд по плотности. Методом компьютерного моделирования Mayer Sampling получены данные для пятого вириального коэффициента. По данным компьютерного моделирования получены аналитические выражения для четвертого и пятого вириальных коэффициентов. В работе были определены аналитические выражения для свободной энергии, давления, химического потенциала и коэффициента градиентной броуновской диффузии. Построенная теория хорошо согласуется с данными компьютерного моделирования в широкой области концентраций феррочастиц (φ≤4) и интенсивности межчастичного диполь-дипольного взаимодействия (λ≤4). / We investigate thermodynamic properties of ferrofluids with low and intermediate dipolar coupling constant. The model of dipolar hard spheres is observed. In this work a new theory is constructed based on the transformation of the virial expansion for the Helmholtz free energy into a logarithmic form. Analytical expressions of the forth and the fifth virial coefficients are obtained on the basis of numerical results of Mayer-sampling calculation. In this work analytical expressions for free energy, chemical potential, pressure and gradient diffusion coefficient were obtained. The comparison between theory and computer simulation shows good agreement for dipolar coupling constant λ≤4 and for particle volume fraction φ≤4.

Une écoulement de cavité ouverte tridimensionnel saturé non-linéairement est étudié par une approche spatio-temporelle utilisant des données expérimentales résolues à la fois en temps et en espace. Ces données ont été acquises dans deux plans longitudinaux, respectivement perpendiculaire et parallèle au fond de la cavité, dans le régime incompressible, en air ou en eau. À l'aide de multiples méthodes de décompositions globales en temps et en espace, les ondes et les structures cohérentes constituant la dynamique dans le régime permanent et pouvant être produites par des mécanismes d'instabilités différents sont identifiées et caractérisées.Tout d'abord, on approfondit la compréhension de l'effet des non-linéarités sur les oscillations auto-entretenues de la couche cisaillée impactante et leurs interactions avec l'écoulement intra-cavitaire. En particulier, l'analyse spectrale d'une portion de l'espace des paramètres permet de mettre en évidence un lien entre l'accrochage des modes d'oscillations auto-entretenues, la modulation d'amplitude au niveau du coin impactant et l'intermittence de ces modes. De plus, l'observation des basses fréquences intéragissant fortement avec les oscillations de la couche de mélange démontre l'existence d'une dynamique tridimensionnelle intrinsèque à l'intérieur de la cavité malgré les perturbations causées par la couche cisaillée instable.Les analyses de stabilité linéaire ont montré que des instabilités centrifuges peuvent résulter de la courbure induite par la recirculation. L'étude de la dynamique après saturation révèle de nombreuses structures cohérentes dont les propriétés sont quantifiées et classées en s'appuyant sur la forme des instabilités sous-jacentes: des ondes transverses progressives ou stationnaires. Enfin, certains comportements des structures saturées suggèrent que les mécanismes non-linéaires gouvernant le développement de l'écoulement une fois sorti du régime linéaire pourraient être étudiés dans le cadre des équations d'amplitude. / A space-time study of a three-dimensional nonlinearly saturated open cavity flow is undertaken using time-resolved space-extended experimental data, acquired in both cross-stream and spanwise planes, in incompressible air and water flows. Through use of multiple modal decompositions in time and space, the waves and coherent structures composing the dynamics in the permanent regime are identified and characterised with respect to the instabilities arising in the flow.Effects of nonlinearities are thoroughly investigated in the impinging shear layer, regarding the self-sustained oscillations and their interactions with the inner-flow. In particular, the analysis conducted throughout the parameter space enlightens a global connection between the selection of locked-on modes and the amplitude modulation at the impingement and the mode switching phenomenon. Furthermore, observations of low frequencies interacting drastically with the shear layer flapping motion underline the existence of intrinsic coherent three-dimensional dynamics inside the cavity in spite of the shear layer disturbances.Linear stability analyses have demonstrated that centrifugal instabilities are at play along the main recirculation. The present investigation of the dynamics after onset of the saturation reveals numerous space-time coherent structures, whose properties are quantified and classified with respect to the underlying instabilities: travelling or standing spanwise waves. Finally, some patterns exhibited by the saturated structures suggest that the nonlinear mechanisms governing the mutations of the flow after the linear regime could gain more insight in the frame of amplitude equations.

Cette thèse est une contribution au calcul des variations et à la théorie du contrôle optimal dans les cadres discret, plus généralement time scale, et fractionnaire. Ces deux domaines ont récemment connu un développement considérable dû pour l’un à son application en informatique et pour l’autre à son essor dans des problèmes physiques de diffusion anormale. Que ce soit dans le cadre time scale ou dans le cadre fractionnaire, nos objectifs sont de : a) développer un calcul des variations et étendre quelques résultats classiques (voir plus bas); b) établir un principe du maximum de Pontryagin (PMP en abrégé) pour des problèmes de contrôle optimal. Dans ce but, nous généralisons plusieurs méthodes variationnelles usuelles, allant du simple calcul des variations au principe variationnel d’Ekeland (couplé avec la technique des variations-aiguilles), en passant par l’étude d’invariances variationnelles par des groupes de transformations. Les démonstrations des PMPs nous amènent également à employer des théorèmes de point fixe et à prendre en considération la technique des multiplicateurs de Lagrange ou encore une méthode basée sur un théorème d’inversion locale conique. Ce manuscrit est donc composé de deux parties : la Partie 1 traite de problèmes variationnels posés sur time scale et la Partie 2 est consacrée à leurs pendants fractionnaires. Dans chacune de ces deux parties, nous suivons l’organisation suivante : 1. détermination de l’équation d’Euler-Lagrange caractérisant les points critiques d’une fonctionnelle Lagrangienne ; 2. énoncé d’un théorème de type Noether assurant l’existence d’une constante de mouvement pour les équations d’Euler-Lagrange admettant une symétrie ; 3. énoncé d’un théorème de type Tonelli assurant l’existence d’un minimiseur pour une fonctionnelle Lagrangienne et donc, par la même occasion, d’une solution pour l’équation d’Euler-Lagrange associée (uniquement en Partie 2) ; 4. énoncé d’un PMP (version forte en Partie 1, version faible en Partie 2) donnant une condition nécessaire pour les trajectoires qui sont solutions de problèmes de contrôle optimal généraux non-linéaires ; 5. détermination d’une condition de type Helmholtz caractérisant les équations provenant d’un calcul des variations (uniquement en Partie 1 et uniquement dans les cas purement continu et purement discret). Des théorèmes de type Cauchy-Lipschitz nécessaires à l’étude de problèmes de contrôle optimal sont démontrés en Annexe. / This dissertation deals with the mathematical fields called calculus of variations and optimal control theory. More precisely, we develop some aspects of these two domains in discrete, more generally time scale, and fractional frameworks. Indeed, these two settings have recently experience a significant development due to its applications in computing for the first one and to its emergence in physical contexts of anomalous diffusion for the second one. In both frameworks, our goals are: a) to develop a calculus of variations and extend some classical results (see below); b) to state a Pontryagin maximum principle (denoted in short PMP) for optimal control problems. Towards these purposes, we generalize several classical variational methods, including the Ekeland’s variational principle (combined with needle-like variations) as well as variational invariances via the action of groups of transformations. Furthermore, the investigations for PMPs lead us to use fixed point theorems and to consider the Lagrange multiplier technique and a method based on a conic implicit function theorem. This manuscript is made up of two parts : Part A deals with variational problems on time scale and Part B is devoted to their fractional analogues. In each of these parts, we follow (with minor differences) the following organization: 1. obtaining of an Euler-Lagrange equation characterizing the critical points of a Lagrangian functional; 2. statement of a Noether-type theorem ensuring the existence of a constant of motion for Euler-Lagrange equations admitting a symmetry;3. statement of a Tonelli-type theorem ensuring the existence of a minimizer for a Lagrangian functional and, consequently, of a solution for the corresponding Euler-Lagrange equation (only in Part B); 4. statement of a PMP (strong version in Part A and weak version in Part B) giving a necessary condition for the solutions of general nonlinear optimal control problems; 5. obtaining of a Helmholtz condition characterizing the equations deriving from a calculus of variations (only in Part A and only in the purely continuous and purely discrete cases). Some Picard-Lindelöf type theorems necessary for the analysis of optimal control problems are obtained in Appendices.