Neste trabalho é proposta uma metodologia numérica como ferramenta para avaliação de ambientes térmicos com manequins. A simulação de CFD do ambiente térmico em simulador comercial é integrada à simulação do sistema térmico do corpo humano realizada em código acadêmico. As soluções dos fluxos de calor e temperaturas nas peles são retroalimentadas e a transferência de informações é realizada via arquivo. A geometria do ambiente térmico é simplificada para minimizar os efeitos de problemas com a malha computacional na simulação de CFD, permitindo uma melhor análise do método interativo proposto. O manequim é separado em 15 segmentos cilíndricos representando cabeça, pescoço, tronco, braços, antebraços, mãos, coxas, pernas e pés. Cada segmento é subdividido em quatro quadrantes, totalizando 60 zonas, para capturar assimetrias térmicas e aerodinâmicas. Foi conduzido estudo de validação geométrica de manequim virtual pela comparação de resultados dos coeficientes térmicos obtidos com CFD e resultados obtidos de ensaios experimentais da literatura. A qualidade da malha e o tratamento de parede são discutidos. Os resultados tornam evidente que uma geometria simplificada do manequim é suficiente para estudos e avaliações de ambiente térmico e de conforto térmico quando se utiliza técnicas numéricas de CFD. Os resultados a partir da integração dos simuladores mostram que o método numérico pode ser instável nos segmentos com baixo metabolismo e baixa vazão de sangue, como nos pés e mãos. Pretende-se introduzir na metodologia proposta algum mecanismo que identifique automaticamente este fenômeno, para evitar a divergência do método e tornar a ferramenta mais robusta. / In this work is proposed a new numerical methodology as a tool for thermal comfort evaluation. This method promotes the interaction of the thermal environment simulation and the thermal system of the human body simulation. The commercial CFD simulator FLUENT R and an academician code for human body simulation are used. The solutions are fed back and the transfer is made by file. The geometry of the room is simplified to minimize the effects of problems with the computational mesh in the CFD simulation, allowing a better analysis of the proposed interactive method. The dummy is separated into 15 cylindrical segments representing head, neck, trunk, arms, forearms, hands, thighs, legs and feet. Each segment is subdivided into four quadrants, totaling 60 zones, to capture asymmetries in the heat flux field and temperature field. Was conducted a geometric validation of virtual dummy by comparing results of heat transfer coefficients from literature and CFD simulation. The mesh quality and near wall treatment are discussed. The results show that a simplified geometry of the dummy is sufficient for thermal environment studies and evaluations in CFD simulations. The results from coupled simulations show that the numerical method can be unstable in the segments with low metabolism and low blood flow, as the feet and hands. So, its intended to introduce mechanisms in the methodology to automatically identify this phenomenon and to avoid the divergence of the method to make more robust this methodology.

Este trabalho lida com o escoamento incompressível ao redor de pares de cilindros rígidos e imóveis, posicionados de forma alinhada em relação a uma corrente uniforme. São dois os objetivos desta pesquisa: o primeiro é estabelecer relações de causalidade entre características físicas do escoamento e as mudanças observadas nas forças e no campo fluido com as variações do número de Reynolds (Re) e do espaçamento entre os corpos; e o segundo consiste em compreender a influência mútua entre tridimensionalidades e interferência no escoamento. Utilizou-se o método de elementos espectrais para realizar simulações bi e tridimensionais do escoamento. Os espaçamentos entre centros (lcc) analisados vão de 1,5 a 8 diâmetros, e eles são comparados com o caso de um cilindro isolado. A faixa de Re vai de 160 a 320, compreendendo a transição na esteira. O foco foi dado nas instabilidades de pequena escala (modos A e B). Dados referentes ao número de Strouhal, coeficiente de arrasto médio, flutuação do coeficiente de sustentação e correlação axial são apresentados. Com auxílio de visualizações do escoamento, são propostos mecanismos que explicam o fenômeno de interferência, refletido no comportamento das curvas. Os resultados mostraram que simulações bidimensionais são insuficientes para se prever a combinação (Re, lcc) de inversão do arrasto. Verificou-se também que, quando o espaçamento é menor do que o crítico, o processo de transição na esteira se dá de forma diferente do observado para um cilindro isolado. / This work deals with the incompressible flow around pairs of rigid and immovable circular cylinders in tandem arrangements. There are two goals in this research: the first one is to find causality relationships between physical characteristics of the flow and the changes that are observed in the forces and in the flow field with the variation of the Reynolds number (Re) and the distance between the bodies; and the second one is to comprehend the mutual influence between three-dimensional structures and interference. The spectral element method was employed to carry out two- and three-dimensional simulations of the flow. The centre-to-centre distance (lcc) of the investigated configurations varies between 1.5 and 8 diameters, and they are compared to the isolated cylinder case. The Re range goes from 160 to 320, covering the transition in the wake. We focused in the small scale instabilities (modes A and B). Data of Strouhal number, mean drag coefficient, RMS of the lift coefficient and axial correlation are presented. With aid of flow visualizations, we propose mechanisms to explain the interference phenomenon, which is reflected in the behaviour of the graphics. The results show that two-dimensional simulations are not sufficient to predict the (Re, lcc) pair correspondent to the drag inversion point. We also verified that, in the cases where lcc is lower than the critical spacing, the transition in the wake happens in a way different from the one observed in the flow around a single cylinder.

