Showing 1–50 of 53 results for author: Hansen, C

Error Field Predictability and Consequences for ITER

Authors: Matthew Pharr, Nikolas Logan, Carlos Paz-Soldan, Jong-Kyu Park, Chistopher Hansen

Abstract: ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils… ▽ More ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scaling laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set. △ Less

Submitted 9 June, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Comments: 11 pages, 13 Figures. Submitted to Nuclear Fusion

MANTA: A Negative-Triangularity NASEM-Compliant Fusion Pilot Plant

Authors: MANTA Collaboration, G. Rutherford, H. S. Wilson, A. Saltzman, D. Arnold, J. L. Ball, S. Benjamin, R. Bielajew, N. de Boucaud, M. Calvo-Carrera, R. Chandra, H. Choudhury, C. Cummings, L. Corsaro, N. DaSilva, R. Diab, A. R. Devitre, S. Ferry, S. J. Frank, C. J. Hansen, J. Jerkins, J. D. Johnson, P. Lunia, J. van de Lindt, S. Mackie , et al. (16 additional authors not shown)

Abstract: The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicin… ▽ More The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Shock Wave Formation in Radiative Plasmas

Authors: F. Garcia-Rubio, V. Tranchant, E. C. Hansen, R. Tabassum, A. Reyes, H. U. Rahman, P. Ney, E. Ruskov, P. Tzeferacos

Abstract: The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thin and thick regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynam… ▽ More The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thin and thick regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thin and thick solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive test-case has been designed to benchmark radiative hydrodynamic codes. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Approved in Phys. Rev. E

Geometric Burn Control For Tokamaks

Authors: J. F. Parisi, J. W. Berkery, A. Sladkomedova, S. Guizzo, M. R. Hardman, J. R. Ball, A. O. Nelson, S. M. Kaye, M. Anastopoulos-Tzanis, S. A. M. McNamara, J. Dominski, S. Janhunen, M. Romanelli, D. Dickinson, A. Diallo, A. Dnestrovskii, W. Guttenfelder, C. Hansen, O. Myatra, H. R. Wilson

Abstract: A new burn control scheme for tokamaks is described where the total fusion power is controlled by adjusting the plasma volume fraction that is packed into power dense regions. In an example spherical tokamak burning plasma, by modifying the plasma edge squareness the total fusion power is doubled at almost constant total plasma volume and fusion power density. Therefore, increased plasma squarenes… ▽ More A new burn control scheme for tokamaks is described where the total fusion power is controlled by adjusting the plasma volume fraction that is packed into power dense regions. In an example spherical tokamak burning plasma, by modifying the plasma edge squareness the total fusion power is doubled at almost constant total plasma volume and fusion power density. Therefore, increased plasma squareness could be extremely beneficial to a fusion reactor and squareness control could be desirable for power load balancing. Experiments have observed the impact of increased edge squareness on modified core plasma volume, highlighting the practical relevance of this approach. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 7 pages, 7 figures

Generative deep learning-enabled ultra-large field-of-view lens-free imaging

Authors: Ronald B. Liu, Zhe Liu, Max G. A. Wolf, Krishna P. Purohit, Gregor Fritz, Yi Feng, Carsten G. Hansen, Pierre O. Bagnaninchi, Xavier Casadevall i Solvas, Yunjie Yang

Abstract: Advancements in high-throughput biomedical applications necessitate real-time, large field-of-view (FOV) imaging capabilities. Conventional lens-free imaging (LFI) systems, while addressing the limitations of physical lenses, have been constrained by dynamic, hard-to-model optical fields, resulting in a limited one-shot FOV of approximately 20 $mm^2$. This restriction has been a major bottleneck i… ▽ More Advancements in high-throughput biomedical applications necessitate real-time, large field-of-view (FOV) imaging capabilities. Conventional lens-free imaging (LFI) systems, while addressing the limitations of physical lenses, have been constrained by dynamic, hard-to-model optical fields, resulting in a limited one-shot FOV of approximately 20 $mm^2$. This restriction has been a major bottleneck in applications like live-cell imaging and automation of microfluidic systems for biomedical research. Here, we present a deep-learning(DL)-based imaging framework - GenLFI - leveraging generative artificial intelligence (AI) for holographic image reconstruction. We demonstrate that GenLFI can achieve a real-time FOV over 550 $mm^2$, surpassing the current LFI system by more than 20-fold, and even larger than the world's largest confocal microscope by 1.76 times. The resolution is at the sub-pixel level of 5.52 $μm$, without the need for a shifting light source. The unsupervised learning-based reconstruction does not require optical field modeling, making imaging dynamic 3D samples (e.g., droplet-based microfluidics and 3D cell models) in complex optical fields possible. This GenLFI framework unlocks the potential of LFI systems, offering a robust tool to tackle new frontiers in high-throughput biomedical applications such as drug discovery. △ Less

Submitted 22 March, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Local stability guarantees for data-driven quadratically nonlinear models

Authors: Mai Peng, Alan Kaptanoglu, Chris Hansen, Jacob Stevens-Haas, Krithika Manohar, Steven L. Brunton

Abstract: The Navier Stokes equations (NSEs) are partial differential equations (PDEs) to describe the nonlinear convective motion of fluids and they are computationally expensive to simulate because of their high nonlinearity and variables being fully coupled. Reduced-order models (ROMs) are simpler models for evolving the flows by capturing only the dominant behaviors of a system and can be used to design… ▽ More The Navier Stokes equations (NSEs) are partial differential equations (PDEs) to describe the nonlinear convective motion of fluids and they are computationally expensive to simulate because of their high nonlinearity and variables being fully coupled. Reduced-order models (ROMs) are simpler models for evolving the flows by capturing only the dominant behaviors of a system and can be used to design controllers for high-dimensional systems. However it is challenging to guarantee the stability of these models either globally or locally. Ensuring the stability of ROMs can improve the interpretability of the behavior of the dynamics and help develop effective system control strategies. For quadratically nonlinear systems that represent many fluid flows, the Schlegel and Noack trapping theorem (JFM, 2015) can be used to check if ROMs are globally stable (long-term bounded). This theorem was subsequently incorporated into system identification techniques that determine models directly from data. In this work, we relax the quadratically energy-preserving constraints in this theorem, and then promote local stability in data-driven models of quadratically nonlinear dynamics. First, we prove a theorem outlining sufficient conditions to ensure local stability in linear-quadratic systems and provide an estimate of the stability radius. Second, we incorporate this theorem into system identification methods and produce a-priori locally stable data-driven models. Several examples are presented to demonstrate the effectiveness and accuracy of the proposed algorithm. △ Less

Submitted 28 August, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Assessment of vertical stability for negative triangularity pilot plants

Authors: S. Guizzo, A. O. Nelson, C. Hansen, F. Logak, C. Paz-Soldan

Abstract: Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When couple… ▽ More Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Low latency optical-based mode tracking with machine learning deployed on FPGAs on a tokamak

Authors: Yumou Wei, Ryan F. Forelli, Chris Hansen, Jeffrey P. Levesque, Nhan Tran, Joshua C. Agar, Giuseppe Di Guglielmo, Michael E. Mauel, Gerald A. Navratil

Abstract: Active feedback control in magnetic confinement fusion devices is desirable to mitigate plasma instabilities and enable robust operation. Optical high-speed cameras provide a powerful, non-invasive diagnostic and can be suitable for these applications. In this study, we process fast camera data, at rates exceeding 100kfps, on $\textit{in situ}$ Field Programmable Gate Array (FPGA) hardware to trac… ▽ More Active feedback control in magnetic confinement fusion devices is desirable to mitigate plasma instabilities and enable robust operation. Optical high-speed cameras provide a powerful, non-invasive diagnostic and can be suitable for these applications. In this study, we process fast camera data, at rates exceeding 100kfps, on $\textit{in situ}$ Field Programmable Gate Array (FPGA) hardware to track magnetohydrodynamic (MHD) mode evolution and generate control signals in real-time. Our system utilizes a convolutional neural network (CNN) model which predicts the $n$=1 MHD mode amplitude and phase using camera images with better accuracy than other tested non-deep-learning-based methods. By implementing this model directly within the standard FPGA readout hardware of the high-speed camera diagnostic, our mode tracking system achieves a total trigger-to-output latency of 17.6$μ$s and a throughput of up to 120kfps. This study at the High Beta Tokamak-Extended Pulse (HBT-EP) experiment demonstrates an FPGA-based high-speed camera data acquisition and processing system, enabling application in real-time machine-learning-based tokamak diagnostic and control as well as potential applications in other scientific domains. △ Less

