The Role of MHD Turbulence in Magnetic Self-Excitation in The

Total Page:16

File Type:pdf, Size:1020Kb

The Role of MHD Turbulence in Magnetic Self-Excitation in The THE ROLE OF MHD TURBULENCE IN MAGNETIC SELF-EXCITATION: A STUDY OF THE MADISON DYNAMO EXPERIMENT by Mark D. Nornberg A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Physics) at the UNIVERSITY OF WISCONSIN–MADISON 2006 °c Copyright by Mark D. Nornberg 2006 All Rights Reserved i For my parents who supported me throughout college and for my wife who supported me throughout graduate school. The rest of my life I dedicate to my daughter Margaret. ii ACKNOWLEDGMENTS I would like to thank my adviser Cary Forest for his guidance and support in the completion of this dissertation. His high expectations and persistence helped drive the work presented in this thesis. I am indebted to him for the many opportunities he provided me to connect with the world-wide dynamo community. I would also like to thank Roch Kendrick for leading the design, construction, and operation of the experiment. He taught me how to do science using nothing but duct tape, Sharpies, and Scotch-Brite. He also raised my appreciation for the artistry of engineer- ing. My thanks also go to the many undergraduate students who assisted in the construction of the experiment, especially Craig Jacobson who performed graduate-level work. My research partner, Erik Spence, deserves particular thanks for his tireless efforts in modeling the experiment. His persnickety emendations were especially appreciated as we entered the publi- cation stage of the experiment. The conversations during our morning commute to the lab will be sorely missed. I never imagined forging such a strong friendship with a colleague, and I hope our families remain close despite great distance. I would also like to thank Jim, Mike, and Julie for both moral support and epicurean sustenance. I would never have endured the struggle of graduate school without the steadfast support of my family. I will always be grateful for the education I received through the sacrifices of my parents. My beautiful wife Julianne has been my comfort and encouragement throughout my graduate career. She checked my discouragement with unwavering optimism and provided the hope that is needed in accomplishing such an endeavor. iii It seems an almost incredible thing to me that an invention of the human intellect and the structure of the universe coincide. — Benedict XVI iv TABLE OF CONTENTS Page List of tables .......................................... vi List of figures .......................................... vii Nomenclature ......................................... x Abstract ............................................. xi 1 Introduction ........................................ 1 1.1 A brief history of the dynamo problem ........................ 1 1.2 Dynamo experiments ................................. 3 1.3 Motivation for the experiment ............................ 6 1.4 Outline of the thesis .................................. 8 2 The kinematic dynamo problem ............................. 10 2.1 The magnetic induction equation ........................... 10 2.2 The Bullard and Gellman formalism ......................... 12 2.3 The kinematic eigenvalue problem .......................... 15 2.4 Simple roll flow dynamos ............................... 17 2.5 The stretch-twist-fold mechanism .......................... 21 2.6 Comparison with simulation ............................. 24 2.7 Anticipated effects of turbulence ........................... 27 2.7.1 Intermittent excitations ............................ 27 2.7.2 Mean-field dynamo .............................. 27 2.7.3 Small-scale dynamo ............................. 30 3 Description of the sodium experiment .......................... 32 3.1 Industrial automation ................................. 36 3.2 Holding tank ..................................... 36 3.3 Experimental vessel .................................. 38 3.4 Expansion tank .................................... 41 3.5 Melt station ...................................... 41 3.6 Pneumatic transfer system .............................. 42 v Page 3.7 Magnetic sensors and data acquisition ........................ 43 3.8 External field coils .................................. 46 3.9 Safety equipment ................................... 47 3.9.1 Protective clothing and gear ......................... 47 3.9.2 Scrubber system ............................... 47 3.9.3 Fume hood for small sodium fires and material tests ............ 48 3.9.4 Sodium cleaning station ........................... 48 4 Hydrodynamic experiments to model the flow ..................... 50 4.1 Description of the water experiment ......................... 50 4.2 Laser Doppler velocimetry .............................. 55 4.3 Model of the mean flow ................................ 56 4.4 Velocity fluctuations ................................. 62 4.5 Kinematic growth rate of the magnetic field ..................... 62 4.6 Kolmogorov turbulence ................................ 66 4.7 Measurements of the velocity spectrum ....................... 70 4.8 Turbulent conductivity ................................ 75 5 Comparison of measurements of the induced field with predictions ......... 80 5.1 The predicted mean induced magnetic field ..................... 80 5.2 Measurements of the magnetic field ......................... 81 5.3 Reconstruction of the mean magnetic field ...................... 85 5.4 Magnetic field amplification and feedback ...................... 87 5.5 The magnetic spectrum due to turbulence ...................... 94 5.6 Measurements of the magnetic spectrum ....................... 98 5.7 Influence of the back reaction on turbulence ..................... 103 6 Observation of an intermittently excited magnetic field ................ 104 7 Summary and discussion ................................. 112 References ........................................... 115 APPENDICES Appendix A: Calculating the Mode Energy ....................... 128 vi LIST OF TABLES Table Page 2.1 Optimized flow parameters ................................ 20 3.1 Experiment parameters .................................. 35 4.1 Turbulence characteristics ................................ 68 6.1 Results from conditional averaging ........................... 109 vii LIST OF FIGURES Figure Page 1.1 Schematic of the sodium experiment ........................... 7 2.1 Dudley and James t2s2 flow ............................... 18 2.2 Growth rate plot of the Dudley and James t2s2 flow .................. 19 2.3 Stretch-twist-fold mechanism .............................. 23 2.4 Structure of the excited magnetic field .......................... 25 2.5 Growth rates for different flows ............................. 26 3.1 Photo of the Madison Dynamo Experiment ....................... 33 3.2 Sodium conductivity ................................... 34 3.3 Automation software display ............................... 37 3.4 Motor power ....................................... 40 3.5 Photo of magnetic sensors ................................ 44 4.1 Photo of the water experiment .............................. 51 4.2 Photo of the impellers .................................. 52 4.3 Schematic of the water experiment ............................ 57 4.4 LDV time series measurements ............................. 59 4.5 Measurements of the mean flow ............................. 60 4.6 Reconstruction of the mean flow ............................. 61 viii Figure Page 4.7 PDF of velocity fluctuations ............................... 63 4.8 Contour plot of velocity fluctuation levels ........................ 64 4.9 Eigenmode growth rate .................................. 65 4.10 Particle arrival time .................................... 71 4.11 Sample and hold interpolation .............................. 72 4.12 Velocity power spectra from LDV ............................ 74 4.13 Velocity correlation functions .............................. 76 4.14 Turbulent conductivity .................................. 78 4.15 Spatial distribution of turbulent conductivity ...................... 79 5.1 Time series measurements of the magnetic field ..................... 83 5.2 Internal magnetic field measurements .......................... 84 5.3 Spherical harmonic modes of the field induced by the flow ............... 88 5.4 Reconstruction of the magnetic field induced by the flow ................ 89 5.5 Predicted gain and growth rate of induced transverse dipole field ............ 92 5.6 Gain and orientation of induced transverse dipole field ................. 93 5.7 Magnetic power spectrum from simulation ....................... 96 5.8 Magnetic power spectra measurements ......................... 99 5.9 Magnetic spectrum varying applied field strength .................... 101 5.10 Dependence of dissipation scale on Rm ......................... 102 6.1 PDF of magnetic field fluctuations ............................ 106 6.2 Surface magnetic field during burst ........................... 107 6.3 Time series of transverse dipole field .......................... 108 ix Appendix Figure Page 6.4 Conditionally averaged magnetic bursts ......................... 108 6.5 Energy PDFs of transverse dipole field .......................... 110 x Nomenclature MHD Magnetohydrodynamics LDV Laser Doppler Velocimetry VKS Von Karm´ an´ Sodium Experiment PDF Probability Distribution Function PSD Power Spectral Density ACF Autocorrelation Function xi Abstract Determining the onset conditions for magnetic field growth in magnetohydrodynamics is fun- damental to understanding how astrophysical dynamos such as the
Recommended publications
  • Chapter 5 Dimensional Analysis and Similarity
    Chapter 5 Dimensional Analysis and Similarity Motivation. In this chapter we discuss the planning, presentation, and interpretation of experimental data. We shall try to convince you that such data are best presented in dimensionless form. Experiments which might result in tables of output, or even mul- tiple volumes of tables, might be reduced to a single set of curves—or even a single curve—when suitably nondimensionalized. The technique for doing this is dimensional analysis. Chapter 3 presented gross control-volume balances of mass, momentum, and en- ergy which led to estimates of global parameters: mass flow, force, torque, total heat transfer. Chapter 4 presented infinitesimal balances which led to the basic partial dif- ferential equations of fluid flow and some particular solutions. These two chapters cov- ered analytical techniques, which are limited to fairly simple geometries and well- defined boundary conditions. Probably one-third of fluid-flow problems can be attacked in this analytical or theoretical manner. The other two-thirds of all fluid problems are too complex, both geometrically and physically, to be solved analytically. They must be tested by experiment. Their behav- ior is reported as experimental data. Such data are much more useful if they are ex- pressed in compact, economic form. Graphs are especially useful, since tabulated data cannot be absorbed, nor can the trends and rates of change be observed, by most en- gineering eyes. These are the motivations for dimensional analysis. The technique is traditional in fluid mechanics and is useful in all engineering and physical sciences, with notable uses also seen in the biological and social sciences.
    [Show full text]
  • Fluid Mechanics of Liquid Metal Batteries
    Fluid Mechanics of Liquid Metal Douglas H. Kelley Batteries Department of Mechanical Engineering, University of Rochester, The design and performance of liquid metal batteries (LMBs), a new technology for grid- Rochester, NY 14627 scale energy storage, depend on fluid mechanics because the battery electrodes and elec- e-mail: [email protected] trolytes are entirely liquid. Here, we review prior and current research on the fluid mechanics of LMBs, pointing out opportunities for future studies. Because the technology Tom Weier in its present form is just a few years old, only a small number of publications have so far Institute of Fluid Dynamics, considered LMBs specifically. We hope to encourage collaboration and conversation by Helmholtz-Zentrum Dresden-Rossendorf, referencing as many of those publications as possible here. Much can also be learned by Bautzner Landstr. 400, linking to extensive prior literature considering phenomena observed or expected in Dresden 01328, Germany LMBs, including thermal convection, magnetoconvection, Marangoni flow, interface e-mail: [email protected] instabilities, the Tayler instability, and electro-vortex flow. We focus on phenomena, materials, length scales, and current densities relevant to the LMB designs currently being commercialized. We try to point out breakthroughs that could lead to design improvements or make new mechanisms important. [DOI: 10.1115/1.4038699] 1 Introduction of their design speed, in order for the grid to function properly. Changes to any one part of the grid affect all parts of the grid. The The story of fluid mechanics research in LMBs begins with implications of this interconnectedness are made more profound one very important application: grid-scale storage.
    [Show full text]
  • Cavitation in Valves
    VM‐CAV/WP White Paper Cavitation in Valves Table of Contents Introduction. 2 Cavitation Analysis. 2 Cavitation Data. 3 Valve Coefficient Data. 4 Example Application. .. 5 Conclusion & Recommendations . 5 References. 6 Val‐Matic Valve & Mfg. Corp. • www.valmatic.com • [email protected] • PH: 630‐941‐7600 Copyright © 2018 Val‐Matic Valve & Mfg. Corp. Cavitation in Valves INTRODUCTION Cavitation can occur in valves when used in throttling or modulating service. Cavitation is the sudden vaporization and violent condensation of a liquid downstream of the valve due to localized low pressure zones. When flow passes through a throttled valve, a localized low pressure zone forms immediately downstream of the valve. If the localized pressure falls below the vapor pressure of the fluid, the liquid vaporizes (boils) and forms a vapor pocket. As the vapor bubbles flow downstream, the pressure recovers, and the bubbles violently implode causing a popping or rumbling sound similar to tumbling rocks in a pipe. The sound of cavitation in a pipeline is unmistakable. The condensation of the bubbles not only produces a ringing sound, but also creates localized stresses in the pipe walls and valve body that can cause severe pitting. FIGURE 1. Cavitation Cavitation is a common occurrence in shutoff valves during the last few degrees of closure when the supply pressure is greater than about 100 psig. Valves can withstand limited durations of cavitation, but when the valve must be throttled or modulated in cavitating conditions for long periods of time, the life of the valve can be drastically reduced. Therefore, an analysis of flow conditions is needed when a valve is used for flow or pressure control.
    [Show full text]
  • Control Valve Sizing Theory, Cavitation, Flashing Noise, Flashing and Cavitation Valve Pressure Recovery Factor
    Control Valve Sizing Theory, Cavitation, Flashing Noise, Flashing and Cavitation Valve Pressure Recovery Factor When a fluid passes through the valve orifice there is a marked increase in velocity. Velocity reaches a maximum and pressure a minimum at the smallest sectional flow area just downstream of the orifice opening. This point of maximum velocity is called the Vena Contracta. Downstream of the Vena Contracta the fluid velocity decelerates and the pressure increases of recovers. The more stream lined valve body designs like butterfly and ball valves exhibit a high degree of pressure recovery where as Globe style valves exhibit a lower degree of pressure recovery because of the Globe geometry the velocity is lower through the vena Contracta. The Valve Pressure Recovery Factor is used to quantify this maximum velocity at the vena Contracta and is derived by testing and published by control valve manufacturers. The Higher the Valve Pressure Recovery Factor number the lower the downstream recovery, so globe style valves have high recovery factors. ISA uses FL to represent the Valve Recovery Factor is valve sizing equations. Flow Through a restriction • As fluid flows through a restriction, the Restriction Vena Contracta fluid’s velocity increases. Flow • The Bernoulli Principle P1 P2 states that as the velocity of a fluid or gas increases, its pressure decreases. Velocity Profile • The Vena Contracta is the point of smallest flow area, highest velocity, and Pressure Profile lowest pressure. Terminology Vapor Pressure Pv The vapor pressure of a fluid is the pressure at which the fluid is in thermodynamic equilibrium with its condensed state.
    [Show full text]
  • Anomalous Viscosity, Resistivity, and Thermal Diffusivity of the Solar
    Anomalous Viscosity, Resistivity, and Thermal Diffusivity of the Solar Wind Plasma Mahendra K. Verma Department of Physics, Indian Institute of Technology, Kanpur 208016, India November 12, 2018 Abstract In this paper we have estimated typical anomalous viscosity, re- sistivity, and thermal difffusivity of the solar wind plasma. Since the solar wind is collsionless plasma, we have assumed that the dissipation in the solar wind occurs at proton gyro radius through wave-particle interactions. Using this dissipation length-scale and the dissipation rates calculated using MHD turbulence phenomenology [Verma et al., 1995a], we estimate the viscosity and proton thermal diffusivity. The resistivity and electron’s thermal diffusivity have also been estimated. We find that all our transport quantities are several orders of mag- nitude higher than those calculated earlier using classical transport theories of Braginskii. In this paper we have also estimated the eddy turbulent viscosity. arXiv:chao-dyn/9509002v1 5 Sep 1995 1 1 Introduction The solar wind is a collisionless plasma; the distance travelled by protons between two consecutive Coulomb collisions is approximately 3 AU [Barnes, 1979]. Therefore, the dissipation in the solar wind involves wave-particle interactions rather than particle-particle collisions. For the observational evidence of the wave-particle interactions in the solar wind refer to the review articles by Gurnett [1991], Marsch [1991] and references therein. Due to these reasons for the calculations of transport coefficients in the solar wind, the scales of wave-particle interactions appear more appropriate than those of particle-particle interactions [Braginskii, 1965]. Note that the viscosity in a turbulent fluid is scale dependent.
    [Show full text]
  • Eulerian–Lagrangian Method for Simulation of Cloud Cavitation ∗ Kazuki Maeda , Tim Colonius
    Journal of Computational Physics 371 (2018) 994–1017 Contents lists available at ScienceDirect Journal of Computational Physics www.elsevier.com/locate/jcp Eulerian–Lagrangian method for simulation of cloud cavitation ∗ Kazuki Maeda , Tim Colonius Division of Engineering and Applied Science, California Institute of Technology, 1200 East California Boulevard, Pasadena, CA 91125, USA a r t i c l e i n f o a b s t r a c t Article history: We present a coupled Eulerian–Lagrangian method to simulate cloud cavitation in a Received 3 December 2017 compressible liquid. The method is designed to capture the strong, volumetric oscillations Received in revised form 10 April 2018 of each bubble and the bubble-scattered acoustics. The dynamics of the bubbly mixture Accepted 16 May 2018 is formulated using volume-averaged equations of motion. The continuous phase is Available online 18 May 2018 discretized on an Eulerian grid and integrated using a high-order, finite-volume weighted Keywords: essentially non-oscillatory (WENO) scheme, while the gas phase is modeled as spherical, Bubble dynamics Lagrangian point-bubbles at the sub-grid scale, each of whose radial evolution is tracked Cavitation by solving the Keller–Miksis equation. The volume of bubbles is mapped onto the Eulerian Eulerian–Lagrangian method grid as the void fraction by using a regularization (smearing) kernel. In the most general Compressible multiphase flows case, where the bubble distribution is arbitrary, three-dimensional Cartesian grids are used Multiscale modeling for spatial discretization. In order to reduce the computational cost for problems possessing Reduced-order modeling translational or rotational homogeneities, we spatially average the governing equations along the direction of symmetry and discretize the continuous phase on two-dimensional or axi-symmetric grids, respectively.
    [Show full text]
  • A Numerical Algorithm for MHD of Free Surface Flows at Low Magnetic Reynolds Numbers
    A Numerical Algorithm for MHD of Free Surface Flows at Low Magnetic Reynolds Numbers Roman Samulyak1, Jian Du2, James Glimm1;2, Zhiliang Xu1 1Computational Science Center, Brookhaven National Laboratory, Upton, NY 11973 2 Department of Applied Mathematics and Statistics, SUNY at Stony Brook, Stony Brook, NY 11794, USA November 7, 2005 Abstract We have developed a numerical algorithm and computational soft- ware for the study of magnetohydrodynamics (MHD) of free surface flows at low magnetic Reynolds numbers. The governing system of equations is a coupled hyperbolic/elliptic system in moving and ge- ometrically complex domains. The numerical algorithm employs the method of front tracking for material interfaces, high resolution hy- perbolic solvers, and the embedded boundary method for the elliptic problem in complex domains. The numerical algorithm has been imple- mented as an MHD extension of FronTier, a hydrodynamic code with free interface support. The code is applicable for numerical simulations of free surface conductive liquids or flows of weakly ionized plasmas. Numerical simulations of the Muon Collider/Neutrino Factory target have been discussed. 1 Introduction Computational magnetohydrodynamics, greatly inspired over the last decades by the magnetic confinement fusion and astrophysics problems, has achieved significant results. However the major research effort has been in the area of highly ionized plasmas. Numerical methods and computational software for MHD of weakly conducting materials such as liquid metals or weakly ionized plasmas have not been developed to such an extent despite the need 1 for fusion research and industrial technologies. Liquid metal MHD, driven by potential applications of flowing liquid metals or electrically conducting liquid salts as coolant in magnetic confinement fusion reactors as well as some industrial problems, has attracted broad theoretical, computational, and experimental studies (see [16, 17, 18] and references therein).
    [Show full text]
  • Rayleigh-Bernard Convection Without Rotation and Magnetic Field
    Nonlinear Convection of Electrically Conducting Fluid in a Rotating Magnetic System H.P. Rani Department of Mathematics, National Institute of Technology, Warangal. India. Y. Rameshwar Department of Mathematics, Osmania University, Hyderabad. India. J. Brestensky and Enrico Filippi Faculty of Mathematics, Physics, Informatics, Comenius University, Slovakia. Session: GD3.1 Core Dynamics Earth's core structure, dynamics and evolution: observations, models, experiments. Online EGU General Assembly May 2-8 2020 1 Abstract • Nonlinear analysis in a rotating Rayleigh-Bernard system of electrical conducting fluid is studied numerically in the presence of externally applied horizontal magnetic field with rigid-rigid boundary conditions. • This research model is also studied for stress free boundary conditions in the absence of Lorentz and Coriolis forces. • This DNS approach is carried near the onset of convection to study the flow behaviour in the limiting case of Prandtl number. • The fluid flow is visualized in terms of streamlines, limiting streamlines and isotherms. The dependence of Nusselt number on the Rayleigh number, Ekman number, Elasser number is examined. 2 Outline • Introduction • Physical model • Governing equations • Methodology • Validation – RBC – 2D – RBC – 3D • Results – RBC – RBC with magnetic field (MC) – Plane layer dynamo (RMC) 3 Introduction • Nonlinear interaction between convection and magnetic fields (Magnetoconvection) may explain certain prominent features on the solar surface. • Yet we are far from a real understanding of the dynamical coupling between convection and magnetic fields in stars and magnetically confined high-temperature plasmas etc. Therefore it is of great importance to understand how energy transport and convection are affected by an imposed magnetic field: i.e., how the Lorentz force affects convection patterns in sunspots and magnetically confined, high-temperature plasmas.
    [Show full text]
  • Toward a Self-Generating Magnetic Dynamo: the Role of Turbulence
    PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 2000 Toward a self-generating magnetic dynamo: The role of turbulence Nicholas L. Peffley, A. B. Cawthorne, and Daniel P. Lathrop* Department of Physics, University of Maryland, College Park, Maryland 20742 ͑Received 6 July 1999͒ Turbulent flow of liquid sodium is driven toward the transition to self-generating magnetic fields. The approach toward the transition is monitored with decay measurements of pulsed magnetic fields. These mea- surements show significant fluctuations due to the underlying turbulent fluid flow field. This paper presents experimental characterizations of the fluctuations in the decay rates and induced magnetic fields. These fluc- tuations imply that the transition to self-generation, which should occur at larger magnetic Reynolds number, will exhibit intermittent bursts of magnetic fields. PACS number͑s͒: 47.27.Ϫi, 47.65.ϩa, 05.45.Ϫa, 91.25.Cw I. INTRODUCTION Reynolds number will be quite large for all flows attempting to self-generate ͑where Re ӷ1 yields Reӷ105)—implying The generation of magnetic fields from flowing liquid m turbulent flow. These turbulent flows will cause the transition metals is being pursued by a number of scientific research to self-generation to be intermittent, showing both growth groups in Europe and North America. Nuclear engineering and decay of magnetic fields irregularly in space and time. has facilitated the safe use of liquid sodium, which has con- This intermittency is not something addressed by kinematic tributed to this new generation of experiments. With the dynamo studies. The analysis in this paper focuses on three highest electrical conductivity of any liquid, sodium retains main points: we quantify the approach to self-generation and distorts magnetic fields maximally before they diffuse with increasing Rem , characterize the turbulence of induced away.
    [Show full text]
  • Numerical Study of Cavitation Within Orifice Flow
    NUMERICAL STUDY OF CAVITATION WITHIN ORIFICE FLOW A Thesis by PENGZE YANG Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Chair of Committee, Robert Handler Committee Members, David Staack Prabir Daripa Head of Department, Andreas Polycarpou December 2015 Major Subject: Mechanical Engineering Copyright 2015 Pengze Yang ABSTRACT Cavitation generally occurs when the pressure at certain location drops to the vapor pressure and the liquid water evaporates as a consequence. For the past several decades, numerous experimental researches have been conducted to investigate this phenomenon due to its degradation effects on hydraulic device structures, such as erosion, noise and vibration. A plate orifice is an important restriction device that is widely used in many industries. It serves functions as restricting flow and measuring flow rate within a pipe. The plate orifice is also subject to intense cavitation at high pressure difference, therefore, the simulation research of the cavitation phenomenon within an orifice flow becomes quite essential for understanding the causes of cavitation and searching for possible preventing methods. In this paper, all researches are simulation-oriented by using ANSYS FLUENT due to its high resolution comparing to experiments. Standard orifice plates based on ASME PTC 19.5-2004 are chosen and modeled in the study with the diameter ratio from 0.2 to 0.75. Steady state studies are conducted for each diameter ratio at the cavitation number roughly from 0.2 to 2.5 to investigate the dependency of discharge coefficient on the cavitation number.
    [Show full text]
  • Turbulent Magnetic Prandtl Number and Magnetic Diffusivity Quenching from Simulations
    A&A 411, 321–327 (2003) Astronomy DOI: 10.1051/0004-6361:20031371 & c ESO 2003 Astrophysics Turbulent magnetic Prandtl number and magnetic diffusivity quenching from simulations T. A. Yousef1, A. Brandenburg2,andG.R¨udiger3 1 Department of Energy and Process Engineering, Norwegian University of Science and Technology, Kolbjørn Hejes vei 2B, 7491 Trondheim, Norway 2 NORDITA, Blegdamsvej 17, 2100 Copenhagen Ø, Denmark 3 Astrophysical Institute Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany Received 21 February 2003 / Accepted 25 August 2003 Abstract. Forced turbulence simulations are used to determine the turbulent kinematic viscosity, νt, from the decay rate of a large scale velocity field. Likewise, the turbulent magnetic diffusivity, ηt, is determined from the decay of a large scale magnetic field. In the kinematic regime, when the field is weak, the turbulent magnetic Prandtl number, νt/ηt, is about unity. When the field is nonhelical, ηt is quenched when magnetic and kinetic energies become comparable. For helical fields the quenching is stronger and can be described by a dynamical quenching formula. Key words. magnetohydrodynamics (MHD) – turbulence 1. Introduction however, this instability is suppressed (R¨udiger & Shalybkov 2002). On the other hand, the Reynolds number of the flow is The concept of turbulent diffusion is often invoked when mod- quite large (105 :::106) and the flow therefore most certainly eling large scale flows and magnetic fields in a turbulent turbulent. This led Noguchi et al. (2002) to invoke a turbulent medium. Turbulent magnetic diffusion is similar to turbulent kinematic viscosity, νt, but to retain the microscopic value of η. thermal diffusion which characterizes the turbulent exchange The resulting effective magnetic Prandtl number they used was 2 of patches of warm and cold gas.
    [Show full text]
  • Excitation and Control of a High-Speed Induction Generator
    Excitation and Control of a High-Speed Induction Generator by Steven Carl Englebretson S.B., Colorado School of Mines (Dec 2002) Submitted to the Department of Electrical Engineering and Computer Science in partial fulfillment of the requirements for the degree of Master of Science at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY September 2005 @ Massachusetts Institute of Technology, MMV. All rights reserved. A uth or .............................. .. ....... I/.. .. ................. Department of Electrical Engineering and Computer Science August 26, 2005 Certified by ....... .... .............. /1 James L. Kirtley Jr. Professor of Electrical Engineering I Thesis Supervisor Accepted by....................... .............. ............ Arthur C. Smith Chairman, Department Committee on Graduate Students MASSACHUSETTS INSIT=EU OF TECHNOLOGY BARKER MAR 2 8"2006 LIBRARIES - 11 - I Excitation and Control of a High-Speed Induction Generator by Steven Carl Englebretson Submitted to the Department of Electrical Engineering and Computer Science on August 26, 2005, in partial fulfillment of the requirements for the degree of Master of Science Abstract This project investigates the use of a high speed, squirrel cage induction generator and power converter for producing DC electrical power onboard ships and submarines. Potential advantages of high speed induction generators include smaller size and weight, increased durability, and decreased cost and maintenance. Unfortunately, induction generators require a "supply of reactive power" to run and suffer from variation in output voltage and frequency with any changes to the input reactive power excitation, mechanical drive speed, and load. A power converter can resolve some of these issues by circulating the changing reactive power demanded by the generator while simultaneously controlling the stator frequency to adjust the machine slip and manage the real output power.
    [Show full text]