Hydrogen fuel cells are devices used to generate electricity from the electrochemical reaction between air and hydrogen gas. An attractive advantage of these devices is that their byproduct is water, which is very safe to the environment. However, hydrogen fuel cells still lack some improvements in terms of increasing their life time and electricity production, decreasing power losses, and optimizing their operating conditions. In this thesis, the cathode part of the hydrogen fuel cell will be considered. This part mainly consists of an air gas channel and a gas diffusion layer. To simulate the fluid dynamics taking place in the cathode, we present two models, a general model and a simple model both based on a set of conservation laws governing the fluid dynamics and chemical reactions. A numerical method to solve these models is presented and verified in terms of accuracy. We also show that both models give similar results and validate the simple model by recovering a polarization curve obtained experimentally. Next, a shape optimization problem is introduced to find an optimal design of the air gas channel. This problem is defined from the simple model and a cost functional, $E$, that measures efficiency factors. The objective of this functional is to maximize the electricity production, uniformize the reaction rate in the catalytic layer and minimize the pressure drop in the gas channel. The impact of the gas channel shape optimization is investigated with a series of test cases in long and short fuel cell geometries. In most instances, the optimal design improves efficiency in on- and off-design operating conditions by shifting the polarization curve vertically and to the right. The second primary goal of the thesis is to analyze mathematical issues related to the introduced shape optimization problem. This involves existence and uniqueness of the solution for the presented model and differentiability of the state variables with respect to the domain of the air channel. The optimization problem is solved using the gradient method, and hence the gradient of $E$ must be found. The gradient of $E$ is obtained by introducing an adjoint system of equations, which is coupled with the state problem, namely the simple model of the fuel cell. The existence and uniqueness of the solution for the adjoint system is shown, and the shape differentiability of the cost functional $E$ is proved.

Hydrogen fuel cells are devices used to generate electricity from the electrochemical reaction between air and hydrogen gas. An attractive advantage of these devices is that their byproduct is water, which is very safe to the environment. However, hydrogen fuel cells still lack some improvements in terms of increasing their life time and electricity production, decreasing power losses, and optimizing their operating conditions. In this thesis, the cathode part of the hydrogen fuel cell will be considered. This part mainly consists of an air gas channel and a gas diffusion layer. To simulate the fluid dynamics taking place in the cathode, we present two models, a general model and a simple model both based on a set of conservation laws governing the fluid dynamics and chemical reactions. A numerical method to solve these models is presented and verified in terms of accuracy. We also show that both models give similar results and validate the simple model by recovering a polarization curve obtained experimentally. Next, a shape optimization problem is introduced to find an optimal design of the air gas channel. This problem is defined from the simple model and a cost functional, $E$, that measures efficiency factors. The objective of this functional is to maximize the electricity production, uniformize the reaction rate in the catalytic layer and minimize the pressure drop in the gas channel. The impact of the gas channel shape optimization is investigated with a series of test cases in long and short fuel cell geometries. In most instances, the optimal design improves efficiency in on- and off-design operating conditions by shifting the polarization curve vertically and to the right. The second primary goal of the thesis is to analyze mathematical issues related to the introduced shape optimization problem. This involves existence and uniqueness of the solution for the presented model and differentiability of the state variables with respect to the domain of the air channel. The optimization problem is solved using the gradient method, and hence the gradient of $E$ must be found. The gradient of $E$ is obtained by introducing an adjoint system of equations, which is coupled with the state problem, namely the simple model of the fuel cell. The existence and uniqueness of the solution for the adjoint system is shown, and the shape differentiability of the cost functional $E$ is proved.

Design optimization of truss structures is a popular topic in aerospace, mechanical, civil, and structural engineering due to benefits to industry. Common design problem for the structures is the weight minimization. Especially in aerospace engineering the minimization of the weight of the total structure gets the highest importance in the design. This study focuses on the design optimization of 2D and 3D truss structures. The objective function is the total mass of the structure which is subjected to stress and nodal displacement constraints. To optimize the design, Genetic Algorithm (GA) is preferred due to its efficiency in dealing with problems with discrete design variables as in the case of truss structures. This technique yields more realistic results than linear programming methods. In the thesis, a finite element code is developed for the analysis of planar and space truss structures. The developed finite element solver is coupled with a genetic algorithm optimization code which is also developed as a part of the thesis study. Different truss optimization case studies are performed to demonstrate the performance of the finite element solver and the genetic algorithm optimization code that are developed. It is shown that with the use of adaptive penalty function employing scaled fitnesses, the arbitrariness issue of the factor multiplying the error term in the augmented fitness function can be resolved. It is also shown that significant weight reduction can v be achieved by employing shape optimization together with size optimization compared to pure size optimization.

