Numerical modelling and simulation are powerful tools for analysis and design, and with the improvement of computational power and numerical methods they are being applied on complex phenomena and systems. This work shows examples of the application of a very sophisticated numerical method, namely the Spectral/hp element method, in the study of the flow inside bioinspired complex geometries. The two topics investigated are fluid dynamic instabilities in a constricted channel and flow inside a paediatric ventricular assist device were studied by means of computational fluid mechanics. The constricted channel is an idealized model of a nasal cavity, which is characterized by complex airway channels, and also bears some resemblance to a human artery in the presence of an atherosclerotic plaques. The paediatric ventricular assist device is an actual device, designed by the Bioengineering research group of the Heart Institute of the Medicine School of the University of São Paulo, which works as a pump that assists the left ventricle of patients waiting for transplantation. Therefore, the aim of this thesis is to contribute in the understanding of biological and bio-inspired geometries flows, using computational tools. Linear and nonlinear stability were carried out for the constricted channel. Three different flow regimes were investigated: symmetric steady flow, which is stable for low Reynolds number, asymmetric steady flow, which rises as a result of the primary bifurcation of the symmetric flow and pulsatile flow. Direct stability analysis was carried out to determine the unstable regions and the critical values for each flow regime. The physical mechanisms behind the transition processes were studied by means of direct numerical simulations to characterize the bifurcations. Since the bifurcations had subcritical behaviour, the relevance of non-normal growth in these flows was assessed. Dependence on phase, Reynolds number and spanwise wavenumber of optimal modes were extensively investigated in stable regions of the three flow regimes. Convective instabilities were also studied in order to comprehend the physical mechanisms which led the optimal modes to their maxima growth, and different convective mechanisms were found. The flow inside the paediatric ventricular assist device was analyzed by means of threedimensional numerical simulations. A computational model based on special boundaries conditions was developed to model the pulsatile flow. In this model, the opening and closure of the mitral valve and diaphragm were represented with the use of specially devised boundary conditions. The driving force and the flow direction of the diaphragm were defined by velocity distribution on the diaphragm wall, and the opening and closure of the mitral valve were performed by a velocity waveform which goes to zero in the systolic period. Flow patterns, velocity fields and time-average wall shear rate were analyzed to evaluate the performance of the device. / Modelagem e simulação numéricas são ferramentas poderosas para análise e design, e com a melhoria do poder computacional e dos métodos numéricos, eles estão sendo aplicados em fenômenos e sistemas complexos. Este trabalho mostra exemplos de aplicações de um método numérico sofisticado, o método dos elementos espectrais/hp, no estudo do escoamento dentro de geometrias complexas bio-inspiradas. Os dois tópicos investigados são: instabilidades em dinâmica de fluido em um canal constrito e o escoamento dentro de um dispositivo de assistência ventricular pediátrica. O canal constrito é um modelo idealizado de uma cavidade nasal, que é caracterizada por canais complexos da via aérea, e também tem semelhança com uma artéria humana na presença de placas ateroscleróticas. O dispositivo de assistência ventricular pediátrica é um dispositivo real, projetado pelo grupo de pesquisa de Bioengenharia do Instituto do Coração da Faculdade de Medicina da Universidade de São Paulo, que funciona como uma bomba que auxilia o ventrículo esquerdo dos pacientes à espera de transplante. Portanto, o objetivo desta tese é contribuir na compreensão de escoamentos em geometrias biológicas e bio-inspiradas, usando ferramentas computacionais. Análises de estabilidade linear e não linear foram feitas para um canal constrito. Três diferentes regimes de escoamento foram empregados: escoamento estacionário simétrico, que é estável para baixo número de Reynolds, escoamento assimétrico, o qual é resultado da primeira bifurcação do escoamento simétrico e escoamento pulsátil. Análise de estabilidade direta foi executada para determinar as regiões instáveis em cada regime de escoamento. Os mecanismos físicos por trás do processo de transição foram estudados por meio de simulação numérica direta para caracterizar as bifurcações. Uma vez que, as bifurcações tiveram um comportamento subcrítico, a relevância do crescimento não normal nestes escoamentos foi avaliado. Assim, dependência com a fase, número de Reynolds e número de onda do modo tridimensional foram extensivamente investigados em regiões estáveis para os três regimes de escoamento. Instabilidades convectivas foram também estudadas a fim de compreender os mecanismos físicos que conduzem os modos ótimos para seus crescimentos máximos, e diferentes mecanismos convectivos foram encontrados. O escoamento dentro do dispositivo de assistência ventricular pediátrico foi analisado por meios de simulações numéricas tridimensionais. Um modelo computacional baseado em condições de contorno especiais foi desenvolvido para modelar o escoamento pulsátil. Neste modelo, a abertura e fechamento da válvula mitral e diafragma foram representados com o uso de condições de contorno especialmente elaboradas. A força motora e o direcionamento do fluxo do diafragma foram definidos por uma distribuição de velocidades na parede do diafragma, e a abertura e fechamento da válvula mitral foram executadas por uma função de onda de velocidade que vai a zero no período sistólico. Padrões do escoamento, campos de velocidade e tensão de cisalhamento no tempo foram analisadas para avaliar o desempenho do dispositivo.

