Kronecker structured representations are used to cope with the state space explosion problem in Markovian modeling and analysis. Currently an open research problem is that of devising strong preconditioners to be used with projection methods for the computation of the stationary vector of Markov chains (MCs) underlying such representations. This paper proposes a block SOR (BSOR) preconditioner for hierarchical Markovian Models (HMMs) that are composed of multiple low level models and a high level model that defines the interaction among low level models. The Kronecker structure of an HMM yields nested block partitionings in its underlying continuous-time MC which may be used in the BSOR preconditioner. The computation of the BSOR preconditioned residual in each iteration of a preconditioned projection method becoms the problem of solving multiple nonsingular linear systems whose coefficient matrices are the diagonal blocks of the chosen partitioning. The proposed BSOR preconditioner solvers these systems using sparse LU or real Schur factors of diagonal blocks. The fill-in of sparse LU factorized diagonal blocks is reduced using the column approximate minimum degree algorithm (COLAMD). A set of numerical experiments are presented to show the merits of the proposed BSOR preconditioner.

Cette thèse présente les résultats d'observations d'étoiles en rotation rapide menées sur le spectro-interféromètre AMBER du VLTI dans ses modes haute et moyenne résolutions spectrales. Les mesures effectuées sont les visibilités estimées sur trois bases simultanées, les phases différentielles en fonction de la longueur d'onde et des phases de clôtures avec, pour certaines nuits une bonne couverture du plan (u,v). Les données utilisées sont issues de plusieurs campagnes d'observation. Ces dernières étaient fortement dégradées par les défauts optiques d'AMBER, et affectés par des bruits classiques d'interférométrie optique à longue base en IR: défauts du détecteur, bruit de lecture, instabilités du suiveur de franges, ...etc. Leur analyse a nécessité la mise au point d'outils numériques de réduction spécifiques pour atteindre les précisions nécessaires à l'interprétation de mesures interférométriques. Pour interpréter ces mesures j'ai développé un modèle semi-analytique chromatique d'étoile en rotation rapide qui m'a permis d'estimer, à partir des phases différentielles; le degré d'aplatissement, le rayon équatorial, la vitesse de rotation, l'angle d'inclinaison, l'angle position de l'axe de rotation de l'étoile sur le ciel, la distribution de la température effective locale et de la gravité à la surface de l'étoile dans le cadre du théorème de von Zeipel. Les résultats concernant 4 étoiles massives de types spectraux B, A et F m'ont permis de les caractériser pour les mécanismes évoqués ci-dessus et d'ouvrir ainsi la perspective d'études plus systématiques d'objets similaires en étendant ultérieurement ces études à la relation photosphère-enveloppe circumstellaire. / This thesis presents the results of rapidly rotating stars observations conducted on the AMBER spectro-interferometer VLTI in its high average spectral modes and resolutions. The measurements are estimated on three simultaneous visibility bases, differential phases depending on the wavelength and closure phases, with good coverage of the (u, v) plane for some nights. The data used are from several observation campaigns. These were highly degraded by the optical defects of AMBER, and assigned by standard optical interferometry long base IR noises: defects of the detector, reading noise, fringes follower instabilities, ... etc. Their analysis required the development of digital reduction of specific tools to achieve the necessary details to the interpretation of interferometric measurements. In interpreting these measures I developed a chromatic semi-analytical model of rapidly rotating star that allowed me to estimate, from the differential phases; the degree of flattening, the equatorial radius, speed of rotation, angle of inclination, the position angle of the star rotation axis in the sky, the local distribution of the actual temperature and the gravity to the surface of the star within the von Zeipel theorem. The results for four massive stars of spectral type B, A and F have allowed me to characterize the mechanisms discussed above and thus open framework for more systematic studies of similar objects subsequently extending these studies to the relationship photosphere circumstellar envelope.

