Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reconstruction algorithms”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

GNET2: an R package for constructing gene regulatory networks from transcriptomic data

Abstract Motivation The Gene Network Estimation Tool (GNET) is designed to build gene regulatory networks (GRNs) from transcriptomic gene expression data with a probabilistic graphical model. The data preprocessing, model construction and visualization modules of the original GNET software were developed on different programming platforms, which were inconvenient for users to deploy and use. Results Here, we present GNET2, an improved implementation of GNET as an integrated R package. GNET2 provides more flexibility for parameter initialization and regulatory module construction based on the core iterative modeling process of the original algorithm. The data exchange interface of GNET2 is handled within an R session automatically. Given the growing demand for regulatory network reconstruction from transcriptomic data, GNET2 offers a convenient option for GRN inference on large datasets. Availability and implementation The source code of GNET2 is available at https://github.com/jianlin-cheng/GNET2. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Closed-Form Approximation of the Total Variation Proximal Operator

Total variation (TV) is a widely used function for regularizing imaging inverse problems that is particularly appropriate for images whose underlying structure is piecewise constant. TV regularized optimization problems are typically solved using proximal methods, but the way in which they are applied is constrained by the absence of a closed-form expression for the proximal operator of the TV function. A closed-form approximation of the TV proximal operator has previously been proposed, but its accuracy was not theoretically explored in detail. Here, we address this gap by making several new theoretical contributions, proving that the approximation leads to a proximal operator of some convex function, it is equivalent to a gradient descent step on a smoothed version of TV, and that its error can be fully characterized and controlled with its scaling parameter. We experimentally validate our theoretical results on image denoising and sparse-view computed tomography (CT) image reconstruction.

97 MATHEMATICS AND COMPUTING↗

Improvement of the NOvA Near Detector Event Reconstruction and Primary Vertexing through the Application of Machine Learning Methods

The purpose of this work is to examine the application of a deep learning model in event reconstruction of neutrino interactions. The challenges faced in event reconstruction include the placement of an accurate primary neutrino interaction vertex which is used to support the particle track and prong algorithms. The result of accurate primary vertex ensures all particles involved in a neutrino interaction are included. We propose a regression-based Convolutional Neural Network (CNN) method to predict the primary vertex of a particle interaction. We show that with raw two-dimensional pixel map views as input, the regression-based CNN can predict the primary vertex in all three coordinates. This work is applied as part of the NOvA (NuMI Off-axis $\nu_e$ Appearance) near detector reconstruction efforts. The primary vertex predicted by the regression-based CNN model shows promising results for future applications. This deep learning method can be extended to secondary vertexing through a Kernel Density Estimate algorithm discussed in this work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of the production of charm jets tagged with D0 mesons in pp collisions at $\sqrt{s}$ = 5.02 and 13 TeV

The measurement of the production of charm jets, identified by the presence of a D 0 meson in the jet constituents, is presented in proton–proton collisions at centre-of-mass energies of $\sqrt{s}$= 5.02 and 13 TeV with the ALICE detector at the CERN LHC. The D 0 mesons were reconstructed from their hadronic decay D 0 → K – π + and the respective charge conjugate. Jets were reconstructed from D 0 -meson candidates and charged particles using the anti-k T algorithm, in the jet transverse momentum range 5 < p T,chjet < 50 GeV/c, pseudorapidity |η jet | < 0.9 – R, and with the jet resolution parameters R = 0.2, 0.4, 0.6. The distribution of the jet momentum fraction carried by a D 0 meson along the jet axis $\left({z}_{\Big\Vert}^{\textrm{ch}}\right)$ was measured in the range 0.4 < ${z}_{\Big\Vert}^{\textrm{ch}}$ < 1.0 in four ranges of the jet transverse momentum. Comparisons of results for different collision energies and jet resolution parameters are also presented. The measurements are compared to predictions from Monte Carlo event generators based on leading-order and next-to-leading-order perturbative quantum chromodynamics calculations. A generally good description of the main features of the data is obtained in spite of a few discrepancies at low p T,chjet . Measurements were also done for R = 0.3 at $\sqrt{s}$= 5.02 and are shown along with their comparisons to theoretical predictions in an appendix to this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigating charm quark energy loss in medium with the nuclear modification factor of D 0 -tagged jets

