No description available.

No description available.

No description available.

Knowledge of the indoor propagation of sound has many important applications including acoustic prediction in homes, office buildings, stores, and schools, and the design of concert halls, auditoriums, classrooms, and factories. At low frequencies, interior acoustics are analyzed with the wave equation, but significant computational expense imposes an upper frequency limit. Thus, energy methods are often sought for high frequency analysis. However, conventional energy methods are significantly limited by vast simplifications or computational costs. Therefore, new improvements are still being sought. The basis of this dissertation is a recently developed mathematical model for interior acoustics known as the acoustic diffusion model. The model extends statistical methods in high frequency acoustics to predict the spatial distribution of acoustic energy in the volume over time as a diffusion process. Previously, solutions to the acoustic diffusion model have been limited to one dimensional (1-D) analytical solutions and to the use of the finite element method (FEM). This dissertation focuses on a new, efficient method for solving the acoustic diffusion model based on a boundary element method (BEM) solution using the Laplace transform. First, a Laplace domain solution to the diffusion model is obtained using the BEM. Then, a numerical inverse Laplace transform is used to efficiently compute the time domain response. The diffusion boundary element-Laplace transform solution (BE-LTS) is validated through comparisons with Sabine theory, ray tracing, and a diffusion FEM solution. All methods demonstrate excellent agreement for three increasingly complex acoustic volumes and the computational efficiency of the BE-LTS is exposed. Structural coupling is then incorporated in the diffusion BE-LTS using two methods. First, a simple transmission coefficient separating two acoustic volumes is implemented. Second, a structural power flow model represents the coupling partition separating acoustic volumes. The validation of these methods is successfully performed in an example through comparisons with statistical theory, a diffusion FEM solution, ray tracing, and experimental data. Finally, the diffusion model and the BE-LTS are shown to possess capabilities beyond that of room acoustics. The acoustic transmission through a heat exchanger, acoustic foam, and mufflers is successfully modeled using the diffusion BE-LTS and compared to experimental data. / Ph. D.

