Métodos analíticos, numéricos e experimentais para o cálculo de ondas de impacto em meios líquidos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Entre Duas Margens. Os Portugueses no Golfo Pérsico (1623-1653)

Dissertação apresentada para cumprimento dos requisitos necessários à obtenção do grau de Mestre em História Moderna e dos Descobrimentos

Modelo do Centro Naval Ensino a Distância – Dicotomia dos sistemas de gestão informação e da aprendizagem

Dissertação apresentada para cumprimento dos requisitos necessários à obtenção do grau de Mestre em Gestão de Sistemas de e-learning

Desenvolvimento de formulações de elementos finitos para barras sujeitas a cargas transitórias

Com o presente trabalho pretende-se formular, implementar e validar duas classes de elementos finitos não-convencionais para problemas elastoestáticos e elastodinâmicos (harmónicos e transitórios) envolvendo barras solicitadas por cargas axiais. O desempenho numérico dos elementos não convencionais é estudado para uma larga gama de situações de interesse prático e comparado com o dos elementos finitos conformes de deslocamento (convencionais). A resolução de problemas transitórios envolve a integração no tempo e no espaço das equações diferenciais governativas, bem como a imposição das respetivas condições iniciais e de fronteira. A metodologia de integração no tempo adotada neste trabalho é baseada no método de Newmark. A resolução de problemas estáticos e harmónicos não carece de integração no tempo, ou a mesma é feita de forma trivial. Concluída a discretização no tempo, a segunda fase da resolução envolve a integração no espaço de cada uma das equações discretizadas, nomeadamente através do método dos elementos finitos. Para esse efeito, apresentam-se as formulações relativas aos elementos finitos convencionais, híbridos e híbridos-Trefftz. As três formulações têm como ponto de partida a forma fraca da equação diferencial de Navier, que é imposta utilizando o método de Galerkin. A principal diferença entre os elementos convencionais e não-convencionais prende-se com a maneira como são impostas as condições de fronteira de Dirichlet e as condições de compatibilidade nas fronteiras interiores. Os elementos não convencionais são implementados numa plataforma computacional desenvolvida de raiz no ambiente Matlab. A implementação é feita de maneira a permitir uma definição muito geral e flexível da estrutura e das respetivas ações, bem como das discretizações no tempo e no espaço e das bases de aproximação, que podem ser diferentes para cada elemento finito. Por fim, efetuam-se testes numéricos com o objetivo de analisar os resultados obtidos com os elementos não convencionais e de os comparar com as respetivas soluções analíticas (caso existam), ou com os resultados obtidos utilizando elementos convencionais. É especialmente focada a convergência das soluções aproximadas sob refinamentos da malha (h), no espaço e no tempo, e das funções de aproximação (p), sendo que o uso simultâneo dos dois refinamentos parece conduzir mais rapidamente a soluções próximas da solução exata. Analisam-se também problemas complexos, envolvendo propagação de ondas de choque, com o fim de se efetuar uma comparação entre os elementos convencionais, disponíveis no programa comercial SAP2000, e os elementos não convencionais fornecidos pela plataforma computacional desenvolvida neste trabalho.

Tratamiento de la litiasis urinaria en pacientes pediátricos mediante litotricia extracorpórea

La patologia litiàsica urinaria infantil afecta a un 5 –10% d’aquesta població. El tractament de la litiasi en edat infantil inclou desde tractament mínimament invasiu com la litotrípsia extracorpòria per ones de xoc (LEOX) a tractament quirúrgic (endourològic o obert). Es va realitzar un estudi retrospectiu amb un total de 82 pacients en edat pediàtrica amb patologia litiàsica tractada mitjançant litotrípsia per ones de xoc. Es van analitzar de forma descriptiva els resultats sobre aquest grup de pacients. Dels 82 pacients estudiats, el 54% van ser nens i el 46% nenes, amb una edat mitjana de 8.6 anys. Es va analitzar la resposta a litotrípsia segons tamany i grups d’edat per estrats observant que no hi havia diferències significatives respecte al número de sessions necessàries por grups d’edat però amb tendència a la significació en el grup de menor edat (de 0 a 3 anys). Tanmateix, aquest grup de pacients presentaven comparativamente litiasi de major tamany que el grup de més edat (de 16 a 20 anys). Les complicacions es van presentar en només 10 pacients (8%) considerant com a tals la presència de fragments obstructius, dolor lumbar o carrer litiàsic. Després d’analitzar els resultats es pot concloure que la litotrícia per ones de xoc constitueix un tractament eficaç pel tractament de la litiasi independentment de la seva localització, no trobant diferències per grups d’edat i aconseguint bones tases de “stone free” inclús per litiasi coraliformes en pacients en edat infantil i amb escasses complicacions.

