Flexible framework for fluid topology optimization with OpenFOAM (R) and finite element-based high-level discrete adjoint method (FEniCS/dolfin-adjoint)

Autor(es)

Alonso, Diego Hayashi ^[1] ; Garcia Rodriguez, Luis Fernando ^[1] ; Silva, Emilio Carlos Nelli ^[1]

Afiliação do(s) autor(es)

^[1] Univ Sao Paulo, Polytech Sch, Dept Mechatron & Mech Syst Engn, Sao Paulo, SP - Brazil

• Tipo de documento: Artigo Científico
• Fonte: STRUCTURAL AND MULTIDISCIPLINARY OPTIMIZATION; v. 64, n. 6 SEP 2021.
• Citações Web of Science: 0

Resumo

In order to implement the topology optimization method, it is necessary to simulate the fluid flow dynamics and also obtain the sensitivities with respect to the design variable (such as through the adjoint method). However, more complex fluid flows, such as turbulent, non-Newtonian, and compressible flows, may turn the implementation of these two aspects difficult and non-intuitive. In order to solve this deadlock, this work proposes the combination of two well-known and established open-source softwares: OpenFOAM (R) and FEniCS/dolfin-adjoint. OpenFOAM (R) already provides efficient implementations for various fluid flow models, while FEniCS, when combined with the dolfin-adjoint library, provides an efficient and automatic high-level discrete adjoint model. There have been various attempts for obtaining the adjoint model directly in OpenFOAM (R) , but they mostly rely on the following: (1) manually deducing the adjoint equations, which may become a hard and cumbersome task for complex models; (2) C++ automatic differentiation tools, which are generally computationally inefficient; and (3) finite differences, which have been developed for shape optimization (not topology optimization, where there are many more design variable values). Nonetheless, these approaches generally do not provide an easy setup, and may be fairly complex to consider. The FEniCS platform does not provide any fluid flow model out of the box, but makes it fairly simple to ``simplistically{''} define them. The main problem of the FEniCS implementation and even implementations ``by hand{''} (such as in C++, Matlab (R) or Python) is the convergence of the simulation, which would possibly require fairly complex adjustments in the implementation in order to reach convergence. Therefore, the combination proposed in this work (OpenFOAM (R) and FEniCS/dolfin-adjoint) is a simpler but efficient approach to consider more complex fluid flows, countering the difficult adjoint model implementation in OpenFOAM (R) and also the convergence issues in FEniCS. The implemented framework, referred as ``FEniCS TopOpt Foam{''}, can perform the coupling between the two softwares. Numerical examples are presented considering laminar and turbulent flows (Spalart-Allmaras model) for 2D, 2D axisymmetric, and 3D domains. (AU)

• Processo FAPESP: 13/24434-0 - Sistemas propulsores eletromagnéticos para coração artificial implantável e dispositivos de suporte circulatório mecânico
• Beneficiário: José Roberto Cardoso
• Modalidade de apoio: Auxílio à Pesquisa - Temático
• Processo FAPESP: 14/50279-4 - Brasil Research Centre for Gas Innovation
• Beneficiário: Julio Romano Meneghini
• Modalidade de apoio: Auxílio à Pesquisa - Programa Centros de Pesquisa em Engenharia
• Processo FAPESP: 17/27049-0 - Método de otimização topológica aplicado ao projeto de rotor e voluta de dispositivos de assistência ventricular baseados no efeito de viscosidade (princípio Tesla)
• Beneficiário: Diego Hayashi Alonso
• Modalidade de apoio: Bolsas no Brasil - Doutorado Direto