[EN] The behavior of the heart is governed by electrical currents generated in the myocardium, and therefore, the study of the cardiac electrical activity is essential for the diagnosis of cardiac diseases. The forward problem of the electrocardiography (FP) entails the calculation of the torso potentials from the electrical activity of the heart and the 3D body model, while the inverse problem (IP) resolution allows the noninvasive reconstruction of the electrical activity of the heart from surface potentials. The IP is of great importance in clinical applications since it allows estimating the electrical activity of the myocardium with only noninvasive recordings. However, IP resolution is still a big challenge in electrocardiography since it is ill-posed, very unstable and has multiple solutions. In this thesis different algorithms and strategies based on the IP resolution were developed and applied in the noninvasive diagnosis of ventricular and atrial arrhythmias and evaluated with mathematical cellular models and clinical data bases. The thesis focuses on the IP resolution for the noninvasive reconstruction of the myocardial electrical activity for different diseases and propagation patterns, implementing a novel system for complex propagation patterns. The obtained results and propagation patterns were evaluated and classified with the corresponding optimal resolution strategy that minimizes the error and increases the stability of the system, proving its advantages and disadvantages depending on the different diseases and their activation pattern. A novel iterative method was implemented for the IP dipolar resolution optimized for representing simple propagation patterns, achieving a high stability and robustness against noise by constraining the solution to a limited number of dipoles. However, propagation patterns not representable by few dipoles need to be computed with the IP in terms of epicardial solutions which provide a more detailed estimation of the myocardial activity. IP resolution in the voltage and phase domains showed a good accuracy for simple and organized propagation patterns. This method allowed the noninvasive diagnosis of the Brugada syndrome or the location of ectopic focus in atrial arrhythmias by performing a parametric analysis of the electrograms morphology or the activation map reconstruction. However, mathematical and patient results presented in this thesis proved that, for complex propagation patterns like atrial fibrillation (AF), inverse solutions in the voltage and phase domains are over-smoothed and over-optimistic, simplifying the complex AF activity, leading to non-physiological results that do not match with the complex intracardiac electrograms recorded in AF patients. In this thesis, we proposed a novel technique for the noninvasive identification and location of high dominant frequency AF sources, based on the assumption that in many cases atrial drivers present the highest activation rate with an intermittent propagation to the rest of the tissue that activates at a slower rate. Although, voltage and phase inverse solutions for AF complex propagation patterns were over smoothed and inaccurate, the noninvasive estimation of frequency maps was significantly more accurate, allowing the identification of the AF frequency gradient and location of high frequency sources. This technique may help in planning ablation procedures, avoiding unnecessary interseptal punctures for right-to-left frequency gradients cases and facilitating the targeting of the AF drivers, reducing risk and time of the clinical procedure. / [ES] El comportamiento del corazón se rige por corrientes eléctricas generadas en el miocardio y, por lo tanto, el estudio de su actividad eléctrica es esencial para el diagnóstico de enfermedades cardíacas. El problema directo (PD) de la electrocardiografía implica el cálculo de los potenciales del torso a partir de la actividad eléctrica del corazón y el modelo 3D del cuerpo, mientras que la resolución del problema inverso (PI) permite la reconstrucción no invasiva de la actividad eléctrica del corazón a partir de los potenciales de superficie, cobrando una gran importancia en la práctica clínica. Sin embargo, sigue siendo un gran desafío para la electrocardiografía ya que está mal planteado, es muy inestable y tiene múltiples soluciones. A lo largo de esta tesis se han desarrollado diferentes estrategias para la resolución del PI, aplicándolas en el diagnóstico no invasivo de arritmias ventriculares y auriculares, verificándolas mediante modelos celulares matemáticos y bases de datos clínicas. La tesis se centra en la resolución del PI para la reconstrucción no invasiva de la actividad eléctrica del miocardio para diferentes enfermedades cardiacas con diferentes patrones de propagación, implementando un novedoso sistema para patrones de propagación complejos. Además, se han validado los resultados obtenidos y se han clasificado los diferentes patrones de propagación con la estrategia de resolución del PI óptima que minimice el error y aumente la estabilidad del sistema. Un nuevo método iterativo fue implementado para la resolución del PI para fuentes dipolares, siendo óptimo para representar patrones de propagación simples, logrando una alta estabilidad e inmunidad al ruido al restringir la solución a un número limitado de dipolos. Sin embargo, los patrones de propagación que no pueden ser representados por un número limitado de dipolos deben calcularse mediante la resolución del PI en términos de potenciales epicárdicos, proporcionando una estimación más detallada de la actividad del miocardio. La resolución del PI en el dominio de la tensión y fase mostró ser muy preciso para patrones de propagación simples y organizados. Este método permite el diagnóstico no invasivo del síndrome de Brugada o la ubicación de focos ectópicos en arritmias auriculares mediante un análisis paramétrico de la morfología de los electrogramas o la reconstrucción de los mapas de activación. Sin embargo, los resultados matemáticos y clínicos presentados en esta tesis demostraron que, para patrones de propagación complejos como la fibrilación auricular (FA), los resultados obtenidos mediante la resolución del PI en el dominio de la tensión y fase son demasiado suaves y optimistas, simplificando enormemente la complejidad de la FA, llevando a resultados no fisiológicos que no coinciden con la actividad compleja de los electrogramas intracardiacos registrados en pacientes con FA. En esta tesis, se ha propuesto una novedosa técnica para la identificación y localización no invasiva de fuentes con una frecuencia dominante alta, basado en la suposición de que en muchos casos las fuentes eléctricas que generan y mantienen la FA presentan una tasa de activación más alta, con una propagación intermitente hacia el resto del tejido auricular cuya frecuencia de activación es más lenta. Aunque las soluciones en el dominio de la tensión y fase para patrones de propagación complejos fueron más suaves y menos precisas, la estimación no invasiva de los mapas de frecuencia fue significativamente más precisa, permitiendo la identificación del gradiente de frecuencia y ubicación de las fuentes de FA de alta frecuencia. Esta técnica puede ser de gran ayuda en la planificación de los procedimientos de ablación, evitando punciones interseptales innecesarias para casos con un gradiente de frecuencia de derecha a izquierda y facilitando la localización de las fuentes de alta frecuencia / [CA] El comportament del cor es regeix per corrents elèctrics generades en el miocardi i, per tant, l'estudi de la seua activitat elèctrica és essencial per al diagnòstic de malalties cardíaques. El problema directe (PD) de l'electrocardiografia implica el càlcul dels potencials del tors a partir de l'activitat elèctrica del cor i el model 3D del cos, mentre que la resolució del problema invers (PI) permet la reconstrucció no invasiva de l'activitat elèctrica del cor a partir de els potencials de superfície. La resolució del PI de l'electrocardiografia té una gran importància en la pràctica clínica atès que fa possible una estimació de l'activitat elèctrica del miocardi únicament a partir de registres no invasius. No obstant això, la resolució del PI segueix sent un gran desafiament per a la electrocardiografia ja que està mal plantejat, és molt inestable i té múltiples solucions. Al llarg d'aquesta tesi s'han desenvolupat diferents estratègies basades en la resolució PI, aplicant-les en el diagnòstic no invasiu d'arítmies ventriculars i auriculars, verificant mitjançant models cel·lulars matemàtics i bases de dades clíniques. La tesi se centra en la resolució del PI per a la reconstrucció no invasiva de l'activitat elèctrica del miocardi per a diferents malalties cardíaques amb diferents patrons de propagació, implementant un nou sistema per a patrons de propagació complexos. A més se han validat els resultats obtinguts i se han classificat els diferents patrons de propagació amb l'estratègia de resolució del PI òptima que minimitze l'error i augmente l'estabilitat del sistema. Un nou mètode iteratiu va ser implementat per a la resolució del PI per fonts dipolars, sent òptim per representar patrons de propagació simples, aconseguint una alta estabilitat i immunitat al soroll en restringir la solució a un nombre limitat de dipols. No obstant això, els patrons de propagació que no poden ser representats per un nombre limitat de dipols s'han de calcular mitjançant la resolució del PI en termes de potencials epicàrdics, proporcionant una estimació més detallada de l'activitat del miocardi. La resolució del PI en el domini de la tensió i fase va mostrar ser molt precís per a patrons de propagació simples i organitzats. Aquest mètode permet el diagnòstic no invasiu de la síndrome de Brugada o la ubicació de focus ectòpics en arítmies auriculars mitjançant una anàlisi paramètric de la morfologia dels electrogrames o la reconstrucció dels mapes d'activació. No obstant això, els resultats matemàtics i clínics presentats en aquesta tesi van demostrar que, per patrons de propagació complexos com la fibril·lació auricular (FA), els resultats obtinguts mitjançant la resolució del PI en el domini de la tensió i fase són massa suaus i optimistes, simplificant enormement la complexitat de la FA, obtenint resultats no fisiològics que no coincideixen amb l'activitat complexa dels electrogrames intracardiacos registrats en pacients amb FA. En aquesta tesi, s'ha proposat una nova tècnica per a la identificació i localització no invasiva de fonts amb una freqüència dominant alta, basat en la suposició que en molts casos les fonts elèctriques que generen i mantenen la FA presenten una taxa d'activació més alta, amb una propagació intermitent cap a la resta del teixit auricular on la freqüència d'activació és més lenta. Encara que, les solucions en el domini de la tensió i fase per patrons de propagació complexos van ser més suaus i menys precises, l'estimació no invasiva dels mapes de freqüència va ser significativament més precisa, permetent la identificació del gradient de freqüència i ubicació de les fonts de FA d'alta freqüència. Aquesta tècnica pot ser de gran ajuda en la planificació dels procediments d'ablació, evitant puncions interseptales innecessaris per a casos amb un gradient de freqüència de dreta a esquerra i facilitant la / Pedrón Torrecilla, J. (2015). Non-invasive Reconstruction of the Myocardial Electrical Activity from Body Surface Potential Recordings [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/58268