Submitted 9 July, 2024; v1 submitted 30 November, 2023; originally announced December 2023.

Comments: This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in Rev. Sci. Instrum. 95, 073509 (2024) and may be found at https://doi.org/10.1063/5.0190354

Report number: FERMILAB-PUB-23-655-CSAID

Journal ref: Rev. Sci. Instrum. 95, 073509 (2024)

Monitoring the evolution of relative product populations at early times during a photochemical reaction

Authors: Joao Pedro Figueira Nunes, Lea Maria Ibele, Shashank Pathak, Andrew R. Attar, Surjendu Bhattacharyya, Rebecca Boll, Kurtis Borne, Martin Centurion, Benjamin Erk, Ming-Fu Lin, Ruaridh J. G. Forbes, Nate Goff, Christopher S. Hansen, Matthias Hoffmann, David M. P. Holland, Rebecca A. Ingle, Duan Luo, Sri Bhavya Muvva, Alex Reid, Arnaud Rouzée, Artem Rudenko, Sajib Kumar Saha, Xiaozhe Shen, Anbu Selvam Venkatachalam, Xijie Wang , et al. (9 additional authors not shown)

Abstract: Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging re… ▽ More Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (~50%) yield of an episulfide isomer containing a strained 3-membered ring within ~1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Feasibility and Performance of the Staged Z-Pinch: A One-dimensional Study with FLASH and MACH2

Authors: Edward C. Hansen, Fernando Garcia-Rubio, Marissa B. P. Adams, Milad Fatenejad, Kasper Moczulski, Paul Ney, Hafiz U. Rahman, Adam C. Reyes, Emil Ruskov, Victor Tranchant, Petros Tzeferacos

Abstract: Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to… ▽ More Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (> 300), but achieve fuel ion temperatures of many keV. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept. △ Less

Submitted 3 April, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Authors: C. Hansen, I. G. Stewart, D. Burgess, M. Pharr, S. Guizzo, F. Logak, A. O. Nelson, C. Paz-Soldan

Abstract: In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like description… ▽ More In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial orders for analytic and cross-code test cases. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Journal ref: Computer Physics Communications 298, 109111 (2024)

Quantum-informed plasmonics for strong coupling: the role of electron spill-out

Authors: Ida Juliane Bundgaard, Christian Nicolaisen Hansen, P. Elli Stamatopoulou, Christos Tserkezis

Abstract: The effect of nonlocality on the optical response of metals lies at the forefront of research in nanoscale physics and, in particular, quantum plasmonics. In alkali metals, nonlocality manifests predominantly as electron density spill-out at the metal boundary, and as surface-enabled Landau damping. For an accurate description of plasmonic modes, these effects need be taken into account in the the… ▽ More The effect of nonlocality on the optical response of metals lies at the forefront of research in nanoscale physics and, in particular, quantum plasmonics. In alkali metals, nonlocality manifests predominantly as electron density spill-out at the metal boundary, and as surface-enabled Landau damping. For an accurate description of plasmonic modes, these effects need be taken into account in the theoretical modelling of the material. The resulting modal frequency shifts and broadening become particularly relevant when dealing with the strong interaction between plasmons and excitons, where hybrid modes emerge and the way they are affected can reflect modifications of the coupling strength. Both nonlocal phenomena can be incorporated in the classical local theory by applying a surface-response formalism embodied by the Feibelman parameters. Here, we implement surface-response corrections in Mie theory to study the optical response of spherical plasmonic--excitonic composites in core--shell configurations. We investigate sodium, a jellium metal dominated by spill-out, for which it has been anticipated that nonlocal corrections should lead to an observable change in the coupling strength, appearing as a modification of the width of the mode splitting. We show that, contrary to expectations, the influence of nonlocality on the anticrossing is minimal, thus validating the accuracy of the local response approximation in strong-coupling photonics. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Design of Passive and Structural Conductors for Tokamaks Using Thin-Wall Eddy Current Modeling

Authors: A. F. Battey, C. Hansen, D. Garnier, D. Weisberg, C. Paz-Soldan, R. Sweeney, R. A. Tinguely, A. J. Creely

Abstract: A new three-dimensional electromagnetic modeling tool ThinCurr has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vess… ▽ More A new three-dimensional electromagnetic modeling tool ThinCurr has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vessel and other conducting structural elements with realistic material resistivities. This model was leveraged to support the design of a passive runaway electron mitigation coil (REMC), studying the effect of various design parameters, including coil resistivity, current quench duration, and plasma vertical position, on the effectiveness of the coil. The REMC is a non-axisymmetric coil designed to passively drive large non-axisymmetric fields during the plasma disruption thereby destroying flux surfaces and deconfining RE seed populations. These studies indicate that current designs should apply substantial 3D fields at the plasma surface during future plasma current disruptions as well as highlight the importance of having the REMC conductors away from the machine midplane in order to ensure they are robust to off-normal disruption scenarios. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 33 pages and 14 figures. to be published in Nuclear Fusion

Extension of the law of the wall exploiting weak similarity of velocity fluctuations in turbulent channels

Authors: Christoffer Hansen, Jens N. Sørensen, Xiang I. A. Yang, Mahdi Abkar

Abstract: This paper explores the similarity of the streamwise velocity fluctuations in a channel. In the analysis, we employ a one-dimensional scalar variant of the proper orthogonal decomposition (POD). This approach naturally motivates the introduction of two different levels of similarity which we will refer to as strong and weak similarity. Strong similarity requires that the two-point correlation, and… ▽ More This paper explores the similarity of the streamwise velocity fluctuations in a channel. In the analysis, we employ a one-dimensional scalar variant of the proper orthogonal decomposition (POD). This approach naturally motivates the introduction of two different levels of similarity which we will refer to as strong and weak similarity. Strong similarity requires that the two-point correlation, and thus, all POD modes, show Reynolds number similarity, while weak similarity only requires that the first few POD modes show similarity. As POD concerns information at more than one location, these similarities are more general than various similarities found in the literature concerning single-point flow statistics. We examine flows at $Re_τ=$180, 540, 1000, and 5200. Strong similarity is observed in the viscous layer and the wake region, and weak similarity is found in both the viscous wall region and the outer part of the logarithmic layer. The presence of weak similarity suggests the existence of an extension to the law of the wall (LoW). We propose such an extension based on the results from the one-dimensional POD analysis. The usefulness of the LoW extension is then assessed by comparing flow reconstructions according to the conventional equilibrium LoW and the extended LoW. We show that the extended LoW provides accurate flow reconstructions in the wall layer, capturing fine-scale motions that are entirely missed by the equilibrium LoW. △ Less

Submitted 29 December, 2023; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 17 pages, 16 figures

Optimization of Graded Metamaterials for Control of Energy Transmission Using a Genetic Algorithm

Authors: Joshua Morris, Weidi Wang, Thomas Plaisted, Christopher J. Hansen, Alireza V. Amirkhizi

Abstract: Optimization of functionally graded metamaterial arrays with a high dimensional and continuous geometric design space is cumbersome and could be accelerated via machine learning tools. Mechanical metamaterials can manipulate acoustic or ultrasonic waves by introducing large dispersive and attenuative effects near their natural frequency. In this work functionally graded structures are designed and… ▽ More Optimization of functionally graded metamaterial arrays with a high dimensional and continuous geometric design space is cumbersome and could be accelerated via machine learning tools. Mechanical metamaterials can manipulate acoustic or ultrasonic waves by introducing large dispersive and attenuative effects near their natural frequency. In this work functionally graded structures are designed and optimized to combine the energy attenuation performance of a number of unit cells with varying frequency responses and to reduce the interlayer mismatch effects. Optimization through genetic algorithm avoids the many local minima related to high dimensionality of the design space but requires many iterations. A reduced order model (ROM) is applied that can reproduce the transmission response that is traditionally calculated with FEM, in a fraction of the time. Pairing GA and the ROM together, an array of 6 unit cells (with a total of 18 independent geometric design variables) is optimized to have stop bands with extended width and sharper boundaries. Symmetric functionally graded structures are determined to be optimal geometric configurations. Measured 3D printed features are projected onto the ROM solutions to quantify the effect of printing uncertainty on array performance. Repeatability error of $\pm~20$ $μ$m is determined to reduce the mean depth of the transmission stop band by a factor of $10^2$ and introduce small shifts in center frequency and band width. Proposed methods to improve the resolution of accessible points in the ROM space, reduce sensitivity to geometric uncertainty, and add design freedom include introducing out-of-plane perforations and varying constituent materials using tunable filled resin systems. △ Less