In this thesis, structural optimization study is performed by using three different methods. In the first method, optimization is performed using MSC.NASTRAN Optimization Module, a commercial structural analysis program. In the second method, optimization is performed using the optimization code prepared in MATLAB and MSC.NASTRAN as the solver. As the third method, optimization is performed by using the optimization code prepared in MATLAB and analytical equations as the solver. All three methods provide certain advantages in the solution of optimization problems. Therefore, within the context of the thesis these methods are demonstrated and the interface codes specific to the programs used in this thesis are explained in detail. In order to compare the results obtained by the methods, the verification study has been performed on a cantilever beam with rectangular cross-section. In the verification study, the height and width of the cross-section of the beam are taken as the two design parameters. This way it has been possible to show the design space on the two dimensional graph, and it becomes easier to trace the progress of the optimization methods during each step. In the last section structural optimization of a multi-element wing torque box has been performed by the MSC.NASTRAN optimization module. In this section geometric property optimization has been performed for constant tip loading and variable loading along the wing span. In addition, within the context of shape optimization optimum rib placement problem has also been solved.

No description available.

No description available.

No description available.

No description available.

Cette thèse porte sur l'optimisation de forme via la méthode des "level sets" pour deux comportements mécaniques induisant des déplacements non différentiables par rapport à la forme: le contact et la plasticité. Pour y remédier, nous utilisons des problèmes approchés issus de méthode de pénalisation et de régularisation.Dans la première partie, nous présentons quelques notions fondamentales d'optimisation de forme (chapitre 1). Puis nous exposons les résultats qui seront utiles à l'analyse des deux problèmes mécaniques considérés et nous illustrons ces résultats.La deuxième partie introduit les modèles statiques de contact (chapitre 3) et le modèle statique de plasticité (chapitre 4) que nous utilisons dans le manuscrit. Pour chacun, nous donnons les bases de la modélisation mécanique, une analyse mathématique des inéquations variationnelles associées et nous expliquons quels solveurs nous avons implémentés.La dernière partie se focalise sur l'optimisation de forme. Dans chacun des chapitres nous donnons les versions pénalisées et régularisées des modèles, prouvons, pour certains, leur convergence vers les modèles exactes, calculons leurs gradients de forme et proposons des exemples 2D et, en contact, 3D. Ainsi, dans le chapitre 5, traitons-nous du contact et considérons deux sortes de problèmes: le premier dans lequel la zone de contact est fixe, le second dans lequel la zone de contact est optimisable. Pour ce dernier, nous introduisons deux méthodes pour résoudre du contact sans discrétiser la zone de contact. Dans le chapitre 6, nous abordons le modèle de Hencky que nous approximons grâce à une pénalisation de Perzyna ainsi que grâce à un modèle de notre crue. / The main purpose of this thesis is to perform shape optimisation, in the framework of the level set method, for two mechanical behaviours inducing displacement which are not shape differentiable: contact and plasticity. To overcome this obstacle, we use approximate problems found by penalisation and regularisation.In the first part, we present some classical notions in optimal design (chapter 1). Then we give the mathematical results needed for the analysis of the two mechanical problems in consideration and illustrate these results.The second part is meant to introduce the five static contact models (chapter 3) and the static plasticity model (chapter 4) we use in the manuscript. For each chapter we provide the basis of the mechanical modeling, a mathematical analysis of the related variational inequations and, finally, explain how we implement the associated solvers.Eventually the last part, consisting of two chapters is devoted to shape optimisation. In each of them, we state the regularised versions of the models, prove, for some of them, the convergence to the exact ones, compute shape gradients and perform some numerical experiments in 2D and, for contact, in 3D. Thus, in chapter 5, we focus on contact and consider two types of optimal design problems: one with a fixed contact zone and another one with a mobile contact zone. For this last type, we introduce two ways to solve frictionless contact without meshing the contact zone. One of them is new and the other one has never been employed in this framework. In chapter 6, we deal with the Hencky model which we approximate thanks to a Perzyna penalised problem as well as a home-made one.