[EN] This Thesis is focused on the development and implementation of efficient numerical methods for the acoustic modelling and design of noise control devices in the exhaust system of combustion engines. Special attention is paid to automotive perforated dissipative silencers, in which significant differences are likely to appear in their acoustic behaviour, depending on the temperature variations within the absorbent material. Also, material heterogeneities can alter the silencer attenuation performance. Therefore, numerical techniques considering all these features are required to guarantee the accuracy of the results. A literature review is carried out, mainly related to one-dimensional models, as well as to acoustic models for absorbent materials and perforated surfaces. However, plane wave model limitations make indispensable using alternative multidimensional methods. In addition, the possibility of using new acoustic elements is explored. These elements have as an objective being a potential alternative to the fibrous absorbent materials, which can have a negative impact on health. The Thesis considers the use of microperforated and sintered surfaces. The latter have, in some cases, a nearly constant acoustic impedance, whose value depends, among others, on the thickness and porosity of the plates. To avoid the limitations of plane wave models, a finite element (FE) approach is proposed for the acoustic analysis of dissipative silencers including a perforated duct with uniform axial mean flow and an outer chamber with a heterogeneous distribution of the absorbent material. On the other hand, property variations can be also produced by temperature gradients. In this case, a hybrid FE model has been derived for perforated dissipative silencers including: (1) Thermal gradients in the central duct and the chamber; (2) A perforated passage carrying non-uniform axial mean flow. A FE approach has been implemented to solve the pressure-based wave equation for a non-moving heterogeneous medium, associated with the chamber. Also, the governing equation in the central duct has been written and solved in terms of an acoustic velocity potential to allow the presence of an axially inhomogeneous flow. The coupling between both regions has been carried out by means of a perforated duct and its acoustic impedance, adapted here to include absorbent material heterogeneities and mean flow effects. It has been found that the presence of non-homogeneities can have a significant influence on the acoustic attenuation of a silencer and should be included in the theoretical models. Optimization techniques for industrial noise control devices are relevant, since they lead to the production of elements with better characteristics. Evolutionary algorithms are emergent techniques able to obtain a solution, even in those problems in which the traditional optimization have difficulties. Optimization techniques are combined with the FE method to achieve the maximum attenuation in the frequency range of interest. A multichamber silencer optimization problem is defined and several analyses are carried out to obtain the most suitable configuration for each application. Under certain assumptions of axial uniformity, several techniques have been considered to reduce the computational effort of a full 3D FE analysis for dissipative silencers with temperature gradients and mean flow. These are based on a decomposition of the acoustic field into transversal and axial modes within each silencer subdomain, and a matching procedure of the modal expansions at the silencer area changes through the continuity conditions of the acoustic fields. The relative computational efficiency and accuracy of predictions for the matching techniques are studied, including point collocation at nodes and Gauss points and also mode-matching with weighted integration. All of them provide accurate predictions of the attenuation and improve the computational cost of a FE calculation / [ES] Esta Tesis se centra en el desarrollo e implementación de métodos numéricos eficientes para el diseño y modelado de componentes de la línea de escape en motores de combustión interna. Merecen especial atención los silenciadores disipativos perforados de automóviles, ya que su comportamiento acústico puede sufrir variaciones importantes debidas a las variaciones de temperatura en el material absorbente, así como a las heterogeneidades de la fibra. Por tanto, se requieren técnicas numéricas que consideren estos casos para garantizar la precisión de los resultados. Se lleva a cabo una revisión bibliográfica que recoge los modelos de onda unidimensionales, así como modelos acústicos de materiales absorbentes y superficies perforadas. Sin embargo, las limitaciones de los primeros hacen indispensable el uso de modelos multidimensionales. Además se explora la posibilidad de usar nuevos elementos acústicos, cuyo objetivo es ser una alternativa potencial a los materiales absorbentes, que pueden tener un efecto negativo sobre la salud. La Tesis considera el uso de superficies microperforadas y sinterizadas. Estas últimas en algunos casos presentan una impedancia casi constante, cuyo valor depende, entre otras cosas, del espesor y la porosidad de las placas. Para evitar las limitaciones de los modelos de onda plana, se propone un enfoque en elementos finitos (EF) para el análisis acústico de silenciadores disipativos que incluyen un conducto con flujo medio axial uniforme y una cámara externa con una distribución heterogénea de material absorbente. Por otro lado, la variación de las propiedades también puede producirse por gradientes térmicos. En este caso, se propone una formulación híbrida de EF para silenciadores disipativos perforados que incluye: (1) Gradientes térmicos en el conducto central y la cámara; (2) Un conducto perforado que canaliza flujo medio axial no uniforme. Se ha implementado una formulación de EF para resolver la ecuación de ondas en términos de presión para el medio estacionario heterogéneo asociado a la cámara. Además, la ecuación asociada al conducto central, expresada en términos de potencial de velocidad acústica, permite la presencia de flujo axial no uniforme. El acoplamiento entre ambas regiones se ha realizado mediante un conducto perforado y su impedancia acústica y se ha adaptado para incluir la citada falta de homogeneidad. Se ha visto que las heterogeneidades pueden influir notablemente en la atenuación acústica de un silenciador, debiéndose incluir en los modelos teóricos. Las técnicas de optimización para componentes industriales de control de ruido son importantes, ya que producen elementos con mejores características. Los algoritmos evolutivos son técnicas emergentes capaces de obtener una solución, incluso cuando la optimización tradicional tiene dificultades. Las técnicas de optimización se combinan con el MEF para conseguir la máxima atenuación posible en el rango de frecuencias de interés. Se ha definido un problema de optimización de un silenciador multicámara y se han llevado a cabo varios análisis para obtener la configuración más adecuada para cada caso. Bajo ciertas hipótesis de uniformidad axial, se han considerado varias técnicas para reducir el coste computacional de un análisis 3D completo para silenciadores disipativos con gradientes de temperatura y flujo medio. Éstas se basan en la descomposición del campo acústico en modos axiales y transversales dentro de cada subdominio, y un procedimiento de acoplamiento de las expansiones modales en los cambios de sección del silenciador mediante las condiciones de continuidad de los campos acústicos. Se estudia la eficiencia computacional y precisión de las predicciones de las técnicas de acoplamiento, incluyendo colocación puntual en nodos y puntos de Gauss, así como ajuste modal. Todos ellos proporcionan predicciones precisas de la atenuación mejorando el coste / [CAT] Aquesta Tesi es centra en el desenvolupament i implementació de mètodes numèrics eficients per al disseny i modelatge de components de la línia d'escapament en motors de combustió interna. Mereixen especial atenció els silenciadors dissipatius perforats d'automòbils, ja que el seu comportament acústic pot patir variacions importants degudes a les variacions de temperatura en el material absorbent, així com a les heterogeneïtats de la fibra. Per tant, es requereixen tècniques numèriques que considerin aquests casos per garantir la precisió dels resultats. Es porta a terme una revisió bibliogràfica que recull els models d'ona unidimensionals, així com models acústics de materials absorbents i superfícies perforades. No obstant això, les limitacions dels primers fan indispensable l'ús de models multidimensionals. A més s'explora la possibilitat d'usar nous elements acústics amb l'objectiu que siguen una alternativa potencial als materials absorbents, que poden tenir un efecte negatiu sobre la salut. La Tesi considera l'ús de superfícies microperforades i sinteritzades. Aquestes últimes en alguns casos presenten una impedància gairebé constant. El seu valor depèn, entre altres coses, del gruix i la porositat de les plaques. Per evitar les limitacions dels models d'ona plana, es proposa un enfocament amb elements finits (EF) per a l'anàlisi acústic de silenciadors dissipatius que inclouen un conducte amb flux mig axial uniforme i una càmera externa amb una distribució heterogènia de material absorbent. D'altra banda, la variació de les propietats també es pot produir per gradients tèrmics. En aquest cas, es proposa una formulació híbrida d'EF per silenciadors dissipatius perforats que inclou: (1) Gradients tèrmics en el conducte central i la càmera; (2) Un conducte perforat que canalitza flux mig axial no uniforme. S'ha implementat una formulació d'EF per resoldre l'equació d'ones en termes de pressió per al medi estacionari heterogeni associat a la càmera. A més, l'equació associada al conducte central, expressada en termes de potencial de velocitat acústica, permet la presència de flux axial no uniforme. L'acoblament entre les dues regions s'ha realitzat mitjançant un conducte perforat i la seva impedància acústica i s'ha adaptat per incloure la esmentada falta d'homogeneïtat. S'ha vist que les heterogeneïtats poden influir notablement en l'atenuació acústica d'un silenciador i s'han d'incloure en els models teòrics. Les tècniques d'optimització per a components industrials de control de soroll són importants, ja que produeixen elements amb millors característiques. Els algoritmes evolutius són tècniques emergents capaces d'obtenir una solució, fins i tot quan l'optimització tradicional té dificultats. Les tècniques d'optimització es combinen amb el mètode d'elements finits (MEF) per aconseguir la màxima atenuació possible en el rang de freqüències d'interès. S'ha definit un problema d'optimització d'un silenciador multicàmera i s'han dut a terme diverses anàlisis per obtenir la configuració més adequada per a cada cas. Sota certes hipòtesis d'uniformitat axial, s'han considerat diverses tècniques per reduir el cost computacional d'una anàlisi 3D complet per silenciadors dissipatius amb gradients de temperatura i flux mig. Aquestes es basen en la descomposició del camp acústic en modes axials i transversals dins de cada subdomini, i un procediment d'acoblament de les expansions modals en els canvis de secció del silenciador mitjançant les condicions de continuïtat dels camps acústics. S'estudia l'eficiència computacional i precisió de les prediccions de les tècniques d'acoblament, incloent col·locació puntual en nodes i punts de Gauss, així com ajust modal. Tots ells proporcionen prediccions precises de l'atenuació millorant el cost computacional d'EF. / Sánchez Orgaz, EM. (2016). Advanced numerical techniques for the acoustic modelling of materials and noise control devices in the exhaust system of internal combustion engines [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/64090 / TESIS