Submitted 18 June, 2023; originally announced June 2023.

Effect of Moisture Absorption on Curing of Wind Blades during Repair

Authors: Sagar P. Shah, Michael N. Olaya, Evgenia Plaka, Joseph McDonald, Christopher J. Hansen, Marianna Maiarù

Abstract: Efficient structural repair of wind turbine blades is essential to limiting global warming and reducing the Levelized Cost of Energy (LCOE). Repairs carried out up-tower are sensitive to environmental conditions whose effect on the material properties during processing needs to be accounted for to accurately predict the repair outcome. This study investigates the effect of moisture content from en… ▽ More Efficient structural repair of wind turbine blades is essential to limiting global warming and reducing the Levelized Cost of Energy (LCOE). Repairs carried out up-tower are sensitive to environmental conditions whose effect on the material properties during processing needs to be accounted for to accurately predict the repair outcome. This study investigates the effect of moisture content from environmental exposure on the cure kinetics of an infusion resin system used in wind turbine blade manufacturing and repair and provides an experimentally validated finite element tool for the analysis of cure cycle repairs as a function of repair geometry and moisture content. Moisture absorption tests on the two-part infusion system reported up to 12% moisture uptake by the curing agent under high temperature and relative humidity conditions. Differential scanning calorimetric measurements of resin in the presence of moisture revealed an accelerated cure behavior. Numerical predictions of a repair model agreed very well with the corresponding lab-scale repair and revealed a substantial temperature lag within the repair patch which resulted in thermal gradients and spatial distribution of the degree of cure. It was shown that the repair geometry and the accelerated-cure kinetics greatly influenced the temperature and cure distribution within the repair. The proposed approach can be used to reduce turbine downtime by minimizing the curing time. △ Less

Submitted 3 April, 2023; originally announced April 2023.

Comments: 29 pages, 12 figures

POD-mode-augmented wall model and its applications to flows at non-equilibrium conditions

Authors: Christoffer Hansen, Xiang IA Yang, Mahdi Abkar

Abstract: Insights gained from modal analysis are invoked for predictive large-eddy simulation (LES) wall modeling. Specifically, we augment the law of the wall (LoW) by an additional mode based on a one-dimensional proper orthogonal decomposition (POD) applied to a 2D turbulent channel. The constructed wall model contains two modes, i.e., the LoW mode and the POD-based mode, and the model matches with the… ▽ More Insights gained from modal analysis are invoked for predictive large-eddy simulation (LES) wall modeling. Specifically, we augment the law of the wall (LoW) by an additional mode based on a one-dimensional proper orthogonal decomposition (POD) applied to a 2D turbulent channel. The constructed wall model contains two modes, i.e., the LoW mode and the POD-based mode, and the model matches with the LES at two, instead of one, off-wall locations. To show that the proposed model captures non-equilibrium effects, we perform a-priori and a-posteriori tests in the context of both equilibrium and non-equilibrium flows. The a-priori tests show that the proposed wall model captures extreme wall-shear stress events better than the equilibrium wall model. The model also captures non-equilibrium effects due to adverse pressure gradients. The a-posteriori tests show that the wall model captures the rapid decrease and the initial decrease of the streamwise wall-shear stress in channels subjected to suddenly imposed adverse and transverse pressure gradients, respectively, both of which are missed by currently available wall models. These results show promise in applying modal analysis for turbulence wall modeling. In particular, the results show that employing multiple modes helps in the modeling of non-equilibrium flows. △ Less

Submitted 16 October, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Alternating north-south brightness ratio of Ganymede's auroral ovals: Hubble Space Telescope observations around the Juno PJ34 flyby

Authors: Joachim Saur, Stefan Duling, Alexandre Wennmacher, Clarissa Willmes, Lorenz Roth, Darrell F. Strobel, Frédéric Allegrini, Fran Bagenal, Scott J. Bolton, Bertrand Bonfond, George Clark, Randy Gladstone, T. K. Greathouse, Denis C. Grodent, Candice J. Hansen, W. S. Kurth, Glenn S. Orton, Kurt D. Retherford, Abigail M. Rymer, Ali H. Sulaiman

Abstract: We report results of Hubble Space Telescope observations from Ganymede's orbitally trailing side which were taken around the flyby of the Juno spacecraft on June 7, 2021. We find that Ganymede's northern and southern auroral ovals alternate in brightness such that the oval facing Jupiter's magnetospheric plasma sheet is brighter than the other one. This suggests that the generator that powers Gany… ▽ More We report results of Hubble Space Telescope observations from Ganymede's orbitally trailing side which were taken around the flyby of the Juno spacecraft on June 7, 2021. We find that Ganymede's northern and southern auroral ovals alternate in brightness such that the oval facing Jupiter's magnetospheric plasma sheet is brighter than the other one. This suggests that the generator that powers Ganymede's aurora is the momentum of the Jovian plasma sheet north and south of Ganymede's magnetosphere. Magnetic coupling of Ganymede to the plasma sheet above and below the moon causes asymmetric magnetic stresses and electromagnetic energy fluxes ultimately powering the auroral acceleration process. No clear statistically significant time variability of the auroral emission on short time scales of 100s could be resolved. We show that electron energy fluxes of several tens of mW m$^{-2}$ are required for its OI 1356 Å$\;$ emission making Ganymede a very poor auroral emitter. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted for publication with Geophys. Res. Let

CUBES, the Cassegrain U-Band Efficient Spectrograph

Authors: S. Cristiani, J. M. Alcalá, S. H. P. Alencar, S. A. Balashev, N. Bastian, B. Barbuy, U. Battino, A. Calcines, G. Calderone, P. Cambianica, R. Carini, B. Carter, S. Cassisi, B. V. Castilho, G. Cescutti, N. Christlieb, R. Cirami, I. Coretti, R. Cooke, S. Covino, G. Cremonese, K. Cunha, G. Cupani, A. R. da Silva, V. De Caprio , et al. (52 additional authors not shown)

Abstract: In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lowe… ▽ More In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: SPIE proceedings, SPIE Astronomical Telescopes + Instrumentation 2022, Montréal, Canada; 20 pages, 13 figures, 2 tables

Long-term quantification and characterisation of wind farm noise amplitude modulation

Authors: Phuc D. Nguyen, Kristy L. Hansen, Peter Catcheside, Colin Hansen, Branko Zajamsek

Abstract: The large-scale expansion of wind farms has prompted community debate regarding adverse impacts of wind farm noise (WFN). One of the most annoying and potentially sleep disturbing components of WFN is amplitude modulation (AM). Here we quantified and characterised AM over one year using acoustical and meteorological data measured at three locations near three wind farms. We found that the diurnal… ▽ More The large-scale expansion of wind farms has prompted community debate regarding adverse impacts of wind farm noise (WFN). One of the most annoying and potentially sleep disturbing components of WFN is amplitude modulation (AM). Here we quantified and characterised AM over one year using acoustical and meteorological data measured at three locations near three wind farms. We found that the diurnal variation of outdoor AM prevalence was substantial, the nighttime prevalence was approximately 2 to 5 times higher than the daytime prevalence. On average, indoor AM occurred during the nighttime from 1.1 to 1.7 times less often than outdoor AM, but the indoor AM depth was higher than that measured outdoors. We observed an association between AM prevalence and sunset and sunrise. AM occurred more often at downwind and crosswind conditions. These findings provide important insights into long term WFN characteristics that will help to inform future WFN assessment guidelines. △ Less

Submitted 29 May, 2022; originally announced June 2022.

Journal ref: Measurement 2021

Multi-input model uncertainty analysis for long-range wind farm noise predictions

Authors: Phuc D. Nguyen, Kristy L. Hansen, Branko Zajamsek, Peter Catcheside, Colin H. Hansen