Numerical modelling and simulation are powerful tools for analysis and design, and with the improvement of computational power and numerical methods they are being applied on complex phenomena and systems. This work shows examples of the application of a very sophisticated numerical method, namely the Spectral/hp element method, in the study of the flow inside bioinspired complex geometries. The two topics investigated are fluid dynamic instabilities in a constricted channel and flow inside a paediatric ventricular assist device were studied by means of computational fluid mechanics. The constricted channel is an idealized model of a nasal cavity, which is characterized by complex airway channels, and also bears some resemblance to a human artery in the presence of an atherosclerotic plaques. The paediatric ventricular assist device is an actual device, designed by the Bioengineering research group of the Heart Institute of the Medicine School of the University of São Paulo, which works as a pump that assists the left ventricle of patients waiting for transplantation. Therefore, the aim of this thesis is to contribute in the understanding of biological and bio-inspired geometries flows, using computational tools. Linear and nonlinear stability were carried out for the constricted channel. Three different flow regimes were investigated: symmetric steady flow, which is stable for low Reynolds number, asymmetric steady flow, which rises as a result of the primary bifurcation of the symmetric flow and pulsatile flow. Direct stability analysis was carried out to determine the unstable regions and the critical values for each flow regime. The physical mechanisms behind the transition processes were studied by means of direct numerical simulations to characterize the bifurcations. Since the bifurcations had subcritical behaviour, the relevance of non-normal growth in these flows was assessed. Dependence on phase, Reynolds number and spanwise wavenumber of optimal modes were extensively investigated in stable regions of the three flow regimes. Convective instabilities were also studied in order to comprehend the physical mechanisms which led the optimal modes to their maxima growth, and different convective mechanisms were found. The flow inside the paediatric ventricular assist device was analyzed by means of threedimensional numerical simulations. A computational model based on special boundaries conditions was developed to model the pulsatile flow. In this model, the opening and closure of the mitral valve and diaphragm were represented with the use of specially devised boundary conditions. The driving force and the flow direction of the diaphragm were defined by velocity distribution on the diaphragm wall, and the opening and closure of the mitral valve were performed by a velocity waveform which goes to zero in the systolic period. Flow patterns, velocity fields and time-average wall shear rate were analyzed to evaluate the performance of the device. / Modelagem e simulação numéricas são ferramentas poderosas para análise e design, e com a melhoria do poder computacional e dos métodos numéricos, eles estão sendo aplicados em fenômenos e sistemas complexos. Este trabalho mostra exemplos de aplicações de um método numérico sofisticado, o método dos elementos espectrais/hp, no estudo do escoamento dentro de geometrias complexas bio-inspiradas. Os dois tópicos investigados são: instabilidades em dinâmica de fluido em um canal constrito e o escoamento dentro de um dispositivo de assistência ventricular pediátrica. O canal constrito é um modelo idealizado de uma cavidade nasal, que é caracterizada por canais complexos da via aérea, e também tem semelhança com uma artéria humana na presença de placas ateroscleróticas. O dispositivo de assistência ventricular pediátrica é um dispositivo real, projetado pelo grupo de pesquisa de Bioengenharia do Instituto do Coração da Faculdade de Medicina da Universidade de São Paulo, que funciona como uma bomba que auxilia o ventrículo esquerdo dos pacientes à espera de transplante. Portanto, o objetivo desta tese é contribuir na compreensão de escoamentos em geometrias biológicas e bio-inspiradas, usando ferramentas computacionais. Análises de estabilidade linear e não linear foram feitas para um canal constrito. Três diferentes regimes de escoamento foram empregados: escoamento estacionário simétrico, que é estável para baixo número de Reynolds, escoamento assimétrico, o qual é resultado da primeira bifurcação do escoamento simétrico e escoamento pulsátil. Análise de estabilidade direta foi executada para determinar as regiões instáveis em cada regime de escoamento. Os mecanismos físicos por trás do processo de transição foram estudados por meio de simulação numérica direta para caracterizar as bifurcações. Uma vez que, as bifurcações tiveram um comportamento subcrítico, a relevância do crescimento não normal nestes escoamentos foi avaliado. Assim, dependência com a fase, número de Reynolds e número de onda do modo tridimensional foram extensivamente investigados em regiões estáveis para os três regimes de escoamento. Instabilidades convectivas foram também estudadas a fim de compreender os mecanismos físicos que conduzem os modos ótimos para seus crescimentos máximos, e diferentes mecanismos convectivos foram encontrados. O escoamento dentro do dispositivo de assistência ventricular pediátrico foi analisado por meios de simulações numéricas tridimensionais. Um modelo computacional baseado em condições de contorno especiais foi desenvolvido para modelar o escoamento pulsátil. Neste modelo, a abertura e fechamento da válvula mitral e diafragma foram representados com o uso de condições de contorno especialmente elaboradas. A força motora e o direcionamento do fluxo do diafragma foram definidos por uma distribuição de velocidades na parede do diafragma, e a abertura e fechamento da válvula mitral foram executadas por uma função de onda de velocidade que vai a zero no período sistólico. Padrões do escoamento, campos de velocidade e tensão de cisalhamento no tempo foram analisadas para avaliar o desempenho do dispositivo.