-
Vorticity and magnetic dynamo from subsonic expansion waves II: Dependence on magnetic Prandtl number, forcing scale, cooling time
Authors:
Albert Elias-López,
Fabio del Sordo,
Daniele Viganò
Abstract:
The amplification of astrophysical magnetic fields takes place via dynamo instability in turbulent environments. The presence of vorticity is crucial for the dynamo to happen. However, the role of vorticity is not yet fully understood. This work is an extension of previous research on the effect of an irrotational subsonic forcing on a magnetized medium in the presence of rotation or a differentia…
▽ More
The amplification of astrophysical magnetic fields takes place via dynamo instability in turbulent environments. The presence of vorticity is crucial for the dynamo to happen. However, the role of vorticity is not yet fully understood. This work is an extension of previous research on the effect of an irrotational subsonic forcing on a magnetized medium in the presence of rotation or a differential velocity profile, aimed at exploring a wider parameter space in terms of Reynolds numbers, magnetic Prandtl number, forcing scale, cooling timescale in a Newtonian cooling. We study the effect of imposing either the acceleration or the velocity forcing function to be curl-free and evaluate the terms responsible for the evolution vorticity. We use Direct Numerical Simulations (DNS) to solve the fully compressible, resistive magnetohydrodynamical (MHD) equations with the Pencil Code. We study both isothermal and non-isothermal regimes and address the relative importance of different vorticity source terms.We report no small-scale dynamo for the models that do not include shear. We find a hydro instability, followed by a magnetic one, when a shearing velocity profile is applied. The vorticity production is found to be numerical in the purely irrotational case. Non-isothermality, rotation, shear or forcing in the form of a velocity curl-free, when included, contribute to increasing vorticity. Consistently with our previous study, we find that turbulence driven by subsonic expansion waves can amplify vorticity and magnetic field only in the presence of a background shearing profile. The presence of a cooling function make the instability happens on a shorter timescale. We estimate critical Reynolds and Magnetic Reynolds Numbes of 40 and 20, respectively.
△ Less
Submitted 18 April, 2024; v1 submitted 16 April, 2024;
originally announced April 2024.
-
Vorticity and magnetic dynamo from subsonic expansion waves
Authors:
Albert Elias-López,
Fabio Del Sordo,
Daniele Viganò
Abstract:
This work concentrates on the effect of an irrotational forcing on a magnetized flow in the presence of rotation, baroclinicity, shear, or a combination of them. By including magnetic field in the model we can evaluate the occurrence of dynamo on both small and large scales. We aim at finding what are the minimum ingredients needed to trigger a dynamo instability and what is the relation between d…
▽ More
This work concentrates on the effect of an irrotational forcing on a magnetized flow in the presence of rotation, baroclinicity, shear, or a combination of them. By including magnetic field in the model we can evaluate the occurrence of dynamo on both small and large scales. We aim at finding what are the minimum ingredients needed to trigger a dynamo instability and what is the relation between dynamo and the growth of vorticity.
We use the Pencil code to run resistive MHD direct numerical simulations. We report no dynamo in all cases where only rotation is included, regardless on the equation of state. Conversely, the inclusion of a background sinusoidal shearing profile leads to an hydrodynamic instability that produces an exponential growth of the vorticity at all scales, starting from small ones. This is know as vorticity dynamo. The onset of this instability occurs after a rather long temporal evolution of several thousand turbulent turnover times. The vorticity dynamo in turn drives an exponential growth of the magnetic field, first at small scales, then also at large one. The instability then saturates and the magnetic field approximately reaches equipartition with the turbulent kinetic energy. During the saturation phase we can observe a winding of the magnetic field in the direction of the shearing flow. By varying the intensity of the shear we see that the growth rates of this instability change. The inclusion of the baroclinic term delays the onset of the vorticity dynamo but leads to a more rapid growth.
We demonstrate how in the presence of shear, even a purely irrotational forcing amplifies the field to equipartition. At the same time, we confirm how this forcing alone does not lead to vorticity nor magnetic field growth, and this picture does not change in the presence of rotation or baroclinicity up to $256^3$ meshpoints.
△ Less
Submitted 26 July, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
QUBIC VII: The feedhorn-switch system of the technological demonstrator
Authors:
F. Cavaliere,
A. Mennella,
M. Zannoni,
P. Battaglia,
E. S. Battistelli,
D. Burke,
G. D'Alessandro,
P. de Bernardis,
M. De Petris,
C. Franceschet,
L. Grandsire,
J. -Ch. Hamilton,
B. Maffei,
E. Manzan,
S. Marnieros,
S. Masi,
C. O'Sullivan,
A. Passerini,
F. Pezzotta,
M. Piat,
A. Tartari,
S. A. Torchinsky,
D. Viganò,
F. Voisin,
P. Ade
, et al. (106 additional authors not shown)
Abstract:
We present the design, manufacturing and performance of the horn-switch system developed for the technological demonstrator of QUBIC (the $Q$\&$U$ Bolometric Interferometer for Cosmology). This system is constituted of 64 back-to-back dual-band (150\,GHz and 220\,GHz) corrugated feed-horns interspersed with mechanical switches used to select desired baselines during the instrument self-calibration…
▽ More
We present the design, manufacturing and performance of the horn-switch system developed for the technological demonstrator of QUBIC (the $Q$\&$U$ Bolometric Interferometer for Cosmology). This system is constituted of 64 back-to-back dual-band (150\,GHz and 220\,GHz) corrugated feed-horns interspersed with mechanical switches used to select desired baselines during the instrument self-calibration. We manufactured the horns in aluminum platelets milled by photo-chemical etching and mechanically tightened with screws. The switches are based on steel blades that open and close the wave-guide between the back-to-back horns and are operated by miniaturized electromagnets. We also show the current development status of the feedhorn-switch system for the QUBIC full instrument, based on an array of 400 horn-switch assemblies.
△ Less
Submitted 1 April, 2022; v1 submitted 28 August, 2020;
originally announced August 2020.
-
QUBIC V: Cryogenic system design and performance
Authors:
S. Masi,
E. S. Battistelli,
P. de Bernardis,
C. Chapron,
F. Columbro,
G. D'Alessandro,
M. De Petris,
L. Grandsire,
J. -Ch. Hamilton,
S. Marnieros,
L. Mele,
A. May,
A. Mennella,
C. O'Sullivan,
A. Paiella,
F. Piacentini,
M. Piat,
L. Piccirillo,
G. Presta,
A. Schillaci,
A. Tartari,
J. -P. Thermeau,
S. A. Torchinsky,
F. Voisin,
M. Zannoni
, et al. (104 additional authors not shown)
Abstract:
Current experiments aimed at measuring the polarization of the Cosmic Microwave Background (CMB) use cryogenic detector arrays and cold optical systems to boost the mapping speed of the sky survey. For these reasons, large volume cryogenic systems, with large optical windows, working continuously for years, are needed. Here we report on the cryogenic system of the QUBIC (Q and U Bolometric Interfe…
▽ More
Current experiments aimed at measuring the polarization of the Cosmic Microwave Background (CMB) use cryogenic detector arrays and cold optical systems to boost the mapping speed of the sky survey. For these reasons, large volume cryogenic systems, with large optical windows, working continuously for years, are needed. Here we report on the cryogenic system of the QUBIC (Q and U Bolometric Interferometer for Cosmology) experiment: we describe its design, fabrication, experimental optimization and validation in the Technological Demonstrator configuration. The QUBIC cryogenic system is based on a large volume cryostat, using two pulse-tube refrigerators to cool at ~3K a large (~1 m^3) volume, heavy (~165kg) instrument, including the cryogenic polarization modulator, the corrugated feedhorns array, and the lower temperature stages; a 4He evaporator cooling at ~1K the interferometer beam combiner; a 3He evaporator cooling at ~0.3K the focal-plane detector arrays. The cryogenic system has been tested and validated for more than 6 months of continuous operation. The detector arrays have reached a stable operating temperature of 0.33K, while the polarization modulator has been operated from a ~10K base temperature. The system has been tilted to cover the boresight elevation range 20 deg -90 deg without significant temperature variations. The instrument is now ready for deployment to the high Argentinean Andes.
△ Less
Submitted 25 August, 2021; v1 submitted 24 August, 2020;
originally announced August 2020.
-
Wideband 67-116 GHz receiver development for ALMA Band 2
Authors:
P. Yagoubov,
T. Mroczkowski,
V. Belitsky,
D. Cuadrado-Calle,
F. Cuttaia,
G. A. Fuller,
J. -D. Gallego,
A. Gonzalez,
K. Kaneko,
P. Mena,
R. Molina,
R. Nesti,
V. Tapia,
F. Villa,
M. Beltran,
F. Cavaliere,
J. Ceru,
G. E. Chesmore,
K. Coughlin,
C. De Breuck,
M. Fredrixon,
D. George,
H. Gibson,
J. Golec,
A. Josaitis
, et al. (21 additional authors not shown)
Abstract:
ALMA has been operating since 2011, but has not yet been populated with the full suite of intended frequency bands. In particular, ALMA Band 2 (67-90 GHz) is the final band in the original ALMA band definition to be approved for production. We aim to produce a wideband, tuneable, sideband-separating receiver with 28 GHz of instantaneous bandwidth per polarisation operating in the sky frequency ran…
▽ More
ALMA has been operating since 2011, but has not yet been populated with the full suite of intended frequency bands. In particular, ALMA Band 2 (67-90 GHz) is the final band in the original ALMA band definition to be approved for production. We aim to produce a wideband, tuneable, sideband-separating receiver with 28 GHz of instantaneous bandwidth per polarisation operating in the sky frequency range 67-116 GHz. Our design anticipates new ALMA requirements following the recommendations in the 2030 ALMA Development Roadmap. The cryogenic cartridge is designed to be compatible with the ALMA Band 2 cartridge slot, where the coldest components -- the feedhorns, orthomode transducers, and cryogenic low noise amplifiers -- operate at a temperature of 15 K. We use multiple simulation methods and tools to optimise our designs for both the passive optics and the active components. The cryogenic cartridge interfaces with a room temperature cartridge hosting the local oscillator (LO) and the downconverter module. This warm cartridge is largely based on GaAs semiconductor technology and is optimised to match the cryogenic receiver bandwidth with the required instantaneous LO tuning range. Our collaboration has designed, fabricated, and tested multiple technical solutions for each of the components, producing a state-of-the-art receiver covering the full ALMA Band 2 & 3 atmospheric window. The receiver is suitable for deployment on ALMA in the coming years, and is capable of dual-polarisation, sideband-separating observations in intermediate frequency bands spanning 4-18 GHz, for a total of 28 GHz on-sky bandwidth per polarisation channel. We conclude that the 67-116 GHz wideband implementation for ALMA Band 2 is now feasible, and this receiver is a compelling instrumental upgrade that will enhance observational capabilities and scientific reach.
△ Less
Submitted 20 February, 2020; v1 submitted 20 December, 2019;
originally announced December 2019.
-
Extension of the sub-grid-scale gradient model for compressible magnetohydrodynamics turbulent instabilities
Authors:
Daniele Viganò,
Ricard Aguilera-Miret,
Carlos Palenzuela
Abstract:
Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evol…
▽ More
Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evolution equation of a compressible ideal MHD fluid with a generic equation of state. We assess the model by considering box simulations of the turbulence triggered across a shear layer by the Kelvin-Helmholtz instability, testing cases where the small-scale dynamics cannot be fully captured by the resolution considered, such that the efficiency of the simulated dynamo effect depends on the resolution employed. This lack of numerical convergence is actually a currently common issue in several astrophysical problems, where the integral and fastest-growing-instability scales are too far apart to be fully covered numerically. We perform a-priori and a-posteriori tests of the extended gradient model. In the former, we find that, for many different initial conditions and resolutions, the gradient model outperforms other commonly used models in terms of correlation with the residuals coming from the filtering of a high-resolution run. In the second test, we show how a low-resolution run with the gradient model is able to quantitatively reproduce the evolution of the magnetic energy (the integrated value and the spectral distribution) coming from higher-resolution runs. This extension is the first step towards the implementation in relativistic magnetohydrodynamics.
△ Less
Submitted 12 September, 2019; v1 submitted 8 April, 2019;
originally announced April 2019.
-
A Simflowny-based high-performance 3D code for the generalized induction equation
Authors:
Daniele Viganò,
David Martínez-Gómez,
José A. Pons,
Carlos Palenzuela,
Federico Carrasco,
Borja Miñano,
Antoni Arbona,
Carles Bona,
Joan Massó
Abstract:
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage…
▽ More
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em SAMRAI} architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources.
△ Less
Submitted 11 December, 2018; v1 submitted 20 November, 2018;
originally announced November 2018.
-
Thermal architecture for the QUBIC cryogenic receiver
Authors:
A. J. May,
C. Chapron,
G. Coppi,
G. D'Alessandro,
P. de Bernardis,
S. Masi,
S. Melhuish,
M. Piat,
L. Piccirillo,
A. Schillaci,
J. -P. Thermeau,
P. Ade,
G. Amico,
D. Auguste,
J. Aumont,
S. Banfi,
G. Barbara,
P. Battaglia,
E. Battistelli,
A. Bau,
B. Belier,
D. Bennett,
L. Berge,
J. -Ph. Bernard,
M. Bersanelli
, et al. (105 additional authors not shown)
Abstract:
QUBIC, the QU Bolometric Interferometer for Cosmology, is a novel forthcoming instrument to measure the B-mode polarization anisotropy of the Cosmic Microwave Background. The detection of the B-mode signal will be extremely challenging; QUBIC has been designed to address this with a novel approach, namely bolometric interferometry. The receiver cryostat is exceptionally large and cools complex opt…
▽ More
QUBIC, the QU Bolometric Interferometer for Cosmology, is a novel forthcoming instrument to measure the B-mode polarization anisotropy of the Cosmic Microwave Background. The detection of the B-mode signal will be extremely challenging; QUBIC has been designed to address this with a novel approach, namely bolometric interferometry. The receiver cryostat is exceptionally large and cools complex optical and detector stages to 40 K, 4 K, 1 K and 350 mK using two pulse tube coolers, a novel 4He sorption cooler and a double-stage 3He/4He sorption cooler. We discuss the thermal and mechanical design of the cryostat, modelling and thermal analysis, and laboratory cryogenic testing.
△ Less
Submitted 6 November, 2018;
originally announced November 2018.
-
A Simflowny-based finite-difference code for high-performance computing in Relativity
Authors:
Carlos Palenzuela,
Borja Miñano,
Daniele Viganò,
Antoni Arbona,
Carles Bona-Casas,
Andreu Rigo,
Miguel Bezares,
Carles Bona,
Joan Massó
Abstract:
The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processor…
▽ More
The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processors and advanced computational strategies are adopted. To achieve this goal we have developed Simflowny, an open platform for scientific dynamical models which automatically generates parallel code for different simulation frameworks, allowing the use of HPC infrastructures to non-specialist scientists. One of these frameworks is SAMRAI, a mature patch-based structured adaptive mesh refinement infrastructure, capable of reaching exascale in some specific problems. Here we present the numerical techniques that we have implemented on this framework by using Simflowny in order to perform fast, efficient, accurate and highly-scalable simulations. These techniques involve high-order schemes for smooth and non-smooth solutions, Adaptive Mesh Refinement with arbitrary resolution ratios and an optimal strategy for the sub-cycling in time. We validate the automatically generated codes for the SAMRAI infrastructure with some simple test examples (i.e., wave equation and Newtonian MHD) and finally with the Einstein equations.
△ Less
Submitted 11 June, 2018;
originally announced June 2018.