In this thesis the impact of inverted cones on a liquid surface is studied. It is known that with the right combination of velocity, geometry, and surface treatment, a cavity of air can be formed behind an impacting body and extended for a considerable distance. Other investigators have shown that the time and depth of the cone when this cavity collapses and seals follows a different power law for flat objects such as disks, then it does for slender objects such as cylinders. Intuitively it can be expected that a more slender body will have less drag and that the streamlined shape will not push the fluid out of it's way at impact to the same extent as a more blunt body, therefore forming a smaller cavity behind it. With a smaller initial cavity, the time and depth of it's eventual collapse can be expected to be less than that of a much more blunt object, such as a flat disk. To study this, a numerical model has been developed to simulate cones with the same base radius but different angles impacting on a liquid surface over a range of velocities, showing how the seal depth, time at cavity seal, and drag forces change. In order to ensure the numerical model is accurate, it is compared with experimental data including high speed video and measurements made of the force with time. It is expected that the results will fall inside the power law exponents reported by other authors for very blunt objects such as disks on one end of the spectrum, and long slender cylinders on the other. Furthermore, we expect that the drag force exerted on the cones will become lower as the L/D of the cone is increased. / Master of Science / In this thesis the impact of inverted cones on a liquid surface is studied. It is known that with the right combination of velocity, geometry, and surface treatment, a cavity of air can be formed behind an impacting body and extended for a considerable distance. Other investigators have shown that the time and depth of the cone when this cavity collapses and seals follows a different power law for flat objects such as disks, then it does for slender objects such as cylinders. Intuitively it can be expected that a more slender body will have less drag and that the streamlined shape will not push the fluid out of it’s way at impact to the same extent as a more blunt body, therefore forming a smaller cavity behind it. With a smaller initial cavity, the time and depth of it’s eventual collapse can be expected to be less than that of a much more blunt object, such as a flat disk. To study this, a numerical model has been developed to simulate cones with the same base radius but different angles impacting on a liquid surface over a range of velocities, showing how the seal depth, time at cavity seal, and drag forces change. In order to ensure the numerical model is accurate, it is compared with experimental data including high speed video and measurements made of the force with time. It is expected that the results will fall inside the power law exponents reported by other authors for very blunt objects such as disks on one end of the spectrum, and long slender cylinders on the other. Furthermore, we expect that the drag force exerted on the cones will become lower as the <i>L/D</i> of the cone is increased.