Abstract: One of the major sources of uncertainty in predictions of wind farm noise (WFN) reflect parametric and model structure uncertainty. The model structure uncertainty is a systematic uncertainty, which relates to uncertainty about the appropriate mathematical structure of the models. Here we quantified the model structure uncertainty in predicting WFN arising from multi-input models, including nine g… ▽ More One of the major sources of uncertainty in predictions of wind farm noise (WFN) reflect parametric and model structure uncertainty. The model structure uncertainty is a systematic uncertainty, which relates to uncertainty about the appropriate mathematical structure of the models. Here we quantified the model structure uncertainty in predicting WFN arising from multi-input models, including nine ground impedance and four wind speed profile models. We used a numerical ray tracing sound propagation model for predicting the noise level at different receivers. We found that variations between different ground impedance models and wind speed profile models were significant sources of uncertainty, and that these sources contributed to predicted noise level differences in excess of 10 dBA at distances greater than 3.5 km. We also found that differences between atmospheric vertical wind speed profile models were the main source of uncertainty in predicting WFN at long-range distances. When predicting WFN, it is important to acknowledge variability associated with different models as this contributes to the uncertainty of the predicted values. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Diagnosing Magnetic Fields in Cylindrical Implosions with Oblique Proton Radiography

Authors: P. V. Heuer, L. S. Leal, J. R. Davies, E. C. Hansen, D. H. Barnak, J. L. Peebles, F. García-Rubio, B. Pollock, J. Moody, A. Birkel, F. H. Seguin

Abstract: Two experiments on the OMEGA Laser System used oblique proton radiography to measure magnetic fields in cylindrical implosions with and without an applied axial magnetic field. Although the goal of both experiments was to measure the magnitude of the compressed axial magnetic field in the core of the implosion, this field was obfuscated by two features in the coronal plasma produced by the compres… ▽ More Two experiments on the OMEGA Laser System used oblique proton radiography to measure magnetic fields in cylindrical implosions with and without an applied axial magnetic field. Although the goal of both experiments was to measure the magnitude of the compressed axial magnetic field in the core of the implosion, this field was obfuscated by two features in the coronal plasma produced by the compression beams: an azimuthal self-generated magnetic field and small length scale, high-amplitude structures attributed to collisionless effects. In order to understand these features, synthetic radiographs are generated using fields produced by 3-D HYDRA simulations. These synthetic radiographs reproduce the features of the experimental radiographs with the exception of the small-scale structures. A direct inversion algorithm is successfully applied to a synthetic radiograph, but is only partially able to invert the experimental radiographs in part because some protons are blocked by the field coils. The origins of the radiograph features and their dependence on various experimental parameters are explored. The results of this analysis should inform future measurements of compressed axial magnetic fields in cylindrical implosions. △ Less

Submitted 22 July, 2022; v1 submitted 22 March, 2022; originally announced March 2022.

Impact of Electron Transport Models on Capillary Discharge Plasmas

Authors: A. Diaw, S. J. Coleman, N. M. Cook, J. Edelen, E. C. Hansen, P. Tzeferacos

Abstract: Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing nex… ▽ More Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al. (ii) Spitzer, and (iii) Epperlein-Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over $20\%$ while predicting a lower azimuthal magnetic field. Moreover, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 9 pages, 5 figures

Iterative image reconstruction for CT with unmatched projection matrices using the generalized minimal residual algorithm

Authors: Emil Y. Sidky, Per Christian Hansen, Jakob S. Jørgensen, Xiaochuan Pan

Abstract: The generalized minimal residual (GMRES) algorithm is applied to image reconstruction using linear computed tomography (CT) models. The GMRES algorithm iteratively solves square, non-symmetric linear systems and it has practical application to CT when using unmatched back-projector/projector pairs and when applying preconditioning. The GMRES algorithm is demonstrated on a 3D CT image reconstructio… ▽ More The generalized minimal residual (GMRES) algorithm is applied to image reconstruction using linear computed tomography (CT) models. The GMRES algorithm iteratively solves square, non-symmetric linear systems and it has practical application to CT when using unmatched back-projector/projector pairs and when applying preconditioning. The GMRES algorithm is demonstrated on a 3D CT image reconstruction problem where it is seen that use of unmatched projection matrices does not prevent convergence, while using an unmatched pair in the related conjugate gradients for least-squares (CGLS) algorithm leads to divergent iteration. Implementation of preconditioning using GMRES is also demonstrated. △ Less

Submitted 3 May, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: Accepted to the 2022 CT Meeting

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

Authors: Joseph M. Levesque, Andy S. Liao, Patrick Hartigan, Rachel P. Young, Matthew Trantham, Sallee Klein, William Gray, Mario Manuel, Gennady Fiksel, Joseph Katz, Chikang Li, Andrew Birkel, Petros Tzeferacos, Edward C. Hansen, Benjamin Khiar, John M. Foster, Carolyn Kuranz

Abstract: The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising… ▽ More The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvénic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counter-propagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. We compare features in the reconstructed fields to two-dimensional MHD simulations of the system. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 36 pages, 12 figures

Journal ref: Published in Physics of Plasmas 29, 012106 (2022)

Hypotheses for Triton's Plumes: New Analyses and Future Remote Sensing Tests

Authors: Jason D. Hofgartner, Samuel P. D. Birch, Julie Castillo, Will M. Grundy, Candice J. Hansen, Alexander G. Hayes, Carly J. A. Howett, Terry A. Hurford, Emily S. Martin, Karl L. Mitchell, Tom A. Nordheim, Michael J. Poston, Louise M. Prockter, Lynnae C. Quick, Paul Schenk, Rebecca N. Schindhelm, Orkan M. Umurhan

Abstract: At least two active plumes were observed on Neptune's moon Triton during the Voyager 2 flyby in 1989. Models for Triton's plumes have previously been grouped into five hypotheses, two of which are primarily atmospheric phenomena and are generally considered unlikely, and three of which include eruptive processes and are plausible. These hypotheses are compared, including new arguments, such as com… ▽ More At least two active plumes were observed on Neptune's moon Triton during the Voyager 2 flyby in 1989. Models for Triton's plumes have previously been grouped into five hypotheses, two of which are primarily atmospheric phenomena and are generally considered unlikely, and three of which include eruptive processes and are plausible. These hypotheses are compared, including new arguments, such as comparisons based on current understanding of Mars, Enceladus, and Pluto. An eruption model based on a solar-powered, solid-state greenhouse effect was previously considered the leading hypothesis for Triton's plumes, in part due to the proximity of the plumes to the subsolar latitude during the Voyager 2 flyby and the distribution of Triton's fans that are putatively deposits from former plumes. The other two eruption hypotheses are powered by internal heat, not solar insolation. Based on new analyses of the ostensible relation between the latitude of the subsolar point on Triton and the geographic locations of the plumes and fans, we argue that neither the locations of the plumes nor fans are strong evidence in favor of the solar-powered hypothesis. We conclude that all three eruption hypotheses should be considered further. Five tests are presented that could be implemented with remote sensing observations from future spacecraft to confidently distinguish among the eruption hypotheses for Triton's plumes. The five tests are based on the: (1) composition and thickness of Triton's southern hemisphere terrains, (2) composition of fan deposits, (3) distribution of active plumes, (4) distribution of fans, and (5) surface temperature at the locations of plumes and/or fans. The tests are independent, but complementary, and implementable with a single flyby mission such as the Trident mission concept. We note that, in the case of the solar-driven hypothesis, the 2030s and 2040s may be the last ... △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Icarus

Promoting global stability in data-driven models of quadratic nonlinear dynamics

Authors: Alan A. Kaptanoglu, Jared L. Callaham, Christopher J. Hansen, Aleksandr Aravkin, Steven L. Brunton

Abstract: Modeling realistic fluid and plasma flows is computationally intensive, motivating the use of reduced-order models for a variety of scientific and engineering tasks. However, it is challenging to characterize, much less guarantee, the global stability (i.e., long-time boundedness) of these models. The seminal work of Schlegel and Noack (JFM, 2015) provided a theorem outlining necessary and suffici… ▽ More Modeling realistic fluid and plasma flows is computationally intensive, motivating the use of reduced-order models for a variety of scientific and engineering tasks. However, it is challenging to characterize, much less guarantee, the global stability (i.e., long-time boundedness) of these models. The seminal work of Schlegel and Noack (JFM, 2015) provided a theorem outlining necessary and sufficient conditions to ensure global stability in systems with energy-preserving, quadratic nonlinearities, with the goal of evaluating the stability of projection-based models. In this work, we incorporate this theorem into modern data-driven models obtained via machine learning. First, we propose that this theorem should be a standard diagnostic for the stability of projection-based and data-driven models, examining the conditions under which it holds. Second, we illustrate how to modify the objective function in machine learning algorithms to promote globally stable models, with implications for the modeling of fluid and plasma flows. Specifically, we introduce a modified "trapping SINDy" algorithm based on the sparse identification of nonlinear dynamics (SINDy) method. This method enables the identification of models that, by construction, only produce bounded trajectories. The effectiveness and accuracy of this approach are demonstrated on a broad set of examples of varying model complexity and physical origin, including the vortex shedding in the wake of a circular cylinder. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Journal ref: Phys. Rev. Fluids 6, 094401 (2021)