The nuclear modification factor R AA of charm jets, identified by the presence of a D 0 meson among the jet constituents, has been measured for the first time in Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}} = 5.02$ TeV with the ALICE detector at the LHC. The D 0 mesons and their charge conjugates are reconstructed from the hadronic decay D 0 → K – π + . Jets are reconstructed from D 0 -meson candidates and charged particles using the anti-k T algorithm with jet resolution parameter R = 0.3, in the jet transverse momentum (p T ) range 5 < $p^{ch}_{T}$$^{jet}$ < 50 GeV/c and pseudorapidity |η ch jet| < 0.6. A hint of reduced suppression in the charm-jet R AA is observed in comparison to inclusive jets in central Pb–Pb collisions with a significance of about 2σ in 20 < $p^{ch}_{T}$$^{jet}$ < 50 GeV/c, suggesting the in-medium energy loss to depend on both the difference between quark and gluon coupling strength (Casimir colour-charge effect) and quark mass (dead-cone effect). The data are compared with model calculations that include mass effects in the in-medium energy loss. Several state-of-the-art models are consistent with the data, with the LIDO model providing the best description of the data in the common kinematic range of inclusive and D 0 -tagged jets, highlighting the role of mass effects in interpreting the results.

Acharya, S. (ORCID:0000000292135329)↗

Deep learning based x-ray spectrometer for high repetition rate characterization of betatron radiation

Betatron radiation produced from a laser-wakefield accelerator is a broadband, hard x-ray (>1 keV) source that has been used in a variety of applications in medicine, engineering, and fundamental science. Further development and optimization of stable, high repetition rate (HRR) (>1 Hz) betatron sources will provide a means to extend their application base to include single-shot dynamical measurements of ultrafast processes or dense materials. Recent advances in laser technology used in such experiments have enabled increases in shot-rate and system stability, providing improved statistical analysis and detailed parameter scans. However, unique challenges exist at high repetition rate, where data throughput and source optimization are now limited by diagnostic acquisition rates and analysis. Here, we present the development of a machine-learning algorithm for the real-time analysis of betatron radiation. We report on the fielding of this deep learning algorithm for online source characterization at the Institut National de la Recherche Scientifique's Advanced Laser Light Source. By fine-tuning an algorithm originally trained on a fully synthetic dataset using a subset of experimental data, the algorithm can predict the betatron critical energy with a percent error of 7.2 % with a reconstruction time of 1.5 ms, providing a valuable tool for real-time, multi-objective optimization at HRR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved Localization Precision and Angular Resolution of a Cylindrical, Time-Encoded Imaging System From Adaptive Detector Movements

To the first order, the localization precision and angular resolution of a cylindrical, time-encoded imaging (c-TEI) system is governed by the geometry of the system. Improving either measure requires increasing the mask radius or decreasing the detector diameter, both of which are undesirable. Here, we propose an alternative option of repositioning the detector within the mask to increase the detector-to-mask distance in the direction of a source, thereby improving the localization precision and angular resolution in that direction. Since the detector-to-mask distance only increases for a small portion of the field of view (FOV), we propose implementing adaptive imaging where one leverages data collected during the measurement to optimize the system configuration. This article utilizes both simulations and experiments to set upper bounds on the potential gain from adaptive detector movements for one and two sources in the FOV. When only one source is present, adaptive detector movements can improve the localization precision and angular resolution by 20% for a source at 90 cm and by 32% for a far-field source. When two sources are present, adaptive detector movements can improve localization precision and angular resolution by up to 50% for sources that are ~10° apart (90 cm from the system). We experimentally verify these results through maximum likelihood estimation of the source position(s) and image reconstruction of point sources that are close together. As a demonstration of an adaptive imaging algorithm, we image a complex arrangement of special nuclear material at the Zero Power Physics Reactor facility at Idaho National Laboratory.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Development of NDE/NDT Tools for High-Volume & High-Speed Inspection of CFRP Structures in Automotive Manufacturing

