No description available.

<p>Ultrarelativistic heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) are thought to have produced a state of matter called the Quark-Gluon-Plasma (QGP). The QGP forms when nuclear matter governed by Quantum Chromodynamics (QCD) reaches a temperature and baryochemical potential necessary to achieve the transition of hadrons (bound states of quarks and gluons) to {it deconfined} quarks and gluons. Such conditions have been achieved at RHIC, and the resulting QGP created exhibits properties of a near perfect fluid. In particular, strong evidence shows that the QGP exhibits a very small shear viscosity to entropy density ratio &eta/s, near the lower bound predicted for that quantity by Anti-deSitter space/Conformal Field Theory (AdS/CFT) methods of &eta/s = $hbar$/ 4 &pi $k_B$ where $hbar$ is Planck's constant and $k_B$ is Boltzmann's constant. As the produced matter expands and cools, it evolves through a phase described by a hadron gas with rapidly increasing $eta/s$.</p><p>This thesis presents robust calculations of $eta/s$ for hadronic and partonic media as a function of temperature using the Green-Kubo formalism. An analysis is performed for the behavior of $eta/s$ to mimic situations of the hadronic media at RHIC evolving out of chemical equilibrium, and systematic uncertainties are assessed for our method. In addition, preliminary results are presented for the bulk viscosity to entropy density ratio $zeta/s$, whose behavior is not well-known in a relativistic heavy ion collisions. The diffusion coefficient for baryon number is investigated, and an algorithm is presented to improve upon the previous work of investigation of heavy quark diffusion in a thermal QGP. </p><p>By combining the results of my investigations for $eta/s$ from our microscopic transport models with what is currently known from the experimental results on elliptic flow from RHIC, I find that the trajectory of $eta/s$ in a heavy ion collision has a rich structure, especially near the deconfinement transition temperature $T_c$. I have helped quantify the viscous hadronic effects to enable investigators to constrain the value of $eta/s$ for the QGP created at RHIC.</p> / Dissertation

[pt] Determinou-se experimentalmente o coeficiente local de troca de calor, para o escoamento turbulento interno, em tubos de seção transversal elíptica. As condições de contorno impostas foram, temperatura constante na parede e entrada abrupta, com canto vivo, em uma placa plana. Foram obtidos resultados na região de entrada e na região desenvolvida, cobrindo a faixa de número de Reynolds de 7 000 a 60 000. Foram investigadas as razões de aspecto 0,12, 0,25 e 0,5 e os resultados, comparados com os encontrados na literatura para placa plana e tubo circular. Uma diferença máxima de 16 por cento, foi encontrada em relação aos resultados para tubo circular. Portanto, a utilização de tubos elípticos em trocadores de calor é, em geral, vantajosa do ponto de vista térmico. / [en] Experiments were performed to determine local heat transfer coefficients for the turbulent flow in a duct with elliptical cross section. The boundary conditions imposed were isothermal wall and abrupt-contraction at the entrance built into a large wall. Both entrance-region and fully-developed results were obtained, whereas the Reynolds number was varied in the overall range 7,000-60,000. The aspect ratios 0.12, 0.25 and 0.5 were investigated and the results compared with the ones for parallel planes and for circular tubes, found in the open literature. The measured heat transfer coefficients displayed a maximum deviation of 16 percent relative to the results for the circular tube. Therefore, using elliptical ducts in heat exchangers is in general advisable as far as heat transfer is concerned.