Au cours de cette thèse, nous avons développé un outil numérique, basé sur une formulation intégrale en éléments de frontière, qui permet une analyse classique du comportement d'une fissure 3D soumise à des sollicitations mécaniques complexes. Cet outil industriel est destiné à être intégré dans un code de calcul à usage industriel. Dans le but d'appréhender l'impact de la microstructure sur le comportement en fissuration fragile, nous nous sommes intéressés aux modèles de comportement non local. Nous avons commencé par adopter le modèle de comportement élastique non local de Eringen, qui permet de décrire plus finement le comportement élastique au voisinage de la fissure en prenant en compte les interactions à longue distance au sein du matériau. Cette modélisation du comportement conduit, contrairement à l'approche classique, à un un champ de contrainte fini sur le front de la fissure et localement maximal en avant du front. Ces résultats montrent qu'il est possible de prévoir la stabilité et la direction de propagation de la fissure à l'aide d'un critère plus simple et plus naturel, basé sur les variations du champ de contrainte au voisinage du front de la fissure. La stratégie numérique adoptée permet de traiter indifféremment des cas de fissure en traction, compression, cisaillement ou sollicitation mixte. L'intérêt de l'approche non-locale étant clairement démontré, nous avons considéré la version améliorée du modèle de Eringen telle que proposée par Polizzotto. Cette modélisation est la plus appropriée pour les milieux finis et requiert une mise en oeuvre numérique particulière. Les bases d'une méthodologie numérique, initiée par R. Kouitat ont été établies. Cette méthode est fondée sur un couplage des éléments de frontière avec une méthode de collocation par points d'équations aux dérivées partielles. Les premiers résultats obtenus dans ce cadre sont très encourageants et montrent qu'il sera effectivement possible de traiter le phénomène irréversible de fissuration de la même façon que les problèmes de plasticité / In this thesis, we have developed a numerical tool, based on a classical boundary elements method, which allows a conventional analysis of a stationary crack in a 3D specimen under complex mechanical loading. In order to assess the impact of the microstructure on the brittle fracture, we were interested in non local models of behavior. First, we have adopted the non local elastic model due to Eringen. This refined constitutive equation allows to account for long range interactions in the description of the elastic behavior in the vicinity of the crack front. Unlike the traditional approach, this type of model leads to a finite stress field at the crack front. Moreover, the stress is locally maximal ahead of the front. These interesting results indicate that it is possible to predict the stability and direction of crack propagation in a simple and more naturel way by using stress based criteria. The implemented numerical strategy can handle cases of crack in tension or compression, under shear stress or mixed loadings. Having clearly highlighted the interest of non local models, we have adopted the improved version of Eringen elastic model as proposed by Polizzotto. This elastic model is applicable to finite domains and requires a specific numerical treatment. The basis of such a numerical strategy initiated by R. Kouitat has been established. The method couples the conventional boundary element method with local point interpolation of a strong form differential equation. Promising results are obtained and show that with such modeling of material behavior, it is possible to describe the irreversible process of fracturing in a similar way as plasticity