Circuitos de convecção natural ou sistemas de circulação natural são empregados em diversas áreas da engenharia. Reatores nucleares refrigerados a água utilizam circuitos de circulação natural como método passivo de seguranca. Em situações críticas, sem qualquer controle externo, o sistema permanece em segurança por suas próprias características de funcionamento (intrinsecamente seguro). O trabalho proposto consiste em estudar numericamente o circuito de circulação natural de água, localizado no Instituto de Pesquisas Energéticas e Nucleares / Comissão Nacional de Energia Nuclear em São Paulo, por meio do uso de modelos matemáticos, objetivando determinar o padrão do escoamento em condições sem mudança de fase líquido-vapor. A comparação dos resultados de temperatura obtidos por cada um dos modelos de turbulência aos pontos instrumentados no circuito experimental, na condição transitória, revelou desvios significativos nas respostas do modelo de zero equação. Desvios intermediário foram observados nos modelos de transporte da viscosidade turbulenta (EVTE), k - ω, SST e SSG e resultados melhores foram vericados nos modelos k - ε e DES (com significativa superioridade do primeiro modelo). / Natural circulation loops apply to many engineering applications such as: water heating solar energy system (thermo-siphons), thermal management of electrical components (voltage converter), geothermal energy, nuclear reactors, etc. In pressurized water nuclear reactors, known as PWR\'s, the natural circulation loops are employed to ensure passive safety. In critical situations, the heat transfer will occur only by natural convection, without any external control or mechanical devices. This feature is desired and has been considered in modern nuclear reactor projects. This work consists of a numerical study of the natural circulation loop, located at the Instituto de Pesquisas Energeticas e Nucleares / Comissão Nacional de Energia Nuclear in São Paulo, Brazil, in order to establish the ow pattern in single phase conditions. The comparison of numerical results to experiments in transient condition revealed significant deviations for the Zero Equation turbulence model. Intermediate deviations for the Eddy Viscosity Turbulence Equation (EVTE), k - ω, SST e SSG models. And the best results are obtained by the k - ε e DES models (with better results for the k - ε model).

Este trabalho tem por objetivo estudar através de simulações numéricas a interação fluido-estrutura de um escoamento ao redor de um cilindro livre para vibrar na direção transversal do escoamento, com o objetivo de compreender melhor o fenômeno muito recorrente em estruturas oceânicas, principalmente os dutos de extração de petróleo, chamados risers. A condição a ser estudada é a de escoamento com número de Reynolds igual a 10000, em um cilindro com parâmetro de massa igual a 3,3 e de amortecimento igual a 0,0026. Nestas condições, diversos modelos serão criados para representar toda a faixa de velocidades reduzidas em que ocorre o fenômeno de sincronização. Estes parâmetros foram escolhidos devido à presença de dados experimentais para o comparativo entre as simulações e os dados experimentais, que será apresentado ao final do trabalho. Resultados como a média do coeficiente de arrasto e o RMS do coeficiente de sustentação também serão comparados com dados experimentais disponíveis na literatura e uma análise dos perfis de escoamento na região da esteira é feita para a identificação de padrões recorrentes nestes tipos de caso. / The objective of this work is to study through numerical simulations the fluid-structure interaction in a cylinder, which is free to oscillate in the cross flow direction, aiming to better comprehend the phenomena that frequently appears in offshore structures, mainly in pipelines for oil extraction, called risers, The studied condition is a flow with Reynolds 10000 around a cylinder with mass ratio of 3.3 and damping factor of 0.0026. Under these conditions, multiples models are created to represent the reduced velocities range where the synchronization phenomena happen. Those parameters were chosen due to the availability of experimental data in the literature, which allows the comparison between the simulations and experiments that is presented in the end of the work. Results such as the drag coefficient mean and the lift coefficient RMS are also compared with the experimental data available in the literature and an analysis of the flow in the wake region is done with the intention of identifying the patterns found in this type of case.

