Objectives: The main objective of the current dissertation is to establish a formal shape optimization procedure for a given bypass grafts end-to-side distal anastomosis (ETSDA). The motivation behind this dissertation is that most of the previous ETSDA shape optimization research activities cited in the literature relied on direct optimization approaches that do not guaranty accurate optimization results. Three different ETSDA models are considered herein: The conventional, the Miller cuff, and the hood models. Materials and Methods: The ETSDA shape optimization is driven by three computational objects: a localized collocation meshless method (LCMM) solver, an automated geometry pre-processor, and a genetic-algorithm-based optimizer. The usage of the LCMM solver is very convenient to set an autonomous optimization mechanism for the ETSDA models. The task of the automated pre-processor is to randomly distribute solution points in the ETSDA geometries. The task of the optimized is the adjust the ETSDA geometries based on mitigation of the abnormal hemodynamics parameters. Results: The results reported in this dissertation entail the stabilization and validation of the LCMM solver in addition to the shape optimization of the considered ETSDA models. The LCMM stabilization results consists validating a custom-designed upwinding scheme on different one-dimensional and two-dimensional test cases. The LCMM validation is done for incompressible steady and unsteady flow applications in the ETSDA models. The ETSDA shape optimization include single-objective optimization results in steady flow situations and bi-objective optimization results in pulsatile flow situations. Conclusions: The LCMM solver provides verifiably accurate resolution of hemodynamics and is demonstrated to be third order accurate in a comparison to a benchmark analytical solution of the Navier-Stokes. The genetic-algorithm-based shape optimization approach proved to be very effective for the conventional and Miller cuff ETSDA models. The shape optimization results for those two models definitely suggest that the graft caliber should be maximized whereas the anastomotic angle and the cuff height (in the Miller cuff model) should be chosen following a compromise between the wall shear stress spatial and temporal gradients. The shape optimization of the hood ETSDA model did not prove to be advantageous, however it could be meaningful with the inclusion of the suture line cut length as an optimization parameter.

No description available.

Advances in numerical optimization have raised the possibility that efficient and novel aircraft configurations may be ``discovered'' by an algorithm. To begin exploring this possibility, a fast and robust set of tools for aerodynamic shape optimization is developed. Parameterization and mesh-movement are integrated to accommodate large changes in the geometry. This integrated approach uses a coarse B-spline control grid to represent the geometry and move the computational mesh; consequently, the mesh-movement algorithm is two to three orders faster than a node-based linear elasticity approach, without compromising mesh quality. Aerodynamic analysis is performed using a flow solver for the Euler equations. The governing equations are discretized using summation-by-parts finite-difference operators and simultaneous approximation terms, which permit nonsmooth mesh continuity at block interfaces. The discretization results in a set of nonlinear algebraic equations, which are solved using an efficient parallel Newton-Krylov-Schur strategy. A gradient-based optimization algorithm is adopted. The gradient is evaluated using adjoint variables for the flow and mesh equations in a sequential approach. The flow adjoint equations are solved using a novel variant of the Krylov solver GCROT. This variant of GCROT is flexible to take advantage of non-stationary preconditioners and is shown to outperform restarted flexible GMRES. The aerodynamic optimizer is applied to several studies of induced-drag minimization. An elliptical lift distribution is recovered by varying spanwise twist, thereby validating the algorithm. Planform optimization based on the Euler equations produces a nonelliptical lift distribution, in contrast with the predictions of lifting-line theory. A study of spanwise vertical shape optimization confirms that a winglet-up configuration is more efficient than a winglet-down configuration. A split-tip geometry is used to explore nonlinear wake-wing interactions: the optimized split-tip demonstrates a significant reduction in induced drag relative to a single-tip wing. Finally, the optimal spanwise loading for a box-wing configuration is investigated.