A hidden Markov model for describing turbostratic disorder applied to carbon blacks and graphene

Authors: A G Hart, T C Hansen, W F Kuhs

Abstract: We present a mathematical framework to represent turbostratic disorder in materials like carbon blacks, smectites, and twisted $n$-layer graphene. In particular, the set of all possible disordered layers, including rotated, shifted, and curved layers form a stochastic sequence governed by a hidden Markov model. The probability distribution over the set of layer types is treated as an element of a… ▽ More We present a mathematical framework to represent turbostratic disorder in materials like carbon blacks, smectites, and twisted $n$-layer graphene. In particular, the set of all possible disordered layers, including rotated, shifted, and curved layers form a stochastic sequence governed by a hidden Markov model. The probability distribution over the set of layer types is treated as an element of a Hilbert space, and using tools of Fourier analysis and functional analysis, we develop expressions for the scattering cross sections of a broad class of disordered materials. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 15 pages 3 figures

Journal ref: May 2019 Acta Crystallographica Section A: Foundations and Advances 75(3) A75, 501-516

A Markov theoretic description of stacking disordered aperiodic crystals including ice and opaline silica

Authors: Allen G Hart, Thomas C Hansen, Werner F Kuhs

Abstract: We review the Markov theoretic description of 1D aperiodic crystals, describing the stacking-faulted crystal polytype as a special case of an aperiodic crystal. Under this description we generalise the centrosymmetric unit cell underlying a topologically centrosymmetric crystal to a reversible Markov chain underlying a reversible aperiodic crystal. We show that for the close-packed structure, almo… ▽ More We review the Markov theoretic description of 1D aperiodic crystals, describing the stacking-faulted crystal polytype as a special case of an aperiodic crystal. Under this description we generalise the centrosymmetric unit cell underlying a topologically centrosymmetric crystal to a reversible Markov chain underlying a reversible aperiodic crystal. We show that for the close-packed structure, almost all stackings are irreversible when the interaction reichweite is greater than 4. Moreover, we present an analytic expression of the scattering cross section of a large class of stacking disordered aperiodic crystals, lacking translational symmetry of their layers, including ice and opaline silica (opal CT). We then relate the observed stackings and their underlying reichweite to the physics of various nucleation and growth processes of disordered ice. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 16 pages 5 figures

Journal ref: July 2018 Acta Crystallographica Section A: Foundations and Advances 74(4):357-372

Stellar astrophysics in the near UV with VLT-CUBES

Authors: H. Ernandes, C. J. Evans, B. Barbuy, B. Castilho, G. Cescutti, N. Christlieb, S. Cristiani, G. Cupani, P. Di Marcantonio, M. Franchini, C. Hansen, A. Quirrenbach, R. Smiljanic

Abstract: Alongside future observations with the new European Extremely Large Telescope (ELT), optimised instruments on the 8-10m generation of telescopes will still be competitive at 'ground UV' wavelengths (3000-4000 A). The near UV provides a wealth of unique information on the nucleosynthesis of iron-peak elements, molecules, and neutron-capture elements. In the context of development of the near-UV CUB… ▽ More Alongside future observations with the new European Extremely Large Telescope (ELT), optimised instruments on the 8-10m generation of telescopes will still be competitive at 'ground UV' wavelengths (3000-4000 A). The near UV provides a wealth of unique information on the nucleosynthesis of iron-peak elements, molecules, and neutron-capture elements. In the context of development of the near-UV CUBES spectrograph for ESO's Very Large Telescope (VLT), we are investigating the impact of spectral resolution on the ability to estimate chemical abundances for beryllium and more than 30 iron-peak and heavy elements. From work ahead of the Phase A conceptual design of CUBES, here we present a comparison of the elements observable at the notional resolving power of CUBES (R~20,000) to those with VLT-UVES (R~40,000). For most of the considered lines signal-to-noise is a more critical factor than resolution. We summarise the elements accessible with CUBES, several of which (e.g. Be, Ge, Hf) are now the focus of quantitative simulations as part of the ongoing Phase A study. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Journal ref: Proceedings of the SPIE, Volume 11447, id. 1144760 11 pp. (2020)

The structure of global conservation laws in Galerkin plasma models

Authors: Alan A. Kaptanoglu, Kyle D. Morgan, Christopher J. Hansen, Steven L. Brunton

Abstract: Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD). Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, can furnish this gap in the lower levels of the model hierarchy. In the present work, we d… ▽ More Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD). Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, can furnish this gap in the lower levels of the model hierarchy. In the present work, we develop low-dimensional Galerkin plasma models which preserve global conservation laws by construction. This additional model structure enables physics-constrained machine learning algorithms that can discover these types of low-dimensional plasma models directly from data. This formulation relies on an energy-based inner product which takes into account all of the dynamic variables. The theoretical results here build a bridge to the extensive Galerkin literature in fluid mechanics, and facilitate the development of physics-constrained reduced-order models from plasma data. △ Less

Submitted 9 January, 2021; originally announced January 2021.

Comments: arXiv admin note: text overlap with arXiv:2004.10389

The upgrade of the ALICE TPC with GEMs and continuous readout

Authors: J. Adolfsson, M. Ahmed, S. Aiola, J. Alme, T. Alt, W. Amend, F. Anastasopoulos, C. Andrei, M. Angelsmark, V. Anguelov, A. Anjam, H. Appelshäuser, V. Aprodu, O. Arnold, M. Arslandok, D. Baitinger, M. Ball, G. G. Barnaföldi, E. Bartsch, P. Becht, R. Bellwied, A. Berdnikova, M. Berger, N. Bialas, P. Bialas , et al. (210 additional authors not shown)

Abstract: The upgrade of the ALICE TPC will allow the experiment to cope with the high interaction rates foreseen for the forthcoming Run 3 and Run 4 at the CERN LHC. In this article, we describe the design of new readout chambers and front-end electronics, which are driven by the goals of the experiment. Gas Electron Multiplier (GEM) detectors arranged in stacks containing four GEMs each, and continuous re… ▽ More The upgrade of the ALICE TPC will allow the experiment to cope with the high interaction rates foreseen for the forthcoming Run 3 and Run 4 at the CERN LHC. In this article, we describe the design of new readout chambers and front-end electronics, which are driven by the goals of the experiment. Gas Electron Multiplier (GEM) detectors arranged in stacks containing four GEMs each, and continuous readout electronics based on the SAMPA chip, an ALICE development, are replacing the previous elements. The construction of these new elements, together with their associated quality control procedures, is explained in detail. Finally, the readout chamber and front-end electronics cards replacement, together with the commissioning of the detector prior to installation in the experimental cavern, are presented. After a nine-year period of R&D, construction, and assembly, the upgrade of the TPC was completed in 2020. △ Less

Submitted 25 March, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 88 pages, 60 figures

Journal ref: JINST 16 (2021) P03022

Physics-constrained, low-dimensional models for MHD: First-principles and data-driven approaches

Authors: Alan A. Kaptanoglu, Kyle D. Morgan, Chris J. Hansen, Steven L. Brunton

Abstract: Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by pro… ▽ More Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally-consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a 3D spheromak experiment. This manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics, and facilitates future principled development of projection-based and data-driven models for plasmas. △ Less

Submitted 15 July, 2021; v1 submitted 22 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. E 104, 015206 (2021)

Advanced modeling for the HIT-SI Experiment

Authors: Alan A. Kaptanoglu, Thomas E. Benedett, Kyle D. Morgan, Chris J. Hansen, Thomas R. Jarboe