As ever more challenging designs are required in present-day industries, the traditional trial-and-error procedure frequently used for designing mechanical parts slows down the design process and yields suboptimal designs, so that new approaches are needed to obtain a competitive advantage. With the ascent of the Finite Element Method (FEM) in the engineering community in the 1970s, structural shape optimization arose as a promising area of application. However, due to the iterative nature of shape optimization processes, the handling of large quantities of numerical models along with the approximated character of numerical methods may even dissuade the use of these techniques (or fail to exploit their full potential) because the development time of new products is becoming ever shorter. This Thesis is concerned with the formulation of a 3D methodology based on the Cartesian-grid Finite Element Method (cgFEM) as a tool for efficient and robust numerical analysis. This methodology belongs to the category of embedded (or fictitious) domain discretization techniques in which the key concept is to extend the structural analysis problem to an easy-to-mesh approximation domain that encloses the physical domain boundary. The use of Cartesian grids provides a natural platform for structural shape optimization because the numerical domain is separated from a physical model, which can easily be changed during the optimization procedure without altering the background discretization. Another advantage is the fact that mesh generation becomes a trivial task since the discretization of the numerical domain and its manipulation, in combination with an efficient hierarchical data structure, can be exploited to save computational effort. However, these advantages are challenged by several numerical issues. Basically, the computational effort has moved from the use of expensive meshing algorithms towards the use of, for example, elaborate numerical integration schemes designed to capture the mismatch between the geometrical domain boundary and the embedding finite element mesh. To do this we used a stabilized formulation to impose boundary conditions and developed novel techniques to be able to capture the exact boundary representation of the models. To complete the implementation of a structural shape optimization method an adjunct formulation is used for the differentiation of the design sensitivities required for gradient-based algorithms. The derivatives are not only the variables required for the process, but also compose a powerful tool for projecting information between different designs, or even projecting the information to create h-adapted meshes without going through a full h-adaptive refinement process. The proposed improvements are reflected in the numerical examples included in this Thesis. These analyses clearly show the improved behavior of the cgFEM technology as regards numerical accuracy and computational efficiency, and consequently the suitability of the cgFEM approach for shape optimization or contact problems. / La competitividad en la industria actual impone la necesidad de generar nuevos y mejores diseños. El tradicional procedimiento de prueba y error, usado a menudo para el diseño de componentes mecánicos, ralentiza el proceso de diseño y produce diseños subóptimos, por lo que se necesitan nuevos enfoques para obtener una ventaja competitiva. Con el desarrollo del Método de los Elementos Finitos (MEF) en el campo de la ingeniería en la década de 1970, la optimización de forma estructural surgió como un área de aplicación prometedora. El entorno industrial cada vez más exigente implica ciclos cada vez más cortos de desarrollo de nuevos productos. Por tanto, la naturaleza iterativa de los procesos de optimización de forma, que supone el análisis de gran cantidad de geometrías (para las se han de usar modelos numéricos de gran tamaño a fin de limitar el efecto de los errores intrínsecamente asociados a las técnicas numéricas), puede incluso disuadir del uso de estas técnicas. Esta Tesis se centra en la formulación de una metodología 3D basada en el Cartesian-grid Finite Element Method (cgFEM) como herramienta para un análisis numérico eficiente y robusto. Esta metodología pertenece a la categoría de técnicas de discretización Immersed Boundary donde el concepto clave es extender el problema de análisis estructural a un dominio de aproximación, que contiene la frontera del dominio físico, cuya discretización (mallado) resulte sencilla. El uso de mallados cartesianos proporciona una plataforma natural para la optimización de forma estructural porque el dominio numérico está separado del modelo físico, que podrá cambiar libremente durante el procedimiento de optimización sin alterar la discretización subyacente. Otro argumento positivo reside en el hecho de que la generación de malla se convierte en una tarea trivial. La discretización del dominio numérico y su manipulación, en coalición con la eficiencia de una estructura jerárquica de datos, pueden ser explotados para ahorrar coste computacional. Sin embargo, estas ventajas pueden ser cuestionadas por varios problemas numéricos. Básicamente, el esfuerzo computacional se ha desplazado. Del uso de costosos algoritmos de mallado nos movemos hacia el uso de, por ejemplo, esquemas de integración numérica elaborados para poder capturar la discrepancia entre la frontera del dominio geométrico y la malla de elementos finitos que lo embebe. Para ello, utilizamos, por un lado, una formulación de estabilización para imponer condiciones de contorno y, por otro lado, hemos desarrollado nuevas técnicas para poder captar la representación exacta de los modelos geométricos. Para completar la implementación de un método de optimización de forma estructural se usa una formulación adjunta para derivar las sensibilidades de diseño requeridas por los algoritmos basados en gradiente. Las derivadas no son sólo variables requeridas para el proceso, sino una poderosa herramienta para poder proyectar información entre diferentes diseños o, incluso, proyectar la información para crear mallas h-adaptadas sin pasar por un proceso completo de refinamiento h-adaptativo. Las mejoras propuestas se reflejan en los ejemplos numéricos presentados en esta Tesis. Estos análisis muestran claramente el comportamiento superior de la tecnología cgFEM en cuanto a precisión numérica y eficiencia computacional. En consecuencia, el enfoque cgFEM se postula como una herramienta adecuada para la optimización de forma. / Actualment, amb la competència existent en la industria, s'imposa la necessitat de generar nous i millors dissenys . El tradicional procediment de prova i error, que amb freqüència es fa servir pel disseny de components mecànics, endarrereix el procés de disseny i produeix dissenys subòptims, pel que es necessiten nous enfocaments per obtindre avantatge competitiu. Amb el desenvolupament del Mètode dels Elements Finits (MEF) en el camp de l'enginyeria en la dècada de 1970, l'optimització de forma estructural va sorgir com un àrea d'aplicació prometedora. No obstant això, a causa de la natura iterativa dels processos d'optimització de forma, la manipulació dels models numèrics en grans quantitats, junt amb l'error de discretització dels mètodes numèrics, pot fins i tot dissuadir de l'ús d'aquestes tècniques (o d'explotar tot el seu potencial), perquè al mateix temps els cicles de desenvolupament de nous productes s'estan acurtant. Esta Tesi se centra en la formulació d'una metodologia 3D basada en el Cartesian-grid Finite Element Method (cgFEM) com a ferramenta per una anàlisi numèrica eficient i sòlida. Esta metodologia pertany a la categoria de tècniques de discretització Immersed Boundary on el concepte clau és expandir el problema d'anàlisi estructural a un domini d'aproximació fàcil de mallar que conté la frontera del domini físic. L'utilització de mallats cartesians proporciona una plataforma natural per l'optimització de forma estructural perquè el domini numèric està separat del model físic, que podria canviar lliurement durant el procediment d'optimització sense alterar la discretització subjacent. A més, un altre argument positiu el trobem en què la generació de malla es converteix en una tasca trivial, ja que la discretització del domini numèric i la seua manipulació, en coalició amb l'eficiència d'una estructura jeràrquica de dades, poden ser explotats per estalviar cost computacional. Tot i això, estos avantatges poden ser qüestionats per diversos problemes numèrics. Bàsicament, l'esforç computacional s'ha desplaçat. De l'ús de costosos algoritmes de mallat ens movem cap a l'ús de, per exemple, esquemes d'integració numèrica elaborats per poder capturar la discrepància entre la frontera del domini geomètric i la malla d'elements finits que ho embeu. Per això, fem ús, d'una banda, d'una formulació d'estabilització per imposar condicions de contorn i, d'un altra, desevolupem noves tècniques per poder captar la representació exacta dels models geomètrics Per completar la implementació d'un mètode d'optimització de forma estructural es fa ús d'una formulació adjunta per derivar les sensibilitats de disseny requerides pels algoritmes basats en gradient. Les derivades no són únicament variables requerides pel procés, sinó una poderosa ferramenta per poder projectar informació entre diferents dissenys o, fins i tot, projectar la informació per crear malles h-adaptades sense passar per un procés complet de refinament h-adaptatiu. Les millores proposades s'evidencien en els exemples numèrics presentats en esta Tesi. Estes anàlisis mostren clarament el comportament superior de la tecnologia cgFEM en tant a precisió numèrica i eficiència computacional. Així, l'enfocament cgFEM es postula com una ferramenta adient per l'optimització de forma. / Marco Alacid, O. (2017). Structural Shape Optimization Based On The Use Of Cartesian Grids [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/86195 / TESIS