We review and extend the theory of ideal, homogeneous, incompressible, magnetohydrodynamic (MHD) turbulence. The theory contains a solution to the ‘dynamo problem’, i.e., the problem of determining ...how a planetary or stellar body produces a global dipole magnetic field. We extend the theory to the case of ideal MHD turbulence with a mean magnetic field that is aligned with a rotation axis. The existing theory is also extended by developing the thermodynamics of ideal MHD turbulence based on entropy. A mathematical model is created by Fourier transforming the MHD equations and dynamical variables, resulting in a dynamical system consisting of the independent Fourier coefficients of the velocity and magnetic fields. This dynamical system has a large but finite-dimensional phase space in which the phase flow is divergenceless in the ideal case. There may be several constants of the motion, in addition to energy, which depend on the presence, or lack thereof, of a mean magnetic field or system rotation or both imposed on the magnetofluid; this leads to five different cases of MHD turbulence that must be considered. The constants of the motion (ideal invariants)—the most important being energy and magnetic helicity—are used to construct canonical probability densities and partition functions that enable ensemble predictions to be made. These predictions are compared with time averages from numerical simulations to test whether or not the system is ergodic. In the cases most pertinent to planets and stars, nonergodicity is observed at the largest length-scales and occurs when the components of the dipole field become quasi-stationary and dipole energy is directly proportional to magnetic helicity. This nonergodicity is evident in the thermodynamics, while dipole alignment with a rotation axis may be seen as the result of dynamical symmetry breaking, i.e., ‘broken ergodicity’. The relevance of ideal theoretical results to real (forced, dissipative) MHD turbulence is shown through numerical simulation. Again, an important result is a statistical solution of the ‘dynamo problem’.
Transition of ideal, homogeneous, incompressible, magnetohydrodynamic (MHD) turbulence to near-equilibrium from non-equilibrium initial conditions is examined through new long-time numerical ...simulations on a 1283 periodic grid. Here, we neglect dissipation because we are primarily concerned with behavior at the largest scale which has been shown to be essentially the same for ideal and real (forced and dissipative) MHD turbulence. A Fourier spectral transform method is used to numerically integrate the dynamical equations forward in time and results from six computer runs are presented with various combinations of imposed rotation and mean magnetic field. There are five separate cases of ideal, homogeneous, incompressible, MHD turbulence: Case I, with no rotation or mean field; Case II, where only rotation is imposed; Case III, which has only a mean magnetic field; Case IV, where rotation vector and mean magnetic field direction are aligned; and Case V, which has nonaligned rotation vector and mean field directions. Dynamic coefficients are predicted by statistical mechanics to be zero-mean random variables, but largest-scale coherent magnetic structures emerge in all cases during transition; this implies dynamo action is inherent in ideal MHD turbulence. These coherent structures are expected to occur in Cases I, II and IV, but not in Cases III and V; future studies will determine whether they persist.
We continue our study of the transition of ideal, homogeneous, incompressible, magnetohydrodynamic (MHD) turbulence from non-equilibrium initial conditions to equilibrium using long-time numerical ...simulations on a 1283 periodic grid. A Fourier spectral transform method is used to numerically integrate the dynamical equations forward in time. The six runs that previously went to near equilibrium are here extended into equilibrium. As before, we neglect dissipation as we are primarily concerned with behavior at the largest scale where this behavior has been shown to be essentially the same for ideal and real (forced and dissipative) MHD turbulence. These six runs have various combinations of imposed rotation and mean magnetic field and represent the five cases of ideal, homogeneous, incompressible, and MHD turbulence: Case I (Run 1), with no rotation or mean field; Case II (Runs 2a and 2b), where only rotation is imposed; Case III (Run 3), which has only a mean magnetic field; Case IV (Run 4), where rotation vector and mean magnetic field direction are aligned; and Case V (Run 5), which has non-aligned rotation vector and mean field directions. Statistical mechanics predicts that dynamic Fourier coefficients are zero-mean random variables, but largest-scale coherent magnetic structures emerge and manifest themselves as Fourier coefficients with very large, quasi-steady, mean values compared to their standard deviations, i.e., there is ‘broken ergodicity.’ These magnetic coherent structures appeared in all cases during transition to near equilibrium. Here, we report that, as the runs were continued, these coherent structures remained quasi-steady and energetic only in Cases I and II, while Case IV maintained its coherent structure but at comparatively low energy. The coherent structures that appeared in transition in Cases III and V were seen to collapse as their associated runs extended into equilibrium. The creation of largest-scale, coherent magnetic structure appears to be a dynamo process inherent in ideal MHD turbulence, particularly in Cases I and II, i.e., those cases most pertinent to planets and stars. Furthermore, the statistical theory of ideal MHD turbulence has proven to apply at the largest scale, even when dissipation and forcing are included. This, along with the discovery and explanation of dynamically broken ergodicity, is essentially a solution to the ‘dynamo problem’.
Solar magnetism is believed to originate through dynamo action in the tachocline. Statistical mechanics, in turn, tells us that dynamo action is an inherent property of magnetohydrodynamic (MHD) ...turbulence, depending essentially on magnetic helicity. Here, we model the tachocline as a rotating, thin spherical shell containing MHD turbulence. Using this model, we find an expression for the entropy and from this develop the thermodynamics of MHD turbulence. This allows us to introduce the macroscopic parameters that affect magnetic self-organization and dynamo action, parameters that include magnetic helicity, as well as tachocline thickness and turbulent energy.
We find the analytical form of inertial waves in an incompressible, rotating fluid constrained by concentric inner and outer spherical surfaces with homogeneous boundary conditions on the normal ...components of velocity and vorticity. These fields are represented by Galerkin expansions whose basis consists of toroidal and poloidal vector functions, i.e., products and curls of products of spherical Bessel functions and vector spherical harmonics. These vector basis functions also satisfy the Helmholtz equation and this has the benefit of providing each basis function with a well-defined wavenumber. Eigenmodes and associated eigenfrequencies are determined for both the ideal and dissipative cases. These eigenmodes are formed from linear combinations of the Galerkin expansion basis functions. The system is truncated to numerically study inertial wave structure, varying the number of eigenmodes. The largest system considered in detail is a 25 eigenmode system and a graphical depiction is presented of the five lowest dissipation eigenmodes, all of which are non-oscillatory. These results may be useful in understanding data produced by numerical simulations of fluid and magnetofluid turbulence in a spherical shell that use a Galerkin, toroidal–poloidal basis as well as qualitative features of liquids confined by a spherical shell.
The Earth’s magnetic field is measured on and above the crust, while the turbulent dynamo in the outer core produces magnetic field values at the core–mantle boundary (CMB). The connection between ...the two sets of values is usually assumed to be independent of the electrical conductivity in the mantle. However, the turbulent magnetofluid in the Earth’s outer core produces a time-varying magnetic field that must induce currents in the lower mantle as it emerges, since the mantle is observed to be electrically conductive. Here, we develop a model to assess the possible effects of mantle electrical conductivity on the magnetic field values at the CMB. This model uses a new method for mapping the geomagnetic field from the Earth’s surface to the CMB. Since numerical and theoretical results suggest that the turbulent magnetic field in the outer core as it approaches the CMB is mostly parallel to this boundary, we assume that this property exists and set the normal component of the model magnetic field to zero at the CMB. This leads to a modification of the Mauersberger–Lowes spectrum at the CMB so that it is no longer flat, i.e., the modified spectrum depends on mantle conductance. We examined several cases in which mantle conductance ranges from low to high in order to gauge how CMB magnetic field strength and mantle ohmic heat generation may vary.
A technique for the kinetic simulation of plasma flow in strong external magnetic fields was developed which captures the compression and expansion of plasma bound to a magnetic flux tube as well as ...forces on magnetized particles within the flux tube. This quasi-one-dimensional (Q1D) method resolves a single spatial dimension while modeling two-dimensional effects. The implementation of this method in a Particle-In-Cell (PIC) code was verified with newly formulated test cases which include two-particle motion and particle dynamics in a magnetic mirror. Results from the Q1D method and fully two dimensional simulations were compared and error analyses performed verifying that the Q1D model reproduces the fully 2D results in the correct regimes. The Q1D method was found to be valid when the hybrid Larmor radius was less than 10% of the magnetic field scale length for magnetic field guided plasma expansions and less than 1% of the magnetic field scale length for a plasma in a converging–diverging magnetic field. The simple and general Q1D method can readily be incorporated in standard 1D PIC codes to capture multi-dimensional effects for plasma flow along magnetic fields in parameter spaces currently inaccessible by fully kinetic methods.
Large Area Picosecond Photodetectors (LAPPD) are a new generation of microchannel plate photomultipliers being manufactured by Incom. These devices feature large sensitive area of 350 cm2, high ...quantum efficiency (∼20%), and tens of picosecond single photon timing resolution. LAPPDs use an anode structure with 28 striplines, allowing for spatial resolution of 1÷3 mm while minimizing the number of readout channels. In this report, we present our design of integrated readout electronics for the LAPPD . These electronics read out all 2×28 channels of a stripline-anode LAPPD . Waveform sampling at up to 5GSPS is performed by 8 DRS4 switched-capacitor array ASICs. All DRS4 channels are digitized in parallel with two 32-channel ADCs. An on-board FPGA coordinates digitization and readout of waveforms, and could further be expanded to include some waveform processing. Data packages built in the FPGA are sent to a data acquisition system (DAQ) via optical fiber, with a baseline Gigabit Ethernet interface implemented entirely on the FPGA . The electronics is designed to accommodate different triggering options: self-triggering using DRS4 transparent mode and external triggering, making event control very flexible. Further flexibility is enhanced with embedded software for an on-FPGA soft-core processor, as well as DAQ readout and control software. The device is plug and play with any existing IP network. An open-source ecosystem is being developed to provide full control of the device operation and an easy way to integrate it to any environment. In the report we describe the status of the electronics development, its firmware and readout software. Also the results of the first tests of the electronics with the LAPPD devices will be presented.
The cross section of the process e + e − → π + π − has been measured in the center-of-mass energy range from 0.32 to 1.2 GeV with the CMD-3 detector at the electron-positron collider VEPP-2000. The ...measurement is based on an integrated luminosity of about 88 pb − 1 , of which 62 pb − 1 represent a complete dataset collected by CMD-3 at center-of-mass energies below 1 GeV. In the dominant region near the ρ resonance a systematic uncertainty of 0.7% was achieved. The implications of the presented results for the evaluation of the hadronic contribution to the anomalous magnetic moment of the muon are discussed. Published by the American Physical Society 2024
We present theoretical and computational results in magnetohydrodynamic turbulence that we feel are essential to understanding the geodynamo. These results are based on a mathematical model that ...focuses on magnetohydrodynamic (MHD) turbulence, but ignores compressibility and thermal effects, as well as imposing model-dependent boundary conditions. A principal finding is that when a turbulent magnetofluid is in quasi-equilibrium, the magnetic energy in the internal dipole component is equal to the magnetic helicity multiplied by the dipole wavenumber. In the case of the Earth, measurement of the exterior magnetic field gives us, through boundary conditions, the internal poloidal magnetic field. The connection between magnetic helicity and dipole field in the liquid core then gives us the toroidal part of the internal dipole field and a model value of 3 mT for the average core dipole magnetic field. Here, we present the theoretical analysis and numerical simulations that lead to these conclusions. We also test an earlier assertion that differential oblateness may be related to dipole alignment, and while there is an effect, rotation appears to be far more important. In addition, the relationship between dipole quasi-stationarity, broken ergodicity and broken symmetry is clarified. Lastly, we discuss how inertial waves in a rotating magnetofluid can affect dipole alignment.
PDF
