Cette thèse adresse des problèmes liés à la modélisation des écosystèmes microbiens dans les sols.Les modèles de transport réactif sont très performants pour modéliser la biogéochimie d'un milieu poreux variablement saturé en tenant compte de plusieurs paramètres physiques et géochimiques mais ils n'intègrent pas (ou sinon de facon très sommaire) les activités microbiologiques dans les sols. Le modèle mathématique du chemostat est couramment utilisé en écologie microbienne. Dans le but d'une représentation spatiale plus simplifiée que les modèles de transport réactif mais suffisamment pertinente pour rendre compte de phénomènes biologiques, cette thèse est une première tentative pour comprendre la fonction entrées-sorties dans un chemostat structuré, et étudier comment une structure spatiale peut modifier cette fonction, et quels sont les paramètres clés. / The aim of this thesis is to study some problems related to the modeling of microbial ecosystems in soil. Considering different physical and biogeochemical parameters, the reactif transport models are very performing to represent the biochemistry of a variably saturated porous media. In contrast, these models don't integrate (or sparsely) the microbiological activities in the soils. The mathematical model of chemostat is frequently used in microbial ecology. For the purpose of a more simplified spatial representation than the reactif transport models, but sufficiently relevant to represent biological phenomena, this thesis is a first attempt to understand the "input-output" function in a structured chemostat, and to study how a special structure can alter this function, and what are the key parameters.

L’ablation partielle du foie est une chirurgie qui intervient dans le traitement des lésions du foie et lors d’une transplantation partielle de foie. Les relations entre l’hémodynamique du foie, son volume et ses fonctions restent à élucider pour mieux comprendre les causes des complications de ces chirurgies. Lors de la chirurgie, l’hémodynamique du foie est altérée suite à l’augmentation de la résistance au flux sanguin de l’organe. La régénération du foie semble dépendante des changements de débit et de pression dans la veine porte. D’autre part, comme le foie reçoit 25% du débit cardiaque, la chirurgie impacte la circulation sanguine globale.
 Dans ce contexte, le premier objectif est de mieux comprendre, grâce à des modèles mathématiques, l’influence de l’hépatectomie sur l’hémodynamique. Le second objectif est l’analyse de la perfusion et de la fonction du foie. Premièrement, la procédure chirurgicale, les conditions expérimentales ainsi que les mesures obtenues sont détaillées.
 Ensuite, les valeurs moyennes mesurées lors de douze chirurgies sont reproduites par un modèle de circulation entière, basé sur des équations différentielles ordinaires. Lors des différentes hépatectomies, des changements de forme de courbe sont observés. Un modèle de circulation entière, basée sur des équations 1D et 0D est proposé pour analyser ces changements. Ce travail pourrait permettre une meilleure compréhension des changements d’architecture du foie induits par l’hépatectomie.
 Puis, le transport dans le sang d’un composé ainsi que son traitement par le foie sont modélisés. Un modèle pharmacocinétique est développé et grâce aux mesures, les paramètres du modèle sont estimés. / Major liver resection is being performed to treat liver lesions or for adult-to-adult living donor liver transplantation. Complications of these surgeries are related to a poor liver function. The links between liver hemodynamics, liver volume and liver function remain unclear and are important to better understand these complications. The surgery increases the resistance to blood flow in the organ, therefore it modifies liver hemodynamics. Large modifications of the portal vein hemodynamics have been associated with poor liver regeneration. Moreover the liver receives 25% of the cardiac outflow, therefore liver surgery may impact the whole blood circulation. In this context, the first goal is to investigate with mathematical models the impact of liver surgery on liver hemodynamics. The second goal is to study the liver perfusion and function with mathematical models. The first part describes the experimental conditions and reports the measurements recorded. Then, the second part focuses on the liver hemodynamics during partial hepatectomy. On one hand, the hemodynamics during several surgeries is quantitatively reproduced and explained by a closed-loop model based on ODE. On the other hand, the change of waveforms observed after different levels of liver resection is reproduced with a model of the global circulation, including 0D and 1D equations. This may contribute to a better understanding of the change of liver architecture induced by hepatectomy. Next, the transport in blood of a compound is studied. And a pharmacokinetics model and its parameter identification are developed to quantitatively analyze indocyanine green fluorescence dynamics in the liver tissue.

