Risers são tubos que transportam fluidos do fundo do mar até as plataformas flutuantes e vice-versa. Diversas configurações e materiais são utilizados, porém apenas os steel catenary risers (SCR) são estudados neste trabalho. Os risers são estruturas extremamente esbeltas e, por isso, grande parte de seu trecho suspenso tem comportamento de cabo. Apenas em duas regiões a rigidez flexional é relevante: no hang-off (topo) e na touch-down zone (TDZ), sendo esta última a região mais complexa para análise devido ao contato unilateral com o solo. Em função dos diversos carregamentos dinâmicos a que o riser é submetido, grandes variações na curvatura ocorrem na TDZ, além de impacto e atrito com o solo, que podem reduzir a vida útil da estrutura ou até mesmo por em risco a sua integridade. Por estas razões, este trabalho visa à elaboração de uma metodologia analítica para a construção de um modelo de ordem reduzida (MOR) capaz de analisar o comportamento dinâmico não linear da TDZ de um SCR. Como na TDZ a rigidez flexional predomina sobre a rigidez geométrica, o riser é modelado como uma viga semi-infinita, tendo uma parte suspensa e outra apoiada sobre solo hipoteticamente elástico com contato unilateral. Na extremidade suspensa são impostos deslocamentos verticais e trações dinâmicas que fazem com que a posição do touch-down point (TDP) também varie com o tempo. Trata-se, portanto, de um problema com condições de contorno móveis. A metodologia adotada para a resolução deste problema foi transformá-lo em um problema de condições de contorno fixas por meio de uma transformação de variáveis. Contudo, paga-se um preço por tal transformação, e fortes não linearidades surgem na equação diferencial de movimento, tornando-a extremamente complexa para uma resolução analítica direta. Para o problema de flexão simples, consegue-se obter os modos de vibração não lineares através do método das múltiplas escalas. De posse destes modos, utiliza-se o método de Galerkin não linear para projetar a equação completa em um modo escolhido, transformando o modelo contínuo em um modelo de ordem reduzida com apenas um grau de liberdade, cuja coordenada generalizada modal é o deslocamento horizontal do TDP. Obtida a equação do MOR, nota-se que existem coeficientes que variam com o tempo, como na clássica equação de Mathieu, indicando a possibilidade de ocorrer ressonância paramétrica. Neste tipo de ressonância, entre outras possibilidades, pode ocorrer que a frequência de excitação seja o dobro da frequência natural trata-se da ressonância paramétrica principal. A equação do MOR é integrada numericamente e suas respostas são comparadas com as respostas obtidas por modelos de elementos finitos elaborados em softwares comerciais, como o Abaqus e o Orcaflex. Por fim, discutem-se as potencialidades e limitações do MOR, sendo uma grande vantagem a possibilidade de processar diversos casos facilmente, variando os parâmetros que influem nas respostas. Com este mapeamento das respostas, é possível estimar as amplitudes dos estados estacionários pós-críticos. / Risers are pipes that convey fluids from the seabed up to the floating platforms and vice-versa. Many configurations and materials are used, but only steel catenary risers (SCR) are studied in this work. Risers are extremely slender structures, and for this reason, most of the suspended part has cable behavior. Only in two regions the bending stiffness is important: at the hang-off and at the touch-down zone (TDZ), which is the most complex region for analysis because of the unilateral contact with the seabed. Due to several dynamic loads that the riser is subjected to, great curvature variations occur at the TDZ, apart from impacts and friction with the soil, which can reduce the life time of the structure or even jeopardize its integrity. For these reasons, this work aims at the development of an analytical methodology for the construction of a reduced-order model (ROM) able to analyze the nonlinear dynamic behavior of the TDZ of a SCR. As at the TDZ the bending stiffness prevails over the geometric stiffness, the riser is modeled as a semi-infinite beam, having a suspended part and another one resting on the elastic soil with unilateral contact. At the end of the suspended part, vertical displacements and dynamic tensions are imposed, that cause the TDPs position to vary with time. It is, therefore, a problem with moving boundary conditions. The methodology adopted for solving this problem was to transform it into a problem with fixed boundary conditions via a variable transformation. However, a price is paid for such a transformation, and strong nonlinearities appear in the differential equation of motion, making it extremely complex to solve analytically. For the simple bending problem, nonlinear vibration modes are obtained via the method of multiple scales. In possession of these modes, the nonlinear Galerkin method is used to project the complete equation into a chosen mode, transforming the continuum model into a reduced-order model (ROM) with only one degree of freedom whose modal generalized coordinate is the horizontal displacement of the TDP. After obtaining the ROM, it is noticed that there are coefficients that vary with time, as in the classic Mathieu equation, indicating the possibility of parametric resonance. In this kind of resonance, among other possibilities, the excitation frequency may be twice the natural frequency it is the so-called principal parametric resonance. The ROMs equation is integrated numerically and the responses are compared to those given by finite-element models studied with the help of commercial softwares, like Abaqus and Orcaflex. Finally, the potentialities and limitations of the ROM are discussed. One of the advantages is the possibility of processing several cases easily, changing the parameters that affect the responses. With this response mapping, it is possible to estimate the post-critical steady-state amplitudes that take place.