For the purpose of numerically simulating continuum mechanical structures, different types of material may be represented by the extreme values {<img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />,1}, where 0&lt;<img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /><img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cll" />1, of a varying coefficient <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> in the governing equations. The paramter <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> is not allowed to vanish in order for the equations to be solvable, which means that the exact conditions are approximated. For example, for linear elasticity problems, presence of material is represented by the value <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = 1, while <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> provides an approximation of void, meaning that material-free regions are approximated with a weak material. For acoustics applications, the value <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = 1 corresponds to air and <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> to an approximation of sound-hard material using a dense fluid. Here we analyze the convergence properties of such material approximations as <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />!0, and we employ this type of approximations to perform design optimization. In Paper I, we carry out boundary shape optimization of an acoustic horn. We suggest a shape parameterization based on a local, discrete curvature combined with a fixed mesh that does not conform to the generated shapes. The values of the coefficient <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" />, which enters in the governing equation, are obtained by projecting the generated shapes onto the underlying computational mesh. The optimized horns are smooth and exhibit good transmission properties. Due to the choice of parameterization, the smoothness of the designs is achieved without imposing severe restrictions on the design variables. In Paper II, we analyze the convergence properties of a linear elasticity problem in which void is approximated by a weak material. We show that the error introduced by the weak material approximation, after a finite element discretization, is bounded by terms that scale as <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> and <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />1/2hs, where h is the mesh size and s depends on the order of the finite element basis functions. In addition, we show that the condition number of the system matrix scales inversely proportional to <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />, and we also construct a left preconditioner that yields a system matrix with a condition number independent of <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />. In Paper III, we observe that the standard sound-hard material approximation with <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> gives rise to ill-conditioned system matrices at certain wavenumbers due to resonances within the approximated sound-hard material. To cure this defect, we propose a stabilization scheme that makes the condition number of the system matrix independent of the wavenumber. In addition, we demonstrate that the stabilized formulation performs well in the context of design optimization of an acoustic waveguide transmission device. In Paper IV, we analyze the convergence properties of a wave propagation problem in which sound-hard material is approximated by a dense fluid. To avoid the occurrence of internal resonances, we generalize the stabilization scheme presented in Paper III. We show that the error between the solution obtained using the stabilized soundhard material approximation and the solution to the problem with exactly modeled sound-hard material is bounded proportionally to <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />.

Ce travail est motivé par une demande croissante d'informations quantitatives sur l'évolution du littoral. Nous avons étudié deux approches pour l'analyse de la dynamique sédimentaire. Les deux techniques aboutissent à la résolution de modèles non-locaux pour le fond. L'étude mathématique a porté sur l'analyse de l'existence et l'unicité de perturbations autour des ondes progressives solutions du modèle de Fowler. Nous avons montré que les solutions constantes de l'équation de Fowler sont instables. Pour la simulation numérique de ce modèle, nous avons dans un premier temps considéré des schémas aux différences finies explicites pour lesquels nous avons obtenu des critères de stabilité numérique. Dans un second temps, nous avons utilisé une approche par splitting de sorte à pouvoir résoudre la convection, puis la diffusion et l'anti-diffusion fractionnaire de façon exacte. Ensuite, il est apparu que nous pouvions utiliser les principes de minimisation pour décrire l'évolution d'un lit érodable sous l'action de l'eau où le fond est considéré comme une structure déformable de faible rigidité s'adaptant en minimisant une certaine fonctionnelle d'énergie. Il est intéressant de constater que cette seconde approche peut être liée à la première car elle débouche aussi sur une équation de type Exner avec un terme non-local. En nous inspirant du modèle morphodynamique non-local de Fowler, nous concluons cette thèse par une application exotique au traitement de signal où nous proposons une nouvelle méthode de filtrage. / This work is motivated by a growing demand for quantitative information on the evolution of the coastline.We have studied two approaches for the analysis of sand morphodynamics.Both techniques lead to the resolution of nonlocal models for the seabottom.The mathematical study focused on the analysis of the existence and uniqueness of perturbations around the travelling-waves solutions of the Fowler model. We have shown that constant solutions of Fowler's equation are unstable.For the numerical simulation of this model, we have first considered explicit finite difference schemes for which we got numerical stability criteria. We have next used an approach by splitting method in order to solve first the convection, then the diffusion/fractional anti-diffusion exactly. We have also used minimization principles to describe the evolution of an erodible bed sheared by a fluid flow where the seabed is considered as a deformable structure with low stiffness whichadapts itself by minimizing a certain energy functional. It is interesting to note that this secondapproach can be linked to the first one because it also leads to a new Exner equation with a nonlocal term for the flux. Inspired by Fowler's morphodynamical model, we conclude this dissertation with an unexpected application to signal processing.

