Orientador: Carlos Roberto Minussi / Coorientador: Geraldo de Freitas Maciel / Banca: Anna Diva Plasencia Lotufo / Banca: Mara Lúcia Martins Lopes / Banca: Luis Miguel Chagas Costa Gil / Banca: Mônica Feijó Naccache / Resumo: Nesta tese propõem-se investigar os escoamentos com superfície livre de fluidos newtonianos e hiperconcentrados de reologia não-newtoniana, com aplicação no âmbito de catástrofes naturais. Esta temática apresenta duas vertentes: uma visão das ciências sociais sobre riscos e catástrofes e seus impactos de ordem econômica e social; uma segunda vertente no plano da engenharia (visão mais mecanicista) que trata do estudo dos escoamentos de fluidos lamosos em corredores de lavas torrenciais, canais e encostas íngremes, com presença de um tipo específico de instabilidade na superfície livre, as roll waves. O fenômeno roll wave é caracterizado por sua configuração periódica, com frentes de ondas que se deslocam com elevada velocidade. Assim, a capacidade destrutiva do escoamento é potencializada com a presença destas instabilidades. Nesta tese procura-se identificar as condições de geração destas ondas e suas propriedades através de modelagem matemática e numérica. Para se estabelecer um modelo matemático que represente o fenômeno em questão, parte-se das equações fundamentais de conservação da massa e quantidade de movimento (equações de Cauchy), levando-se em consideração as hipóteses de águas rasas, com introdução do modelo reológico de Herschel Bulkley na parte viscosa do tensor de tensões. A técnica de análise de estabilidade linear é usada no equacionamento em variáveis adimensionais, tendo sido estabelecido dois critérios para a geração das roll waves: o primeiro está relacionado com número de Froude > , em que é o número de Froude mínimo necessário para geração de roll waves e o segundo é determinado em função da frequência de perturbação < , onde é a frequência máxima de perturbação possível à geração destas ondas. Identificadas tais condições, um... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: This thesis follows the investigation of free surface flows of Newtonian and hyperconcentrated (non-Newtonian) fluids on natural disasters context. This research presents two approaches, being the first approach a vision of natural hazards and disasters and their social consequences; and the second, the engineering plan (mechanical point of view) which is the study of muddy fluids in corridors, steep slopes and open channels with the presence of a specific type of instability on the free surface known as the roll waves. Roll wave phenomenon is characterized by its periodic configuration with moving wavefronts at high speed. Thus, the destructive capability of flow is enhanced by the presence of roll waves. This thesis attempts to identify generation conditions for these waves and also their properties, both through mathematical and numerical modeling. In order to establish a mathematical model that represents the phenomenon, the system of equations composed by mass and momentum conservation (Cauchy's equations) is used, taking into account the assumptions of shallow water, and introducing the Herschel-Bulkley rheological model on the viscous stress tensor. The technique of linear stability analysis is used on the system on dimensionless variables. Two criteria are established for the generation of roll waves: the first is related to Froude number > , where is the minimum Froude number for the generation of roll waves; and the second is determined by the disturbance < , where is the maximum disturbance frequency possible to the generation of these waves. Identified these conditions, a mathematical model (1D) determines the main properties of roll waves, such as amplitude, length and propagation velocity. This model is established by rewriting the equations in a moving coordinate system. As this model i... (Complete abstract click electronic access below) / Doutor

A presente dissertação investiga o comportamento do escoamento em um túnel de cavitação, através de simulações fluido-dinâmicas computacionais, excluindo a bomba de circulação. Para tanto um extenso estudo buscou selecionar a abordagem numérica e configurações de simulação mais adequadas, de uma forma a reproduzir as características hidrodinâmica inerentes à operação do túnel. Portanto, as informações referentes à perda de carga, uniformidade no perfil de velocidade, tendência de cavitação, descolamento e altura da camada limite foram apreciadas. Esse estudo foi direcionado por dados provenientes de métodos empíricos da literatura, resultados experimentais do próprio túnel pesquisado e de outros túneis. Em geral, os resultados foram satisfatórios, pois a perda de carga estimada foi semelhante ao obtido experimentalmente. Além de serem identificados os trechos com efeito desfavorável na uniformidade do escoamento, foram determinadas a distribuição de pressão nas aletas e o perfil de velocidade incidente na bomba. / This study deals with a cavitation tunnel flow (except the pump region) through computational fluid dynamic simulations. Therefore, previously one made a literature survey to aid the selection of the most appropriate techniques and simulation settings to simulate the hydrodynamic characteristics of the tunnel. The pressure drop, uniformity of the velocity profile, existence of cavitation, detachment and height of the boundary layer were evaluated. This study used empirical methods and also experimental results. The results were satisfactory, as long as the estimated pressure drop was close to the one obtained experimentally. In addition, it was forecasted low flow uniformity in some parts of the tunnel, the pressure distribution on the fins and the incident velocity profile at the pump.

No presente trabalho, foi desenvolvido um modelo para o escoamento em reator de membrana do tipo contator ativo, através da abordagem de dinâmica dos fluidos computacional (CFD), utilizando o código comercial ANSYS FLUENT. O modelo incluiu todo o módulo de membrana constituído por uma membrana tubular e um casco metálico. A reação modelo estudada foi a hidrogenação parcial de 1,5-ciclooctadieno, realizada pelo bombeamento da mistura reacional, dissolvida em n-heptano, através da membrana, a partir das extremidades do tubo. Como catalisador, considerou-se a presença de nanopartículas de Pd impregnadas na membrana. O meio poroso foi aproximado por leito granular representado pela equação de Ergun, tendo como parâmetros a porosidade e o tamanho de grão da membrana de alfa-Al2O3 . O valor para o tamanho de grão foi adotado como equivalente ao diâmetro de partícula determinado com uso do código aberto de estereologia ImageJ, do instituto Nacional de Saúde dos Estados Unidos. O modelo de turbulência adotado foi o RNG k-epsilon. Um estudo de sensibilidade incluiu simulações comparando escoamento desprezando reações como escoamento reativo, variação da velocidade, alteração da saída do fluxo e ativação de modelo de turbulência no meio poroso. Foram realizadas simulações de defeitos estruturais na membrana, correspondendo a regiões de porosidade alterada, com e sem perda de sua uniformidade azimutal. Conclui-se que a presença de defeitos estruturais que afetem a uniformidade azimutal da membrana pode resultar em sensível alteração do escoamento em CMRs. / This study focused on the development of a model for the flow in a reactor membrane of the type active contactor, approached through computational fluid dynamics (CFD), using the commercial code ANSYS FLUENT. The model included the entire membrane module, consisted of a tubular membrane and a metal shell. The model reaction studied was the partial hydrogenation of 1,5-cyclooctadiene initiated by the pumping of the reaction mixture, dissolved in n-heptane, through the membrane, from the ends of the tube. As a catalyst, the study considered the presence of impregnated Pd nanoparticles in the membrane. The porous medium was approximated by a granular bed as represented by the Ergun equation, having as parameters the porosity and the grain size of the alfa-Al2O3 membrane. The value for the grain size was adopted as equivalent to particle diameter determined through the open source stereology software ImageJ, of the National Institute of Health USA. The turbulence model used was the RNG k-epsilon. A sensitivity study included simulations of flow neglecting and including reactions, speed variation, change the flow outlet and activation of turbulence model in the porous media. Simulations of structural defects in the membrane were performed, defining regions of porosity changes with and without loss of azimuthal uniformity. The conclusion was that the presence of structural defects that affect the azimuthal uniformity of the membrane can result in marked alteration of the flow regime in CMRs.

The fow around bluff bodies is an essential topic in fuid dynamics. This fow is characterized by large vortical fow regions separating from the surface of the bluff body, and they cause oscillating drag and lift forces on the structure. The fow around an infinite long cylinder is a well-known case being studied in the literature. However, a cylinder with low aspect-ratio piercing the free surface has not been studied much although such an arrangement can be found in many floating offshore structures. In this thesis the results of CFD calculations are presented for a fixed, free surface piercing cylinder with an aspect-ratio L/D equals to 2. The Reynolds number was equal to 4:3 x 104 indicating that the flow is in the sub-critical turbulent regime. An extensive methodology of verication and validation is followed to investigate the reliability of the results. To investigate the effect of the free surface on the calculated hydrodynamic loads, two approaches have been considered: a double-body symmetry condition and an interface capturing Volume-of-Fluid (VoF) method. Additionally, two turbulence models are investigated: a two-equation turbulence model; a non-linear Explicit Algebraic Reynolds Stress Model (EARSM); and the Improved Delayed Detached Eddy Simulation (IDDES) turbulence model. The results are presented in terms of integral results (drag and lift coefficients) and flow visualizations. Based on the results of the cases in which the free surface was modeled as a double body symmetry boundary condition, it is concluded that the model is not suitable for this type of flow as the model damps out the flow dynamics due to over-production of eddy-viscosity. Hence, the characteristic oscillating lift forces are not captured using this turbulence model. However, this turbulence model showed good agreements regarding the flow fields in comparison with experimental PIV measurements. Results of the case modeled with EARSM turbulence model shows better agreement with the experimental results compared with the turbulence model. In the cases where the free-surface is considered, results with the EARSM turbulence model show similar results for the drag forces whereas the lift uctuations were one order of magnitude smaller, compared with the double body case. Lastly, the results using the IDDES turbulence model and free-surface VoF modeling are shown to produce the best comparison with the experimental results, regarding both, integral values and flow field results. / O escoamento ao redor de corpos rombudos é um tópico essencial na dinâmica de fluidos. O escoamento é caracterizado por regiões com grande vorticidade que se separam do corpo e causam oscilações das forças de arrasto e sustentação sobre a estrutura. O escoamento ao redor de cilindros longos é um tema que tem sido extensivamente estudado com muitos trabalhos encontrados na literatura. Entretanto, o cilindro com baixa razão de aspecto perfurante à superfície livre é um caso pouco estudado, apesar desta estrutura ser encontrada em várias estruturas oceânicas flutuantes. Esta dissertação apresenta cálculos numéricos para o escoamento ao redor de um cilindro fixo, que trespassa a superfície livre com razão de aspecto L/D igual a 2. O problema é estudado em um regime subcrítico de turbulência, com número de Reynolds igual a 4:3 x 104. Uma vasta metodologia de verificação e validação foi seguida para avaliar a confiabilidade dos resultados obtidos numericamente. Para investigar os efeitos da superfície livre nas cargas hidrodinâmicas, duas abordagens s~ao consideradas: condição de simetria de duplo corpo e um método de captura de interface Volume of Fluid. Além disso, dois modelos de turbulência foram investigados: o modelo não linear Explicit Algebraic Reynolds Stress Model (EARSM), e o modelo de turbulência Improved Delayed Detached Eddy Simulation (IDDES). Os resultados relacionados aos coeficientes de arrasto e sustentação são apresentados a partir de análise estatística, complementados através de ilustrações que permitem visualizar os campos de escoamento e pressão. Com base nos resultados de casos em que a superfície livre é modelada com uma condição de contorno de simetria, conclui-se que o modelo de turbulência não é adequado para este tipo de escoamento, pois o modelo amortece a dinâmica do escoamento devido à superprodução de viscosidade turbulenta. Consequentemente, as oscilações na força de sustentação não são capturadas usando este modelo. Entretanto, resultados dos campos médios do escoamento mostram-se concordantes com imagens experimentais obtidas com técnicas de PIV - Particle Image Velocimetry. Resultados do caso modelado com o modelo de turbulência EARSM mostram melhores concordâncias na comparação dos parâmetros estatísticos com experimentos do que o modelado com o modelo EARSM. Nos casos em que a superfície livre é modelada com o método VoF, o modelo de turbulência EARSM mostra resultados semelhantes para o arrasto, enquanto as flutuações da sustentação apresentam-se uma ordem de grandeza menores, quando comparadas ao caso de duplo corpo. Resultados usando o modelo de turbulência IDDES e VoF apresentam melhores comparações aos resultados experimentais, tanto para os parâmetros estatísticos quando para as visualizações do escoamento.

A simulação numérica de escoamentos reativos, como a combustão, tem um caráter altamente não-linear devido a presença de diversas reações químicas que acontecem entre as espécies que descrevem o processo de oxidação do combustível. Além disso, tais processos ocorrem a nível molecular, tornando o sistema de equações governantes rígido, o que implica na necessidade de esquemas numéricos de alta ordem bem como malhas finas e passo de tempo pequeno, aumentando consideravelmente o custo computacional. Neste sentido, o uso de mecanismos de oxidação detalhados na simulação numérica é proibitivo, e técnicas de redução química são necessárias de modo a desenvolver modelos reduzidos com menos variáveis e rigidez moderado, mantendo a precisão e abrangência do modelo detalhado. O objetivo do presente trabalho é obter uma comparação dos resultados obtidos para duas técnicas de redução química diferentes, Directed Relation Graph - DRG, baseada no desenvolvimento de mecanismos esqueletos, e a Reaction Diffusion Manifolds - REDIM, baseada na separação das escalas de tempo. Como validação dos modelos propostos, simulações numéricas 1D de chamas pré-misturadas e não pré-misturadas, bem como de reatores homogêneos, são desenvolvidas. Além disso, uma estratégia que une as duas técnicas de redução é apresentada, com o objetivo de ser aplicada em mecanismos cinéticos grandes. / Numerical simulation of reactive flows, such as combustion, has a highly non-linear character due to the presence of several chemical reactions that occur among the chemical species that describe the process of fuel’s oxidation. Besides, such processes occur at a molecular level, making the system of governing equations stiff, which implies in the need of high order numerical schemes as well as fine meshes and small time step, enhancing considerably the computational cost. In this sense, the use of detailed oxidation mechanisms in the numerical simulation is prohibitive, and chemical reduction techniques are needed in order to develop reduced models with less variables and moderate stiffness, while maintaining the accuracy and comprehensiveness of the detailed model. The objective of the present works if to obtain a comparison between two chemical reduction techniques, the Directed Relation Graph - DRG, based on the skeletal mechanisms generation, and the Reaction Diffusion Manifolds - REDIM, based on the separation of time scales. As validation of the proposed models, one-dimensional numerical simulations of premixed and non-premixed flames, as well as homogeneous reactors, are carry out. Besides, a coupled methodology between DRG and REDIM is presented, that will provide a useful tool for simulation of fuels with very large detailed kinetic mechanisms.

Na matriz elétrica brasileira, as termelétricas vêm sendo utilizadas como fontes de energia complementar. Para minimizar os impactos ambientais na geração de energia a partir do carvão mineral, diversos estudos apontam para o uso da tecnologia de leito fluidizado como uma alternativa flexível, eficiente e limpa. Neste tipo de reator a fluidodinâmica determina a distribuição das partículas e do gás, o que influencia na mistura entre as fases, e consequentemente, na eficiência do processo, porém, não é totalmente compreendida. Baseando-se nisto, esta tese tem como objetivo caracterizar a fluidodinâmica do leito fluidizado borbulhante, a partir de modelos matemáticos e análise experimental. Busca-se desenvolver uma maior compreensão da influência fluidodinâmica no desempenho da planta em processos de combustão do carvão mineral brasileiro, o qual possui alto teor de cinzas. Para isso, os modelos matemáticos foram validados na planta de bancada a temperatura ambiente, e posteriormente, aplicados para a planta piloto a 850 °C. O desempenho de cada operação foi avaliado a partir do rendimento de conversão e teor de carbono residual nas cinzas para diferentes alturas de leito e velocidade superficial do gás. Os resultados obtidos mostram que para a obtenção de um maior tempo de residência, a velocidade superficial do gás mais adequada depende da altura do leito fixo. O teor de carbono residual nas cinzas foi menor para o leito mais alto, o que é explicado pelas características do combustível utilizado e pelo ponto de alimentação do carvão. Não foi observada relação direta entre os coeficientes de troca gasosa para a transferência de massa no leito, com o teor de carbono residual nas cinzas, ou com o rendimento de conversão. Para uma análise mais completa da influência fluidodinâmica no desempenho da planta, a região de freeboard deve ser incluída na modelagem. A análise em planta de bancada a temperatura ambiente foi fundamental para compreender como ocorre a movimentação das partículas e das bolhas. As variações na altura do leito e velocidade do gás geraram diferentes variações nos coeficientes de troca gasosa para a transferência de massa, no tempo de residência e no teor de carbono nas cinzas. / In the Brazilian electric matrix, thermoelectric power plants are used as complementary energy sources. In order to minimize the environmental impacts of coal power generation, several studies point to the use of fluidized bed technology as a flexible, efficient and clean alternative. In this type of reactor, the fluid dynamics determines the distribution of particles and the gas, which influences in the mixing between the phases, and consequently, in the efficiency of the process, however, is not fully understood. Based on this information, this thesis aims to characterize the fluid dynamics of the bubbling fluidized bed, from mathematical models and experimental analysis. Aim to develop a better understanding of the fluid dynamics influence on the plant performance in the combustion processes of the Brazilian coal, which has a high ash content. For this purpose, the mathematical models were validated in the bench plant at room temperature, and later, applied to the pilot plant at 850 °C. The performance of each operation was evaluated from the combustion efficiency and residual carbon content in the ashes for different bed heights and gassurface velocity. The obtained results show that to attain a longer residence time the most suitable gas-surface velocity depends on the height of the fixed bed. The residual carbon content in the ashes was lower in the upper bed, which is explained by the characteristics of the fuel used and the point of coal feeding. Besides, was not observed a direct relation between the gas exchange coefficients for mass transfer in the bed with the residual carbon content in the ashes or with the combustion efficiency. Toward a more complete analysis of fluid dynamics influence on the plant performance, the freeboard region should be included in the modeling. The analysis in bench plant at room temperature was fundamental to understand how the movement of particles and bubbles occurs. The variations in bed height and gas velocity generated different variations in the gas exchange coefficients for mass transfer, residence time and carbon content in the ashes.

O presente trabalho apresenta uma investigação experimental sobre a influência da presença de uma zona de recirculação interna (ZRI) no processo de combustão de gás natural. Basicamente os experimentos foram conduzidos para verificar a influência da ZRI no limite de "blow out", balanço de energia e emissões de CO e NO. Um queimador com configuração especial foi projetado para permitir mudanças na intensidade da ZRI e os experimentos foram conduzidos em uma câmara de combustão horizontal refrigerada e de escala laboratorial. A intensidade da ZRI foi caracterizada pelo número de "swirl" e o número de Reynolds do jato de gás natural, ambos calculados pela geometria do queimador e condições dos escoamentos dos reagentes. Os resultados mostraram que aumentando o número de "swirl" é possível estender o limite de "blow out", principalmente para baixos números de Reynolds. Outra observação importante é a tendência da redução do NO quando a intensidade da ZRI aumenta e a emissão desse poluente não é muito influenciada pelo número de Reynolds do jato. Contudo, o número de Reynolds tem uma importante influência na emissão de CO, maior números de Reynolds aumentam o processo de mistura entre o ar e o jato de combustível por efeito de carreamento. Mas aumentando a intensidade de ZRI para minimizar a emissão de NO pode produzir altos níveis de CO. Desta forma, os resultados apontaram que a ZRI ideal para conciliar as emissões de CO e NO é uma solução de compromisso.

A presente proposta tem como objetivo investigar numericamente a influência de uma inserção porosa em um tubo com contração súbita, no qual duas contrações abruptas diferentes foram empregadas com razão entre as seções de 0,285 e 0,10. Casos com escoamento em regime laminar e turbulento são estudados. Na simulação numérica, são empregadas as equações macroscópicas de conservação de massa, quantidade de movimento e de grandezas turbulentas baseadas no conceito da dupla decomposição, as quais são válidas em todo o domínio computacional. Para o fechamento das equações, modelos de turbulência macroscópico k-e lineares de Alto e Baixo Reynolds e não linear de Alto Reynolds são empregados. Primeiramente, os resultados da simulação, sem a inserção porosa, são comparados com dados obtidos de experimentos existentes na literatura, onde o principal parâmetro de comparação é o coeficiente de perda de carga localizada (kc) devido às súbitas mudanças de seção. Em seguida, o escoamento no duto é simulado numericamente com uma inserção porosa. Parâmetros como porosidade, permeabilidade e espessura são variados, onde, para cada caso, as perdas de carga no duto e o comportamento do escoamento na região de contração são analisados e comparados com os respectivos casos sem inserção porosa. No escoamento em dutos com contração súbita, observa-se uma bolha de recirculação logo após a contração que causa uma diminuição da área útil de passagem de fluido, fenômeno conhecido como vena contracta. A contração adicional (vena contracta) gera uma perda de carga ainda mais elevada, assim uma inserção porosa colocada logo após a contração pode diminuir ou eliminar essa recirculação reduzindo seu efeito no aumento da perda de carga. Porém, por outro lado, a própria inserção porosa gera uma perda de carga adicional. Além disso, em muitas situações na engenharia é desejável eliminar ou reduzir a recirculação, bem como uniformizar o escoamento à jusante da contração no duto. No entanto, dependendo de certos parâmetros, tais como porosidade, espessura e permeabilidade da inserção porosa, a perda de carga gerada pode ser muito elevada. Logo, este trabalho visa investigar o problema de forma a obter um melhor entendimento da interdependência dos parâmetros que compõem o meio poroso no comportamento global do escoamento. No cálculo da perda de carga localizada no duto para os casos sem inserção porosa em escoamento turbulento, os resultados obtidos utilizando o modelo de turbulência não linear apresentam uma melhor concordância com os resultados experimentais encontrados na literatura se comparados aos obtidos através dos modelos lineares. Para os casos com inserção porosa, os resultados mostram que apesar da redução ou eliminação da recirculação, as perdas de carga são sempre mais elevadas se comparadas aos casos sem inserção porosa. Além disso, este trabalho mostra o quanto a bolha de recirculação é amortecida e, ao mesmo tempo, o quanto de perda de carga adicional é gerada dependendo de cada inserção porosa considerada, o que pode ser de extrema utilidade em muitos projetos termo-mecânicos, onde, por exemplo, se deseja reduzir ou eliminar a bolha de recirculação e aumentar a troca de calor.

Procurar uma turbina axial eficiente é evitar desperdício de energia, já que o aumento da eficiência dos componentes de uma turbina a gás refletem em uma máquina mais eficiente, que fornece trabalho com um menor custo de combustível. Com a dinâmica dos fluidos computacional é possivel obter informações quantitativas do desempenho da turbina axial pela simulação do escoamento entre suas pás. Permite também visualizar detalhadamente os pontos onde há separação, ondas de choque, gradientes elevados; fenômenos diretamente relacionados à geração de perdas. Este é o estudo da análise do desempenho de turbinas axiais de alto desempenho utilizando uma ferramenta para simulação computacional do seu escoamento tridimensional e viscoso, em regime permanente. Um software comercial foi empregado para resolver as equações de Navier-Stokes com média de Reynolds com modelo de turbulência Spalart-Allmaras. A simulação em diversos pontos de operação permitiu a construção de um mapa de desempenho e sua comparação foi feita com o gerado usando uma metodologia de previsão de desempenho que contabiliza as perdas por meio das correlações de Ainley-Mathieson e Kacker-Okapuu.

Injetores supersônicos são instalados em túneis de vento transônicos com a finalidade de estender o envelope operacional baseado no número de Reynolds, sem demandar maior potência do compressor principal. O objetivo deste trabalho é estudar numericamente o processo de injeção do Túnel Transônico Piloto do CTA. A câmara de injeção contém cinco bicos injetores localizados no piso e cinco no teto do túnel, os quais operam com número de Mach 1,9, fornecendo quantidade de movimento adicional à corrente principal do circuito do túnel. Tal câmara na verdade é a própria seção de transição do túnel, a qual, por sua vez, localiza-se imediatamente antes do difusor. Devido à grande diferença entre as dimensões dos injetores e da seção transversal do túnel, onde os mesmos encontram-se montados, o tratamento tem que ser necessariamente tridimensional. Para tanto foi desenvolvido um código numérico baseado nas equações de Navier-Stokes com média de Reynolds, seguindo-se o princípio do algoritmo diagonal em diferenças finitas, e efeitos de turbulência foram previstos por meio do esquema de Spalart e Allmaras. Várias hipóteses simplificadoras foram introduzidas para tornar o problema factível e a integração numérica é realizada dividindo o domínio de cálculo em sub-regiões. Mesmo assim, devido ao grande número de pontos, foi aplicada uma técnica de malhas seqüenciais para economizar o tempo computacional. O escoamento na região de mistura foi simulado com sucesso e os resultados foram todos consistentes. Com os valores dos parâmetros conhecidos em todo o domínio computacional, foram calculados o coeficiente de perda de carga, o ganho e a eficiência do processo de injeção - parâmetros estes de grande importância na engenharia de túneis de vento. Entre diversos aspectos físicos muito interessantes, destaca-se a formação de domos de choque e expansão ao longo da direção longitudinal no processo de mistura.