Este trabalho apresenta um estudo sobre o escoamento de ar em dutos com múltiplas saídas que são identificados como caixas sopradoras, quando instalados em máquinas desaguadoras de celulose. Estas caixas auxiliam o processo de secagem conduzindo ar aquecido até um conjunto de orifícios que produzem jatos que são impingidos sobre as folhas de celulose fabricadas. É realizada revisão bibliográfica abordando processos e equipamentos de secagem, com foco no escoamento de ar, métodos para solução numérica de escoamentos, seus modelos e parâmetros, e abordagem semi-empírica para fluxos de ar na saída de dutos. O escoamento do ar no interior da caixa sopradora e nos jatos de saída são estudados através de simulação numérica, e de modelo semi-empírico. São propostas duas alternativas de geometria para estes dutos de modo que se garanta a maior uniformidade possível na vazão e na velocidade média do ar que saí através dos orifícios Na primeira alternativa considera-se que a seção transversal da caixa sopradora se mantém constante. Na segunda se considera que a seção transversal seja reduzida de modo uniforme ao longo da caixa sopradora. Um estudo experimental sobre o comportamento do escoamento de ar para caixas sopradoras em escala reduzida apresenta resultados que possibilitam validar um conjunto de valores numéricos obtidos na simulação e também estabelecer características para o escoamento de modo a atender o processo de secagem. Para cada alternativa são adotadas duas condições de vazão de ar. Os resultados obtidos pela da simulação numérica e modelagem semi-empírica validam um perfil geométrico adequado para que a caixa sopradora produza um conjunto de jatos nas suas saídas com valores de velocidade média. Um dos principais resultados foi obtido para a caixa com 275 orificios, na qual as velocidades médias estiveram dentro do intervalo de 32,5 m/s e 33,1 m/s, que representou um ótimo resultado em termos de balanceamento e uniformidade das vazões de ar. / This study presents an analysis about airflow in ducts with multiple outlets that are identified as blow boxes, when installed in pulp dewatering machines. These boxes work in the drying process leading heated air to a set of nozzles that produce air jets that are impinged against the pulp sheet. It is performed literature review about drying processes and equipment, with a focus on air flow, methods for numerical solution of airflows, their models and parameters, and semi-empirical approach about airflow at the outlets of the duct. The airflow inside the blow box and the output jets are studied by numerical simulation, and semi-empirical model. It is proposed two alternative geometries for these ducts in order to guarantee the greatest possible uniformity in flow rate, and in the average velocity of air jets. In the first alternative is considered that the cross section of the blow box remains constant. In the second one it is considered that the cross section is reduced uniformly along the blow box. An experimental study about the behavior of the airflow in blow boxes, using scaled models, presents results that can validate a set of numerical values obtained in the simulation and also establish characteristics for the drying process. For each alternative two conditions of airflow are adopted. The results obtained through numerical simulation and semi-empirical modeling validates an appropriate geometrical profile for the blow box in order to set the values of average speed within the range required by the drying process. One of the principal results, obtained for the case of boxes with 275 holes, was the average speed within the range of 32,5 m/s and 33,1 m/s, this range represents a great result in tterms of balance and uniformity of air flow.

