La reconstruction des sources de l'activité cérébrale à partir des mesures de potentiel fournies par un électroencéphalographie (EEG) nécessite de résoudre le problème connu sous le nom de « problème inverse de l'EEG ». La solution de ce problème dépend de la solution du « problème direct de l'EEG », qui fournit à partir de sources de courant connues, le potentiel mesuré au niveau des électrodes. Pour des modèles de tête réels, ce problème ne peut être résolut que de manière numérique. En particulier, les équations intégrales de surfaces requièrent uniquement la discrétisation des interfaces entre les différents compartiments constituant le milieu cérébral. Cependant, les formulations intégrales existant actuellement ne prennent pas en comptent l'anisotropie du milieu. Le travail présenté dans cette thèse introduit deux nouvelles formulations intégrales permettant de palier à cette faiblesse. Une formulation indirecte capable de prendre en compte l'anisotropie du cerveau est proposée. Elle est discrétisée à l'aide de fonctions conformes aux propriétés spectrales des opérateurs impliqués. L'effet de cette discrétisation de type mixe lors de la reconstruction des sources cérébrales est aussi étudié. La seconde formulation se concentre sur l'anisotropie due à la matière blanche. Calculer rapidement la solution du système numérique obtenu est aussi très désirable. Le travail est ainsi complémenté d'une preuve de l'applicabilité des stratégies de préconditionnement de type Calderon pour les milieux multicouches. Le théorème proposé est appliqué dans le contexte de la résolution du problème direct de l'EEG. Un préconditionneur de type Calderon est aussi introduit pour l'équation intégrale du champ électrique (EFIE) dans le cas de structures unidimensionnelles. Finalement, des résultats préliminaires sur l'impact d'un solveur rapide direct lors de la résolution rapide du problème direct de l'EEG sont présentés. / Electroencephalography (EEG) is a very useful tool for characterizing epileptic sources. Brain source imaging with EEG necessitates to solve the so-called EEG inverse problem. Its solution depends on the solution of the EEG forward problem that provides from known current sources the potential measured at the electrodes positions. For realistic head shapes, this problem can be solved with different numerical techniques. In particular surface integral equations necessitates to discretize only the interfaces between the brain compartments. However, the existing formulations do not take into account the anisotropy of the media. The work presented in this thesis introduces two new integral formulations to tackle this weakness. An indirect formulation that can handle brain anisotropies is proposed. It is discretized with basis functions conform to the mapping properties of the involved operators. The effect of this mixed discretization on brain source reconstruction is also studied. The second formulation focuses on the white matter fiber anisotropy. Obtaining the solution to the obtained numerical system rapidly is also highly desirable. The work is hence complemented with a proof of the preconditioning effect of Calderon strategies for multilayered media. The proposed theorem is applied in the context of solving the EEG forward problem. A Calderon preconditioner is also introduced for the wire electric field integral equation. Finally, preliminary results on the impact of a fast direct solver in solving the EEG forward problem are presented.

Cette thèse est consacrée à l'étude numérique des interactions non-linéaires entre des vagues et un corps rigide perçant la surface libre. La méthode développée repose sur un modèle d'éléments de frontière qui réduit la dimensionnalité du problème d'une dimension. Dans un premier temps, un modèle2D est appliqué à des géométries simples et permet de démontrer la pertinence de l'approche envisagée pour la prédiction des mouvements d'une structure flottante soumise à des vagues monochromatiques régulières. Dans un second temps, en nous inspirant d'un modèle potentiel non-linéaire 3D développé par Grilli textit{et al.}~cite{grilli2001fully}, nous proposons une généralisation de la méthode pour des maillages triangulaires non-structurés de surfaces 3D. Le modèle développé permet de traiter des configurations arbitraires de plusieurs cylindres verticaux en interaction avec les vagues. Nous présentons des cas de validation de nature académique qui permettent d'apprécier le comportement du modèle numérique. Puis nous nous tournons vers l'application visée par EDF R&D, qui concerne le dimensionnement d'éoliennes off-shore flottantes. Un flotteur de type semi-submersible est évalué à l'aide du modèle non-linéaire / This PhD work is devoted to the study of nonlinear interactions between waves and floating rigid structures. The developed model relies on a boundary element method which reduces the dimensionality of the problem by one. First, a 2D model is applied to basic geometries and allows us to demonstrate the validity of the method for predicting the motion of a floating structrure subject to incoming monochromatic regular waves. Secondly, getting inspired by the 3D fully nonlinear potential flow model of Grilli textit{et al.}~cite{grilli2001fully}, we propose a novel model which generalizes the method for unstructured triangular meshes of 3D surfaces. The proposed model is able to deal with arbitrary configurations of multiple vertical cylinders interacting with the waves. We present academic validation test cases which show how the model works and behaves. Finally, we study situations of interest for EDF R&D related to floating off-shore wind turbines. A semi-submersible platform is evaluated with the nonlinear model