Abstract: A two-temperature magnetohydrodynamic (MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet code and used to model plasma dynamics in the HIT-SI experiment. When compared with single-temperature Hall-MHD, the two-temperature Hall-MHD model demonstrates improved qualitative agreement with experimental measurements, including: far-infrared interferome… ▽ More A two-temperature magnetohydrodynamic (MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet code and used to model plasma dynamics in the HIT-SI experiment. When compared with single-temperature Hall-MHD, the two-temperature Hall-MHD model demonstrates improved qualitative agreement with experimental measurements, including: far-infrared interferometry, ion Doppler spectroscopy, Thomson scattering, and magnetic probe measurements. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5-68.5 kHz range. At all frequencies the two-temperature models result in increased toroidal current, lower chord-averaged density, and symmetrization of the current centroid, relative to single-temperature simulations. Both codes produce higher average temperatures and toroidal currents as the injector frequency is increased. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and indicate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for the artificial diffusivity, and Dirichlet wall temperature and density. Artificial diffusivity and the density boundary condition both significantly modify the plasma density profiles, leading to larger average temperatures, higher toroidal current, and increased relative density fluctuations at low diffusivity and low wall density. High power, low density simulations at 14.5 kHz achieve sufficiently high gain (G = 5) to generate significant volumes of closed flux lasting 1-2 injector periods. △ Less

Submitted 1 March, 2020; originally announced March 2020.

Tracking the Ultraviolet Photochemistry of Thiophenone During and Beyond the Initial Ultrafast Ring Opening

Authors: Shashank Pathak, Lea M. Ibele, Rebecca Boll, Carlo Callegari, Alexander Demidovich, Benjamin Erk, Raimund Feifel, Ruaridh Forbes, Michele Di Fraia, Luca Giannessi, Christopher S. Hansen, David M. P. Holland, Rebecca A. Ingle, Robert Mason, Oksana Plekan, Kevin C. Prince, Arnaud Rouzée, Richard J. Squibb, Jan Tross, Michael N. R. Ashfold, Basile F. E. Curchod, Daniel Rolles

Abstract: Photoinduced isomerization reactions, including ring-opening reactions, lie at the heart of many processes in nature. The mechanisms of such reactions are determined by a delicate interplay of coupled electronic and nuclear dynamics unfolding on the femtosecond scale, followed by the slower redistribution of energy into different vibrational degrees of freedom. Here we apply time-resolved photoele… ▽ More Photoinduced isomerization reactions, including ring-opening reactions, lie at the heart of many processes in nature. The mechanisms of such reactions are determined by a delicate interplay of coupled electronic and nuclear dynamics unfolding on the femtosecond scale, followed by the slower redistribution of energy into different vibrational degrees of freedom. Here we apply time-resolved photoelectron spectroscopy with a seeded extreme ultraviolet free electron laser to trace the ultrafast ring opening of gas phase thiophenone molecules following photoexcitation at 265 nm. When combined with cutting edge ab initio electronic structure and molecular dynamics calculations of both the excited and ground state molecules, the results provide unprecedented insights into both electronic and nuclear dynamics of this fundamental class of reactions. The initial ring opening and non-adiabatic coupling to the electronic ground state is shown to be driven by ballistic SC bond extension and to be complete within 350 femtoseconds. Theory and experiment also allow clear visualization of the rich ground-state dynamics involving formation of, and interconversion between, several ring opened isomers and the reformed cyclic structure, and fragmentation (CO loss) over much longer timescales. △ Less

Submitted 14 March, 2020; v1 submitted 1 December, 2019; originally announced December 2019.

Comments: 40 pages, 21 figures Changes from the previous version: 1) Added theoretical calculations for explaining long timescale changes in shown in Figure 4(a). 2) Reworked on fitting (modelling) of the experimental data in Figure 2(b). Added another panel i.e. Figure 2(c). 3) Other minor changes and rewording in response to questions and suggestions by the referees

Characterizing Magnetized Plasmas with Dynamic Mode Decomposition

Authors: Alan A. Kaptanoglu, Kyle D. Morgan, Chris J. Hansen, Steven L. Brunton

Abstract: Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear, computationally intensive, and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable r… ▽ More Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear, computationally intensive, and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the HIT-SI experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas. △ Less

Submitted 15 November, 2019; originally announced November 2019.

4D MRI: Robust sorting of free breathing MRI slices for use in interventional settings

Authors: Gino Gulamhussene, Fabian Joeres, Marko Rak, Maciej Pech, Christian Hansen

Abstract: Purpose: We aim to develop a robust 4D MRI method for large FOVs enabling the extraction of irregular respiratory motion that is readily usable with all MRI machines and thus applicable to support a wide range of interventional settings. Method: We propose a 4D MRI reconstruction method to capture an arbitrary number of breathing states. It uses template updates in navigator slices and search re… ▽ More Purpose: We aim to develop a robust 4D MRI method for large FOVs enabling the extraction of irregular respiratory motion that is readily usable with all MRI machines and thus applicable to support a wide range of interventional settings. Method: We propose a 4D MRI reconstruction method to capture an arbitrary number of breathing states. It uses template updates in navigator slices and search regions for fast and robust vessel cross-section tracking. It captures FOVs of 255 mm x 320 mm x 228 mm at a spatial resolution of 1.82 mm x 1.82 mm x 4mm and temporal resolution of 200ms. A total of 37 4D MRIs of 13 healthy subjects were reconstructed to validate the method. A quantitative evaluation of the reconstruction rate and speed of both the new and baseline method was performed. Additionally, a study with ten radiologists was conducted to assess the subjective reconstruction quality of both methods. Results: Our results indicate improved mean reconstruction rates compared to the baseline method (79.4\% vs. 45.5\%) and improved mean reconstruction times (24s vs. 73s) per subject. Interventional radiologists perceive the reconstruction quality of our method as higher compared to the baseline (262.5 points vs. 217.5 points, p=0.02). Conclusions: Template updates are an effective and efficient way to increase 4D MRI reconstruction rates and to achieve better reconstruction quality. Search regions reduce reconstruction time. These improvements increase the applicability of 4D MRI as a base for seamless support of interventional image guidance in percutaneous interventions. △ Less

Submitted 6 July, 2020; v1 submitted 4 October, 2019; originally announced October 2019.

Comments: 16 pages, 11 figures

Mobile Phone Data for Children on the Move: Challenges and Opportunities

Authors: Vedran Sekara, Elisa Omodei, Laura Healy, Jan Beise, Claus Hansen, Danzhen You, Saskia Blume, Manuel Garcia-Herranz

Abstract: Today, 95% of the global population has 2G mobile phone coverage and the number of individuals who own a mobile phone is at an all time high. Mobile phones generate rich data on billions of people across different societal contexts and have in the last decade helped redefine how we do research and build tools to understand society. As such, mobile phone data has the potential to revolutionize how… ▽ More Today, 95% of the global population has 2G mobile phone coverage and the number of individuals who own a mobile phone is at an all time high. Mobile phones generate rich data on billions of people across different societal contexts and have in the last decade helped redefine how we do research and build tools to understand society. As such, mobile phone data has the potential to revolutionize how we tackle humanitarian problems, such as the many suffered by refugees all over the world. While promising, mobile phone data and the new computational approaches bring both opportunities and challenges. Mobile phone traces contain detailed information regarding people's whereabouts, social life, and even financial standing. Therefore, developing and adopting strategies that open data up to the wider humanitarian and international development community for analysis and research while simultaneously protecting the privacy of individuals is of paramount importance. Here we outline the challenging situation of children on the move and actions UNICEF is pushing in helping displaced children and youth globally, and discuss opportunities where mobile phone data can be used. We identify three key challenges: data access, data and algorithmic bias, and operationalization of research, which need to be addressed if mobile phone data is to be successfully applied in humanitarian contexts. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 13 pages, book chapter

Compressed Sensing MRI With Variable Density Averaging (CS-VDA) Outperforms Full Sampling At Low SNR

Authors: Jasper Schoormans, Gustav J. Strijkers, Anders C. Hansen, Aart J. Nederveen, Bram F. Coolen

Abstract: We investigated whether a combination of k-space undersampling and variable density averaging enhances image quality for low-SNR MRI acquisitions. We implemented 3D Cartesian k-space prospective undersampling with a variable number of averages for each k-line. The performance of this compressed sensing with variable-density averaging (CS-VDA) method was evaluated in retrospective analysis of fully… ▽ More We investigated whether a combination of k-space undersampling and variable density averaging enhances image quality for low-SNR MRI acquisitions. We implemented 3D Cartesian k-space prospective undersampling with a variable number of averages for each k-line. The performance of this compressed sensing with variable-density averaging (CS-VDA) method was evaluated in retrospective analysis of fully sampled phantom MRI measurements, as well as for prospectively accelerated in vivo 3D brain and knee MRI scans. Both phantom and in vivo results showed that acquisitions using the CS-VDA approach resulted in better image quality as compared to full sampling of k-space in the same scan time. Specifically, CS-VDA with a higher number of averages in the center of k-space resulted in the best image quality, apparent from increased anatomical detail with preserved soft-tissue contrast. This novel approach will facilitate improved image quality of inherently low SNR data, such as those with high-resolution or specific contrast-weightings with low SNR efficiency. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 11 pages, 9 figures