Main advantages of the air-coupled ultrasound testing (ACUT) and electromagnetic testing (EMT) techniques for NDE of CFRP composites were non-contact sensing, scalability for high-speed inspection, cost-effectiveness, and non-hazardous operation. Despite these advantages, no systems that would satisfy the project requirements were commercially available. Hence, one of the major efforts of the Michigan State University (MSU) team at the initial stage of the project was to close this technological gap by developing, optimizing, and validating array sensors that would provide sufficient sensitivity, spatial coverage, and resolution for robust defect detection. Optimization of the ACUT and EMT sensor designs was performed using experimentally validated finite element models. Initial experiments using array probes were conducted on relatively flat CFRP samples. In parallel, the MSU team designed and assembled a portable platform with two robotic arms. The robots were equipped with newly designed sensors that enabled high-speed NDE of curved CFRP parts. Presently, the developed robotic platform can be used as a demo/template NDE system, which is easily adaptable to manufacturing environments and in-line NDE. The ACUT NDE system developed by the MSU team used a high-power 4-channel pulser receiver for parallel data acquisition. The array probes were designed by stacking commercially available ACUT transducers, which operated in the frequency range between 100 kHz and 500 kHz. MSU optimized the excitation procedure and developed wave focusing cones so as to reduce the crosstalk between the transducers and to provide higher pulse repletion frequency (PRF). The through-transmission (TT) and single-side access (SSA) inspection modes were successfully implemented. In the TT-ACUT, structural defects in CFRP were detected by passing ultrasonic waves through the test part. Hence, the ACUT transmitters and receivers needed to be placed on the opposite sides of the test part. In the SSA-ACUT, guided waves (GW) were excited in the test part using the transmitters and were sensed by the receivers from the same side. Multi-channel TT-ACUT and SSA-ACUT provided high-speed NDE, and were successfully validated on CFRP test samples with interlaminar delaminations and other embedded defects The EM techniques developed by the MSU team included: 1) eddy current testing (ECT), 2) capacitive imaging (CI) and hybrid dual-mode imaging. In ECT, structural damage was detected in CFRP using coils sensor arrays. In ECT, the excitation magnetic field is generated by passing an alternating current through a coil, which is placed above the test sample. The excitation field penetrates the conductive sample and induces the eddy currents in its transect. In turn, the eddy currents generate the reaction field, which affects the total field sensed by a coil. Hence, the presence of structural flaws will alter the eddy current flow and the picked-up signal. ECT is mostly sensitive to local changes of the electric conductivity of the test sample, and CFRPs are mostly conductive in the direction of carbon fibers. Hence, ECT was well suited for the detection of fiber damage/fiber irregularities. The MSU team developed printed circuit boards (PCB) with coil sensor arrays optimized for NDE of CFRP. Unlike most commercial probes designed for ECT of metallic structures, the MSU array probes were designed for operation in [1-10] MHz frequency range, which was optimal for low-conductive CFRP. Multiple sensing topologies (coil groups excitation/sensing arrangements) were implemented and successfully validated. Capacitive Imaging (CI) technique developed by MSU was complementary to ECT. In contrast to ECT, which was sensitive to local changes of the electrical conductivity, the CI was sensitive to local changes of the dielectric constant. Therefore, CI could provide information about matrix damage/matrix irregularities in CFRP. The MSU CI sensor arrays were made of multiple circular or rectangular open-plate capacitors printed on PCB. Sensors of this type are not commercially available. In addition to ECT and CI, the MSU team developed a hybrid (dual-mode) inductive/capacitive measurement technique that synergistically combined the benefits of inductive and capacitive sensing for rapid NDE of fiber reinforced polymer (FRP) composite structures. Fiber damage and fiber irregularities in FRPs were detected by configuring hybrid sensors as coil sensors. Similarly, matrix damage, matrix irregularities and interlaminar delaminations were detected by configuring hybrid sensors as capacitive sensors. ECT and CI were performed sequentially by means of electronic switching. Hence, eliminating the need for mounting two separate sensor arrays on the probe. Portable robotic platform was developed by MSU for multi-technique high-speed NDE of CFRP test parts. The platform had two 6-axis robots, which enabled inspection of curved parts in approximately a 6×6×6 ft 3 active scan area. On the software side, the MSU team integrated scripts for NDE hardware control with scripts for robot motion control. MSU also implemented automated path planning for the robots, reconstruction of part’s surfaces via stereovision, 3D rendering of inspection data, and image processing algorithms for enhanced defect detection. Automotive composite parts manufactured by Plasan Composites from Phase I were used to validate the ACUT and EMT techniques on representative testbeds. Among those parts were three X-braces for a Dodge Viper, one composite calibration plaque with known defects at known locations, and four other test sections, including sections from a front splitter, a corner section from a composite hood, and a high-pressure RTM panel made using non crimp fabric. Other test samples included CFRP and GFRP calibration plates with fiber/matrix defects fabricated at MSU/CVRC.