Orientador: Paulo Sollero / Tese (doutorado) - Universidade Estaudal de Campinas. Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-07T18:05:59Z (GMT). No. of bitstreams: 1 Paiva_WilliamPortilhode_D.pdf: 26643195 bytes, checksum: 744132c2c47025a4ed4eec9de86c860b (MD5) Previous issue date: 2005 / Resumo: Esta tese apresenta o desenvolvimento, a implementação e aplicações de uma formulação de elementos de contorno para análise de problemas estáticos e dinâmicos em placas anisotrópicas submetidas a carregamentos de flexão. O trabalho se baseia na teoria de placas finas de Kirchhoff Uma vez que a formulação dinâmica de elementos de contorno não foi encontrada na literatura consultada, as análises são realizadas utilizando-se uma formulação desenvolvida a partir das soluções fundamentais para elastostática. As integrais de domínio, provenientes dos carregamentos distribuídos ou dos termos de inércia, são transformadas em integrais de contorno usando-se o método da reciprocidade dual. Esta formulação é aplicada no cálculo da deflexão de placas isotrópicas e anisotrópicas, na análise dinâmica de estruturas ortotrópicas, anisotrópicas e estruturas de materiais compósitos laminados, submetidas a cargas harmônicas. É analisada a influência dos nós internos e da malha na precisão dos resultados. Os resultados numéricos são apresentados para problemas quaseisotrópicos, ortotrópicos e totalmente anisotrópicos. Os resultados obtidos são comparados com resultados presentes na literatura e mostram boa concordância / Abstract: Here, static and dynamic analysis of a number of anisotropic problems using the boundary element method are presented. The formulation used is developed fram the elastostatic fundamental solutions and the domain integraIs due to inertial terms are transformed into boundary integraIs by using the dual recipracity method. Application is made to the calculation of deflections in anisotrapic plates, to the dynamic analysis of orthotrapic plates and to the analysis of anisotropic and composite laminar plates under harmonic loads. Numerical results are presented for quasi-isotropic, orthotropic and anisotrapic problems. The influence of number of internal nodes and number of boundary elements on the accuracy of the results is analyzed. The results are compared with other results fram the literature showing good agreement / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Orientadores: Paulo Sollero, Eder Lima de Albuquerque / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-28T12:08:28Z (GMT). No. of bitstreams: 1 Shaterzadeh-Yazdi_MohammadHossein_M.pdf: 5411191 bytes, checksum: 83da697ff892a31af99059f3e88bd338 (MD5) Previous issue date: 2015 / Resumo: Em problemas de contato mecânico entre dois corpos elásticos, o cálculo de tensões e deformações dos componentes é de grande importância. Em casos particulares os corpos estão sujeitos a cargas normal e tangencial na presença de atrito, o qual aumenta a complexidade do problema. O estudo do fenômeno e a modelagem do problema, empregando o método dos elementos de contorno (MEC), é apresentado neste trabalho. Devido à presença de atrito e restrições de contato, esse problema torna-se um caso não linear. A não linearidade do problema foi contornada com a aplicação incremental de carga e o uso de um método de resolução de sistemas não lineares. A zona de contato é uma das variáveis do problema e pode conter estados de adesão e escorregamento, simultaneamente. Esses estados dependem dos esforços normais e tangenciais no componente e podem variar durante o processo de aplicação de carga. Dessa forma, cada incremento de carga pode perturbar em relação ao estado anterior. Portanto, o cálculo de variáveis e a atualização do sistema de equações em cada iteração é indispensável. Por este motivo, um algoritmo robusto para definição dos estados de contato é proposto. Como o sistema de equações obtido é não linear, o uso de um método numérico adequado é exigido. Para a solução deste sistema, o método de Newton foi aplicado, o qual permite a verificação do estado de contato em cada incremento. A análise é feita com o uso de elementos quadráticos contínuos, apresentando resultados contínuos e sem oscilação. A comparação dos resultados com as soluções analíticas de Hertz e Mindlin-Cattaneo mostram boa concordância / Abstract: The computation of stresses and strains on the components is of great importance, when the contact mechanics problems between two elastic bodies are analyzed. In particular cases, bodies are subjected to normal and shear loading in the presence of friction, which increases the complexity of the problem. The study of the phenomenon and modeling of the problem, using the boundary element method (BEM), are presented in this work. Due to the presence of friction and natural restrictions, this problem becomes non-linear. The non-linearity of the problem was solved with an incremental applied load and with the use of solvers to non linear systems. The contact zone can contain stick and slip states, simultaneously. These states are dependent on the normal and shear forces on the component and can vary during the application load process. Thus, each load increment can violate the previous state and therefore, the evaluation of variables and the updating of the system of equations after each iteration is indispensable For this reason, a robust algorithm for contact state definition is suggested. Since a non linear system of equations is obtained, an appropriate numerical method is required. To solve this system, Newton¿s method is applied, which allows the verification of the state of contact at each increment. The analysis is done with the use of quadratic continuous elements and provides continuous and non-oscillatory results. Comparisons of the results with the analytical solutions of Hertz and Mindlin-Cattaneo show good agreement / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