Esse trabalho trata da análise de corpos multi-fraturados utilizando o método dos elementos de contorno. Este método numérico é conhecido por ser robusto e preciso neste tipo de problema e também por requerer pequeno esforço computacional na criação da malha de elementos para o crescimento das fissuras. Duas metodologias para a análise do comportamento das fissuras são consideradas. A primeira é a já consagrada metodologia dual. Por meio desta técnica equações integrais distintas são aplicadas às faces da fissura. Estas equações integrais são escritas em termos de deslocamentos e forças de superfície. A segunda metodologia é a que emprega a formulação singular onde a fissura é considerada como um vazio no domínio sendo as faces da fissura separadas por uma pequena distância. No tocante ao crescimento das fissuras foi desenvolvido um procedimento especial para a determinação da direção de crescimento das fissuras o qual mostrou-se muito eficiente levando a resultados precisos. O crescimento das fissuras é efetuado considerando o fator de intensidade de tensão atuante na extremidade de cada fissura. Dessa forma, as fissuras mais solicitadas apresentam maior comprimento de propagação tornando a análise mais realista. Os fatores de intensidade de tensão são calculados por meio de duas técnicas. A primeira é a já conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Uma técnica alternativa é também utilizada a qual emprega o campo de tensões presente na extremidade da fissura. Após a determinação dos fatores de intensidade de tensão quatro diferentes teorias de interação de modos podem ser utilizadas para a determinação do ângulo de propagação. Foram analisadas estruturas sendo os resultados comparados aos previstos analiticamente e também numericamente. As respostas obtidas foram satisfatórias validando assim a metodologia proposta neste trabalho / This work deals with analysis of multi-fractured bodies using boundary element method. This numerical method is known to be robust and accurate in this kind of problem and by small computational effort to create elements mesh of crack growth. Two methodologies to analyze of crack behavior are considerate. The first is consecrated dual methodology. Through this technique different kind of integral equations are applied to crack boundaries. These integrals equations are written in displacements and traction variables. Second methodology is singular formulation. Through this technique crack is represented like a hole in body and the crack boundaries is separated by a small gap. For crack growth was created a special proceeding to determination crack growth direction. This method is very efficient and your results are accurate. Crack growth is made through the stress intensity factor performed in crack tip. Then the cracks more requested going to propagate with a larger length’s growth turning this model very realistic. The stress intensity factors are calculated through two techniques in this work. First is the known correlation displacement technique which related displacement in crack boundaries. An alternative technique is also used which consider stress field in crack tip. After determination of stress intensity factors four different theories are used to calculate the crack growth angle. In this work were analyzed structures with results are compared with analytical and numerical answers. The results obtained went very satisfactory validating the methodology proposed