Tratamento cirúrgico da litíase vesical: revisão de literatura

Cálculos vesicais são raros e a maioria dos casos ocorre em homens adultos com obstrução infravesical. Atualmente, existem poucos dados sobre o melhor tratamento desta doença. O objetivo desta revisão foi discutir alguns aspectos da patogênese e abordar o tratamento da litíase vesical. Uma ampla pesquisa na base de dados da "National Library of Medicine"/Pubmed foi realizada com os seguintes unitermos e descritores: vesical ou bexiga associados a cálculo, pedra ou litíase, e cistolitotripsia. Cento e setenta e um artigos foram identificados. Os artigos foram avaliados independentemente por dois revisores com experiência em urolitíase. Foram incluídos quando os resultados, complicações e seguimento foram claramente reportados. No final, 32 estudos preencheram os critérios de inclusão. Nota-se que diversas opções para o tratamento da litíase vesical estão disponíveis, porém não há estudos randomizados comparando-as. Diferentes taxas de pacientes livres de cálculo são descritas, sendo: litotripsia extracorpórea com ondas de choque (75-100%), cistolitotripsia transureteroscópica (63-100%), cistolitotripsia percutânea (89-100%) e cirurgia aberta (100%). O acesso percutâneo apresenta menor morbidade com resultados semelhantes ao tratamento transuretral, enquanto a litotrispsia extracorpórea apresenta as menores de taxas de eliminação de cálculos e fica reservada aos pacientes de alto risco cirúrgico.

¿Qué factores inciden en la hospitalización de pacientes con cálculos ureterales menores a 10 mm en Fundación Santa Fe?

Debido a que el 12% de la población tendrá un cálculo en uréter hacia la mitad de su vida y a que las tasas de recurrencia en los que ya lo presentan son del 50% es necesario estudiar esta patología para aproximarse a un manejo adecuado en el servicio de urgencias. La literatura identifica un conjunto de factores que pueden contribuir a un cambio en el manejo médico. Objetivo: El objetivo de este estudio fue determinar los factores demográficos y clínicos asociados a manejo hospitalario en los pacientes con diagnóstico de cálculo ureteral menor de 10 mm. Métodos: Se diseñó un estudio de casos y controles no emparejados. Un caso fue definido como un paciente de 18 o más años con diagnóstico de urolitiasis con cálculo menor a 10 mm realizado por urotac que consultó (por primera vez para ese episodio) al servicio de urgencias de la Fundación Santa Fe de Bogotá entre el 1 de marzo de 2007 y 30 de abril de 2012. Se indagaron factores como edad, sexo, tamaño y localización del cálculo, respuesta a los analgésicos, evidencia de obstrucción e infección urinaria, además de otros antecedentes medicamentosos y clínicos. Se utilizó regresión logística no condicional bivariada y multivariada para evaluar la asociación entre tipo de manejo (hospitalario o ambulatorio) y las variables recolectadas, calculando odds ratio (OR) e intervalos de confianza al 95% (IC95%). Resultados: El riesgo de hospitalización se incrementó con: 1. La localización del cálculo en tercio superior o medio (OR=1.49; IC95%: 0.751-2.966) al comparar con el inferior, 2. El aumento del tamaño del cálculo (OR=1.49; IC95%: 0.751-2.966, por cada milímetro de incremento), y 3: Por la evidencia de obstrucción o infección urinaria y elevación de azoados. Por el contrario, hubo menos riesgo de hospitalización en aquellos pacientes con una respuesta analgésica apropiada en urgencias.

Estudo de cargas impulsivas com ênfase em explosões : estratégias de análises utilizando métodos numéricos