In this thesis we use the holographic gauge/gravity duality to study two different aspects of strongly coupled non-Abelian plasmas. In the first topic we study the effects of strong (Abelian) magnetic fields on the transport coefficients of a strongly coupled non-Abelian plasma. Due to the spatial anisotropy created by the magnetic field, the most general viscosity tensor of a magnetized plasma has 5 shear viscosity coefficients and 2 bulk viscosities. We use the holographic correspondence for a strongly coupled N=4 Supersymmetric Yang-Mills (SYM) plasma to evaluate the shear viscosity perpendicular to the magnetic field and the shear viscosity parallel to the field. In the presence of a magnetic field, the shear viscosity perpendicular to the field saturates the Kovtun-Son-Starinets viscosity bound while in the direction parallel to the field the bound is violated. The second topic investigated in this thesis is motivated by the study of the near equilibrium behavior of strongly interacting non-Abelian plasmas that display a critical point in their phase diagram. We focus on the spectra of non-hydrodynamic quasinormal modes of a strongly coupled N=4 SYM plasma in the presence of a chemical potential, which displays a critical point in equilibrium. Except close to the critical point, we observe that by increasing the chemical potential one generally increases the damping rate of the quasinormal modes, which leads to a reduction of the characteristic equilibration times in the dual strongly coupled plasma. However, as one approaches the critical point the typical equilibration time increases though its derivative with respect to the chemical potential diverges with an exponent equal to -1/2. We also find a purely imaginary non-hydrodynamical mode in the vector diffusion channel at nonzero chemical potential which dictates the equilibration time in this channel near the critical point. / Nessa tese usamos a dualidade holográfica calibre/gravidade para estudar dois aspectos diferentes de plasmas não-Abelianos fortemente acoplados. No primeiro tópico estudamos os efeitos de campos magnéticos (Abelianos) intensos sobre o coeficientes de transporte de um plasma não-Abeliano fortemente acoplado. Devido à anisotropia espacial criada pelo campo magnético, o tensor de viscosidade mais geral de um plasma magnetizado deve possuir 5 coeficientes de viscosidade de cisalhamento e 2 de viscosidade volumétrica. Usamos a correspondência holográfica para um plasma N=4 Supersimétrico de Yang-Mills (SYM) fortemente acoplado para calcular a viscosidade de cisalhamento perpendicular ao campo magnético e a viscosidade de cisalhamento paralela ao campo. Na presença do campo magnético, a viscosidade de cisalhamento perpendicular ao campo satura o limite viscoso de Kovtun-Son-Starinets enquanto que na direção paralela ao campo o limite é violado. O segundo tópico investigado nessa tese é motivado pelo estudo do comportamento próximo ao equilíbrio de plasmas não-Abelianos fortemente interagentes que exibem um ponto crítico em seus diagramas de fase. Focamos no espectro dos modos quasinormais não-hidrodinâmicos de um plasma N=4 SYM fortemente acoplado na presença de um potencial químico, que exibe um ponto crítico no equilíbrio. Exceto próximo ao ponto crítico, observamos que ao aumentar o potencial químico geralmente se intensifica a taxa de amortecimento dos modos quasinormais, que levam à redução dos tempos de equilibração característicos do plasma dual fortemente acoplado. Entretanto, aproximando-se do ponto crítico o tempo de equilibração típico aumenta embora sua derivada em relação ao potencial químico diverge com um expoente igual à -1/2. Encontramos também um modo não-hidrodinâmico puramente imaginário no canal de difusão vetorial com potencial químico não-nulo que dita o tempo de equilibração neste canal próximo do ponto crítico.

Desenvolvemos urna teoria cinética para urna mistura de gases ionizados em presença de campos elétricos e magnéticos. As leis de Ohm, Fourier e Navier-Stokes são obtidas por dois métodos distintos que se baseiam na equação de Boltzmann. Verificamos que o emprego de teoremas de representação torna o método de Chapman-Enskog mais direto. Entretanto o método combinado mostrou-se extremamente simples, onde os coeficientes de transporte são determinados através da inversão de tensores de segunda e quarta ordens. Calculamos também a integral de colisão para as possíveis interações em gases ionizados tais como, entre partículas carregadas, partícula carregada e partícula neutra e entre partículas neutras. Como uma aplicação do método combinado, determinamos os coeficientes de condutividade elétrica, condutividade térmica, coeficiente termo-elétrico e o coeficiente de viscosidade cisalhante para um gás totalmente ionizado. Obtemos seus respectivos gráficos, considerando então um gás ionizado formado a partir do gás de hélio. / We develop a kinetic theory for ionized gases mixtures under the presence of electric and magnetic fields. The laws of Ohm, Fourier and Navier-Stokes are obtained by two different methods based on the Boltzmann equation. We verify that the use of representation theorems makes the Chapman-Enskog method more direct. However the combined method shows up as extremely simple where the transport coefficients are determined through inversion of second-order and fourth order tensors. We calculate also the collision integrals for possible interactions in ionized gases like: between charged particles, between charged particles and neutral particles and between neutral particles. As an application of the combined method, we determine the electrical and thermal conductivity coefficients, thermo-electric and shear viscosity coefficients for a completely ionized gas. We obtain their respective graphics considering an ionized gas of helium.

