[ca] Aquesta tesi cobreix el desenvolupament de mètodes numèrics a gran escala per a la simulació d'equacions amb derivades parcials en geometries arbitràriament complexes. L'aplicació d'aquesta tesi està orientada a la simulació estructural d'edificis i infraestructures civils, en les quals les demandes estètiques i de lleugeresa augmenten la complexitat de les geometries. En aquestes aplicacions, l'experimentació empírica no és factible durant el procés de disseny. Per tant, es confia en simulacions numèriques per a predir el comportament d'aquestes estructures sotmeses a càrregues realistes, e.g., efectes del vent que es compliquen amb la sofisticació de la geometria. Les eines de simulació actuals es basen en malles no estructurades adaptades a la geometria. La creació de malles no estructurades és un procés demandant que requereix la intervenció humana. A més, els mètodes de malles adaptades a la geometria no poden explotar els recursos computacionals d'alt rendiment moderns de manera eficient. L'objectiu principal d'aquesta tesi és el disseny de noves eines per a simular ràpidament, acurada i automàtica equacions amb derivades parcials sobre geometries definides per disseny assistit per ordinador. Per aquest motiu, es cerca un sistema que combina (1) un procés automàtic des del disseny assistit per ordinador fins a l'anàlisi d'elements finits, (2) una nova formulació espai-temps per a geometries en moviment, i (3) una implementació escalable per explotar recursos computacionals d'alt rendiment. Els desenvolupaments estan accessibles mitjançant programari de codi obert dins de l'ecosistema Gridap i del paquet FEMPAR (escrits en Julia i Fortran, respectivament). Les contribucions d'aquesta tesi incrementen la funcionalitat dels mètodes d'elements finits embeguts ("immersed", "embeddedd" o "unfitted" en anglès) respecte a l'estat de l'art. Aquests mètodes utilitzen malles de fons estructurades per resoldre equacions amb derivades parcials en dominis complexos. En la literatura, aquests dominis es representen implícitament per conjunts de nivell. Per a abordar aquesta limitació, (1) s'ha desenvolupat un algorisme robust que resol problemes en dominis descrits per representacions de contorn lineals. Aquest algorisme es basa en algorismes robustos per retallar poliedres. S'ha testejat l'algorisme amb totes les geometries (fitxers STL) analitzables del recurs Thingi10k (gairebé 5.000). S'ha ampliat l'algorisme per a representacions de contorn d'alt ordre. En aquesta extensió, es fan servir les bases de Bernstein-Bézier i algorismes de cerca d'arrels multivariables. S'ha validat el mètode mitjançant referències analítiques i geometries reals creades amb eines de disseny assistit per ordinador. A continuació, (2) s'ha formulat un marc d'elements finits embegut i espai-temps per a geometries explícites que es mouen. En aquesta formulació, s'utilitzen malles espacials per a evitar la necessitat d'algorismes geomètrics en 4D. En canvi, s'ha desenvolupat un mètode de transferència de valors inicials a les llesques temporals. Els resultats coincideixen amb l'anàlisi matemàtica i amb els experiments numèrics externs. Així mateix, s'ha demostrat l'aplicabilitat a problemes de fluids amb geometries complexes (2D i 3D) en rotació. Finalment, (3) s'ha proposat l'acceleració dels mètodes d'aquesta tesi a través d'algorismes altament escalables. Aquests algorismes aborden colls d'ampolla de la paral·lelització dels algorismes d'intersecció. S'ha demostrat l'escalabilitat d'aquests algorismes amb més de mil milions de cel·les i 12.000 processadors. Per afegiment, es poden combinar aquests algorismes amb tècniques d'adaptació de malla per reduir, encara més, el cost computacional. Aquestes eines proporcionen els mitjans per a accelerar el procés de disseny a simulació tot incrementant la fidelitat dels resultats.

[en] This thesis covers the development of large-scale numerical methods for the simulations of partial differential equations on arbitrarily complex geometries. The target application of this thesis is the structural simulation of buildings and civil infrastructures, in which lightweight and aesthetical demands usually increase the complexity of their geometries. In these applications, empirical experimentation is often not feasible during the design loop. Thus, we rely on numerical simulations to predict the performance of these shapes under realistic loads, e.g., wind loads that increase the complexity with the sophistication of the geometry. Current simulation tools are based on unstructured body-fitted meshes. The generation of unstructured meshes is time-consuming and involves human intervention. Furthermore, body-fitted methods cannot efficiently exploit modern high performance computing resources. The main goal of this thesis is to design novel simulation tools for rapid, accurate, and automated solutions of partial differential equations on geometries described by computer-aided design. Thus, we aim for a framework that combines (1) an automated pipeline from computer-aided design to finite element analysis, (2) a novel space-time formulation for moving geometries, and (3) a scalable implementation for high performance computing resources. Our developments are accessible through open-source software within the Gridap ecosystem and FEMPAR packages (written in Julia and Fortran, respectively). The contributions of this thesis increase the functionality of state-of-the-art unfitted (or immersed or embedded) finite element methods. These methods utilize structured background meshes to solve partial differential equations on complex domains. In the literature, these domains are implicitly represented by level sets. To address this limitation, (1) we developed a robust algorithm that solves problems on domains described by linear boundary representations. This algorithm is based on robust polyhedra clipping algorithms. We have tested the algorithm against all the analysis-ready geometries (STL files) in the Thingi10k dataset (almost 5,000). We have extended this algorithm to high-order boundary representations. In this extension, we utilize Bernstein-Bézier basis and multi-variate root-finding algorithms. We validate the resulting method with analytical benchmarks and real-world geometries from computer-aided design files. Then, (2) we formulated an unfitted space-time finite element framework for moving explicit geometries. In this formulation, we utilized space-only meshes, circumventing the need for 4D geometrical algorithms. In turn, we developed a transfer method for evaluating the initial values at each time slab. The results matched with analytical analysis and external numerical experiments. Furthermore, we have demonstrated the applicability of fluid problems on rotating complex (2D and 3D) geometries. Finally, (3) we proposed the acceleration of the methods of this thesis through highly scalable algorithms. These algorithms tackle the bottlenecks of parallelization of the intersection algorithms. We have demonstrated the scalability of these algorithms over one billion cells and 12,000 cores. Furthermore, we can combine these algorithms with adaptive mesh refinement techniques to reduce the computational cost further. These tools provide the means to significantly accelerate the design-to-simulation pipeline while increasing the fidelity of the results.