A avaliação da solicitação produzida por explosões, assim como da resposta de estruturas, são temas de muito interesse na engenharia atualmente, tanto pela quantidade crescente de acidentes relacionados com explosões quanto pelas ações terroristas muitas vezes vinculadas a estes tipos de ações. Neste contexto, o presente trabalho tem por objetivo explorar técnicas de análise tanto na modelagem da excitação quanto na resposta de estruturas consideradas como alvos. Para isto, são utilizadas metodologias de diferentes tipos: implementações baseadas em sistema de elementos finitos comerciais como Ansys [2000] e LS-Dyna [2003] e técnicas simplificativas que permitem realizar uma avaliação preliminar. As aplicações consideradas são indicadas a seguir: • Análise da Resposta de Estruturas Laminares Submetidas à Ação de Cargas Explosivas: determina-se a pressão produzida por explosivos sólidos a certa distância do epicentro, através de métodos simplificados, determinando a resposta esperada em placas retangulares; • Efeito da Pressão Interna em Vasos de Pressão (Extintores de Incêndio): comparando resultados numéricos e experimentais verifica-se a influência da pressão interna nas propriedades dinâmicas do sistema; • Estudo de Um Vaso Esférico de GLP Sob Ação de Uma Carga Explosiva: aplica-se a ação de uma onda explosiva produzida por um gás inflamável pesado sobre uma estrutura de vaso de pressão esférico com fluido e gás em seu interior, determinando sua resposta, avaliando também a influência de diferentes quantidades de líquido e pressão interna na resposta da estrutura; • Modelamento de uma Camada de Solo / Propagação das Ondas: verifica-se o comportamento da propagação de ondas em um meio elástico, comparando valores encontrados em testes experimentais. Analisa-se numericamente o efeito da inserção de uma valeta na atenuação de tais ondas; • Simulação Numérica Completa de uma Explosão: modela-se um semi-espaço submetido à ação de um explosivo sólido sobre sua superfície, avaliando os campos de pressão gerados. Ao final de cada aplicação são apresentadas conclusões parciais obtidas e as possibilidades de trabalhos futuros vislumbrados. Finalmente, conclui-se que as técnicas empregadas para as simulações são extremamente eficientes, considerando conhecidos todos os parâmetros envolvidos em cada modelo. Neste ponto é fundamental o trabalho do engenheiro, utilizando-se de seus conhecimentos técnicos para transformação do evento real em um modelo numérico, considerando e selecionando as simplificações necessárias.

Abordagem matemática de roll waves em escoamentos hiperconcentrados com superfície livre

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise de riscos envolvendo reações exotérmicas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A 2D BEM-FEM approach for time harmonic fluid-structure interaction analysis of thin elastic bodies.

[EN]This paper deals with two-dimensional time harmonic fluid-structure interaction problems when the fluid is at rest, and the elastic bodies have small thicknesses. A BEM-FEM numerical approach is used, where the BEM is applied to the fluid, and the structural FEM is applied to the thin elastic bodies.

Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.

Comportamiento de hormigones de alta resistencia con fibras de acero frente al impacto de proyectiles

Este trabajo presenta los resultados de la investigación llevada a cabo por los autores sobre el comportamiento de hormigón de 80 MPa de resistencia característica a compresión reforzado con diferentes cuantías de fibras de acero de alto contenido en carbono sometido al impacto de proyectiles de distintos calibres, determinando el espesor de muros de este tipo de hormigón que sería preciso disponer para impedir su perforación por dichos proyectiles, así como los valores máximos de penetración, para que en el caso de no producirse perforación y sólo penetración, no se genera cráter, “scabbing”, en el trasdós de los mismos. Previamente a los ensayos balísticos fue preciso diseñar los hormigones para que, presentaran determinadas características mecánicas, especialmente las relacionadas con la ductilidad, dado que estos hormigones especiales deben absorber la elevada energía que le transmiten los proyectiles y las ondas de choque que los acompañan.

Mejora y optimización de un código de simulación hidrodinámico con transporte de radiación para el estudio de sistemas en condiciones de alta densidad de energía