This work is focused on computer simulations of fractal capacitors. The geometry of capacitors and its influence is investigated. Simulations are realized in programs Matlab, SolidWorks and Comsol Multiphysics. There are also several specific examples of different geometrics of capacitors their comparisons and assessment.

La méso-combustion est le régime de combustion où la taille caractéristique du domaine est juste supérieure à la distance de coincement de la flamme , typiquement de l'ordre du centimètre. La difficile réalisation de systèmes de combustion fonctionnant en ce régime de flamme particulier suscite l'intérêt : il devient alors possible de tirer parti de la haute densité énergétique des hydrocarbures pour concevoir des systèmes de production d'énergie plus compacts. Nous nous intéressons à la réalisation de calculs DNS compressibles d'une chambre de combustion cubique de 8 x 10 x 8 mm3. Ce travail présente autant la mise en œuvre des calcules, en particulier la problématique de la condition frontière mur, que les résultats obtenus. Ces derniers nous permettent d'analyser la phénoménologie complexe de cet écoulement réactif confiné et serviront de base à des modélisations futures. / Meso-combustion can be defined as the combustion regime where the involved lenghts scales are close but slightly larger than the quenching distance of the flame, tipically smaller than a cm. By taking advantage of the high energetic density of liquid hydrocarbons, it would become possible to build small-sized combustion-based long-lived lighter electrical power systems. However combsution phenomena at these meso-scales have their own shortcomings. Indeed, by decreasing the system size, the usual phenomenological balance betwenne chemical reactions, mixing, turbulence and heat transfer is changed. In the present work, we focus on the DNS calculation of a cubic meso-combsution chamber of 8 x 10 x 8 mm3. This works presents the implementation of the numerical strategy used, with a specific attention to the no-slip wall compressible boundary condition. We then present an analysis of this particular reactive flow. The results are useful for future modeling of such a combustor.

Orientador: José Manoel Balthazar / Resumo: A microscopia da força atômica (AFM) com medidas em nanoescala, vêm da nanotecnologia e constitui atualmente um campo de pesquisa multidisciplinar. A proposta de pesquisa apresentada tem como objetivo contribuir a pesquisa cientifica sobre AFM considerando que sistema está operando no modo intermitente e em meio líquido, e que a dinâmica do amortecimento do meio líquido é determinado por derivadas de ordem fracionária. Através de ferramentas de simulação numérica serão realizadas a análise dinâmica do sistema para os casos em que o sistema tenha comportamento caótico. Para supressão deste comportamento caótico serão utilizadas e analisadas duas estratégias de controle, o controle SDRE (Estados Dependentes da Equação de Riccati) e o controle OLFC (Controle Linear por Realimentação Ótimo). / Abstract: The atomic force microscope (AFM) with nanoscale measurements comes from nanotechnology and is currently a multidisciplinary field of research. The present research proposal aims to contribute to scientific research on AFM considering that the system is operating in the intermittent mode and in liquid medium, and that the damping dynamics of the liquid medium is represented by fractional order derivatives. Through numerical simulation tools, the dynamic analysis of the system will be performed in cases where the system has chaotic behavior. To suppress this chaotic behavior we will use and analyze two control strategies, the SDRE (Riccati Equation Dependent States) and OLFC (Linear Control for Optimum Feedback). / Doutor