In this thesis we use the holographic gauge/gravity duality to study two different aspects of strongly coupled non-Abelian plasmas. In the first topic we study the effects of strong (Abelian) magnetic fields on the transport coefficients of a strongly coupled non-Abelian plasma. Due to the spatial anisotropy created by the magnetic field, the most general viscosity tensor of a magnetized plasma has 5 shear viscosity coefficients and 2 bulk viscosities. We use the holographic correspondence for a strongly coupled N=4 Supersymmetric Yang-Mills (SYM) plasma to evaluate the shear viscosity perpendicular to the magnetic field and the shear viscosity parallel to the field. In the presence of a magnetic field, the shear viscosity perpendicular to the field saturates the Kovtun-Son-Starinets viscosity bound while in the direction parallel to the field the bound is violated. The second topic investigated in this thesis is motivated by the study of the near equilibrium behavior of strongly interacting non-Abelian plasmas that display a critical point in their phase diagram. We focus on the spectra of non-hydrodynamic quasinormal modes of a strongly coupled N=4 SYM plasma in the presence of a chemical potential, which displays a critical point in equilibrium. Except close to the critical point, we observe that by increasing the chemical potential one generally increases the damping rate of the quasinormal modes, which leads to a reduction of the characteristic equilibration times in the dual strongly coupled plasma. However, as one approaches the critical point the typical equilibration time increases though its derivative with respect to the chemical potential diverges with an exponent equal to -1/2. We also find a purely imaginary non-hydrodynamical mode in the vector diffusion channel at nonzero chemical potential which dictates the equilibration time in this channel near the critical point. / Nessa tese usamos a dualidade holográfica calibre/gravidade para estudar dois aspectos diferentes de plasmas não-Abelianos fortemente acoplados. No primeiro tópico estudamos os efeitos de campos magnéticos (Abelianos) intensos sobre o coeficientes de transporte de um plasma não-Abeliano fortemente acoplado. Devido à anisotropia espacial criada pelo campo magnético, o tensor de viscosidade mais geral de um plasma magnetizado deve possuir 5 coeficientes de viscosidade de cisalhamento e 2 de viscosidade volumétrica. Usamos a correspondência holográfica para um plasma N=4 Supersimétrico de Yang-Mills (SYM) fortemente acoplado para calcular a viscosidade de cisalhamento perpendicular ao campo magnético e a viscosidade de cisalhamento paralela ao campo. Na presença do campo magnético, a viscosidade de cisalhamento perpendicular ao campo satura o limite viscoso de Kovtun-Son-Starinets enquanto que na direção paralela ao campo o limite é violado. O segundo tópico investigado nessa tese é motivado pelo estudo do comportamento próximo ao equilíbrio de plasmas não-Abelianos fortemente interagentes que exibem um ponto crítico em seus diagramas de fase. Focamos no espectro dos modos quasinormais não-hidrodinâmicos de um plasma N=4 SYM fortemente acoplado na presença de um potencial químico, que exibe um ponto crítico no equilíbrio. Exceto próximo ao ponto crítico, observamos que ao aumentar o potencial químico geralmente se intensifica a taxa de amortecimento dos modos quasinormais, que levam à redução dos tempos de equilibração característicos do plasma dual fortemente acoplado. Entretanto, aproximando-se do ponto crítico o tempo de equilibração típico aumenta embora sua derivada em relação ao potencial químico diverge com um expoente igual à -1/2. Encontramos também um modo não-hidrodinâmico puramente imaginário no canal de difusão vetorial com potencial químico não-nulo que dita o tempo de equilibração neste canal próximo do ponto crítico.

Desenvolvemos urna teoria cinética para urna mistura de gases ionizados em presença de campos elétricos e magnéticos. As leis de Ohm, Fourier e Navier-Stokes são obtidas por dois métodos distintos que se baseiam na equação de Boltzmann. Verificamos que o emprego de teoremas de representação torna o método de Chapman-Enskog mais direto. Entretanto o método combinado mostrou-se extremamente simples, onde os coeficientes de transporte são determinados através da inversão de tensores de segunda e quarta ordens. Calculamos também a integral de colisão para as possíveis interações em gases ionizados tais como, entre partículas carregadas, partícula carregada e partícula neutra e entre partículas neutras. Como uma aplicação do método combinado, determinamos os coeficientes de condutividade elétrica, condutividade térmica, coeficiente termo-elétrico e o coeficiente de viscosidade cisalhante para um gás totalmente ionizado. Obtemos seus respectivos gráficos, considerando então um gás ionizado formado a partir do gás de hélio. / We develop a kinetic theory for ionized gases mixtures under the presence of electric and magnetic fields. The laws of Ohm, Fourier and Navier-Stokes are obtained by two different methods based on the Boltzmann equation. We verify that the use of representation theorems makes the Chapman-Enskog method more direct. However the combined method shows up as extremely simple where the transport coefficients are determined through inversion of second-order and fourth order tensors. We calculate also the collision integrals for possible interactions in ionized gases like: between charged particles, between charged particles and neutral particles and between neutral particles. As an application of the combined method, we determine the electrical and thermal conductivity coefficients, thermo-electric and shear viscosity coefficients for a completely ionized gas. We obtain their respective graphics considering an ionized gas of helium.