Dans cette thèse nous abordons des aspects topologiques de la matière condensée. Les états topologiques sont insensibles à un large spectre des perturbations externes et au désordre – une propriété indispensable dans le domaine d'information quantique. L’effet des interactions dans des systèmes topologiques est pourtant loin d’être bien maîtrisé à ce jour. Dans ce travail, nous étudions la corrélation entre la description topologique et l'effet des interactions. Afin d'accomplir notre but, nous utilisons des méthodes analytiques et numériques. Nous nous intéressons aussi à des sondes expérimentales qui peuvent être utilisées pour vérifier nos prédictions théoriques. Tout d’abord, nous étudions la version bosonique en interactions du modèle de Haldane – le modèle célèbre qui décrit l’effet Hall anomal. Nous proposons son implémentation expérimentale dans des circuits quantiques, basée sur l’application de perturbation périodique dépendantes du temps – méthodologie qui s’appelle l’ingénierie de Floquet. En poursuivant ces idées, nous étudions la version bosonique du modèle de Kane-Mele d’un isolant topologique. Ce modèle possède un diagramme de phase très riche. En particulier, lorsque les interactions sont fortes, nous observons l’apparition d’un modèle de magnétisme frustrée présentant une variété d'états exotiques. La mise en œuvre de ces modèles dans des réseaux d'atomes ultra-froids ou des circuits quantiques permettra de sonder expérimentalement les propriétés exotiques que nous avons observées. Ensuite, nous abordons d’une manière plus détaillée la réalisation expérimentale des modèles topologiques dans des circuits quantiques, en considérant le cas particulier du modèle de Su-Schrieffer-Heeger en couplage fort. Nous testons aussi des nouvelles sondes qui peuvent être utilisées afin de mesurer la phase de Zak et en déduire la topologie du système. Finalement, nous nous intéressons aux sondes hors d’équilibre et des méthodes pour tester les propriétés spectrales de systèmes quantiques, en utilisant l’approche de purification, pertinent pour le numérique et les expériences. / In this thesis we study the topological aspects of condensed matter physics, that received a revolutionary development in the last decades. Topological states of matter are protected against perturbations and disorder, making them very promising in the context of quantum information. The interplay between topology and interactions in such systems is however far from being well understood, while the experimental realization is challenging. Thus, in this work we investigate analytically such strongly correlated states of matter and explore new protocols to probe experimentally their properties. In order to do this, we use various analytical and numerical techniques. First, we analyze the properties of an interacting bosonic version of the celebrated Haldane model – the model for the quantum anomalous Hall effect. We propose its quantum circuit implementation based on the application of periodic time-dependent perturbations – Floquet engineering. Continuing these ideas, we study the interacting bosonic version of the Kane-Mele model – the first model of a topological insulator. This model has a very rich phase diagram with an emergence of an effective frustrated magnetic model and a variety of symmetry broken spin states in the strongly interacting regime. Ultra-cold atoms or quantum circuits implementation of both Haldane and Kane-Mele bosonic models would allow for experimental probes of the exotic states we observed. Second, in order to deepen the perspectives of quantum circuit simulations of topological phases we analyze the strong coupling limit of the Su-Schrieffer-Heeger model and we test new experimental probes of its topology associated with the Zak phase. We also work on the out-of-equilibrium protocols to study bulk spectral properties of quantum systems and quantum phase transitions using a purification scheme which could be implemented both numerically and experimentally.