De nombreux modèles numériques ont été proposés pour mieux comprendre comment le courant électrique est répartit lors d’une stimulation électrique par implant cochléaire. Ceci permet à terme d’optimiser la géométrie des électrodes et les stratégies de stimulation. Les modèles précédemment proposés modélisent les modèles d'interface électrochimique de façon très basique, et ne prennent généralement compte que de l'intensité du courant sur les électrodes. Par conséquent, il leur est difficile de simuler la dynamique temporelle de la stimulation ou de modéliser la répartition du courant en fonction de différents modes de stimulation contrôlés en tension, tels que le mode de de retour commun (Common Ground), ou de retour multiple (Multi-Mode Grounding). Dans cette thèse, nous avons développé un nouveau modèle surfacique de la cochlée. Le modèle géométrique dépend d'un ensemble de paramètres permettant d'ajuster la forme de la cochlée, en utilisant par exemple des données histologiques, des scans CT, ou encore des maillages de surface. Un modèle paramétrique nous a permis de comparer les courants générés par les modèles d'électrodes actuellement disponibles et par un nouveau type d'électrode - faisceau transmodiolaire. Le modèle peut prendre en compte des courants ou des tensions en entrée à chaque électrode, ce qui permet de simuler le mode de retour commun ou multiple. Afin de valider les résultats de simulation et calibrer les paramètres du modèle, nous avons créé un système permettant d'acquérir des mesures in-situ et in-vitro. Les données enregistrées ont permis de valider le modèle combinant le modèle d'interface électrochimique et le modèle tridimensionnel de cochlée. / Computational models have been widely used to improve the knowledge of the current distribution behavior in cochlear implant stimulations, optimizing electrode designs and stimulation strategies. The existing models employed no or simple electrochemical interface models and took current intensity on the electrodes as input. Therefore they have difficulties in making time domain simulations and simulating the stimulation modes that have voltage constraints, such as the Common Ground and the Multi-Mode Grounding modes. In this PhD work, a new parametric surface mesh model of the cochlea has been developed. The shape of the model is controlled by a set of input parameters which can be tuned to fit the cochlear shape acquired from histological images, CT scans or existing cochlear mesh models. The symmetric boundary element method, which was implemented in OpenMEEG, has been applied on the model to simulate the current distribution of the cochlear implant stimulation. Using the parametric model, comparisons on the current field has been made between the existing electrode layouts and a new transmodiolar electrodes. The new model can take either current or voltage as input for each electrode to simulate the common ground and multi-mode grounding modes. By coupling the surface model with lumped capacitor and constant phase element models, time domain simulation of the stimulation waveform has also been achieved. To validate the simulation results and calibrate the parameters of the model, in-situ and in-vitro measurements have been carried out with self-made devices. The recorded data proved the effectiveness of combining lumped components with the surface model.

In this dissertation, we introduce and analyse a scheme for the fast evaluation of integrals stemming from boundary element methods including discretisations of the classical single and double layer potential operators. Our method is based on the parametrisation of boundary elements in terms of a d-dimensional parameter tuple. We interpret the integral as a real-valued function f depending on d parameters and show that f is smooth in a d-dimensional box. A standard interpolation of f by polynomials leads to a d-dimensional tensor which is given by the values of f at the interpolation points. This tensor may be approximated in a low rank tensor format like the canonical format or the hierarchical format. The tensor approximation has to be done only once and allows us to evaluate interpolants in O(dr(m+1)) operations in the canonical format, or O(dk³ + dk(m + 1)) operations in the hierarchical format, where m denotes the interpolation order and the ranks r, k are small integers. In particular, we apply an efficient black box scheme in the hierarchical tensor format in order to adaptively approximate tensors even in high dimensions d with a prescribed (but heuristic) target accuracy. By means of detailed numerical experiments, we demonstrate that highly accurate integral values can be obtained at very moderate costs.