Visando atender requisitos cada vez mais rigorosos de projeto e exigência dos consumidores, é necessário extrair o máximo desempenho de uma dada estrutura de produto, buscando sempre propriedades superiores as atuais. Para obter-se tais propriedades dinâmicas superiores (resistência, rigidez, peso) existem vários métodos de otimização estrutural, entre os quais a otimização de parâmetros, utilizada em ajustes finos de projeto; a otimização topológica, mais complexa e condicionada pelos processos de fabricação disponíveis e a otimização de forma utilizada em chapas estruturais. Dentre os métodos recentes de otimização de forma merece destaque o método de Padrões de Bossas por Formas Modais desenvolvido por Fredö e Hedlung (2004), que permite grandes ganhos de rigidez estrutural com pequenas deformações no formato das chapas. Entretanto, tal método tem sua aplicabilidade restrita, pois depende de fatores de ponderação cujo critério de escolha não foi explorado pelos autores. O presente trabalho analisa teoricamente o método desenvolvido por Fredö e Hedlung (2004), utilizando para tal uma chapa metálica em condições controladas para determinar um critério coerente de definição dos parâmetros de ponderação do método via otimização computacional suportado por uma análise modal via método dos elementos finitos. Com os resultados dessa análise pode-se criar um programa para implementação do método de Padrões de Bossas por Formas Modais em aplicações industriais, com ganhos significativos nas características estruturais de produtos sem impactos no custo final. / To meet increasingly higher design requirements and consumer demands for design, it is necessary to extract the maximum performance of a given product structure, always seeking superior properties versus current design. To obtain such superior dynamic properties (strength, stiffness, weight) there are several methods of structural optimization, including the parameters optimization used to fine-tune design; the topology optimization, more complex and conditioned by manufacturing processes and the shape optimization as used in structural plates. Among the newer methods for shape optimization, it is worth mentioning the Panel embossing pattern optimization method developed by Frëdo and Hedlung (2004), which allows large structural rigidity gains with small deformations in the plate shape. However, this method has a limited applicability because it depends on weighting factors whose selection criterion was unexplored by the authors. This work theoretically analyzes the method developed by Frëdo and Hedlung (2004), using for such a metal sheet under controlled conditions to determine a coherent criterion for the weighting parameters definition process using computational optimization supported by a modal analysis via finite element method. With the results of this analysis, it was possible to create a program to implement the method embossing pattern optimization method in industrial applications, with significant gains in structural characteristics of products without affecting the final cost.