La modélisation d'un disjoncteur haute tension par un modèle MHD nécessite des banques de données de propriétés du plasma afin de simuler le comportement de l'arc électrique. L'hypothèse d'équilibre thermodynamique local est souvent utilisée afin de calculer ces propriétés. Cependant les conditions d'équilibre ne sont pas remplies près des parois ou des électrodes, ou lors du passage par zéro du courant. Les écarts à l'équilibre (thermique notamment) modifient alors considérablement la chimie et les propriétés du plasma. L'étude de ces propriétés dans ces conditions nécessite alors de supposer que le plasma est à deux températures (2T). Le calcul des propriétés du plasma à 2T fait l'objet de cette thèse. La première partie de cette thèse présente le contexte industriel à l'origine de cette étude. Les hypothèses de base qui sous-tendent l'hypothèse 2T sont ensuite rappelées. Dans cette partie une attention spéciale est portée aux températures caractéristiques de peuplement des modes d'énergies internes. Deux jeux d'hypothèses portant sur ces températures sont utilisés dans ces travaux, et le choix de ces hypothèses est discuté. La deuxième partie de ce travail est consacrée au calcul de composition d'un plasma de SF6-C2F4 à 2T. Ce calcul sera fait avec deux méthodes différentes : la première repose sur la loi d'action de masse étendue à 2T, et la seconde sur un calcul collisionnel-radiatif. Des exemples de composition obtenus avec ces deux méthodes sont présentés. Les hypothèses portant sur les températures de peuplement des niveaux d'énergies internes sont discutées au vu de ces résultats. Le troisième chapitre de cet ouvrage aborde le calcul des propriétés thermodynamiques à 2T du plasma. Les formulations théoriques de chacune de ces propriétés sont d'abord rappelés, et les résultats issus de ces expressions sont ensuite présentés et discutés, pour les deux méthodes de calcul de la composition. Le quatrième chapitre est dédié au calcul des coefficients de transports d'un plasma de SF6-C2F4 à 2T. Cette partie s'appuie sur une étude bibliographique des méthodes de calcul déjà existantes et des données indispensables à l'obtention de ces propriétés (intégrales de collision). Pour chaque propriété (viscosité, conductivité électrique et conductivité thermique) les différentes méthodes de calcul recensées dans la littérature sont comparées. Le choix de la technique de calcul la plus appropriée est réalisé par confrontation de résultats à l'ETL. Une attention toute particulière est portée au calcul de la partie réactive de la conductivité thermique, et une formulation adaptée aux besoins de ce travail est proposée. Les résultats issus de ces expressions sont présentés et discutés suivant la même logique quand dans le chapitre précédent. / Modelling a high voltage circuit breaker using a MHD model needs plasma properties databanks to simulate the electric arc behavior. The local thermodynamic equilibrium hypothesis is often used to calculate these properties. However, the equilibrium conditions are not satisfied near the walls or the electrodes, or during the zero crossing of the current. The thermal nonequilibrium considerably modify the chemistry and properties of the plasma. The study of these properties for a 2t plasma is the subject of this thesis. The first part of this thesis presents the industrial context at the origin of this study. Basics assumptions for the 2t hypothesis are then explained. In this section, special attention is given to the temperatures characteristic of the internal energy modes. Two sets of hypotheses concerning these temperatures are used in this work, and the choice of these hypotheses is discussed. The second part of this work is dedicated to the calculation of sf6-c2f4 plasma composition at 2t. This calculation will be done with two different methods: the first is based on the mass action law extended to 2t, and the second on a collisional-radiative calculation. Examples of compositions obtained with these two methods are presented. The hypotheses concerning the temperatures of populating of the internal energy levels are discussed in the light of these results. The third chapter of this thesis deals with the calculation of the thermodynamic properties at 2t of the plasma. The theoretical formulations of each of these properties are first recalled, and the results from these expressions are then presented and discussed, for the two methods of calculating the composition. The fourth chapter is dedicated to calculating the transport coefficients of a sf6-c2f4 plasma at 2t. This part is based on a bibliographic study of the already existing methods of calculation and the essential data to obtain these properties (collision integrals). For each property (viscosity, electrical conductivity and thermal conductivity) the various calculation methods identified in the literature are compared. The choice of the most appropriate calculation technique is made by comparing the results to the ETL. Particular attention is paid to the calculation of the reactive part of the thermal conductivity, and a formulation adapted to the needs of this work is proposed. The results from these expressions are presented and discussed following the same logic when in the previous chapter.