36 MATERIALS SCIENCE↗

Structures via Reasoning - Applying AI to Cryo Electron Microscopy to Reveal Structural Variability

There have been breakthroughs of latest cryo electron microscopy (cryo-EM) data analysis algorithms to classify cryo-EM image data. However, most of these cryo-EM reconstruction methods have focused on classifying distinctly different biomolecule structures. Here, we present our approaches of deep learning to differentiate homologous structures that are distinguishable only with inner morphological differences. We succeeded supervised classification of these subtly different homologues. However, we could not differentiate them with unsupervised methods. Here we discuss what further approaches are likely needed for successful unsupervised classification.

47 OTHER INSTRUMENTATION↗

Calibrating LArPix for TinyTPC

LArTPCs provide sensitivity to GeV signals, such as accelerator neutrinos and part of the supernova neutrino spectrum. TinyTPC is a LArTPC test stand for R&D of LAr doping to expand the reach of LArTPCs down to the 1-10 MeV range, which would substantially enhance the flagship analyses of experiments like DUNE, while enabling low energy analyses. We aim to dope LAr with Xe and pho- tosensitive dopants to expand the LArTPC range by converting hard-to-detect scintillation light to efficiently detected ionization charge. A critical element of the data analysis in TinyTPC is calibrating the readout. This poster will cover the calibration of TinyTPC, a pixelated liquid argon detector, where we find the distance a muon travels through each pixel. We then calculate the energy loss of muons traveling through the detector from cosmic data. We can reconstruct the path of the particles through the TPC using a density-based clustering algorithm designed to sort straight cosmic muon tracks from low energy radioactive decay curled paths and electronic noise

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimizing the hit finding algorithm for liquid argon TPC neutrino detectors using parallel architectures

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming even more difficult as the detectors increase in size to reach their physics goals. Liquid argon time projection chamber (LArTPC) neutrino experiments are expected to grow in the next decade to have 100 times more wires than in currently operating experiments, and modernization of LArTPC reconstruction code, including parallelization both at data- and instruction-level, will help to mitigate this challenge. The LArTPC hit finding algorithm is used across multiple experiments through a common software framework. In this paper we discuss a parallel implementation of this algorithm. Using a standalone setup we find speedup factors of two times from vectorization and 30–100 times from multi-threading on Intel architectures. The new version has been incorporated back into the framework so that it can be used by experiments. On a serial execution, the integrated version is about 10 times faster than the previous one and, once parallelization is enabled, more speedups comparable to the standalone program are achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Different CT slice thickness and contrast‐enhancement phase in radiomics models on the differential performance of lung adenocarcinoma

Abstract Background To investigate the effects of computed tomography (CT) reconstruction slice thickness and contrast‐enhancement phase on the differential diagnosis performance of radiomic signature in lung adenocarcinoma. Methods A total of 187 patients who had been pathologically confirmed with lung adenocarcinoma and nonadenocarcinoma were divided into a training cohort ( n = 149) and validation cohort ( n = 38). All the patients underwent contrast‐enhanced CT and the images were reconstructed with different slice thickness. The radiomic features were extracted from different slice thickness and scan phase. The logistic regression (LR) algorithm was used to build a machine learning model for each group. The area under the curve (AUC) obtained from the receiver operating characteristic (ROC) curve and DeLong test was used to evaluate its discriminating performance. Results Finally, 34 image features and five semantic features were selected to establish a radiomics model. Based on the three contrast‐enhanced CT phases and four reconstruction slice thickness, 12 groups of radiomics models showed good discrimination ability with the AUCs range from 0.9287 to 0.9631, sensitivity range from 0.8349 to 0.9083, specificity range from 0.825 to 0.925 in the training group. Similar results were observed in the validation group. However, there was no statistical significance between the different CT scan phase groups and different slice thickness ( p > 0.05). Conclusions The radiomic analysis of contrast‐enhanced CT can be used for the differential diagnosis of lung adenocarcinoma. Moreover, different slice thickness and contrast‐enhanced scan phase did not affect the discriminating ability in the radiomics models.