Dispositivos de reação miniaturizados tendem a ganhar espaço na indústria de processos químicos por elevarem o transporte de massa e de calor e a segurança dos processos. Para compreender o papel de cada elemento constituinte de um microrreator sobre seu campo de velocidades e fenômenos de mistura, foram simulados um conjunto de dois canais, uma junção em T, 30 canais em formato de serpentina e um microchip completo através da fluidodinâmica computacional. A seção transversal destes microdispositivos têm dimensões entre 100 e 300 µm, enquanto o comprimento dos canais varia de 3000 a 25190 µm. Os modelos computacionais foram discretizados por malhas hexaédricas e os campos de velocidade em estado estacionário foram calculados para vazões de alimentação entre 12,5 e 2000 µL min-1, considerando regime laminar. A mistura foi avaliada pela injeção de traçadores não-reativos e distribuição das respectivas frações mássicas. As simulações foram validadas usando microvelocimetria por imagens de partículas. Os campos de velocidade possuem magnitudes significativas apesar das dimensões reduzidas e baixas vazões de operação dos sistemas. As imagens experimentais do escoamento evidenciaram o formato parabólico do campo de velocidades e o deslocamento de seu ponto máximo nas regiões curvas causado pela força centrífuga, como estimado pelo modelo computacional. Tal força, em conjunto com as forças viscosas na parede, gera fluxos secundários no escoamento. A distribuição de traçadores não-reativos evidenciou a importância dos fluxos secundários para promover mistura na direção ortogonal ao escoamento principal, ocorrendo sob o regime estacionário nas vazões analisadas. O estudo aqui realizado evidencia o emprego da fluidodinâmica computacional como ferramenta para melhor compreensão da fluidodinâmica e como apoio ao projeto de microdispositivos. / Miniaturized reaction vessels are drawing attention of chemical industries because they promote better mass and heat transfer and also enhance process safety. To understand the relevance of each element of a microreactor on the velocity field of the equipment and the corresponding mixing processes, several microdevices were simulated using computational fluid dynamics: an assembly of two channels, a T-junction, 30 channels in a serpentine assembly and a full microreactor. The cross section of the devices is 100 - 300 µm wide and the length of the channels varies between 3000 and 25190 µm. Computational domains were discretized using hexahedral meshes and steady-state velocity fields were computed considering laminar flow for flow rates between 12,5 and 2000 µL min-1. Mixing was evaluated by injecting inert tracers and monitoring their distribution. Simulations were validated against experimental micro particle image velocimetry data. Velocities throughout the devices are relatively high despite the small dimensions of the cross sections and small flow rates. Experimental images of the flow elucidated the parabolic shape of the velocity profile and its distortion on curved segments caused by centrifugal forces, matching predictions of the computational model. Tracer maps indicated secondary flows play an important role in mixing stream perpendicular to the main flow direction. This study emphasizes the use of computational fluid dynamics as a tool for understating flow throughout microdevices and supporting their design.

Circuitos de convecção natural ou sistemas de circulação natural são empregados em diversas áreas da engenharia. Reatores nucleares refrigerados a água utilizam circuitos de circulação natural como método passivo de seguranca. Em situações críticas, sem qualquer controle externo, o sistema permanece em segurança por suas próprias características de funcionamento (intrinsecamente seguro). O trabalho proposto consiste em estudar numericamente o circuito de circulação natural de água, localizado no Instituto de Pesquisas Energéticas e Nucleares / Comissão Nacional de Energia Nuclear em São Paulo, por meio do uso de modelos matemáticos, objetivando determinar o padrão do escoamento em condições sem mudança de fase líquido-vapor. A comparação dos resultados de temperatura obtidos por cada um dos modelos de turbulência aos pontos instrumentados no circuito experimental, na condição transitória, revelou desvios significativos nas respostas do modelo de zero equação. Desvios intermediário foram observados nos modelos de transporte da viscosidade turbulenta (EVTE), k - ω, SST e SSG e resultados melhores foram vericados nos modelos k - ε e DES (com significativa superioridade do primeiro modelo). / Natural circulation loops apply to many engineering applications such as: water heating solar energy system (thermo-siphons), thermal management of electrical components (voltage converter), geothermal energy, nuclear reactors, etc. In pressurized water nuclear reactors, known as PWR\'s, the natural circulation loops are employed to ensure passive safety. In critical situations, the heat transfer will occur only by natural convection, without any external control or mechanical devices. This feature is desired and has been considered in modern nuclear reactor projects. This work consists of a numerical study of the natural circulation loop, located at the Instituto de Pesquisas Energeticas e Nucleares / Comissão Nacional de Energia Nuclear in São Paulo, Brazil, in order to establish the ow pattern in single phase conditions. The comparison of numerical results to experiments in transient condition revealed significant deviations for the Zero Equation turbulence model. Intermediate deviations for the Eddy Viscosity Turbulence Equation (EVTE), k - ω, SST e SSG models. And the best results are obtained by the k - ε e DES models (with better results for the k - ε model).