O estudo da dinâmica estrutural de risers oceânicos apresenta instigantes desafios aos pesquisadores da área da engenharia de estruturas, uma vez que os meios tradicionais de análises dinâmicas lineares nem sempre se ajustam às suas complexas particularidades. No atual estágio do desenvolvimento científico da área de engenharia de estruturas, a aplicação de técnicas de análise dinâmica não linear, dentro de determinadas hipóteses, mostra-se como uma das alternativas possíveis e viáveis à tradicional análise dinâmica linear. Com vistas a uma nova abordagem do problema, o presente trabalho adota uma metodologia de análise não linear dinâmica de risers oceânicos em configuração de lançamento de catenária, conjugada a uma técnica de processamento de Modelos de Ordem Reduzida para o estudo dos fenômenos dinâmicos manifestados por risers. Trata-se de um método de modelagem local, restrito à região de contato unilateral do riser com o solo, considerado este último um meio viscoelástico. Os resultados da aplicação desta metodologia são demonstrados nos estudos de caso apresentados com comparações com modelos numéricos (Método dos Elementos Finitos) e modelos físicos. / The dynamic study of offshore risers still demands large efforts from structural engineering researchers, since these systems may behave in a way that is not well modeled and understood using simply linear dynamic theories. Nevertheless, the current development stage of non linear dynamic theories gives hope that their use for the analyses of such systems can be of great value, even though, this must be carefully done specially by the analyst. The present work refers to a non linear dynamic methodology application to offshore risers, particularly steel catenary risers, by a technique known as reduced-order modeling, in the study of dynamic phenomena that these structures may present. The model is local, which means that it represents the touch-down zone of the riser-soil system. The soil modeling was presumed to be viscoelastic. The results obtained in case studies are compared with those from numerical (Finite Element Method) and small scale physical models.

