-
Direct Comparison of Isobaric and Isochoric Vitrification of Two Aqueous Solutions with Photon Counting X-Ray Computed Tomography
Authors:
Jason T. Parker,
Anthony N. Consiglio,
Boris Rubinsky,
Simo A. Mäkiharju
Abstract:
Vitrification is a promising approach for ice-free cryopreservation of biological material, but progress is hindered by the limited set of experimental tools for studying processes in the interior of the vitrified matter. Isochoric cryopreservation chambers are often metallic, and their opacity prevents direct visual observation. In this study, we introduce photon counting X-ray computed tomograph…
▽ More
Vitrification is a promising approach for ice-free cryopreservation of biological material, but progress is hindered by the limited set of experimental tools for studying processes in the interior of the vitrified matter. Isochoric cryopreservation chambers are often metallic, and their opacity prevents direct visual observation. In this study, we introduce photon counting X-ray computed tomography (CT) to compare the effects of rigid isochoric and unconfined isobaric conditions on vitrification and ice formation during cooling of two aqueous solutions: 50 wt% DMSO and a coral vitrification solution, CVS1. Previous studies have only compared vitrification in isochoric systems with isobaric systems that have an exposed air-liquid interface. We use a movable piston to replicate the surface and thermal boundary conditions of the isochoric system yet maintain isobaric conditions. When controlling for the boundary conditions we find that similar ice and vapor volume fractions form during cooling in isochoric and isobaric conditions. Interestingly, we observe distinct ice and vapor cavity morphology in the isochoric systems, possibly due to vapor outgassing or cavitation as rapid cooling causes the pressure to drop in the confined systems. These observations highlight the array of thermo-fluidic processes that occur during vitrification in confined aqueous systems and motivate the further application of imaging techniques such as photon counting X-ray CT in fundamental studies of vitrification.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
In-Lab X-ray Particle Velocimetry for Multiphase Flows: Design Principles and Demonstration of $O$(1 kHz) XPV
Authors:
Jason T. Parker,
Till Dreier,
Daniel Nilsson,
Simo A. Mäkiharju
Abstract:
We combine X-ray-specific tracer particles, a photon counting detector, and a liquid metal jet anode X-ray source to achieve $O$(1 kHz) X-ray imaging speeds in the laboratory, 15$\times$ faster than previous comparable studies with $O$(50 $μ$m) tracers. To examine the limits of this measurement technique we conduct three experiments: 2D and 3D X-ray particle velocimetry (XPV) of Poiseuille pipe fl…
▽ More
We combine X-ray-specific tracer particles, a photon counting detector, and a liquid metal jet anode X-ray source to achieve $O$(1 kHz) X-ray imaging speeds in the laboratory, 15$\times$ faster than previous comparable studies with $O$(50 $μ$m) tracers. To examine the limits of this measurement technique we conduct three experiments: 2D and 3D X-ray particle velocimetry (XPV) of Poiseuille pipe flow, 3D XPV of flow around a Taylor bubble, and 3D scalar mixing with a laminar jet. These experiments demonstrate the performance improvement achievable by combining the aforementioned elements, the applicability to multiphase flows and deforming systems, and the potential to capture scalar and vector quantities simultaneously. Most importantly, these experiments are conducted with a laboratory-scale system, showing that in-lab X-ray particle velocimetry techniques are now becoming usable for a wider range of flows of interest. Furthermore, the design of XPV experiments is discussed to clarify the trade offs between achievable imaging speed, domain size and spatiotemporal resolution.
△ Less
Submitted 13 October, 2023; v1 submitted 25 April, 2023;
originally announced April 2023.
-
Enhanced Laboratory X-ray Particle Tracking Velocimetry With Newly Developed Tungsten-Coated O(50 $μ$m) Tracers
Authors:
Jason T. Parker,
Jessica DeBerardinis,
Simo A. Mäkiharju
Abstract:
Tracer particles designed specifically for X-ray particle tracking and imaging velocimetry (XPTV and XPIV) are necessary to widen the range of flows that can be studied with these techniques. In this study, we demonstrate in-lab XPTV using new, custom-designed $O$(50 $μ$m) diameter tungsten-coated hollow carbon spheres and a single energy threshold photon counting detector. To explore the measurem…
▽ More
Tracer particles designed specifically for X-ray particle tracking and imaging velocimetry (XPTV and XPIV) are necessary to widen the range of flows that can be studied with these techniques. In this study, we demonstrate in-lab XPTV using new, custom-designed $O$(50 $μ$m) diameter tungsten-coated hollow carbon spheres and a single energy threshold photon counting detector. To explore the measurement quality enhancement enabled by the new tracer particles and photon counting detector, a well understood Poiseulle pipe flow is measured. The data show agreement with the analytical solution for the depth-averaged velocity profile. The experiment also shows that the tungsten-coated particles achieve higher contrast and are better localized than previously available silver-coated particles, making faster and more precise measurements attainable. The particles are manufactured with a readily scalable chemical vapor deposition process.
We further show that laboratory XPTV is practical with currently available energy-resolving photon counting detectors (PCDs), despite their presently lower spatiotemporal resolution compared to scintillating detectors. This finding suggests that energy-thresholding identification of different classes of tracers is feasible, further motivating the exploration of the X-ray tracer particle design space. The latest generation of PCDs are incorporating multiple energy thresholds, and have higher count rate limits. In the near future one could potentially expand on the work presented and track multiple tracer species and scalar fields simultaneously.
△ Less
Submitted 15 July, 2022;
originally announced July 2022.
-
Experimentally Validated Simulations of 50 μm X-ray PIV Tracer Particles
Authors:
Jason T. Parker,
Simo A. Mäkiharju
Abstract:
We evaluate Beer-Lambert (BL) ray-tracing and Monte Carlo N-Particle (MCNP) photon tracking simulations for prediction and comparison of X-ray imaging system performance. These simulation tools can aid the methodical design of laboratory-scale X-ray particle image velocimetry (XPIV) experiments and tracer particles by predicting image quality. Particle image signal-to-noise ratio (SNR) is used as…
▽ More
We evaluate Beer-Lambert (BL) ray-tracing and Monte Carlo N-Particle (MCNP) photon tracking simulations for prediction and comparison of X-ray imaging system performance. These simulation tools can aid the methodical design of laboratory-scale X-ray particle image velocimetry (XPIV) experiments and tracer particles by predicting image quality. Particle image signal-to-noise ratio (SNR) is used as the metric of system performance. Simulated and experiment data of hollow, silver-coated, glass sphere tracer particles (AGSF-33) are compared. As predicted by the simulations, the AGSF-33 particles are visible with a SNR greater than unity in 100~ms exposure time images, demonstrating their potential as X-ray PIV or particle tracking velocimetry (XPTV) tracers. The BL approach predicts the image contrast, is computationally inexpensive, and enables the exploration of a vast parameter space for system design. MCNP simulations, on the other hand, predict experiment images slightly more accurately, but are more than an order of magnitude more computationally expensive than BL simulations. For most practical XPIV system design applications, the higher computational expense of MCNP is likely not justified by the modest accuracy improvement compared to BL.
△ Less
Submitted 17 January, 2022; v1 submitted 27 September, 2021;
originally announced September 2021.
-
Testing nonlocal models of electron thermal conduction for magnetic and inertial confinement fusion applications
Authors:
Jonathan Peter Brodrick,
Robert J. Kingham,
Michael M. Marinak,
Mehul V. Patel,
Alex V. Chankin,
John Omotani,
Maxim Umansky,
Dario Del Sorbo,
Ben Dudson,
Joseph Thomas Parker,
Gary D. Kerbel,
Mark Sherlock,
Christopher P Ridgers
Abstract:
Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (…
▽ More
Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (NFLF) model of Dimits, Joseph and Umansky; and (iii) Schurtz, Nicolaï and Busquet's multigroup diffusion model (SNB). We find that while the EIC and NFLF models accurately predict the damping rate of a small-amplitude temperature perturbation (within 10% at moderate collisionalities), they overestimate the peak heat flow by as much as 35% and do not predict preheat in the more relevant case where there is a large temperature difference. The SNB model, however, agrees better with VFP results for the latter problem if care is taken with the definition of the mean free path. Additionally, we present for the first time a comparison of the SNB model against a VFP code for a hohlraum-relevant problem with inhomogeneous ionisation and show that the model overestimates the heat flow in the helium gas-fill by a factor of ~2 despite predicting the peak heat flux to within 16%.
△ Less
Submitted 6 September, 2017; v1 submitted 28 April, 2017;
originally announced April 2017.
-
Suppression of phase mixing in drift-kinetic plasma turbulence
Authors:
J. T. Parker,
E. G. Highcock,
A. A. Schekochihin,
P. J. Dellar
Abstract:
Transfer of free energy from large to small velocity-space scales by phase mixing leads to Landau damping in a linear plasma. In a turbulent drift-kinetic plasma, this transfer is statistically nearly canceled by an inverse transfer from small to large velocity-space scales due to "anti-phase-mixing" modes excited by a stochastic form of plasma echo. Fluid moments (density, velocity, temperature)…
▽ More
Transfer of free energy from large to small velocity-space scales by phase mixing leads to Landau damping in a linear plasma. In a turbulent drift-kinetic plasma, this transfer is statistically nearly canceled by an inverse transfer from small to large velocity-space scales due to "anti-phase-mixing" modes excited by a stochastic form of plasma echo. Fluid moments (density, velocity, temperature) are thus approximately energetically isolated from the higher moments of the distribution function, so phase mixing is ineffective as a dissipation mechanism when the plasma collisionality is small.
△ Less
Submitted 22 March, 2016;
originally announced March 2016.
-
Gyrokinetic simulations of fusion plasmas using a spectral velocity space representation
Authors:
Joseph Thomas Parker
Abstract:
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the g…
▽ More
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms.
We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code.
Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.
△ Less
Submitted 15 March, 2016;
originally announced March 2016.
-
Parametric Bilinear Generalized Approximate Message Passing
Authors:
Jason T. Parker,
Philip Schniter
Abstract:
We propose a scheme to estimate the parameters $b_i$ and $c_j$ of the bilinear form $z_m=\sum_{i,j} b_i z_m^{(i,j)} c_j$ from noisy measurements $\{y_m\}_{m=1}^M$, where $y_m$ and $z_m$ are related through an arbitrary likelihood function and $z_m^{(i,j)}$ are known. Our scheme is based on generalized approximate message passing (G-AMP): it treats $b_i$ and $c_j$ as random variables and…
▽ More
We propose a scheme to estimate the parameters $b_i$ and $c_j$ of the bilinear form $z_m=\sum_{i,j} b_i z_m^{(i,j)} c_j$ from noisy measurements $\{y_m\}_{m=1}^M$, where $y_m$ and $z_m$ are related through an arbitrary likelihood function and $z_m^{(i,j)}$ are known. Our scheme is based on generalized approximate message passing (G-AMP): it treats $b_i$ and $c_j$ as random variables and $z_m^{(i,j)}$ as an i.i.d.\ Gaussian 3-way tensor in order to derive a tractable simplification of the sum-product algorithm in the large-system limit. It generalizes previous instances of bilinear G-AMP, such as those that estimate matrices $\boldsymbol{B}$ and $\boldsymbol{C}$ from a noisy measurement of $\boldsymbol{Z}=\boldsymbol{BC}$, allowing the application of AMP methods to problems such as self-calibration, blind deconvolution, and matrix compressive sensing. Numerical experiments confirm the accuracy and computational efficiency of the proposed approach.
△ Less
Submitted 11 December, 2015; v1 submitted 30 August, 2015;
originally announced August 2015.
-
Phase mixing vs. nonlinear advection in drift-kinetic plasma turbulence
Authors:
A. A. Schekochihin,
J. T. Parker,
E. G. Highcock,
P. J. Dellar,
W. Dorland,
G. W. Hammett
Abstract:
A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to…
▽ More
A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the "anti-phase-mixing" effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the "critical balance" between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).
△ Less
Submitted 25 April, 2016; v1 submitted 24 August, 2015;
originally announced August 2015.
-
Fourier-Hermite spectral representation for the Vlasov-Poisson system in the weakly collisional limit
Authors:
J. T. Parker,
P. J. Dellar
Abstract:
We study Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation. We describe the propagation of free energy in phase space using forwards and backwards propagating Hermite modes recently developed for gyrokinetics [Schekochihin et al. (2014)]. The change in the electric field corresponds to the net Hermite flux via a free energy evolution equation. In line…
▽ More
We study Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation. We describe the propagation of free energy in phase space using forwards and backwards propagating Hermite modes recently developed for gyrokinetics [Schekochihin et al. (2014)]. The change in the electric field corresponds to the net Hermite flux via a free energy evolution equation. In linear Landau damping, decay in the electric field corresponds to forward propagating Hermite modes; in nonlinear damping, the initial decay is followed by a growth phase characterised by the generation of backwards propagating Hermite modes by the nonlinear term. The free energy content of the backwards propagating modes increases exponentially until balancing that of the forward propagating modes. Thereafter there is no systematic net Hermite flux, so the electric field cannot decay and the nonlinearity effectively suppresses Landau damping. These simulations are performed using the fully-spectral 5D gyrokinetics code SpectroGK [Parker et al. 2014], modified to solve the 1+1D Vlasov-Poisson system. This captures Landau damping via an iterated Lénard-Bernstein collision operator or via Hou-Li filtering in velocity space. Therefore the code is applicable even in regimes where phase-mixing and filamentation are dominant.
△ Less
Submitted 7 July, 2014;
originally announced July 2014.
-
Irreversible energy flow in forced Vlasov dynamics
Authors:
G. G. Plunk,
J. T. Parker
Abstract:
A recent paper [Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation…
▽ More
A recent paper [Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation exhibits asymptotically zero (irreversible) dissipation to all orders under a detuning of the forcing frequency and the characteristic frequency associated with particle streaming. We first prove this by direct calculation, tracking energy flow in terms of certain exact conservation laws of the linear (collisionless) Vlasov equation. Then we analyze the steady-state solutions in detail using a weakly collisional Hermite-moment formulation, and compare with numerical solution. This leads to a detailed description of the Hermite energy spectrum, and a proof of no dissipation at all orders, complementing the collisionless Vlasov result.
△ Less
Submitted 7 January, 2015; v1 submitted 28 February, 2014;
originally announced February 2014.
-
Bilinear Generalized Approximate Message Passing
Authors:
Jason T. Parker,
Philip Schniter,
Volkan Cevher
Abstract:
We extend the generalized approximate message passing (G-AMP) approach, originally proposed for high-dimensional generalized-linear regression in the context of compressive sensing, to the generalized-bilinear case, which enables its application to matrix completion, robust PCA, dictionary learning, and related matrix-factorization problems. In the first part of the paper, we derive our Bilinear G…
▽ More
We extend the generalized approximate message passing (G-AMP) approach, originally proposed for high-dimensional generalized-linear regression in the context of compressive sensing, to the generalized-bilinear case, which enables its application to matrix completion, robust PCA, dictionary learning, and related matrix-factorization problems. In the first part of the paper, we derive our Bilinear G-AMP (BiG-AMP) algorithm as an approximation of the sum-product belief propagation algorithm in the high-dimensional limit, where central-limit theorem arguments and Taylor-series approximations apply, and under the assumption of statistically independent matrix entries with known priors. In addition, we propose an adaptive damping mechanism that aids convergence under finite problem sizes, an expectation-maximization (EM)-based method to automatically tune the parameters of the assumed priors, and two rank-selection strategies. In the second part of the paper, we discuss the specializations of EM-BiG-AMP to the problems of matrix completion, robust PCA, and dictionary learning, and present the results of an extensive empirical study comparing EM-BiG-AMP to state-of-the-art algorithms on each problem. Our numerical results, using both synthetic and real-world datasets, demonstrate that EM-BiG-AMP yields excellent reconstruction accuracy (often best in class) while maintaining competitive runtimes and avoiding the need to tune algorithmic parameters.
△ Less
Submitted 5 June, 2014; v1 submitted 9 October, 2013;
originally announced October 2013.