En biologie, lorsqu'une quantité importante de phospholipides est insérée dans un milieu aqueux, ceux-Ci s'assemblent alors par paires pour former une bicouche, plus communément appelée vésicule. En 1973, Helfrich a proposé un modèle simple pour décrire la forme prise par une vésicule. Imposant la surface de la bicouche et le volume de fluide qu'elle contient, leur forme minimise une énergie élastique faisant intervenir des quantités géométriques comme la courbure, ainsi qu'une courbure spontanée mesurant l'asymétrie entre les deux couches. Les globules rouges sont des exemples de vésicules sur lesquels sont fixés un réseau de protéines jouant le rôle de squelette au sein de la membrane. Un des principaux travaux de la thèse fut d'introduire et étudier une condition de boule uniforme, notamment pour modéliser l'effet du squelette. Dans un premier temps, on cherche à minimiser l'énergie de Helfrich sans contrainte puis sous contrainte d'aire. Le cas d'une courbure spontanée nulle est connu sous le nom d'énergie de Willmore. Comme la sphère est un minimiseur global de l'énergie de Willmore, c'est un bon candidat pour être un minimiseur de l'énergie de Helfrich parmi les surfaces d'aire fixée. Notre première contribution dans cette thèse a été d'étudier son optimalité. On montre qu'en dehors d'un certain intervalle de paramètres, la sphère n'est plus un minimum global, ni même un minimum local. Par contre, elle est toujours un point critique. Ensuite, dans le cas de membranes à courbure spontanée négative, on se demande si la minimisation de l'énergie de Helfrich sous contrainte d'aire peut être effectuée en minimisant individuellement chaque terme. Cela nous conduit à minimiser la courbure moyenne totale sous contrainte d'aire et à déterminer si la sphère est la solution de ce problème. On montre que c'est le cas dans la classe des surfaces axisymétriques axiconvexes mais que ce n'est pas vrai en général.Enfin, lorsqu'une contrainte d'aire et de volume sont considérées simultanément, le minimiseur ne peut pas être une sphère qui n'est alors plus admissible. En utilisant le point de vue de l'optimisation de formes, la troisième et plus importante contribution de cette thèse est d'introduire une classe plus raisonnable de surfaces, pour laquelle l'existence d'un minimiseur suffisamment régulier est assurée pour des fonctionnelles et des contraintes générales faisant intervenir les propriétés d'ordre un et deux des surfaces. En s'inspirant de ce que fit Chenais en 1975 quand elle a considéré la propriété de cône uniforme, on considère les surfaces satisfaisant une condition de boule uniforme. On étudie d'abord des fonctionnelles purement géométriques puis nous autorisons la dépendance à travers la solution de problèmes aux limites elliptiques d'ordre deux posés sur le domaine intérieur à la surface / In biology, when a large amount of phospholipids is inserted in aqueous media, they immediatly gather in pairs to form bilayers also called vesicles. In 1973, Helfrich suggested a simple model to characterize the shapes of vesicles. Imposing the area of the bilayer and the volume of fluid it contains, their shape is minimizing a free-Bending energy involving geometric quantities like curvature, and also a spontanuous curvature measuring the asymmetry between the two layers. Red blood cells are typical examples of vesicles on which is fixed a network of proteins playing the role of a skeleton inside the membrane. One of the main work of this thesis is to introduce and study a uniform ball condition, in particular to model the effects of the skeleton. First, we minimize the Helfrich energy without constraint then with an area constraint. The case of zero spontaneous curvature is known as the Willmore energy. Since the sphere is the global minimizer of the Willmore energy, it is a good candidate to be a minimizer of the Helfrich energy among surfaces of prescribed area. Our first main contribution in this thesis was to study its optimality. We show that apart from a specific interval of parameters, the sphere is no more a global minimizer, neither a local minimizer. However, it is always a critical point. Then, in the specific case of membranes with negative spontaneous curvature, one can wonder whether the minimization of the Helfrich energy with an area constraint can be done by minimizing individually each term. This leads us to minimize total mean curvature with prescribed area and to determine if the sphere is a solution to this problem. We show that it is the case in the class of axisymmetric axiconvex surfaces but that it does not hold true in the general case. Finally, considering both area and volume constraints, the minimizer cannot be the sphere, which is no more admissible. Using the shape optimization point of view, the third main and most important contribution of this thesis is to introduce a more reasonable class of surfaces, in which the existence of an enough regular minimizer is ensured for general functionals and constraints involving the first- and second-Order geometric properties of surfaces. Inspired by what Chenais did in 1975 when she considered the uniform cone property, we consider surfaces satisfying a uniform ball condition. We first study purely geometric functionals then we allow a dependence through the solution of some second-Order elliptic boundary value problems posed on the inner domain enclosed by the shape