Wang, Yang↗

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗

3D track reconstruction of low-energy electrons in the MIGDAL low pressure optical time projection chamber

Here, we demonstrate three-dimensional track reconstruction of electrons in a low pressure (50 Torr) optical TPC consisting of two glass GEMs with an ITO strip readout in CF 4 and CF 4 /Ar mixtures. The reconstructed tracks show a variety of event topologies, including short tracks from photoelectrons induced by 55 Fe 5.9 keV X-rays and long tracks from gamma ray interactions and beta decays. Algorithms for event identification and track ridge detection are discussed as well as multiple methods for integrating information from the camera image and ITO waveforms with the goal of full 3D reconstruction of the track.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Robust Resilient Signal Reconstruction under Adversarial Attacks

We consider the problem of signal reconstruction for a system under sparse signal corruption by a malicious agent. The reconstruction problem follows the standard error coding problem that has been studied extensively in the literature. We include a new challenge of robust estimation of the attack support. The problem is then cast as a constrained optimization problem merging promising techniques in the area of deep learning and estimation theory. A pruning algorithm is developed to reduce the "false positive" uncertainty of data-driven attack localization results, thereby improving the probability of correct signal reconstruction. Sufficient conditions for the correct reconstruction and the associated reconstruction error bounds are obtained for both exact and inexact attack support estimation. Moreover, a simulation of a water distribution system is presented to validate the proposed techniques.

Robust, Signal reconstruction, Resilient estimator↗

Fast wide-field upconversion luminescence lifetime thermometry enabled by single-shot compressed ultrahigh-speed imaging

Photoluminescence lifetime imaging of upconverting nanoparticles is increasingly featured in recent progress in optical thermometry. Despite remarkable advances in photoluminescent temperature indicators, existing optical instruments lack the ability of wide-field photoluminescence lifetime imaging in real time, thus falling short in dynamic temperature mapping. Here, we report video-rate upconversion temperature sensing in wide field using single-shot photoluminescence lifetime imaging thermometry (SPLIT). Developed from a compressed-sensing ultrahigh-speed imaging paradigm, SPLIT first records wide-field luminescence intensity decay compressively in two views in a single exposure. Then, an algorithm, built upon the plug-and-play alternating direction method of multipliers, is used to reconstruct the video, from which the extracted lifetime distribution is converted to a temperature map. Using the core/shell NaGdF 4 :Er 3+ ,Yb 3+ /NaGdF 4 upconverting nanoparticles as the lifetime-based temperature indicators, we apply SPLIT in longitudinal wide-field temperature monitoring beneath a thin scattering medium. SPLIT also enables video-rate temperature mapping of a moving biological sample at single-cell resolution.

47 OTHER INSTRUMENTATION↗

Determination of jet calibration and energy resolution in proton–proton collisions at $\sqrt{s} = 8~\hbox {TeV}$ using the ATLAS detector

The jet energy scale, jet energy resolution, and their systematic uncertainties are measured for jets reconstructed with the ATLAS detector in 2012 using proton–proton data produced at a centre-of-mass energy of 8 TeV with an integrated luminosity of 20fb -1 . Jets are reconstructed from clusters of energy depositions in the ATLAS calorimeters using the anti-k t algorithm. A jet calibration scheme is applied in multiple steps, each addressing specific effects including mitigation of contributions from additional proton–proton collisions, loss of energy in dead material, calorimeter non-compensation, angular biases and other global jet effects. The final calibration step uses several in situ techniques and corrects for residual effects not captured by the initial calibration. Furthermore, these analyses measure both the jet energy scale and resolution by exploiting the transverse momentum balance in γ + jet, Z + jet, dijet, and multijet events. A statistical combination of these measurements is performed. In the central detector region, the derived calibration has a precision better than 1% for jets with transverse momentum 150 GeV < p T < 1500 GeV, and the relative energy resolution is (8.4 ± 0.6)% for p T =100 GeV and (23 ± 2)% for p T =20 GeV. The calibration scheme for jets with radius parameter R = 1.0, for which jets receive a dedicated calibration of the jet mass, is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