Difusores de ar é um tema de particular interesse na indústria dos sistemas de ar condicionado e climatização. O difusor swirl (ou tangencial) é um tipo de difusor já utilizado em alguns ambientes climatizados. O seu comportamento é mais conhecido em aplicações no campo da combustão, mas, em aplicações de sistemas de ar condicionado o Número de Reynolds é bem menor, não há a combustão e as condições de contorno são diferentes. Além disso, têm-se poucos estudos voltados para estes difusores num domínio 3D. Com esta motivação, o presente trabalho apresenta as simulações de um difusor tangencial em vazões típicas de aplicações de conforto térmico personalizado, utilizando a modelo de turbulência Reynolds Stress em um domínio 3D. Algumas simulações em um domínio 2D são realizadas a fim de se obter algumas características essenciais do escoamento, como abertura e comprimento do jato. Porém, comparados ao experimento, os resultados 2D precisam ser melhorados. Esquemas de discretização de maior ordem são utilizados para se avaliar o desempenho. Nas simulações no domínio 3D, verifica-se que um resultado melhor é alcançado quando se refina a malha na região central do jato, logo abaixo do difusor. Dois métodos de especificação da condição de contorno de entrada são estudados: o primeiro consiste em utilizar os dados experimentais obtidos na saída do difusor para simular o escoamento sem a geometria do difusor e o segundo simula o difusor completo, aplicando a magnitude da velocidade perpendicularmente à superfície de entrada com base na vazão calculada pelos dados do experimento do PIV (Particle Image Velocimetry). Os resultados numéricos são comparados com o experimento. Verifica-se que o método de simulação sem o difusor apresenta resultados mais precisos com relação ao experimento e apresenta maiores vantagens na simulação numérica. / Air diffusers are a topic of particular interest in the industry of acclimatization and air conditioning systems. The swirl (or tangential) diffuser is a type of device already used in some air conditioned environments. Their behavior is best known in combustion applications, but in air conditioning systems applications, the Reynolds number is much lower, there is no combustion and the boundary conditions are different. In addition, there have been few studies on these diffusers in 3D domain. With this motivation, this work presents simulations of a tangential flow diffuser for applications of thermal comfort. The numerical study uses the Reynolds stress turbulence model in a 3D domain. Some simulations in a 2D domain are performed in order to obtain some essential features of the flow, as the width and length of the jet. However, compared to the experiment, the 2D results need to be improved. Higher order discretization schemes are used to evaluate performance. In 3D domain simulations, it is verified that a better result is achieved when the mesh is refined in the jets central region, just below the diffuser. Two methods of the inlet boundary conditions are studied: the first consists of using the experimental data obtained at the exit of the diffuser to simulate the flow without the geometry of the diffuser and the second method simulates the diffuser completely, applying the velocity magnitude perpendicular to the inlet surface based on the calculated flow rate with experimental data of PIV (Particle Image Velocimetry). The numerical results are compared with experiment. It is noted that a simulation method without the geometry of the diffuser provides more accurate results with the experiment and has major advantages in the numerical simulation.

Orientador: José Luiz Gasche / Resumo: Originalmente, stents foram projetados para agir como barreiras mecânicas, impedindo a herniação de coils para a artéria-mãe. Recentemente, estudos mostraram que a atual geração de stents intracranianos auto-expansíveis altera a geometria local das artérias: um fenômeno com efeitos hemodinâmicos em parte incompreendidos. Nós realizamos simulações numéricas para avaliar a influência da remodelagem arterial induzida por stent sobre a hemodinâmica em aneurismas intracranianos. As simulações foram realizadas utilizando o software open-source OpenFOAM. O sangue foi modelado como fluido Newtoniano incompressível e as paredes arteriais foram consideradas rígidas. Para quantificar as alterações hemodinâmicas, avaliamos os parâmetros wall shear stress, WSS, e oscillatory shear index, OSI. Quatro geometrias reais de aneurismas intracranianos em bifurcações foram utilizadas. Em um aneurisma na bifurcação da artéria comunicante anterior (ACoA), um stent foi implantando, levando ao endireitamento das artérias que o receberam. Após o procedimento, os níveis de WSS e OSI aumentaram aproximadamente 60% e 25%, respectivamente. Em dois aneurismas em bifurcações da artéria cerebral média (MCA), dois stents foram implantados em uma configuração em “Y”, resultando em um endireitamento de ambas as artérias-filhas. O WSS máximo na superfície do aneurisma aumentou aproximadamente 5% em um dos casos e 22% no outro. Em outro aneurisma em uma bifurcação da MCA, um stent foi implantado, resultando no en... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Stents were first designed to act as mechanical barriers, preventing coil herniation into the parent artery. The current generation of self-expanding intracranial stents has recently been shown to change the local vascular geometry, a phenomenon with unclear hemodynamic effects. We carried out numerical simulations to assess the role of stent-induced vascular remodeling in modifying intraaneurysmal hemodynamics. Simulations were performed using the open-source software OpenFOAM. Blood was assumed to behave as an incompressible Newtonian fluid; vessel walls were assumed to be rigid. Wall shear stress, WSS, and oscillatory shear index, OSI, were evaluated to quantify the hemodynamic changes in the aneurysm sac. Four pre- and post-stent patient-specific geometries of intracranial bifurcation aneurysm were used. In one aneurysm at the anterior communicating artery (ACoA) bifurcation, a single stent was deployed, resulting in straightening of the host vessels. After stenting, WSS and OSI increased by approximately 60% and 25%, respectively. In two aneurysms at middle cerebral artery (MCA) bifurcations, two stents in a “Y” configuration were deployed, resulting in straightening of both daughter arteries. The maximum WSS on the aneurysm surface increased by approximately 5% in one case and 22% in the other. In another aneurysm at a bifurcation of the MCA, a single stent was deployed, resulting in straightening of the host vessels. After stenting, WSS and OSI reduced by approximately... (Complete abstract click electronic access below) / Mestre