Se presentan las mejoras introducidas en un código de transporte de radiación acoplada a la hidrodinámica llamado ARWEN para el estudio de sistemas en el rango de física de alta densidad de energía (High Energy Density Physics). Los desarrollos introducidos se basan en las siguientes áreas: ít>,~ Ecuaciones de estado: se desarrolla una nueva metodología mediante la cual es posible ajustar los resultados de un modelo simple de ecuaciones de estado como QEOS a datos experimentales y resultados de AIMD. Esta metodología tiene carácter general para poder ser aplicada en gran cantidad de materuales de interés y amplia la flexibilidad de ajuste de los métodos de los que ha partido como base este trabajo. En segundo lugar, se ha desarrollado una librería para la gestión de tablas de datos de ecuaciones de estado que también incluye la gestión de tablas con datos de opacidades y de ionización. Esta nueva librería extiende las capacidades de la anterior al tener llamadas más específicas que aceleran los cálculos, y posibilidad de uso de varias tablas a la vez. Solver de difusión: se ha desarrollado un nuevo paquete para resolver la ecuación de difusión que se aplicará a la conducción de calor dentro del plasma. El método anterior no podía ser ejecutado en paralelo y producía resultados dependientes de la resolución de la malla, mientras que este método es paralelizable y además obtiene una solución con mejor convergencia, lo que supone una solución que no depende del refinamiento del mallado. Revisión del paquete de radiación: en primer lugar se ha realizado una revisión de la implementación del modelo de radiación descubriendo varios errores que han sido depurados. También se ha incluido la nueva librería de gestión de tablas de opacidades que permiten la obtención de las propiedades ópticas del plasma en multigrupos de energía. Por otra parte se ha extendido el cálculo de los coeficientes de transporte al esquema multimaterial que ha introducido David Portillo García en el paquete hidrodinámico del código de simulación. Por último se ha revisado el esquema de resolución del transporte y se ha modificado para hacerlo paralelizable. • Se ha implementado un paquete de trazado de rayos para deposición láser que extiende la utilidad del anterior al ser en 3D y poder utilizar entonces diferentes configuraciones. • Una vez realizadas todas estas tareas se ha aplicado el código ARWEN al estudio de la astrofísica de laboratorio simulando los experimentos llevados a cabo en la instalación PALS por Chantal Stehlé acerca de ondas de choque radiativas. Se han comparado los resultados experimentales frente a las predicciones del código ARWEN obteniéndose una gran concordancia en la velocidad de la onda de choque generada y en las dimensiones del precursor. El código de simulación sobre el que se ha trabajado, junto con los desarrollos aportados por otros investigadores durante la realización de esta tesis, ha permitido participar en colaboraciones con laboratorios de Francia o Japón y se han producido resultados científicos publicados basados en el trabajo descrito en esta tesis. ABSTRACT Improvements in radiation hydrodynamic code ARWEN for the study of systems in the range of physics high energy density (High Energy Density Physics) are presented. The developments introduced are based on the following áreas: • Equations of state: a new methodology was developed to adjust the results of a simple Equation of State model like QEOS to experimental data and results of AIMD. This methodology can be applied to a large amount of materials and it increases the flexibility and range of the previous methods used as basis for this work. Also a new computer library has been developed to manage data tables of thermodynamic properties as well as includes the management of opacity and ionization data tables. This new library extends the capabilities of the previous one with more specific routines, and the possibility of using múltiple tables for several materials. • Diffusion solver: a new package has been developed to solve the diffusion equation applied to the heat conduction of the plasma. The previous method is not parallelizable and it produced mesh dependent results, while this new package can be executed in parallel and achieves a more converged solution that does not depend on the refinement of the mesh. • Radiation package: the check of the radiation model rose several bugs in the implementation that had been removed. The new computer library for EOS managing includes capabilities to store opacity tables for multigroups of energy. Moreover the transport coefficients calculations have been extended for the new multimaterial hydrodynamic package developed by David Portillo García. Also the solving methodology for the transport equation has been modified to make the code run in parallel. • A new ray tracing package has been introduced to extend the previous one to 3D. Once all these tasks has been implemented, the ARWEN code has been applied to study laboratory astrophysics systems. Simulations have been done in order to reproduce the results of the experiments carried out in PALS facility by Chantal Stehlé in radiative shock production. Experimental results are in cióse agreement to the ARWEN estimations of the speed of the shock wave and the length of the precursor. The simulation code used in this thesis, including the work done in ARWEN by other colleagues at the time of this research, allowed the collaboration with other research institution in France and Japan and some of the results presented in this thesis have been published in scientific journals.

Análisis de sistemas de protección y seguridad de muros de vidrio frente a cargas explosivas

El objetivo principal de este proyecto es el estudio y análisis de dos sistemas de seguridad y protección para muros de vidrio en terminales de transporte, frente a ondas de choque ocasionadas por la detonación de cargas explosivas. Se ha ensayado en primer lugar un adhesivo de seguridad y protección ("Ultra S400") colocado en la cara expuesta a la detonación y un sistema de bandejas metálicas ("kalzip"), tanto en el exterior como en el interior, en el segundo ensayo. Los ensayos se han instrumentado con captadores de presión, acelerómetros, sismógrafos y con cámaras convencionales y de alta velocidad, ésta última nos ha permitido calcular la velocidad de choque en el primero de los ensayos. El análisis de las presiones ha permitido caracterizar la explosión y compararla con los valores de referencia para el TNT del manual norteamericano UFC 3-340. El diseño de los ensayos se ha hecho conforme a las normas ISO y GSA de modo que los resultados han permitido clasificar las diferentes configuraciones de muro cortina ensayadas. ABSTRACT The main objective of this thesis is the analysis and evaluation of two safety and protective systems for glass walls at transportation terminals, against impact waves caused by the detonation of an explosive charge. It has been tested in first place, a safety and protective adhesive ("Ultra S400"), located at the exposed side of the explosion, and a system consisting of some metallic trays ("kalzip"), both inside and outside, for the second experiment. The tests have been monitored by pressure captures, accelerometer, seismographs, conventional cameras, and high-speed cameras, which allow us to determine the impact speed in the first experiment. The analysis of the pressure enables the characterization of the explosion and its comparison against typical TNT reference values from the UFC 3-340 North American manual. The design of the tests has been made according to ISO and GSA standards. Therefore, the results allow us to classify different patterns of curtain walls tested.