In this thesis, we present results of a numerical modelling study focused on the thermal evolution of the Earth and terrestrial planets. We focus particularly on two problems: I) constraining the internal structure of Venus and Mercury using their geoid and surface topography data and II) evaluating the effects of a rhe- ologically distinct post-perovskite on the secular cooling of the Earth. In part I, we performed simulations in a broad group of models of the Venusian man- tle, characterised by different rheological descriptions, and we compared spectra of their geoid and their surface topography with the observed quantities. Our analysis suggested that the geoid and the surface topography of Venus are con- sistent with a radially symmetric viscosity model with a strong 200 km thick lithosphere, without an asthenosphere and with a gradual viscosity increase in the underlying mantle. In the case of Mercury, none of our models was able to predict observed data, thus suggesting other than a dynamic origin of observed geoid and topography. In part II, we investigated style of Earth's mantle con- vection and its long-term evolution in the models that take into account a weak post-perovskite. We conclude that the presence of the weak post-perovskite en- hances the core cooling. This effect is comparable in...

The development in structural design and construction of high-speed craft has been extensive during the last decades. Environmental and economical issues have increased the need to develop more optimized structures, using new material concepts, to reduce weight and increase performance efficiency. However, both lack of, and limitations in design methodology, makes this a difficult task. In this thesis a methodological framework which enables detailed studies of the slamming loads and associated responses for high-speed planing craft in irregular waves is established. The slamming loads can either be formulated based on numerical simulations, or on experimental measurements and pressure distribution reconstruction. Structure responses are derived in the time-domain using finite element analysis. Statistical methods are used to determine design loads and lifetime extreme responses. The framework is applied to perform phenomenological studies of the slamming loading conditions for high-speed craft, and used to highlight and quantify the limitations in the prevailing semi-empirical method for design load determination with respect to slamming. A number of clarifications regarding the original derivation and the applicability of the prevailing semi-empirical method are presented. Finally, several potential improvements to the method are presented and the associated implications discussed. The long-term goal of the research project is to establish a method for direct calculation of loads and response for high-speed planing craft, which can enable design of truly efficient craft structures. The methodology and the results presented in this thesis are concluded to be important stepping-stones towards this goal. / <p>In page VII, Paper B is wrong title. The correct title is "Experiental Evaluation of Slamming Pressure Models Used in Structural Design of High-Speed Craft". QC 20130228</p>

Automotive parts may be the cause of very annoying friction-induced noise and the source of many customer complaints. Indeed, when a wiper operates on a windshield, vibratory phenomena may appear due to flutter instabilities and may generate squeal noise. As squeal noise generated by wiper system is a random and complex phenomenon, there are only few studies dealing with the wiper noise. The complexity of this phenomenon is due to the cinematic of the movement and to the various environmental parameters which have an influence on the appearance of the noise. This master thesis is a research and development project and presents a numerical simulation methodology used in the aim to reduce and eradicate squeal noise of wiper systems.  In the first part, the finite element model representing a wiper system and the numerical simulation methodology will be presented in detail. In the second part, stability analysis will be carried out in nominal studies and in designs of experiments. Parametric studies will also be achieved to understand the behavior and the influence of each considered input parameters. Two wiper blades, with the same geometry but with different material, will be considered for the different studies. These two wiper blades will be examined to figure out when squeal noises appear.