D'un point de vue industriel, la mise en place de nouvelles architectures de systèmes mécaniques nécessite un long processus de conception permettant de définir et d'anticiper le comportement. Dans le cas particulier des systèmes aéronautiques tels que les moteurs d'avions, un certain nombre de pièces sont particulièrement sensibles car elles doivent répondre à des impératifs stricts en termes d'encombrement, de performance et de tenue mécanique. Dans ce contexte, la prévision du comportement vibratoire revêt une importance particulière puisqu'elle permet d'évaluer le niveau des sollicitations cycliques appliquées sur le système et guide ainsi la détection en amont d'éventuels problèmes de fatigue des matériaux. La plupart du temps, des modèles numériques sont utilisés pour représenter les structures, et le comportement est simulé en résolvant un ensemble d'équations. Pour atteindre un niveau de détail répondant au besoin industriel, ces modèles peuvent être particulièrement gros, et la résolution des équations associées demande des ressources et des temps de calcul considérables. De plus, pour rendre compte au mieux des comportements observés expérimentalement, il est souvent nécessaire de prendre en compte des phénomènes non-linéaires, ce qui augmente encore la difficulté. Les travaux présentés dans ce manuscrit concernent cette problématique du comportement vibratoire des structures non-linéaires et s'orientent autour de deux axes : la réduction de modèle et le calcul des solutions multiples. L'objectif du premier axe est de contribuer à la construction de modèles numériques non linéaires réduits utilisables en conception de systèmes industriels et de proposer des outils d'exploitation et d'interprétation de ces modèles. En particulier, on considère le cas des méthodes de projection de Galerkin et on montre qu'elles sont à même de construire des modèles réduits réalistes. Des méthodes complémentaires de réduction de modèles sont également présentées dans le cas particulier de la recherche de solutions par la méthode de la balance harmonique (HBM) : on s'intéressera en particulier à des méthodes de sélection d'harmoniques. Après avoir comparé les différentes méthodes proposées sur un exemple simple de poutre non-linéaire, elles sont appliquées à un modèle de structure industrielle représentant une aube d'hélice d'open rotor. Le second axe de ces travaux concerne le calcul de solutions multiples pour les systèmes dynamiques non-linéaires. Une particularité de ces systèmes est en effet de présenter plusieurs configurations stables pour un état de sollicitation donné. Il s'agira ici de proposer des méthodes de calcul permettant de dresser la liste exhaustive des solutions possibles. Le travail présenté se concentre sur la recherche de solutions périodiques par la méthode de la balance harmonique pour des systèmes possédant des non-linéarités polynomiales. Ces restrictions conduisent à la résolution de systèmes polynomiaux pour lesquels il existe des méthodes permettant de calculer l'ensemble des solutions. En particulier, on propose l'utilisation originale de méthodes basées sur le calcul de bases de Groebner pour la résolution de systèmes polynomiaux issus de la mécanique. Les différentes méthodes présentées sont illustrées et comparées sur des exemples simples. Les résultats montrent que même pour des systèmes simples, le comportement dynamique peut être très complexe. / In an industrial context, the design of new mechanical systems requires long design processes in order to define and to anticipate the behavior of all the constitutive parts. In the particular case of aeronautical structures such as plane engines, design is especially critical since they have to meet various and strict needs (life duration, performances . . .). Then, anticipating vibratory behavior is very important as this provides information about cyclic solicitations and fatigue. Most often, numerical models are used to mimic the structure and mechanical behavior is simulated by solving a set of differential equations. In the case of industrial structures, such models can be quite large and their resolution very time-consuming. Moreover, in order to model experimental behavior realistically, it is often necessary to take nonlinear phenomena into account and thus increase the required computational effort. The work presented in this PhD deals with the study of mechanical nonlinear systems. It focuses on two principal directions : model reduction and multiple solutions computation. The goal of the first direction is to contribute to the building of numerical reduced order models usable in industrial context and to propose tools to exploit an interpret them. Particularly, Galerkin projection methods are investigated in the context of nonlinear systems reduction, showing that those methods are, under certain conditions, able to give a reliable picture of full system behavior. In the case of the harmonic balance method, complementary methods are also proposed to reduce the size of the algebraic equations system by using harmonic selection techniques. The presented methods are firstly illustrated and compared on a simple nonlinear beam example ; they are then applied to an industrial model of open rotor blade. The second direction of this work deals with the computation of multiple solutions arising in nonlinear dynamical systems. Indeed, it has been shown that such systems can present different stable configurations for a given solicitation. The objective here is to provide tools for computing such multiple solutions. We only consider the case of periodic solutions for systems with polynomial nonlinearities, treated with harmonic balance method. These hypotheses enable one to search for multiple states as solutions of polynomial algebraic systems of equations, for which some methods exist to compute the entire set of solutions. In particular, we propose to use methods relying on Groebner basis computation, in order to compute the whole set of solutions. The proposed methods are illustrated and compared on simple examples, showing that even such simple systems can present very complex dynamical behavior.

In this thesis, the method of Proper Orthogonal Decomposition (POD) is implemented to construct approximate, reduced order models (ROM) of mesh movement methods. Three mesh movement algorithms are implemented and comparatively evaluated, namely radial basis function interpolation, mesh optimization and elastic deformation. POD models of the mesh movement algorithms are constructed using a series of system observations, or snapshots of a given mesh for a set of boundary deformations. The scalar expansion coefficients for the POD basis modes are computed in three different ways, through coefficient optimization, Galerkin projection of the governing set of equations and coefficient interpolation. It is found that using only coefficient interpolation yields mesh movement models that accurately approximates the full order mesh movement, with CPU cost savings in excess of 99%. We further introduce a novel training procedure whereby the POD models are generated in a fully automated fashion. The technology is applicable to any mesh movement method and enables potential reductions of up to four orders of magnitude in mesh movement related costs. The proposed model can be implemented without having to pre-train the POD model, to any fluid-structure interaction code with an existing mesh movement scheme. Copyright / Dissertation (MEng)--University of Pretoria, 2010. / Mechanical and Aeronautical Engineering / unrestricted