Afin de faire face aux nouveaux challenges de l’industrie automobile, les ingénieurs souhaitent appliquer des méthodes d’optimisation à chaque étape du processus de conception. En élargissant l’espace de conception aux paramètres de forme, en augmentant leur nombre et en étendant les plages de variation, de nouveaux verrous sont apparus. C’est le cas de la résistance aux chocs. Avec les temps de calcul long, la non-linéarité, l’instabilité et la dispersion numérique de ce problème de dynamique rapide, la méthode usuellement employée, l’optimisation par plan d’expériences et surfaces de réponse, devient trop coûteuse pour être utilisée industriellement. Se pose alors la problématique suivante : Comment faire de l’optimisation de forme en dynamique rapide avec un nombre élevé de paramètres ?. Pour y répondre, les méthodes d’optimisation par gradient s’avèrent être les plus judicieuses. Le nombre de paramètres a une influence réduite sur le coût de l’optimisation. Elles permettent donc l’optimisation de problèmes ayant de nombreux paramètres. Cependant, les méthodes classiques de calcul du gradient sont peu pertinentes en dynamique rapide : le coût en nombre de simulations et le bruit empêchent l’utilisation des différences finies et le calcul du gradient en dérivant les équations de dynamique rapide n’est pas encore disponible et serait très intrusif vis-à-vis des logiciels. Au lieu de déterminer le gradient, au sens classique du terme, des problèmes de crash, nous avons cherché à l’estimer. L’Equivalent Static Loads Method est une méthode permettant l’optimisation à moindre coût basée sur la construction d’un problème statique linéaire équivalent au problème de dynamique rapide. En utilisant la dérivée du problème équivalent comme estimation du gradient, il nous a été possible d’optimiser des problèmes de dynamique rapide ayant des épaisseurs comme variables d’optimisation. De plus, si l’on construit les équations du problème équivalent avec la matrice de rigidité sécante, l’approximation du gradient n’en est que meilleure. De cette manière, il est aussi possible d’estimer le gradient par rapport à la position des nœuds du modèle de calcul. Comme il est plus courant de travailler avec des paramètres CAO, il faut déterminer la dérivée de la position des nœuds par rapport à ces paramètres. Nous pouvons le faire de manière analytique si nous utilisons une surface paramétrique pour définir la forme et ses points de contrôle comme variables d’optimisation. Grâce à l’estimation du gradient et à ce lien entre nœuds et paramètres de forme, l’optimisation de forme avec un nombre important de paramètres est désormais possible à moindre coût. La méthode a été développée pour deux familles de critères issues du crash automobile. La première est liée au déplacement d’un nœud, objectif important lorsqu’il faut préserver l’intégrité de l’habitacle du véhicule. La seconde est liée à l’énergie de déformation. Elle permet d’assurer un bon comportement de la structure lors du choc. / In order to face their new industrial challenges, automotive constructors wish to apply optimization methods in every step of the design process. By including shape parameters in the design space, increasing their number and their variation range, new problematics appeared. It is the case of crashworthiness. With the high computational time, the nonlinearity, the instability and the numerical dispersion of this rapid dynamics problem, metamodeling techniques become to heavy for the standardization of those optimization methods. We face this problematic: ”How can we carry out shape optimization in rapid dynamics with a high number of parameters ?”. Gradient methods are the most likely to solve this problematic. Because the number of parameters has a reduced effect on the optimization cost, they allow optimization with a high number of parameters. However, conventional methods used to calculate gradients are ineffective: the computation cost and the numerical noise prevent the use of finite differences and the calculation of a gradient by deriving the rapid dynamics equations is not currently available and would be really intrusive towards the software. Instead of determining the real gradient, we decided to estimate it. The Equivalent Static Loads Method is an optimization method based on the construction of a linear static problem equivalent to the rapid dynamic problem. By using the sensitivity of the equivalent problem as the estimated gradient, we have optimized rapid dynamic problems with thickness parameters. It is also possible to approximate the derivative with respect to the position of the nodes of the CAE model. But it is more common to use CAD parameters in shape optimization studies. So it is needed to have the sensitivity of the nodes position with these CAD parameters. It is possible to obtain it analytically by using parametric surface for the shape and its poles as parameters. With this link between nodes and CAD parameters, we can do shape optimization studies with a large number of parameters and this with a low optimization cost. The method has been developed for two kinds of crashworthiness objective functions. The first family of criterions is linked to a nodal displacement. This category contains objectives like the minimization of the intrusion inside the passenger compartment. The second one is linked to the absorbed energy. It is used to ensure a good behavior of the structure during the crash.

