-
A transition from boundary- to bulk-driven acoustic streaming due to nonlinear thermoviscous effects at high acoustic energy densities
Authors:
Jonas Helboe Joergensen,
Wei Qiu,
Henrik Bruus
Abstract:
Acoustic streaming is studied in a rectangular microfluidic channel. It is demonstrated theoretically, numerically, and experimentally with good agreement, frictional heating can alter the streaming pattern qualitatively at high acoustic energy densities E_ac above 500 J/m^3. The study shows, how as a function of increasing E_ac at fixed frequency, the traditional boundary-driven four streaming ro…
▽ More
Acoustic streaming is studied in a rectangular microfluidic channel. It is demonstrated theoretically, numerically, and experimentally with good agreement, frictional heating can alter the streaming pattern qualitatively at high acoustic energy densities E_ac above 500 J/m^3. The study shows, how as a function of increasing E_ac at fixed frequency, the traditional boundary-driven four streaming rolls created at a half-wave standing-wave resonance, transition into two large streaming rolls. This nonlinear transition occurs because friction heats up the fluid resulting in a temperature gradient, which spawns an acoustic body force in the bulk that drives thermoacoustic streaming.
△ Less
Submitted 21 December, 2021;
originally announced December 2021.
-
Theory and modeling of nonperturbative effects at high acoustic energy densities in thermoviscous acoustofluidics
Authors:
Jonas Helboe Joergensen,
Henrik Bruus
Abstract:
A theoretical model of thermal boundary layers and acoustic heating in microscale acoustofluidic devices is presented. It includes effective boundary conditions allowing for simulations in three dimensions. The model is extended by an iterative scheme to incorporate nonlinear thermoviscous effects not captured by standard perturbation theory. The model predicts that the dominant nonperturbative ef…
▽ More
A theoretical model of thermal boundary layers and acoustic heating in microscale acoustofluidic devices is presented. It includes effective boundary conditions allowing for simulations in three dimensions. The model is extended by an iterative scheme to incorporate nonlinear thermoviscous effects not captured by standard perturbation theory. The model predicts that the dominant nonperturbative effects in these devices are due to the dependency of thermoacoustic streaming on gradients in the steady temperature induced by a combination of internal frictional heating, external heating, and thermal convection. The model enables simulations in a nonperturbative regime relevant for design and fabrication of high-throughput acoustofluidic devices.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
Age and helium content of the open cluster NGC 6791 from multiple eclipsing binary members. III. Constraints from a subgiant
Authors:
K. Brogaard,
F. Grundahl,
E. L. Sandquist,
D. Slumstrup,
M. L. Jensen,
J. B. Thomsen,
J. H. Jørgensen,
J. R. Larsen,
S. T. Bjørn,
C. T. G. Sørensen,
H. Bruntt,
T. Arentoft,
S. Frandsen,
J. Jessen-Hansen,
J. A. Orosz,
R. Mathieu,
A. Geller,
N. Ryde,
D. Stello,
S. Meibom,
I. Platais
Abstract:
Models of stellar structure and evolution can be constrained using accurate measurements of the parameters of eclipsing binary members of open clusters. Multiple binary stars provide the means to tighten the constraints and, in turn, to improve the precision and accuracy of the age estimate of the host cluster. In the previous two papers of this series, we have demonstrated the use of measurements…
▽ More
Models of stellar structure and evolution can be constrained using accurate measurements of the parameters of eclipsing binary members of open clusters. Multiple binary stars provide the means to tighten the constraints and, in turn, to improve the precision and accuracy of the age estimate of the host cluster. In the previous two papers of this series, we have demonstrated the use of measurements of multiple eclipsing binaries in the old open cluster NGC6791 to set tighter constraints on the properties of stellar models than was previously possible, thereby improving both the accuracy and precision of the cluster age. We identify and measure the properties of a non-eclipsing cluster member, V56, in NGC\,6791 and demonstrate how this provides additional model constraints that support and strengthen our previous findings. We analyse multi-epoch spectra of V56 from FLAMES in conjunction with the existing photometry and measurements of eclipsing binaries in NGC6971. The parameters of the V56 components are found to be $M_{\rm p}=1.103\pm 0.008 M_{\odot}$ and $M_{\rm s}=0.974\pm 0.007 M_{\odot}$, $R_{\rm p}=1.764\pm0.099 R_{\odot}$ and $R_{\rm s}=1.045\pm0.057 R_{\odot}$, $T_{\rm eff,p}=5447\pm125$ K and $T_{\rm eff,s}=5552\pm125$ K, and surface [Fe/H]=$+0.29\pm0.06$ assuming that they have the same abundance. The derived properties strengthen our previous best estimate of the cluster age of $8.3\pm0.3$ Gyr and the mass of stars on the lower red giant branch (RGB), which is $M_{\rm RGB} = 1.15\pm0.02M_{\odot}$ for NGC6791. These numbers therefore continue to serve as verification points for other methods of age and mass measures, such as asteroseismology.
△ Less
Submitted 29 April, 2021;
originally announced April 2021.
-
Fast microscale acoustic streaming driven by a temperature-gradient-induced non-dissipative acoustic body force
Authors:
Wei Qiu,
Jonas Helboe Joergensen,
Enrico Corato,
Henrik Bruus,
Per Augustsson
Abstract:
We study acoustic streaming in liquids driven by a non-dissipative acoustic body force created by light-induced temperature gradients. This thermoacoustic streaming produces a velocity amplitude approximately 50 times higher than boundary-driven Rayleigh streaming and 90 times higher than Rayleigh-Benard convection at a temperature gradient of 10 K/mm in the channel. Further, Rayleigh streaming is…
▽ More
We study acoustic streaming in liquids driven by a non-dissipative acoustic body force created by light-induced temperature gradients. This thermoacoustic streaming produces a velocity amplitude approximately 50 times higher than boundary-driven Rayleigh streaming and 90 times higher than Rayleigh-Benard convection at a temperature gradient of 10 K/mm in the channel. Further, Rayleigh streaming is altered by the acoustic body force at a temperature gradient of only 0.5 K/mm. The thermoacoustic streaming allows for modular flow control and enhanced heat transfer at the microscale. Our study provides the groundwork for studying microscale acoustic streaming coupled with temperature fields.
△ Less
Submitted 18 March, 2021;
originally announced March 2021.
-
Theory of pressure acoustics with thermoviscous boundary layers and streaming in elastic cavities
Authors:
Jonas H. Joergensen,
Henrik Bruus
Abstract:
We present an effective thermoviscous theory of acoustofluidics including pressure acoustics, thermoviscous boundary layers, and streaming for fluids embedded in elastic cavities. By including thermal fields, we thus extend the effective viscous theory by Bach and Bruus, J. Acoust. Soc. Am. 144, 766 (2018). The acoustic temperature field and the thermoviscous boundary layers are incorporated analy…
▽ More
We present an effective thermoviscous theory of acoustofluidics including pressure acoustics, thermoviscous boundary layers, and streaming for fluids embedded in elastic cavities. By including thermal fields, we thus extend the effective viscous theory by Bach and Bruus, J. Acoust. Soc. Am. 144, 766 (2018). The acoustic temperature field and the thermoviscous boundary layers are incorporated analytically as effective boundary conditions and time-averaged body forces on the thermoacoustic bulk fields. Because it avoids resolving the thin boundary layers, the effective model allows for numerical simulation of both thermoviscous acoustic and time-averaged fields in 3D models of acoustofluidic systems. We show how the acoustic streaming depends strongly on steady and oscillating thermal fields through the temperature dependency of the material parameters, in particular the viscosity and the compressibility, affecting both the boundary conditions and spawning additional body forces in the bulk. We also show how even small steady temperature gradients (1 K/mm) induce gradients in compressibility and density that may result in very high streaming velocities (1 mm/s) for moderate acoustic energy densities (100 J/m^3).
△ Less
Submitted 14 December, 2020;
originally announced December 2020.
-
Probabilistic PARAFAC2
Authors:
Philip J. H. Jørgensen,
Søren F. V. Nielsen,
Jesper L. Hinrich,
Mikkel N. Schmidt,
Kristoffer H. Madsen,
Morten Mørup
Abstract:
The PARAFAC2 is a multimodal factor analysis model suitable for analyzing multi-way data when one of the modes has incomparable observation units, for example because of differences in signal sampling or batch sizes. A fully probabilistic treatment of the PARAFAC2 is desirable in order to improve robustness to noise and provide a well founded principle for determining the number of factors, but ch…
▽ More
The PARAFAC2 is a multimodal factor analysis model suitable for analyzing multi-way data when one of the modes has incomparable observation units, for example because of differences in signal sampling or batch sizes. A fully probabilistic treatment of the PARAFAC2 is desirable in order to improve robustness to noise and provide a well founded principle for determining the number of factors, but challenging because the factor loadings are constrained to be orthogonal. We develop two probabilistic formulations of the PARAFAC2 along with variational procedures for inference: In the one approach, the mean values of the factor loadings are orthogonal leading to closed form variational updates, and in the other, the factor loadings themselves are orthogonal using a matrix Von Mises-Fisher distribution. We contrast our probabilistic formulation to the conventional direct fitting algorithm based on maximum likelihood. On simulated data and real fluorescence spectroscopy and gas chromatography-mass spectrometry data, we compare our approach to the conventional PARAFAC2 model estimation and find that the probabilistic formulation is more robust to noise and model order misspecification. The probabilistic PARAFAC2 thus forms a promising framework for modeling multi-way data accounting for uncertainty.
△ Less
Submitted 21 June, 2018;
originally announced June 2018.
-
Linear scale bounds on dark matter--dark radiation interactions and connection with the small scale crisis of cold dark matter
Authors:
Maria Archidiacono,
Sebastian Bohr,
Steen Hannestad,
Jonas Helboe Jørgensen,
Julien Lesgourgues
Abstract:
One of the open questions in modern cosmology is the small scale crisis of the cold dark matter paradigm. Increasing attention has recently been devoted to self-interacting dark matter models as a possible answer. However, solving the so-called "missing satellites" problem requires in addition the presence of an extra relativistic particle (dubbed dark radiation) scattering with dark matter in the…
▽ More
One of the open questions in modern cosmology is the small scale crisis of the cold dark matter paradigm. Increasing attention has recently been devoted to self-interacting dark matter models as a possible answer. However, solving the so-called "missing satellites" problem requires in addition the presence of an extra relativistic particle (dubbed dark radiation) scattering with dark matter in the early universe. Here we investigate the impact of different theoretical models devising dark matter dark radiation interactions on large scale cosmological observables. We use cosmic microwave background data to put constraints on the dark radiation component and its coupling to dark matter. We find that the values of the coupling allowed by the data imply a cut-off scale of the halo mass function consistent with the one required to match the observations of satellites in the Milky Way.
△ Less
Submitted 26 June, 2017; v1 submitted 21 June, 2017;
originally announced June 2017.
-
Iodine versus Bromine Functionalization for Bottom-Up Graphene Nanoribbon Growth: Role of Diffusion
Authors:
Christopher Bronner,
Tomas Marangoni,
Daniel J. Rizzo,
Rebecca Durr,
Jakob Holm Jørgensen,
Felix R. Fischer,
Michael F. Crommie
Abstract:
Deterministic bottom-up approaches for synthesizing atomically well-defined graphene nanoribbons (GNRs) largely rely on the surface-catalyzed activation of selected labile bonds in a molecular precursor followed by step growth polymerization and cyclodehydrogenation. While the majority of successful GNR precursors rely on the homolytic cleavage of thermally labile C-Br bonds, the introduction of w…
▽ More
Deterministic bottom-up approaches for synthesizing atomically well-defined graphene nanoribbons (GNRs) largely rely on the surface-catalyzed activation of selected labile bonds in a molecular precursor followed by step growth polymerization and cyclodehydrogenation. While the majority of successful GNR precursors rely on the homolytic cleavage of thermally labile C-Br bonds, the introduction of weaker C-I bonds provides access to monomers that can be polymerized at significantly lower temperatures, thus helping to increase the flexibility of the GNR synthesis process. Scanning tunneling microscopy (STM) imaging of molecular precursors, activated intermediates, and polymers resulting from stepwise thermal annealing of both Br and I substituted precursors for chevron GNRs reveals that the polymerization of both precursors proceeds at similar temperatures on Au(111). This observation is consistent with diffusion-limited polymerization of the surface-stabilized radical intermediates that emerge from homolytic cleavage of either the C-Br or the C-I bonds.
△ Less
Submitted 27 March, 2017;
originally announced March 2017.
-
Probing the spin of the central black hole in the galactic centre with secondary images
Authors:
Jonas Helboe Jørgensen,
Ole Ole Eggers Bjælde,
Steen Hannestad
Abstract:
This paper explores the possibility of determining the spin of the supermassive black hole (SMBH) in Sgr A*, by using secondary images of stars orbiting the SMBH. The photons propagate close to the SMBH and their trajectories probe the space time in a region where the spin of the SMBH is important. We find the appearance of spikes in the secondary image, which depends on the angular momentum and s…
▽ More
This paper explores the possibility of determining the spin of the supermassive black hole (SMBH) in Sgr A*, by using secondary images of stars orbiting the SMBH. The photons propagate close to the SMBH and their trajectories probe the space time in a region where the spin of the SMBH is important. We find the appearance of spikes in the secondary image, which depends on the angular momentum and spin axis of the SMBH and study the specific case of the star S2 in detail. The spikes has a magnitude of $\sim 29$ in the K-band and the required angular resolution is of order 15-20 $μ$as. The combination of these two requirements poses an extreme observational challenge, but might be possible with interferometric observations in the sub-mm regime. The next possible time frame for observing this effect on the star S2 is in the late 2017 and then it repeats with the period of the star.
△ Less
Submitted 18 March, 2016;
originally announced March 2016.
-
Few-view single photon emission computed tomography (SPECT) reconstruction based on a blurred piecewise constant object model
Authors:
Paul A Wolf,
Jakob H Jørgensen,
Taly G Schmidt,
Emil Y Sidky
Abstract:
A sparsity-exploiting algorithm intended for few-view Single Photon Emission Computed Tomography (SPECT) reconstruction is proposed and characterized. The algorithm models the object as piecewise constant subject to a blurring operation. To validate that the algorithm closely approximates the true object in the noiseless case, projection data were generated from an object assuming this model and u…
▽ More
A sparsity-exploiting algorithm intended for few-view Single Photon Emission Computed Tomography (SPECT) reconstruction is proposed and characterized. The algorithm models the object as piecewise constant subject to a blurring operation. To validate that the algorithm closely approximates the true object in the noiseless case, projection data were generated from an object assuming this model and using the system matrix. Monte Carlo simulations were performed to provide more realistic data of a phantom with varying smoothness across the field of view. Reconstructions were performed across a sweep of two primary design parameters. The results demonstrate that the algorithm recovers the object in a noiseless simulation case. While the algorithm assumes a specific blurring model, the results suggest that the algorithm may provide high reconstruction accuracy even when the object does not match the assumed blurring model. Generally, increased values of the blurring parameter and TV weighting parameters reduced noise and streaking artifacts, while decreasing spatial resolution. As the number of views decreased from 60 to 9 the accuracy of images reconstructed using the proposed algorithm varied by less than 3%. Overall, the results demonstrate preliminary feasibility of a sparsity-exploiting reconstruction algorithm which may be beneficial for few-view SPECT.
△ Less
Submitted 4 December, 2012;
originally announced December 2012.
-
Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle-Pock algorithm
Authors:
Emil Y. Sidky,
Jakob H. Jørgensen,
Xiaochuan Pan
Abstract:
The primal-dual optimization algorithm developed in Chambolle and Pock (CP), 2011 is applied to various convex optimization problems of interest in computed tomography (CT) image reconstruction. This algorithm allows for rapid prototyping of optimization problems for the purpose of designing iterative image reconstruction algorithms for CT. The primal-dual algorithm is briefly summarized in the ar…
▽ More
The primal-dual optimization algorithm developed in Chambolle and Pock (CP), 2011 is applied to various convex optimization problems of interest in computed tomography (CT) image reconstruction. This algorithm allows for rapid prototyping of optimization problems for the purpose of designing iterative image reconstruction algorithms for CT. The primal-dual algorithm is briefly summarized in the article, and its potential for prototyping is demonstrated by explicitly deriving CP algorithm instances for many optimization problems relevant to CT. An example application modeling breast CT with low-intensity X-ray illumination is presented.
△ Less
Submitted 8 March, 2012; v1 submitted 23 November, 2011;
originally announced November 2011.
-
Ensuring convergence in total-variation-based reconstruction for accurate microcalcification imaging in breast X-ray CT
Authors:
Jakob H. Jørgensen,
Emil Y. Sidky,
Xiaochuan Pan
Abstract:
Breast X-ray CT imaging is being considered in screening as an extension to mammography. As a large fraction of the population will be exposed to radiation, low-dose imaging is essential. Iterative image reconstruction based on solving an optimization problem, such as Total-Variation minimization, shows potential for reconstruction from sparse-view data. For iterative methods it is important to en…
▽ More
Breast X-ray CT imaging is being considered in screening as an extension to mammography. As a large fraction of the population will be exposed to radiation, low-dose imaging is essential. Iterative image reconstruction based on solving an optimization problem, such as Total-Variation minimization, shows potential for reconstruction from sparse-view data. For iterative methods it is important to ensure convergence to an accurate solution, since important image features, such as presence of microcalcifications indicating breast cancer, may not be visible in a non-converged reconstruction, and this can have clinical significance. To prevent excessively long computational times, which is a practical concern for the large image arrays in CT, it is desirable to keep the number of iterations low, while still ensuring a sufficiently accurate reconstruction for the specific imaging task. This motivates the study of accurate convergence criteria for iterative image reconstruction. In simulation studies with a realistic breast phantom with microcalcifications we compare different convergence criteria for reliable reconstruction. Our results show that it can be challenging to ensure a sufficiently accurate microcalcification reconstruction, when using standard convergence criteria. In particular, the gray level of the small microcalcifications may not have converged long after the background tissue is reconstructed uniformly. We propose the use of the individual objective function gradient components to better monitor possible regions of non-converged variables. For microcalcifications we find empirically a large correlation between nonzero gradient components and non-converged variables, which occur precisely within the microcalcifications. This supports our claim that gradient components can be used to ensure convergence to a sufficiently accurate reconstruction.
△ Less
Submitted 10 November, 2011;
originally announced November 2011.
-
Quantifying admissible undersampling for sparsity-exploiting iterative image reconstruction in X-ray CT
Authors:
Jakob H. Jørgensen,
Emil Y. Sidky,
Xiaochuan Pan
Abstract:
Iterative image reconstruction (IIR) with sparsity-exploiting methods, such as total variation (TV) minimization, investigated in compressive sensing (CS) claim potentially large reductions in sampling requirements. Quantifying this claim for computed tomography (CT) is non-trivial, because both full sampling in the discrete-to-discrete imaging model and the reduction in sampling admitted by spars…
▽ More
Iterative image reconstruction (IIR) with sparsity-exploiting methods, such as total variation (TV) minimization, investigated in compressive sensing (CS) claim potentially large reductions in sampling requirements. Quantifying this claim for computed tomography (CT) is non-trivial, because both full sampling in the discrete-to-discrete imaging model and the reduction in sampling admitted by sparsity-exploiting methods are ill-defined. The present article proposes definitions of full sampling by introducing four sufficient-sampling conditions (SSCs). The SSCs are based on the condition number of the system matrix of a linear imaging model and address invertibility and stability. In the example application of breast CT, the SSCs are used as reference points of full sampling for quantifying the undersampling admitted by reconstruction through TV-minimization. In numerical simulations, factors affecting admissible undersampling are studied. Differences between few-view and few-detector bin reconstruction as well as a relation between object sparsity and admitted undersampling are quantified.
△ Less
Submitted 17 August, 2012; v1 submitted 3 September, 2011;
originally announced September 2011.
-
Accelerated gradient methods for total-variation-based CT image reconstruction
Authors:
Jakob Heide Jørgensen,
Tobias Lindstrøm Jensen,
Per Christian Hansen,
Søren Holdt Jensen,
Emil Y. Sidky,
Xiaochuan Pan
Abstract:
Total-variation (TV)-based Computed Tomography (CT) image reconstruction has shown experimentally to be capable of producing accurate reconstructions from sparse-view data. In particular TV-based reconstruction is very well suited for images with piecewise nearly constant regions. Computationally, however, TV-based reconstruction is much more demanding, especially for 3D imaging, and the reconstru…
▽ More
Total-variation (TV)-based Computed Tomography (CT) image reconstruction has shown experimentally to be capable of producing accurate reconstructions from sparse-view data. In particular TV-based reconstruction is very well suited for images with piecewise nearly constant regions. Computationally, however, TV-based reconstruction is much more demanding, especially for 3D imaging, and the reconstruction from clinical data sets is far from being close to real-time. This is undesirable from a clinical perspective, and thus there is an incentive to accelerate the solution of the underlying optimization problem. The TV reconstruction can in principle be found by any optimization method, but in practice the large-scale systems arising in CT image reconstruction preclude the use of memory-demanding methods such as Newton's method. The simple gradient method has much lower memory requirements, but exhibits slow convergence. In the present work we consider the use of two accelerated gradient-based methods, GPBB and UPN, for reducing the number of gradient method iterations needed to achieve a high-accuracy TV solution in CT image reconstruction. The former incorporates several heuristics from the optimization literature such as Barzilai-Borwein (BB) step size selection and nonmonotone line search. The latter uses a cleverly chosen sequence of auxiliary points to achieve a better convergence rate. The methods are memory efficient and equipped with a stopping criterion to ensure that the TV reconstruction has indeed been found. An implementation of the methods (in C with interface to Matlab) is available for download from http://www2.imm.dtu.dk/~pch/TVReg/. We compare the proposed methods with the standard gradient method, applied to a 3D test problem with synthetic few-view data. We find experimentally that for realistic parameters the proposed methods significantly outperform the gradient method.
△ Less
Submitted 19 May, 2011;
originally announced May 2011.
-
Implementation of an Optimal First-Order Method for Strongly Convex Total Variation Regularization
Authors:
Tobias Lindstrøm Jensen,
Jakob Heide Jørgensen,
Per Christian Hansen,
Søren Holdt Jensen
Abstract:
We present a practical implementation of an optimal first-order method, due to Nesterov, for large-scale total variation regularization in tomographic reconstruction, image deblurring, etc. The algorithm applies to $μ$-strongly convex objective functions with $L$-Lipschitz continuous gradient. In the framework of Nesterov both $μ$ and $L$ are assumed known -- an assumption that is seldom satisfied…
▽ More
We present a practical implementation of an optimal first-order method, due to Nesterov, for large-scale total variation regularization in tomographic reconstruction, image deblurring, etc. The algorithm applies to $μ$-strongly convex objective functions with $L$-Lipschitz continuous gradient. In the framework of Nesterov both $μ$ and $L$ are assumed known -- an assumption that is seldom satisfied in practice. We propose to incorporate mechanisms to estimate locally sufficient $μ$ and $L$ during the iterations. The mechanisms also allow for the application to non-strongly convex functions. We discuss the iteration complexity of several first-order methods, including the proposed algorithm, and we use a 3D tomography problem to compare the performance of these methods. The results show that for ill-conditioned problems solved to high accuracy, the proposed method significantly outperforms state-of-the-art first-order methods, as also suggested by theoretical results.
△ Less
Submitted 18 May, 2011;
originally announced May 2011.
-
Toward optimal X-ray flux utilization in breast CT
Authors:
Jakob H. Joergensen,
Per Christian Hansen,
Emil Y. Sidky,
Ingrid S. Reiser,
Xiaochuan Pan
Abstract:
A realistic computer-simulation of a breast computed tomography (CT) system and subject is constructed. The model is used to investigate the optimal number of views for the scan given a fixed total X-ray fluence. The reconstruction algorithm is based on accurate solution to a constrained, TV-minimization problem, which has received much interest recently for sparse-view CT data.
A realistic computer-simulation of a breast computed tomography (CT) system and subject is constructed. The model is used to investigate the optimal number of views for the scan given a fixed total X-ray fluence. The reconstruction algorithm is based on accurate solution to a constrained, TV-minimization problem, which has received much interest recently for sparse-view CT data.
△ Less
Submitted 8 April, 2011;
originally announced April 2011.