skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A semi-analytical solution and AI-based reconstruction algorithms for magnetic particle tracking
Magnetic particle tracking is a recently developed technology that can measure the translation and rotation of a particle in an opaque environment like a turbidity flow and fluidized-bed flow. The trajectory reconstruction usually relies on numerical optimization or filtering, which involve artificial parameters or thresholds. Existing analytical reconstruction algorithms have certain limitations and usually depend on the gradient of the magnetic field, which is not easy to measure accurately in many applications. This paper discusses a new semi-analytical solution and the related reconstruction algorithm. The new method can be used for an arbitrary sensor arrangement. To reduce the measurement uncertainty in practical applications, deep neural network (DNN)-based models are developed to denoise the reconstructed trajectory. Compared to traditional approaches such as wavelet-based filtering, the DNN-based denoisers are more accurate in the position reconstruction. However, they often over-smooth the velocity signal, and a hybrid method that combines the wavelet and DNN model provides a more accurate velocity reconstruction. All the DNN-based and wavelet methods perform well in the orientation reconstruction.  more » « less
Award ID(s):
1934300 1944187
PAR ID:
10311900
Author(s) / Creator(s):
; ; ;
Editor(s):
Chen, Chi-Hua
Date Published:
Journal Name:
PLOS ONE
Volume:
16
Issue:
7
ISSN:
1932-6203
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Magnetic particle tracking (MPT) is a recently developed non-invasive measurement technique that has gained popularity for studying dense particulate or granular flows. This method involves tracking the trajectory of a magnetically labeled particle, the field of which is modeled as a dipole. The nature of this method allows it to be used in opaque environments, which can be highly beneficial for the measurement of dense particle dynamics. However, since the magnetic field of the particle used is weak, the signal-to-noise ratio is usually low. The noise from the measuring devices contaminates the reconstruction of the magnetic tracer’s trajectory. A filter is then needed to reduce the noise in the final trajectory results. In this work, we present a neural network-based framework for MPT trajectory reconstruction and filtering, which yields accurate results and operates at very high speed. The reconstruction derived from this framework is compared to the state-of-the-art extended Kalman filter-based reconstruction. 
    more » « less
  2. A bstract The MicroBooNE liquid argon time projection chamber located at Fermilab is a neutrino experiment dedicated to the study of short-baseline oscillations, the measurements of neutrino cross sections in liquid argon, and to the research and development of this novel detector technology. Accurate and precise measurements of calorimetry are essential to the event reconstruction and are achieved by leveraging the TPC to measure deposited energy per unit length along the particle trajectory, with mm resolution. We describe the non-uniform calorimetric reconstruction performance in the detector, showing dependence on the angle of the particle trajectory. Such non-uniform reconstruction directly affects the performance of the particle identification algorithms which infer particle type from calorimetric measurements. This work presents a new particle identification method which accounts for and effectively addresses such non-uniformity. The newly developed method shows improved performance compared to previous algorithms, illustrated by a 93.7% proton selection efficiency and a 10% muon mis-identification rate, with a fairly loose selection of tracks performed on beam data. The performance is further demonstrated by identifying exclusive final states in ν μ CC interactions. While developed using MicroBooNE data and simulation, this method is easily applicable to future LArTPC experiments, such as SBND, ICARUS, and DUNE. 
    more » « less
  3. A new method for fluid–structure interaction (FSI) diagnostics to simultaneously capture time-resolved three-dimensional, three-component (3D3C) velocity fields and structural deformations using a single light field camera is presented. A light field camera encodes both spatial and angular information of light rays collected by a conventional imaging lens that allows for the 3D reconstruction of a scene from a single image. Building upon this capability, a light field fluid–structure interaction (LF FSI) methodology is developed with a focus on experimental scenarios with low optical access. Proper orthogonal decomposition (POD) is used to separate particle and surface information contained in the same image. A correlation-based depth estimation technique is introduced to reconstruct instantaneous surface positions from the disparity between angular perspectives and conventional particle image velocimetry (PIV) is used for flow field reconstruction. Validation of the methodology is achieved using synthetic images of simultaneously moving flat plates and a vortex ring with a small increase in uncertainty under ~0.5 microlenses observed in both flow and structure measurement compared to independent measurements. The method is experimentally verified using a flat plate translating along the camera’s optical axis in a flow field with varying particle concentrations. Finally, simultaneous reconstructions of the flow field and surface shape around a flexible membrane are presented, with the surface reconstruction further validated using simultaneously captured stereo images. The findings indicate that the LF FSI methodology provides a new capability to simultaneously measure large-scale flow characteristics and structural deformations using a single camera. 
    more » « less
  4. Abstract High fidelity near-wall velocity measurements in wall bounded fluid flows continue to pose a challenge and the resulting limitations on available experimental data cloud our understanding of the near-wall velocity behavior in turbulent boundary layers. One of the challenges is the spatial averaging and limited spatial resolution inherent to cross-correlation-based particle image velocimetry (PIV) methods. To circumvent this difficulty, we implement an explicit no-slip boundary condition in a wavelet-based optical flow velocimetry (wOFV) method. It is found that the no-slip boundary condition on the velocity field can be implemented in wOFV by transforming the constraint to the wavelet domain through a series of algebraic linear transformations, which are formulated in terms of the known wavelet filter matrices, and then satisfying the resulting constraint on the wavelet coefficients using constrained optimization for the optical flow functional minimization. The developed method is then used to study the classical problem of a turbulent channel flow using synthetic data from a direct numerical simulation (DNS) and experimental particle image data from a zero pressure gradient, high Reynolds number turbulent boundary layer. The results obtained by successfully implementing the no-slip boundary condition are compared to velocity measurements from wOFV without the no-slip condition and to a commercial PIV code, using the velocity from the DNS as ground truth. It is found that wOFV with the no-slip condition successfully resolves the near-wall profile with enhanced accuracy compared to the other velocimetry methods, as well as other derived quantities such as wall shear and turbulent intensity, without sacrificing accuracy away from the wall, leading to state of the art measurements in the y + < 1 region of the turbulent boundary layer when applied to experimental particle images. 
    more » « less
  5. Abstract. Numerical models are a powerful tool for investigating the dynamic processes in the interior of the Earth and other planets, but the reliability and predictive power of these discretized models depends on the numerical method as well as an accurate representation of material properties in space and time. In the specific context of geodynamic models, particle methods have been applied extensively because of their suitability for advection-dominated processes and have been used in applications such as tracking the composition of solid rock and melt in the Earth's mantle, fluids in lithospheric- and crustal-scale models, light elements in the liquid core, and deformation properties like accumulated finite strain or mineral grain size, along with many applications outside the Earth sciences. There have been significant benchmarking efforts to measure the accuracy and convergence behavior of particle methods, but these efforts have largely been limited to instantaneous solutions, or time-dependent models without analytical solutions. As a consequence, there is little understanding about the interplay of particle advection errors and errors introduced in the solution of the underlying transient, nonlinear flow equations. To address these limitations, we present two new dynamic benchmarks for transient Stokes flow with analytical solutions that allow us to quantify the accuracy of various advection methods in nonlinear flow. We use these benchmarks to measure the accuracy of our particle algorithm as implemented in the ASPECT geodynamic modeling software against commonly employed field methods and analytical solutions. In particular, we quantify if an algorithm that is higher-order accurate in time will allow for better overall model accuracy and verify that our algorithm reaches its intended optimal convergence rate. We then document that the observed increased accuracy of higher-order algorithms matters for geodynamic applications with an example of modeling small-scale convection underneath an oceanic plate and show that the predicted place and time of onset of small-scale convection depends significantly on the chosen particle advection method. Descriptions and implementations of our benchmarks are openly available and can be used to verify other advection algorithms. The availability of accurate, scalable, and efficient particle methods as part of the widely used open-source code ASPECT will allow geodynamicists to investigate complex time-dependent geodynamic processes such as elastic deformation, anisotropic fabric development, melt generation and migration, and grain damage. 
    more » « less