Dispositivos de reação miniaturizados tendem a ganhar espaço na indústria de processos químicos por elevarem o transporte de massa e de calor e a segurança dos processos. Para compreender o papel de cada elemento constituinte de um microrreator sobre seu campo de velocidades e fenômenos de mistura, foram simulados um conjunto de dois canais, uma junção em T, 30 canais em formato de serpentina e um microchip completo através da fluidodinâmica computacional. A seção transversal destes microdispositivos têm dimensões entre 100 e 300 µm, enquanto o comprimento dos canais varia de 3000 a 25190 µm. Os modelos computacionais foram discretizados por malhas hexaédricas e os campos de velocidade em estado estacionário foram calculados para vazões de alimentação entre 12,5 e 2000 µL min-1, considerando regime laminar. A mistura foi avaliada pela injeção de traçadores não-reativos e distribuição das respectivas frações mássicas. As simulações foram validadas usando microvelocimetria por imagens de partículas. Os campos de velocidade possuem magnitudes significativas apesar das dimensões reduzidas e baixas vazões de operação dos sistemas. As imagens experimentais do escoamento evidenciaram o formato parabólico do campo de velocidades e o deslocamento de seu ponto máximo nas regiões curvas causado pela força centrífuga, como estimado pelo modelo computacional. Tal força, em conjunto com as forças viscosas na parede, gera fluxos secundários no escoamento. A distribuição de traçadores não-reativos evidenciou a importância dos fluxos secundários para promover mistura na direção ortogonal ao escoamento principal, ocorrendo sob o regime estacionário nas vazões analisadas. O estudo aqui realizado evidencia o emprego da fluidodinâmica computacional como ferramenta para melhor compreensão da fluidodinâmica e como apoio ao projeto de microdispositivos. / Miniaturized reaction vessels are drawing attention of chemical industries because they promote better mass and heat transfer and also enhance process safety. To understand the relevance of each element of a microreactor on the velocity field of the equipment and the corresponding mixing processes, several microdevices were simulated using computational fluid dynamics: an assembly of two channels, a T-junction, 30 channels in a serpentine assembly and a full microreactor. The cross section of the devices is 100 - 300 µm wide and the length of the channels varies between 3000 and 25190 µm. Computational domains were discretized using hexahedral meshes and steady-state velocity fields were computed considering laminar flow for flow rates between 12,5 and 2000 µL min-1. Mixing was evaluated by injecting inert tracers and monitoring their distribution. Simulations were validated against experimental micro particle image velocimetry data. Velocities throughout the devices are relatively high despite the small dimensions of the cross sections and small flow rates. Experimental images of the flow elucidated the parabolic shape of the velocity profile and its distortion on curved segments caused by centrifugal forces, matching predictions of the computational model. Tracer maps indicated secondary flows play an important role in mixing stream perpendicular to the main flow direction. This study emphasizes the use of computational fluid dynamics as a tool for understating flow throughout microdevices and supporting their design.