<p>U tezi je predstavljena implementacija metode teorija funkcionala gustine (DFT) bazirana na metodi za sklapanje naelektrisanja (CPM) koja koristi bazise gausijanskih funkcija. Metod je baziran na pretpostavci da se elektronska gustina naelektrisanja velikih sistema, može predstaviti kao suma doprinosa pojedinačnih atoma, takozvanih motiva gustine naelektrisanja, koji se dobijaju računanjem malog prototip sistema. Talasna funkcija,<br />kao i gustina naelektrisanja, se u na&scaron;oj implementaciji reprezentuju uz pomoć bazise gausijanskih funkcija, dok se motivi opisuju kori&scaron;ćenjem prostornih koordinata. Uz pomoć procedure za minimizaciju se iz motiva opisanih koordinatama, dobija gustina naelektrisanja predstavljena u bazisu Gausijana. Implementacija serijskog programa pokazuje značajno pobolj&scaron;anje u performansama u odnosu na prethodne implementacije bazirane na ravnim talasima. Ova implementacija re&scaron;ava sistem od približno 1000 atoma na jednom procesorskom jezgru za svega nekoliko sati. Paralelna implementacija uz pomoć naprednih metoda paralelizacije i distribucije podataka omogućava re&scaron;avanje sistema od vi&scaron;e desetina hiljada atoma. Najveći testirani sistem ima približno<br />20000 atoma i testiran je na 256 paralelnih procesa.</p> / <p>We present the implementation of the density functional theory (DFT) based charge patching method (CPM) using the basis of Gaussian functions. The method is based on the assumption that the electronic charge density of a large system is the sum of contributions of individual atoms, so called charge density motifs, that are obtained from calculations of small prototype systems.In our implementation wave functions and electronic charge density are represented using the basis of Gaussian functions, while charge density motifs are represented using a real space grid. A constrained minimization procedure is used to obtain Gaussian basis representation of charge density from real space representation of motifs. The code based on this&nbsp; implementation exhibits superior performance in comparison to previous implementation of the charge patching method using the basis of plane waves. It enables calculations of electronic structure of systems with around 1000 atoms on a single CPU core with computational time of just several hours. The parallel implementation enables calculations for the system with more than ten thousand atoms. The largest system tested has around 20000 atoms and was computed on 256 parallel processes.</p>

Many systems designed by electrical & computer engineers rely on electromagnetic (EM) signals to transmit, receive, and extract either information or energy. In many cases, these systems are large and complex. Their accurate, cost-effective design requires high-fidelity computer modeling of the underlying EM field/material interaction problem in order to find a design with acceptable system performance. This modeling is accomplished by projecting the governing Maxwell equations onto finite dimensional subspaces, which results in a large matrix equation representation (Zx = b) of the EM problem. In the case of integral equation-based formulations of EM problems, the M-by-N system matrix, Z, is generally dense. For this reason, when treating large problems, it is necessary to use compression methods to store and manipulate Z. One such sparse representation is provided by so-called H^2 matrices. At low-to-moderate frequencies, H^2 matrices provide a controllably accurate data-sparse representation of Z. The scale at which problems in EM are considered ``large'' is continuously being redefined to be larger. This growth of problem scale is not only happening in EM, but respectively across all other sub-fields of computational science as well. The pursuit of increasingly large problems is unwavering in all these sub-fields, and this drive has long outpaced the rate of advancements in processing and storage capabilities in computing. This has caused computational science communities to now face the computational limitations of standard linear algebraic methods that have been relied upon for decades to run quickly and efficiently on modern computing hardware. This common set of algorithms can only produce reliable results quickly and efficiently for small to mid-sized matrices that fit into the memory of the host computer. Therefore, the drive to pursue larger problems has even began to outpace the reasonable capabilities of these common numerical algorithms; the deterministic numerical linear algebra algorithms that have gotten matrix computation this far have proven to be inadequate for many problems of current interest. This has computational science communities focusing on improvements in their mathematical and software approaches in order to push further advancement. Randomized numerical linear algebra (RandNLA) is an emerging area that both academia and industry believe to be strong candidates to assist in overcoming the limitations faced when solving massive and computationally expensive problems. This thesis presents results of recent work that uses a random sampling method (RSM) to implement algebraic operations involving multiple H^2 matrices. Significantly, this work is done in a manner that is non-invasive to an existing H^2 code base for filling and factoring H^2 matrices. The work presented thus expands the existing code's capabilities with minimal impact on existing (and well-tested) applications. In addition to this work with randomized H^2 algebra, improvements in sparse factorization methods for the compressed H^2 data structure will also be presented. The reported developments in filling and factoring H^2 data structures assist in, and allow for, the further pursuit of large and complex problems in computational EM (CEM) within simulation code bases that utilize the H^2 data structure.