O projeto de estruturas e máquinas deve considerar as restrições impostas pelas condições de contorno. Tais condições podem ser de natureza dinâmica, limitando assim as faixas de frequência às quais a estrutura ou máquina pode operar. Dentre as diferentes ferramentas disponíveis para trabalhar com restrições dinâmicas, a otimização de forma se mostra como uma interessante alternativa para afastar as frequências naturais das faixas problemáticas. Um modelo de elementos finitos de malha de 4x5 elementos é correlacionado com os resultados de uma análise modal experimental, e a otimização é realizada utilizando-se o software Nastran. Após usinar a placa com a espessura otimizada, boa concordância é atingida entre os resultados experimentais e os previstos numericamente. Apesar dos bons resultados, obter a placa com 4x5 elementos, cada qual com sua espessura, foi difícil por conta das dimensões envolvidas. É mais apropriado fabricar uma superfície contínua na placa com uma geometria conhecida. Para isso, modelos com malhas mais finas são necessários no procedimento de otimização, e tal quantidade de variáveis nem sempre converge a uma solução. De fato, um modelo com 48x60 elementos para a placa estudada não convergiu a uma solução para as mesmas frequências desejadas. A contribuição principal deste trabalho é mostrar que uma interpolação (linear ou cúbica) dos resultados de uma otimização de forma de malha grossa leva à solução do problema de otimização de uma malha mais refinada. Em outras palavras, é possível obter uma geometria de superfície contínua que otimiza frequências naturais em placas a partir de modelos de elementos finitos de malha mais grossa, sendo assim desnecessários modelos de malhas muito refinadas e altos custos computacionais. / The design of structures and machines must consider the restrictions imposed by the boundary conditions. Such conditions can be of dynamic nature, thus limiting the frequency ranges that the structure/machine can operate. Among the different design tools available for dealing with dynamics restrictions, shape optimization is an interesting way of deviating natural frequencies from problematic ranges. In this work, one presents the shape optimization of a cantilever plate aiming at desired natural frequencies. A 4x5 elements grid finite element model is correlated to results from experimental modal analysis, and the optimization is done with help of Nastran software. After machining the plate with optimized thickness, good agreement is achieved between experimental and numerically predicted results. Despite successful results, machining the plate in 4x5 discrete elements with individual (stepped) thickness showed to be cumbersome. It is more appropriated to machine a smooth surface on the plate with known geometry. For that, finer grid element models are necessary in the optimization procedure, and such amount of design variables not always converge to a solution. In fact, a model with 48x60 elements for the plate in study did not converge to a solution for the same target frequencies. The main contribution of this work is showing that an interpolation (linear or cubic) of the coarser grid shape optimization results leads to the solution of a finer grid optimization problem. In other words, it is possible to obtain the smooth surface geometry that optimizes natural frequencies in plates from coarser grid finite element models, thus not requiring fine grid models and high computational costs.

No description available.