A fundamental challenge is to understand nonequilibrium statistical mechanics starting from microscopic chaos in the equations of motion of a many-particle system. In this thesis we summarize recent theoretical advances along these lines. We focus on two different approaches to nonequilibrium transport: One considers Hamiltonian dynamical systems under nonequilibrium boundary conditions, another one suggests a non-Hamiltonian approach to nonequilibrium situations created by external electric fields and by temperature or velocity gradients. A surprising result related to the former approach is that in simple low-dimensional periodic models the deterministic transport coefficients are typically fractal functions of control parameters. These fractal transport coefficients yield the first central theme of this thesis. We exemplify this phenomenon by deterministic diffusion in a simple chaotic map. We then construct an arsenal of analytical and numerical methods for computing further transport coefficients such as electrical conductivities andchemical reaction rates. These methods are applied to hierarchies of chaotic dynamical systems that are successively getting more complex, starting from abstract one-dimensional maps generalizing a simple random walk on the line up to particle billiards that should be directly accessible in experiments. In all cases, the resulting transport coefficients turn out to be either strictly fractal, or at least to be profoundly irregular. The impact of random perturbations on these quantities is also investigated. We furthermore provide some access roads towards a physical understanding of these fractalities. The second central theme is formed by a critical assessment of the non-Hamiltonian approach to nonequilibrium transport. Here we consider situations where the nonequilibrium constraints pump energy into a system, hence there must be some thermal reservoir that prevents the system from heating up. For this purpose a deterministic and time-reversible modeling of thermal reservoirs was proposed in form of Gaussian and Nose-Hoover thermostats. This approach yielded simple relations between fundamental quantities of nonequilibrium statistical mechanics and of dynamical systems theory. Our goal is to critically assesses the universality of these results. As a vehicle of demonstration we employ the driven periodic Lorentz gas, a toy model for the classical dynamics of an electron in a metal under application of an electric field. Applying different types of thermal reservoirs to this system we compare the resulting nonequilibrium steady states with each other. Along the same lines we discuss an interacting many-particle system under shear and heat. Finally, we outline an unexpected relationship between deterministic thermostats and active Brownian particles modeling biophysical cell motility.