This report presents the work done within the framework of my master thesis in the program Infrastructure Engineering at KTH Royal Institute of Technology, Stockholm. This project has been proposed and sponsored by the French company Setec TPI, part of the Setec group, located in Paris. The overall goal of this study is to investigate fluid-structure interaction and particularly sloshing in liquid-containing structures subjected to seismic or other dynamic action. After a brief introduction, the report is composed of three main chapters. The first one presents and explains fluid-structure interaction equations. Fluid-structure interaction problems obey a general flow equation and several boundary conditions, given some basic assumptions. The purpose of the two following chapters is to solve the corresponding system of equations. The first approach proposes an analytical solution: the problem is solved for 2D rectangular tanks. Different models are considered and compared in order to analyze and describe sloshing phenomenon. Liquid can be decomposed in two parts: the lower part that moves in unison with the structure is modeled as an impulsive added mass; the upper part that sloshes is modeled as a convective added mass. Each of these two added mass creates hydrodynamic pressures and simple formulas are given in order to compute them. The second approach proposes a numerical solution: the goal is to be able to solve the problem for any kind of geometry. The differential problem is resolved using a singularity method and Gauss functions. It is stated as a boundary integral equation and solved by means of the Boundary Element Method. The linear system obtained is then implemented on Matlab. Scripts and results are presented. Matlab programs are run to solve fluid-structure interaction problems in the case of rectangular tanks: the results concur with the analytical solution which justifies the numerical solution. This report gives a good introduction to sloshing phenomenon and gathers several analytical solutions found in the literature. Besides, it provides a Matlab program able to model effects of sloshing in any liquid-containing structures.

A methodology is formulated for the solution of the inverse problem concerned with the reconstruction of multi-dimensional heat fluxes for film cooling applications. The motivation for this study is the characterization of complex thermal conditions in industrial applications such as those encountered in film cooled turbomachinery components. The heat conduction problem in the metal endwall/shroud is solved using the boundary element method (bem), and the inverse problem is solved using a genetic algorithm (ga). Thermal conditions are overspecified at exposed surfaces amenable to measurement, while the temperature and surface heat flux distributions are unknown at the film cooling hole/slot walls. The latter are determined in an iterative process by developing two approaches. The first approach, developed for 2d applications, solves an inverse problem whose objective is to adjust the film cooling hole/slot wall temperatures and heat fluxes until the temperature and heat flux at the measurement surfaces are matched in an overall heat conduction solution. The second approach, developed for 2d and 3d applications, is to distribute a set of singularities (sinks) at the vicinity of the cooling slots/holes surface inside a fictitious extension of the physical domain or along cooling hole centerline with a given initial strength distribution. The inverse problem iteratively alters the strength distribution of the singularities (sinks) until the measuring surfaces heat fluxes are matched. The heat flux distributions are determined in a post-processing stage after the inverse problem is solved. The second approach provides a tremendous advantage in solving the inverse problem, particularly in 3d applications, and it is recommended as the method of choice for this class of problems. It can be noted that the ga reconstructed heat flux distributions are robust, yielding accurate results to both exact and error-laden inputs. In all cases in this study, results from experiments are simulated using a full conjugate heat transfer (cht) finite volume models which incorporate the interactions of the external convection in the hot turbulent gas, internal convection within the cooling plena, and the heat conduction in the metal endwall/shroud region. Extensive numerical investigations are undertaken to demonstrate the significant importance of conjugate heat transfer in film cooling applications and to identify the implications of various turbulence models in the prediction of accurate and more realistic surface temperatures and heat fluxes in the cht simulations. These, in turn, are used to provide numerical inputs to the inverse problem. Single and multiple cooling slots, cylindrical cooling holes, and fan-shaped cooling holes are considered in this study. The turbulence closure is modeled using several two-equation approach, the four-equation turbulence model, as well as five and seven moment reynolds stress models. The predicted results, by the different turbulence models, for the cases of adiabatic and conjugate models, are compared to experimental data reported in the open literature. Results show the significant effects of conjugate heat transfer on the temperature field in the film cooling hole region, and the additional heating up of the cooling jet itself. Moreover, results from the detailed numerical studies presented in this study validate the inverse problem approaches and reveal good agreement between the bem/ga reconstructed heat fluxes and the cht simulated heat fluxes along the inaccessible cooling slot/hole walls

No description available.

No description available.