Journal ref: Physics in Medicine & Biology 65.4 (2020): 045004

Enceladus and Titan: Emerging Worlds of the Solar System (ESA Voyage 2050 White Paper)

Authors: Ali Sulaiman, Nicholas Achilleos, Sushil Atreya, Cesar Bertucci, Andrew Coates, Michele Dougherty, Lina Hadid, Candice Hansen, Mika Holmberg, Hsiang-Wen Hsu, Tomoki Kimura, William Kurth, Alice Le Gall, James McKevitt, Michiko Morooka, Go Murakami, Leonardo Regoli, Elias Roussos, Joachim Saur, Oleg Shebanits, Anezina Solomonidou, Jan-Erik Wahlund, J. Hunter Waite

Abstract: Some of the major discoveries of the recent Cassini-Huygens mission have put Titan and Enceladus firmly on the Solar System map. The mission has revolutionised our view of Solar System satellites, arguably matching their scientific importance with that of their planet. While Cassini-Huygens has made big surprises in revealing Titan's organically rich environment and Enceladus' cryovolcanism, the m… ▽ More Some of the major discoveries of the recent Cassini-Huygens mission have put Titan and Enceladus firmly on the Solar System map. The mission has revolutionised our view of Solar System satellites, arguably matching their scientific importance with that of their planet. While Cassini-Huygens has made big surprises in revealing Titan's organically rich environment and Enceladus' cryovolcanism, the mission's success naturally leads us to further probe these findings. We advocate the acknowledgement of Titan and Enceladus science as highly relevant to ESA's long-term roadmap, as logical follow-on to Cassini-Huygens. In this white paper, we will outline important science questions regarding these satellites and identify the pertinent science themes we recommend ESA cover during the Voyage 2050 planning cycle. Addressing these science themes would make major advancements to the present knowledge we have about the Solar System, its formation, evolution and likelihood that other habitable environments exist outside the Earth's biosphere. △ Less

Submitted 24 October, 2019; v1 submitted 5 August, 2019; originally announced August 2019.

Comments: White paper submitted in response to Voyage 2050 long-term plan in the ESA Science Programme

A comparison between the two lobes of comet 67P/Churyumov-Gerasimenko based on D/H ratios in H2O measured with the Rosetta/ROSINA DFMS

Authors: Isaac R. H. G. Schroeder I, Kathrin Altwegg, Hans Balsiger, Jean-Jacques Berthelier, Michael R. Combi, Johan De Keyser, Björn Fiethe, Stephen A. Fuselier, Tamas I. Gombosi, Kenneth C. Hansen, Martin Rubin, Yinsi Shou, Valeriy M. Tenishev, Thierry Sémon, Susanne F. Wampfler, Peter Wurz

Abstract: The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its… ▽ More The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: 10 pages, 6 figures. Submitted to MNRAS

MSC Class: 85-05 ACM Class: J.2

Hall Effect in the coma of 67P/Churyumov-Gerasimenko

Authors: Z. Huang, G. Toth, T. I. Gombosi, X. Jia, M. R. Combi, K. C. Hansen, N. Fougere, Y. Shou, V. Tenishev, K. Altwegg, M. Rubin

Abstract: Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium… ▽ More Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium observed regions of very weak magnetic field outside the predicted diamagnetic cavity. In this paper, we simulate the inner coma with the Hall magnetohydrodynamics equations and show that the Hall effect is important in the inner coma environment. The magnetic field topology becomes complex and magnetic reconnection occurs on the dayside when the Hall effect is taken into account. The magnetic reconnection on the dayside can generate weak magnetic filed regions outside the global diamagnetic cavity, which may explain the Rosetta Plasma Consortium observations. We conclude that the substantial change in the inner coma environment is due to the fact that the ion inertial length (or gyro radius) is not much smaller than the size of the diamagnetic cavity. △ Less

Submitted 11 January, 2018; originally announced January 2018.

Comments: 23 pages, 6 figurs

Journal ref: Published in MNRAS in 2018

Equilibrium reconstruction with 3D eddy currents in the Lithium Tokamak eXperiment

Authors: C. Hansen, D. P. Boyle, J. C. Schmitt, R. Majeski

Abstract: Axisymmetric free-boundary equilibrium reconstructions of tokamak plasmas in the Lithium Tokamak eXperiment (LTX) are performed using the PSI-Tri equilibrium code. Reconstructions in LTX are complicated by the presence of long-lived non-axisymmetric eddy currents generated by vacuum vessel and first wall structures. To account for this effect, reconstructions are performed with additional toroidal… ▽ More Axisymmetric free-boundary equilibrium reconstructions of tokamak plasmas in the Lithium Tokamak eXperiment (LTX) are performed using the PSI-Tri equilibrium code. Reconstructions in LTX are complicated by the presence of long-lived non-axisymmetric eddy currents generated by vacuum vessel and first wall structures. To account for this effect, reconstructions are performed with additional toroidal current sources in these conducting regions. The source distributions are fixed poloidally, but their scale is adjusted as part of the full reconstruction. Eddy distributions are computed by toroidally averaging currents, generated by coupling to vacuum field coils, from a simplified 3D filament model of important conducting structures. The full 3D eddy current fields are also used to enable the inclusion of local magnetic field measurements, which have strong 3D eddy current pick-up, as reconstruction constraints. Using this method, equilibrium reconstruction yields good agreement with all available diagnostic signals. An accompanying field perturbation produced by 3D eddy currents on the plasma surface with primarily n=2, m=1 character is also predicted for these equilibria. △ Less

Submitted 20 January, 2017; originally announced January 2017.

Comments: 13 pages, 6 figures

Journal ref: Physics of Plasmas 24, 042513 (2017)

Lattice constants and expansivities of gas hydrates from 10K up to the stability limit

Authors: Thomas C. Hansen, Andrzej Falenty, M. Mangir Murshed, Werner F. Kuhs

Abstract: In a combination of neutron and synchrotron diffraction the lattice constants and expansivities of hydrogenated and deuterated CH4-, CO2-, Xe- (structure type I) and N2-hydrate (structure type II) from 10 K up to the stability limit under pressure were established. Some important results emerge from our analysis: (1) Despite the larger guest-size of CO2 as compared to methane, CO2- hydrate has the… ▽ More In a combination of neutron and synchrotron diffraction the lattice constants and expansivities of hydrogenated and deuterated CH4-, CO2-, Xe- (structure type I) and N2-hydrate (structure type II) from 10 K up to the stability limit under pressure were established. Some important results emerge from our analysis: (1) Despite the larger guest-size of CO2 as compared to methane, CO2- hydrate has the smaller lattice constants at low temperatures which we ascribe to the larger attractive guest-host interaction of the CO2-water system. (2) The expansivity of CO2-hydrate is larger than for CH4-hydrate which leads to larger lattice constants for the former at temperatures above ~ 150 K; this is likely due to the higher motional degrees of freedom of the CO2 guest molecules. (3) The cage filling does not affect significantly the lattice constants in CH4- and CO2-hydrate in contrast to Xe-hydrate for which the effect is quantitatively established. (4) Similar to ice Ih, the deuterated compounds have slightly larger lattice constants for all investigated systems which can be ascribed to the somewhat weaker H-bonding; the isotopic difference is smallest for the Xesystem,in which the large Xe atoms lead to an increase of averaged H-bond distances. (5) Compared to ice Ih the high temperature expansivities are about 50% larger; in contrast to ice Ih, there is no negative thermal expansion at low temperature. △ Less

Submitted 23 October, 2015; originally announced October 2015.

Comments: Unpublished contribution to the 8th International Conference on Gas Hydrates (ICGH-8), Beijing, China, 28 July - 1 August 2015 (was only available to the conference participants)

The EUROnu Project

Authors: T. R. Edgecock, O. Caretta, T. Davenne, C. Densham, M. Fitton, D. Kelliher, P. Loveridge, S. Machida, C. Prior, C. Rogers, M. Rooney, J. Thomason, D. Wilcox, E. Wildner, I. Efthymiopoulos, R. Garoby, S. Gilardoni, C. Hansen, E. Benedetto, E. Jensen, A. Kosmicki, M. Martini, J. Osborne, G. Prior, T. Stora , et al. (146 additional authors not shown)