Kronecker structured representations are used to cope with the state space explosion problem in Markovian modeling and analysis. Currently an open research problem is that of devising strong preconditioners to be used with projection methods for the computation of the stationary vector of Markov chains (MCs) underlying such representations. This paper proposes a block SOR (BSOR) preconditioner for hierarchical Markovian Models (HMMs) that are composed of multiple low level models and a high level model that defines the interaction among low level models. The Kronecker structure of an HMM yields nested block partitionings in its underlying continuous-time MC which may be used in the BSOR preconditioner. The computation of the BSOR preconditioned residual in each iteration of a preconditioned projection method becoms the problem of solving multiple nonsingular linear systems whose coefficient matrices are the diagonal blocks of the chosen partitioning. The proposed BSOR preconditioner solvers these systems using sparse LU or real Schur factors of diagonal blocks. The fill-in of sparse LU factorized diagonal blocks is reduced using the column approximate minimum degree algorithm (COLAMD). A set of numerical experiments are presented to show the merits of the proposed BSOR preconditioner.

Hierarchical Markovian Models (HMMs) are composed of multiple low level models (LLMs) and high level model (HLM) that defines the interaction among LLMs. The essence of the HMM approach is to model the system at hand in the form of interacting components so that its (larger) underlying continous-time Markov chain (CTMC) is not generated but implicitly represented as a sum of Kronecker products of (smaller) component matrices. The Kronecker structure of an HMM induces nested block partitionings in its underlying CTMC. These partitionings may be used in block versions of classical iterative methods based on splittings, such as block SOR (BSOR), to solve the underlying CTMC for its stationary vector. Therein the problem becomes that of solving multiple nonsingular linear systems whose coefficient matrices are the diagonal blocks of a particular partitioning. This paper shows that in each HLM state there may be diagonal blocks with identical off-diagonal parts and diagonals differing from each other by a multiple of the identity matrix. Such diagonal blocks are named candidate blocks. The paper explains how candidate blocks can be detected and how the can mutually benefit from a single real Schur factorization. It gives sufficient conditions for the existence of diagonal blocks with real eigenvalues and shows how these conditions can be checked using component matrices. It describes how the sparse real Schur factors of candidate blocks satisfying these conditions can be constructed from component matrices and their real Schur factors. It also demonstrates how fill in- of LU factorized (non-candidate) diagonal blocks can be reduced by using the column approximate minimum degree algorithm (COLAMD). Then it presents a three-level BSOR solver in which the diagonal blocks at the first level are solved using block Gauss-Seidel (BGS) at the second and the methods of real Schur and LU factorizations at the third level. Finally, on a set of numerical experiments it shows how these ideas can be used to reduce the storage required by the factors of the diagonal blocks at the third level and to improve the solution time compared to an all LU factorization implementation of the three-level BSOR solver.