Avec les récents développements obtenus grâce aux modèles globaux magnétohydrodynamiques en trois dimensions de la convection solaire, il est désormais possible de simuler des champs magnétiques structurés à grande échelle et présentant des inversions de polarité bien synchronisées dans chaque hémisphère. Ces modèles qui n'incluent, pour la plupart, aucune modélisation de la surface du Soleil génèrent donc leur dynamo avec l'action de la force électromotrice turbulente (FEM) et de la rotation différentielle uniquement. À partir de cette FEM, différentes techniques peuvent être utilisées pour extraire les coefficients de transport turbulent. Notamment, différents auteurs ont obtenu un tenseur alpha (coefficient du premier ordre) dont les 9 composantes présentent des amplitudes du même ordre, remettant en doute l'approximation faite dans le cas des modèles dynamo de type alphaOmega qui ne tient en compte qu’une de composante du tenseur. À partir d'un code d'analyse par décomposition en valeurs singulières pour évaluer les coefficients du tenseur alpha, nous avons généralisé la procédure pour extraire 18 des composantes du tenseur de deuxième ordre (tenseur beta). Les tenseurs alpha et beta obtenus par cette nouvelle procédure tel qu'appliquée aux sorties du modèle global EULAG-MHD, sont similaires aux tenseurs alpha et beta équivalant obtenus en utilisant l'approximation « Second Order Correlation Approximation ». À l'aide des coefficients de transport turbulent du premier ordre introduit dans un modèle dynamo en champ moyen, nous avons ensuite étudié certaines solutions magnétiques présentant des doubles dynamos. Cette analyse avait pour but de comparer les résultats obtenus par ce modèle simplifié caractérisé par la FEM provenant de EULAG-MHD aux résultats de EULAG-MHD directement. Cette preuve de concept nous a permis de démontrer que l'oscillation observée dans le champ magnétique en surface de EULAG-MHD pouvait provenir de l'action inductive d'une seconde dynamo. Une oscillation biennale est également observée dans plusieurs indices d'activité solaire dont l'origine n'est toujours pas établie. Il est évident que les deux modèles décrits ci-haut et le Soleil opèrent dans des régimes physiques différents. Toutefois, malgré leurs différences, le fait qu'il soit relativement facile de produire une seconde dynamo dans EULAG-MHD et dans le modèle en champ moyen indique que l'action inductive de la FEM peut facilement générer deux dynamos. Finalement, dans le but d'étudier les périodes de grands minima, phénomène encore non reproduit par les modèles globaux, nous avons ajouté une rétroaction magnétique sur l'écoulement azimutal au modèle dynamo cinématique en champ moyen décrit ci-haut. En analysant les solutions de ce modèle dynamo de type alpha2Omega non cinématique, nous avons pu reproduire la tendance observée jusqu'ici uniquement dans les modèles de type alphaOmega selon laquelle le nombre de nombre de Prandtl magnétique contrôle le rapport des périodes générées. De plus, en analysant une solution sur 50 000 ans présentant des périodes de grands minima et maxima non périodiques, nous avons obtenu une distribution de temps de séparation des grands minima presque exponentielle, caractéristique observée dans les reconstructions de l'activité solaire. La rotation différentielle associée à ces périodes de grands minima présente un niveau de fluctuation de 1% par rapport au profil moyen. Ce niveau de fluctuation est d'ailleurs comparable avec les reconstructions historiques de la rotation différentielle en surface obtenues lors du grand minimum de Maunder. / The recent developments achieved by tri-dimensionals magnetohydrodynamic (3D-MHD) global simulations of solar convection allow us to generate an organized large-scale magnetic fields with well-synchronized hemispheric polarity reversal. Because the vast majority of these simulations do not include a modelization of the Sun's surface layer, the generation of their dynamo is thus solely due to the action of the turbulent electromotive force (EMF) in conjunction with differential rotation. From this EMF, different methods can be used to extract the turbulent transport coefficients. In particular, various authors found a full 9 component alpha-tensor (first order coefficients) where all the components are of the same order of magnitude. This finding calls into question the alphaOmega approximation made by the vast majority of mean field dynamo models. We generalized a first order (alpha-tensor) singular value decomposition (SVD) analysis procedure to extract the 18 additional components of the second order tensor (beta-tensor). The alpha and beta tensors obtained by this new procedure as applied to the EULAG-MHD outputs, are similar to the equivalent alpha and beta tensors obtained using the second order correlation approximation (SOCA). By introducing the first order turbulent transport coefficients in a mean field dynamo model, we study the magnetic solutions where double dynamo modes were observed. This analysis allows us to compare the mean field dynamo solutions produced with the EMF, as extracted from EULAG-MHD, with the real magnetic output of EULAG-MHD. This proof of concept demonstrated that the quasi-biennal oscillation observed in the surface toroidal magnetic field in EULAG-MHD can be produced by the inductive action of a secondary dynamo. A similar quasi-biennal oscillation signal is also observed in multiple proxies of the solar activity whose origin is still not confirmed. Although the physical set of properties under which the two numerical models described above operate are different from the Sun, the fact that both models can reproduce a secondary dynamo shows us that the inductive action of the EMF can easily produce two dynamos. Finally, in order to study epochs of grand minima that still cannot be reproduced in global 3D-MHD simulations of convection, we added a magnetic feedback on the mean azimutal flow in our kinematic mean field model. This non-kinematic alpha2Omega model was able to reproduce the tendency of the Prandtl number (Pm) to control the ratio of the modulation period. More specifically, we found an inverse relation between Pm and the ratio of the main magnetic cycle period to the grand minima occurrence period. Moreover, by analyzing a simulation of a length of 50,000 years, where aperiodic periods of grand minima and maxima are observed, we found a waiting time distribution (WTD) of the grand minima close to an exponential, a characteristic also observed in the reconstruction of the solar activity. Finally, the level of fluctuation in the surface differential rotation associated with epochs of grand minima is ~1%. This level of fluctuation was also observed in historical reconstructions of the surface differential rotation during the Maunder minimum.