La structure de l’oeil permet d’observer la microcirculation, grâce aux caméras de fond d’oeil. Ces appareils sont bon marché et couramment utilisés dans la pratique clinique, permettant le dépistage de maladies oculaires. La capacité des vaisseaux à adapter leur diamètre (autorégulation) afin de réguler le débit sanguin est importante dans la microcirculation. L’hémodynamique de l’oeil est impactée par la pression à l’intérieur du globe oculaire (IOP), qui est à son tour influencée par le flux sanguin oculaire. Les altérations de l’autorégulation et l’IOP jouent un rôle dans les maladies oculaires. La modélisation mathématique peut aider à interpréter l’interaction entre ces phénomènes et à mieux exploiter les données médicales disponibles. Dans la première partie, nous présentons un modèle simplifié d’interaction fluidestructure qui inclut l’autorégulation, appliqué à un reseau 3D obtenu par imagerie médicale. Les cellules musculaires lisses regulant le diamètre du vaisseau sont modélisés dans la structure. Ensuite, nous utilisons des équations de poroélasticité pour décrire le flux sanguin dans la choroïde, dans un modèle multi-compartiments de l’oeil. Cette approche permet de rendre compte de la transmission de la pulsatilité de la choroïde à la chambre antérieure, où l’IOP est mesurée. Nous présentons des résultats préliminaires sur la choroïde, l’humeur aqueuse et sur la choroïde couplée avec la vitrée. Enfin, nous présentons un modèle d’ordre réduit pour accélérer des simulations multi-physique. Des modèles de haute précision sont utilisés pour les compartiments d’intérêt et une représentation réduite de l’opérateur de Steklov-Poincaré est utilisée pour les autres compartiments. / The structure of the eye offers a unique opportunity to directly observe the microcirculation, by means, for instance, of fundus camera, which are cheap devices commonly used in the clinical practice. This can facilitate the screening of systemic deseases such as diabetes and hypertension, or eye diseases such as glaucoma. A key phenomenon in the microcirculation is the autoregulation, which is the ability of certain vessels to adapt their diameter to regulate the blood flow rate in response to changes in the systemic pressure or metabolic needs. Impairments in autoregulation are strongly correlated with pathological states. The hemodynamics in the eye is influenced by the intraocular pressure (IOP), the pressure inside the eye globe, which is in turn influenced by the ocular blood flow. The interest in the IOP stems from the fact that it plays a role in several eye-diseases, such as glaucoma. Mathematical modelling can help in interpreting the interplay between these phenomena and better exploit the available data. In the first part of the thesis we present a simplified fluid-structure interaction model that includes autoregulation. A layer of fibers in the vessel wall models the smooth muscle cells that regulate the diameter of the vessel. The model is applied to a 3D image-based network of retinal arterioles. In the second part, we propose a multi-compartments model of the eye. We use the equations of poroelasticity to model the blood flow in the choroid. The model includes other compartments that transmit the pulsatility from the choroid to the anterior chamber, where the measurements of the IOP are actually performed. We present some preliminary results on the choroid, the aqueous humor and on the choroid coupled with the vitreous. Finally, we present a reduced order modelling technique to speed up multiphysics simulations. We use high fidelity models for the compartments of particular interest from the modelling point of view. The other compartments are instead replaced by a reduced representation of the corresponding Steklov-Poincaré operator.

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.