Abstract: The EUROnu project has studied three possible options for future, high intensity neutrino oscillation facilities in Europe. The first is a Super Beam, in which the neutrinos come from the decay of pions created by bombarding targets with a 4 MW proton beam from the CERN High Power Superconducting Proton Linac. The far detector for this facility is the 500 kt MEMPHYS water Cherenkov, located in the… ▽ More The EUROnu project has studied three possible options for future, high intensity neutrino oscillation facilities in Europe. The first is a Super Beam, in which the neutrinos come from the decay of pions created by bombarding targets with a 4 MW proton beam from the CERN High Power Superconducting Proton Linac. The far detector for this facility is the 500 kt MEMPHYS water Cherenkov, located in the Fréjus tunnel. The second facility is the Neutrino Factory, in which the neutrinos come from the decay of μ+ and μ- beams in a storage ring. The far detector in this case is a 100 kt Magnetised Iron Neutrino Detector at a baseline of 2000 km. The third option is a Beta Beam, in which the neutrinos come from the decay of beta emitting isotopes, in particular 6He and 18Ne, also stored in a ring. The far detector is also the MEMPHYS detector in the Fréjus tunnel. EUROnu has undertaken conceptual designs of these facilities and studied the performance of the detectors. Based on this, it has determined the physics reach of each facility, in particular for the measurement of CP violation in the lepton sector, and estimated the cost of construction. These have demonstrated that the best facility to build is the Neutrino Factory. However, if a powerful proton driver is constructed for another purpose or if the MEMPHYS detector is built for astroparticle physics, the Super Beam also becomes very attractive. △ Less

Submitted 17 May, 2013; originally announced May 2013.

Comments: Results from the Framework Programme 7 project EUROnu, which studied three possible accelerator facilities for future high intensity neutrino oscillation facilities in Europe

Journal ref: Phys. Rev. ST Accel. Beams 16 021002 (2013)

Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector

Authors: K. Abe, N. Abgrall, Y. Ajima, H. Aihara, J. B. Albert, C. Andreopoulos, B. Andrieu, M. D. Anerella, S. Aoki, O. Araoka, J. Argyriades, A. Ariga, T. Ariga, S. Assylbekov, D. Autiero, A. Badertscher, M. Barbi, G. J. Barker, G. Barr, M. Bass, M. Batkiewicz, F. Bay, S. Bentham, V. Berardi, B. E. Berger , et al. (407 additional authors not shown)

Abstract: Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure o… ▽ More Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure of iron target plates and scintillator trackers. INGRID directly monitors the muon neutrino beam profile center and intensity using the number of observed neutrino events in each module. The neutrino beam direction is measured with accuracy better than 0.4 mrad from the measured profile center. The normalized event rate is measured with 4% precision. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 32 pages, 27 figures, submitted to Nucl. Instr. and Meth. A

Beta Beams Implementation at CERN

Authors: Christian Hansen

Abstract: Beta Beam,the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive ions beta decay in a storage ring, called Decay Ring (DR), is the base of one of the proposed next generation neutrino oscillation facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown neutrino oscillation parameters depend on… ▽ More Beta Beam,the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive ions beta decay in a storage ring, called Decay Ring (DR), is the base of one of the proposed next generation neutrino oscillation facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown neutrino oscillation parameters depend on the Decay Ring's ion intensity and of it's duty factor (the filled ratio of the ring). Therefore efficient ion production, stripping, bunching, acceleration and storing are crucial sub-projects under study and development within the Beta Beam collaboration. Specifically the feasibility of these tasks as parts of a Beta Beam implementation at CERN will be discussed in this report. The positive impact of the large θ13 indications from T2K on the Beta Beam performance will also be discussed. △ Less

Submitted 9 September, 2011; originally announced September 2011.

Comments: NUFACT11, 8 pages, 11 figures

Report number: EUROnu-WP4-030

Simulations of Bunch Merging in a Beta Beam Decay Ring

Authors: Daniel Christopher Heinrich, Christian Hansen, Antoine Chance

Abstract: To further study neutrino oscillation properties a Beta Beam facility has been proposed. Beta decaying ions with high kinetic energy are stored in a storage ring ("Decay Ring") with straight sections to create pure focused (anti) electron neutrino beams. However to reach high sensitivity to neutrino oscillation parameters in the experiment the bunched beam intensity and duty cycle in the DR have t… ▽ More To further study neutrino oscillation properties a Beta Beam facility has been proposed. Beta decaying ions with high kinetic energy are stored in a storage ring ("Decay Ring") with straight sections to create pure focused (anti) electron neutrino beams. However to reach high sensitivity to neutrino oscillation parameters in the experiment the bunched beam intensity and duty cycle in the DR have to be optimized. The first CERN-based scenario, using 6He and 18Ne as neutrino sources, has been studied using a bunch merging RF scheme. Two RF cavities at different frequencies are used to capture newly injected bunches and then merge them into the stored bunches. It was shown that this scheme could satisfy the requirements on intensity and duty cycle set by the experiment. This merging scheme has now been revised with new simulation software providing new results for 6He and 18Ne. Furthermore bunch merging has been studied for the second CERN-based scenario using 8Li and 8B. △ Less

Submitted 8 September, 2011; originally announced September 2011.

Comments: NUFACT10, 3 pages, 3 figures

Report number: EUROnu-WP4-012

Collective Effect Studies of a Beta Beam Decay Ring

Authors: Christian Hansen, Giovanni Rumolo

Abstract: The Beta Beam, the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive ions beta decay in a storage ring called the Decay Ring (DR), is the basis of one of the proposed next generation neutrino oscillation facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown neutrino oscillation parameters… ▽ More The Beta Beam, the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive ions beta decay in a storage ring called the Decay Ring (DR), is the basis of one of the proposed next generation neutrino oscillation facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown neutrino oscillation parameters depend on the DR's ion intensity and of its duty factor (the filled ratio of the ring). Different methods, including analytical calculations and multiparticle tracking simulations, were used to estimate the DR's potential to contain enough ions in as small a part of the ring as needed for the sensitivities. Studies of transverse blow up of the beams due to resonance wake fields show that a very challenging upper limit of the transverse broadband impedance is required to avoid instabilities and beam loss. △ Less

Submitted 8 September, 2011; originally announced September 2011.

Comments: NUFACT2010, 3 pages, 2 figures

Report number: EUROnu-WP4-013

The T2K Experiment

Authors: T2K Collaboration, K. Abe, N. Abgrall, H. Aihara, Y. Ajima, J. B. Albert, D. Allan, P. -A. Amaudruz, C. Andreopoulos, B. Andrieu, M. D. Anerella, C. Angelsen, S. Aoki, O. Araoka, J. Argyriades, A. Ariga, T. Ariga, S. Assylbekov, J. P. A. M. de André, D. Autiero, A. Badertscher, O. Ballester, M. Barbi, G. J. Barker, P. Baron , et al. (499 additional authors not shown)

Abstract: The T2K experiment is a long-baseline neutrino oscillation experiment. Its main goal is to measure the last unknown lepton sector mixing angle θ_{13} by observing ν_e appearance in a ν_μ beam. It also aims to make a precision measurement of the known oscillation parameters, Δm^{2}_{23} and sin^{2} 2θ_{23}, via ν_μ disappearance studies. Other goals of the experiment include various neutrino cross… ▽ More The T2K experiment is a long-baseline neutrino oscillation experiment. Its main goal is to measure the last unknown lepton sector mixing angle θ_{13} by observing ν_e appearance in a ν_μ beam. It also aims to make a precision measurement of the known oscillation parameters, Δm^{2}_{23} and sin^{2} 2θ_{23}, via ν_μ disappearance studies. Other goals of the experiment include various neutrino cross section measurements and sterile neutrino searches. The experiment uses an intense proton beam generated by the J-PARC accelerator in Tokai, Japan, and is composed of a neutrino beamline, a near detector complex (ND280), and a far detector (Super-Kamiokande) located 295 km away from J-PARC. This paper provides a comprehensive review of the instrumentation aspect of the T2K experiment and a summary of the vital information for each subsystem. △ Less

Submitted 8 June, 2011; v1 submitted 6 June, 2011; originally announced June 2011.

Comments: 33 pages, 32 figures, Submitted and accepted by NIM A. Editor: Prof. Chang Kee Jung, Department of Physics and Astronomy, SUNY Stony Brook, chang.jung@sunysb.edu, 631-632-8108 Submit Edited to remove line numbers