A otimização topológica de estruturas está relacionada à concepção de projetos que executem suas funções com nível de segurança adequado empregando a quantidade mínima de material. Neste trabalho, determina-se a geometria ótima de estruturas planas por meio do acoplamento do Método dos Elementos de Contorno (MEC) ao Método Level Set (MLS). O algoritmo é composto por 3 etapas: problema mecânico, otimização topológica e reconstrução da estrutura. O problema mecânico é resolvido pelas equações algébricas do MEC. A otimização topológica é determinada pelo MLS, este representa a geometria do corpo e suas evoluções por meio da função Level Set (LS) avaliada em seu nível zero. Na reconstrução realiza-se o remalhamento, pois a cada iteração a estrutura é modificada. O acoplamento proposto resulta na geometria ótima da estrutura sem a necessidade da aplicação de filtros. Os exemplos analisados mostram que algoritmo desenvolvido capta adequadamente a geometria ótima das estruturas. Com esse trabalho, avança-se no campo das aplicações do acoplamento MEC-MLS e no desenvolvimento de soluções inovadoras para problemas complexos de engenharia. / In general, the topological optimization of structures is related to design projects that perform their functions with appropriate security levels using the minimum amount of material. This research determines the optimal geometry of 2D structures by coupling the Boundary Blement Method (BEM) to Level Set Method (LSM). The algorithm consists of 3 steps: mechanical model, topology optimization and structure reconstruction. The mechanical model is solved by BEM algebraic equations. The topology optimization is determined using the MLS, the geometry of the body is determined by the Level Set (LS) function evaluated at the zero level. The reconstruction achieves the remeshing, because for each iteration of the structure is modified. The proposed coupling results in the optimal geometry of the structure without the filters application. The examples show that the algorithm developed captures adequately the optimal geometry of the structures. With this dissertation, it is possible advance in the field of applications of the BEM - LSM and develop innovative solutions to complex engineering problems.

O presente trabalho trata da análise do comportamento mecânico da madeira, material anisotrópico, submetida ao crescimento de fissuras e aos processos mecânicos que causam sua ruptura, utilizando os conceitos de Mecânica da Fratura. A avaliação do colapso de peças de madeira, pela da mecânica da fratura, foi realizada por meio de análises experimentais. Para a determinação das propriedades de fratura da madeira foram determinadas propriedades mecânicas da madeira, tais com: módulo de ruptura a flexão estática, módulo de elasticidade à flexão estática, resistência à tração paralela às fibras, resistência ao cisalhamento, resistência ao fendilhamento paralela às fibras e resistência à tração normal às fibras. Foram efetuados ainda ensaios experimentais referentes à mecânica da fratura em corpos de prova tipo SENB (Single Edge Notched Beam). A modelagem numérica foi realizada utilizando-se um código computacional baseado no método de elementos de contorno. Para as simulações foram utilizados os valores das propriedades físicas, de resistência e de elasticidade obtidas na caracterização experimental dos elementos de madeira. Os resultados experimentais foram comparados com as análises numéricas para determinação de uma formulação específica para situações de colapso de elementos de madeira. / The present work deals with the analysis of the mechanical behavior of wood, anisotropic material, submitted to the growth of cracks and the mechanical processes that cause its rupture, using the concepts of Fracture Mechanics. The evaluation of the collapse of pieces of wood, by the fracture mechanics, was carried out by means of experimental analyzes. To determine the fracture properties of the wood, mechanical properties of the wood were determined, such as: modulus of static bending, modulus of elasticity at static bending, tension resistance parallel to grain, shearing resistance, cleavage strength parallel to grain and tension resistance normal to grain. Experimental tests were performed on fracture mechanic in Single Edge Notched beams (SENB). Numerical modeling was performed using a computational code based on the boundary element method. For the simulations, the values of the physical properties, strength and elasticity obtained in the experimental characterization of the wood elements were used. The experimental results were compared with the numerical analyzes to determine a specific formulation for situations of collapse of wood elements.