The number of electrical machines used in modern road-vehicles is continuously increasing to meet regulatory requirements regarding safety and efficiency, as well as consumer expectations in terms of comfort. For auxiliary applications, such as cooling fan or pumps, permanent-magnet synchronous machines (PMSMs) are extensively used owing to their high power density. This thesis focuses on the modelling aspects of PMSMs, with a particular focus on finite-element and reduced-order models to be used in system-level simulations. 2-D and 3-D parametric finite-element (FE) models are developed, allowing to compute irreversible demagnetization in addition to the standard quantities such as torque, back electromotive force and flux-linkages. The effects of magnet overhang on the performance of an interior PMSM is briefly discussed. Using the FE model, a reduced-order lookup-table (LUT) based electromagnetic model, having similar accuracy as FE analysis, is then developed. Coupled to a mechanical state-space representation obtained from a modal FE analysis, the final model allows to compute electromagnetic induced vibrations under pulse width modulation supply. The validation of the complete workflow is carried out on a 12slot-10pole PMSM used to drive a cooling fan. After fitting the damping coefficient in the structural state-space model, the results are in agreement with the experimental results. Due to the usage of LUTs, the simulation time is low compared to a pure FE analysis. This allows the model to be used to optimize low noise control strategies. To conclude this thesis, the parametric FE model is used in an optimization routine to minimize the cost and vibrations of the motor, whilst satisfying the working points. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

京都大学 / 新制・課程博士 / 博士(工学) / 甲第23511号 / 工博第4923号 / 新制||工||1769(附属図書館) / 京都大学大学院工学研究科電気工学専攻 / (主査)教授 松尾 哲司, 教授 雨宮 尚之, 准教授 久門 尚史 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Membrane Distillation (MD) is an emerging sustainable desalination technique. While MD has many advantages and can be powered by solar thermal energy, its main drawback is the low water production rate. However, the MD process has not been fully optimized in terms of its manipulated and controlled variables. This is largely due to the lack of adequate dynamic models to study and simulate the process. In addition, MD is prone to membrane fouling, which is a fault that degrades the performance of the MD process. This work has three contributions to address these challenges. First, we derive a mathematical model of Direct Contact Membrane Distillation (DCMD), which is the building block for the next parts. Then, the proposed model is extended to account for membrane fouling and an observer-based fouling detection method is developed. Finally, various control strategies are implemented to optimize the performance of the DCMD solar-powered process. In part one, a reduced-order dynamic model of DCMD is developed based on lumped capacitance method and electrical analogy to thermal systems. The result is an electrical equivalent thermal network to the DCMD process, which is modeled by a system of nonlinear differential algebraic equations (DAEs). This model predicts the water-vapor flux and the temperature distribution along the module length. Experimental data is collected to validate the steady-state and dynamic responses of the proposed model, with great agreement demonstrated in both. The second part proposes an extension of the model to account for membrane fouling. An adaptive observer for DAE systems is developed and convergence proof is presented. A method for membrane fouling detection is then proposed based on adaptive observers. Simulation results demonstrate the performance of the membrane fouling detection method. Finally, an optimization problem is formulated to maximize the process efficiency of a solar-powered DCMD. The adapted method is known as Extremum Seeking (ES). A Newton-based ES is designed and the proposed model is used to accurately forecast the distilled water flux. Although good results are obtained with this method, a practical modification to the ES scheme is proposed to enhance the practical stability.

Thermal energy storage (TES) modules are specifically designed to respond to transient thermal loading. Their dynamic response depends on the overall structure of the module, including module geometry and dimensions, the internal spatial distribution of phase change material (PCM) and conductive heat-spreading elements, and the thermophysical properties of the different materials composing the module. However, due to the complexity of analyzing a system’s dynamic thermal response to transient input signals, optimal design of a TES module for a particular application is challenging. Conventional design approaches are limited by (1) the computational cost associated with high fidelity simulation of heat transfer in nonlinear systems undergoing a phase transition and (2) the lack of model integration with robust optimization tools. To overcome these challenges, I derive reduced-order dynamic models of two different metal-PCM composite TES modules and validate them against a high fidelity CFD model. Through simulation and validation of both turbulent and laminar flow cases, I demonstrate the accuracy of the reduced-order models in predicting, both spatially and temporally, the evolution of the dynamic model states and other system variables of interest, such as PCM melt fraction. The validated models are used to conduct univariate and bivariate parametric studies to understand the effects of various design parameters on different performance metrics. Finally, a case study is presented in which the models are used to conduct detailed design optimization for the two HX geometries.