Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel methods”

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 1,513 records · Page 84

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MO-C-17A-10: Comparison of Dose Deformable Accumulation by Using Parallel and Serial Approaches

Purpose: The uncertainty of dose accumulation over multiple CT datasets with deformable fusion may have significant impact on clinical decisions. In this study, we investigate the difference of two dose summation approaches involving deformable fusion. Methods: Five patients, four external beam and one brachytherapy(BT), were chosen for the study. The BT patient was treated with CT-based HDR. The CT image sets acquired in the imageguidance process (8-11 CTs/patient) were used to determine the dose delivered to the four external beam patients. (prostate, pelvis, lung and head and neck). For the HDR patient (cervix), five CT image sets and the corresponding BT plans were used. In total 44 CT datasets and RT dose/plans were imported into the image fusion software MiM (6.0.4) for analysis.For each of the five clinical cases, the dose from each fraction was accumulated into the primary CT dataset by using both Parallel and Serial approaches. The dose-volume histogram (DVH) for CTV and selected organs-at-risks (OAR) were generated. The D95(CTV), OAR(mean) and OAR(max) for the four external beam cases the D90(CTV), and the max dose to bladder and rectum for the BT case were compared. Results: For the four external beam patients, the difference in D95(CTV) were <1.2% PD between the parallel and the serial approaches. The differences of the OAR(mean) and the OAR(max ) range from 0 to 3.7% and <1% PD respectively. For the HDR patient, the dose difference for D90 is 11% PD while that of the max dose to bladder and rectum were 11.5% and 23.3% respectively. Conclusion: For external beam treatments, the parallel and serial approaches have <5% difference probably because tumor volume and OAR have less changes from fraction to fraction. For the brachytherapy case, >10% dose difference between the two approaches was observed as significant volume changes of tumor and OAR were observed among treatment fractions.

60 APPLIED LIFE SCIENCES↗

CO 2 -induced chemo-mechanical alteration in reservoir rocks assessed via batch reaction experiments and scratch testing

Here, the injection of carbon dioxide (CO 2 ) into geological formations results in a chemical re-equilibration between the mineral assemblage and the pore fluid, with ensuing mineral dissolution and re-precipitation. Hence, target rock formations may exhibit changes of mechanical and petrophysical properties due to CO 2 exposure. We conducted batch reaction experiments with Entrada Sandstone and Summerville Siltstone exposed to de-ionized water and synthetic brine under reservoir pressure (9–10 MPa) and temperature (80°C) for up to four weeks. Samples originate from the Crystal Geyser field site, where a naturally occurring CO 2 seepage alters portions of these geologic formations. We conducted micro-scratch tests on rock samples without alteration, altered under laboratory conditions, and naturally altered over geologic time. Scratch toughness and hardness decrease as a function of exposure time and water salinity up to 52% in the case of Entrada and 87% in the case of Summerville after CO 2 -induced alteration in the laboratory. Imaging of altered cores with SEM-EDS and X-ray microCT methods show dissolution of carbonate and silica cements and matrix accompanied by minor dissolution of Fe-oxides, clays, and other silicates. Parallel experiments using powdered samples confirm that dissolution of carbonate and silica are the primary reactions. The batch reaction experiments in the autoclave utilize a high fluid to rock volume ratio and represent an end member of possible alteration associated with CO 2 storage systems. These types of tests serve as a pre-screening tool to identify the susceptibility of rock facies to CO 2 -related chemical-mechanical alteration during long-term CO 2 storage.

54 ENVIRONMENTAL SCIENCES↗

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Parallel Implementation of Nonadditive Gaussian Process Potentials for Monte Carlo Simulations

A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deciphering chemical order/disorder and material properties at the single-atom level

Perfect crystals are rare in nature. Real materials often contain crystal defects and chemical order/disorder such as grain boundaries, dislocations, interfaces, surface reconstructions and point defects. Such disruption in periodicity strongly affects material properties and functionality. Despite rapid development of quantitative material characterization methods, correlating three-dimensional (3D) atomic arrangements of chemical order/disorder and crystal defects with material properties remains a challenge. On a parallel front, quantum mechanics calculations such as density functional theory (DFT) have progressed from the modelling of ideal bulk systems to modelling ‘real’ materials with dopants, dislocations, grain boundaries and interfaces; but these calculations rely heavily on average atomic models extracted from crystallography. To improve the predictive power of first-principles calculations, there is a pressing need to use atomic coordinates of real systems beyond average crystallographic measurements. Here we determine the 3D coordinates of 6,569 iron and 16,627 platinum atoms in an iron-platinum nanoparticle, and correlate chemical order/disorder and crystal defects with material properties at the single-atom level. We identify rich structural variety with unprecedented 3D detail including atomic composition, grain boundaries, anti-phase boundaries, anti-site point defects and swap defects. We show that the experimentally measured coordinates and chemical species with 22 picometre precision can be used as direct input for DFT calculations of material properties such as atomic spin and orbital magnetic moments and local magnetocrystalline anisotropy. The work presented here combines 3D atomic structure determination of crystal defects with DFT calculations, which is expected to advance our understanding of structure–property relationships at the fundamental level.

36 MATERIALS SCIENCE↗

Application of gas-coupled laser acoustic detection to gelatins and underwater sensing

Gas-coupled Laser Acoustic Detection (GCLAD) has been used as a method to sense ultrasound waves in materials without contact of the material surface. To sense the waveform, a laser beam is directed parallel to the material surface and displaced or deflected when the radiated waveform traverses the beam. We present recent tests that demonstrate the potential of using this technique for detecting ultrasound in gelatin phantoms and in water. As opposed to interferometric detection, GCLAD operates independently of the optical surface properties of the material. This allows the technique to be used in cases where the material is transparent or semi-transparent. We present results on sensing ultrasound in gelatin phantoms that are used to mimic biological materials. As with air-coupled transducers, the frequency response of GCLAD at high frequencies is limited by the high attenuation of ultrasound in air. In contrast, water has a much lower attenuation. Here we demonstrate the use of a GCLAD-like system in water, measuring the directivity response at 1 MHz and sensing waveforms with higher frequency content.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Did the GPU obfuscate the load imbalance in my MPI simulation?

The current proliferation of GPU-based HPC systems necessitates a method for assessing the performance of simulations on heterogeneous machines. The addition of GPUs to a system adds multiple hierarchical levels of parallelism to the node architecture. In this paper, we demonstrate that the traditional load imbalance metric is insufficient for capturing the load imbalance on GPU-based machines, since it treats the GPU as a monolithic entity and ignores the internal parallelism. We propose a new hierarchical metric that improves the correlation of measured performance and application workload by up to 20.61%. Using our metric for determining application load instead of the traditional metric as the input for the load balancing algorithm reduces the residual load imbalance by up to 4× in our application.

Eberius, David↗

SU-F-BRE-05: Development and Evaluation of a Real-Time Robotic 6D Quality Assurance Phantom

Purpose: A 6 degree-of-freedom robotic phantom capable of reproducing dynamic tumor motion in 6D was designed to more effectively match solid tumor movements throughout pre-treatment scanning and radiation therapy. With the abundance of optical and x-ray 6D real-time tumor tracking methodologies clinically available, and the substantial dosimetric consequences of failing to consider tumor rotation as well as translation, this work presents the development and evaluation of a 6D instrument with the facility to improve quality assurance. Methods: An in-house designed and built 6D robotic motion phantom was constructed following the so-called Stewart-Gough parallel kinematics platform archetype. The device was then controlled using an inverse kinematics formulation, and precise movements in all six degrees of freedom (X, Y, Z, pitch, roll, and yaw) as well as previously obtained cranial motion, were effectively executed. The robotic phantom movements were verified using a 15 fps 6D infrared marker tracking system (Polaris, NDI), and quantitatively compared to the input trajectory. Thus, the accuracy and repeatability of 6D motion was investigated and the phantom performance was characterized. Results: Evaluation of the 6D platform demonstrated translational RMSE values of 0.196 mm, 0.260 mm, and 0.101 mm over 20 mm in X and Y and 10 mm in Z, respectively, and rotational RMSE values of 0.068 degrees, 0.0611 degrees, and 0.095 degrees over 10 degrees of pitch, roll, and yaw, respectively. The robotic stage also effectively performed controlled 6D motions, as well as reproduced cranial trajectories over 15 minutes, with a maximal RMSE of 0.044 mm translationally and 0.036 degrees rotationally. Conclusion: This 6D robotic phantom has proven to be accurate under clinical standards and capable of reproducing tumor motion in 6D. Consequently, such a robotics device has the potential to serve as a more effective system for IGRT QA that involves both translational and rotational dimensions. Research was partially funded by NIH Grant T32 EB002103-21 from NIBIB. Contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIBIB or NIH.

60 APPLIED LIFE SCIENCES↗

Dose enhancement in radiotherapy of small lung tumors using inline magnetic fields: A Monte Carlo based planning study

Purpose: To report on significant dose enhancement effects caused by magnetic fields aligned parallel to 6 MV photon beam radiotherapy of small lung tumors. Findings are applicable to future inline MRI-guided radiotherapy systems. Methods: A total of eight clinical lung tumor cases were recalculated using Monte Carlo methods, and external magnetic fields of 0.5, 1.0, and 3 T were included to observe the impact on dose to the planning target volume (PTV) and gross tumor volume (GTV). Three plans were 6 MV 3D-CRT plans while 6 were 6 MV IMRT. The GTV’s ranged from 0.8 to 16 cm{sup 3}, while the PTV’s ranged from 1 to 59 cm{sup 3}. In addition, the dose changes in a 30 cm diameter cylindrical water phantom were investigated for small beams. The central 20 cm of this phantom contained either water or lung density insert. Results: For single beams, an inline magnetic field of 1 T has a small impact in lung dose distributions by reducing the lateral scatter of secondary electrons, resulting in a small dose increase along the beam. Superposition of multiple small beams leads to significant dose enhancements. Clinically, this process occurs in the lung tissue typically surrounding the GTV, resulting in increases to the D{sub 98%} (PTV). Two isolated tumors with very small PTVs (3 and 6 cm{sup 3}) showed increases in D{sub 98%} of 23% and 22%. Larger PTVs of 13, 26, and 59 cm{sup 3} had increases of 9%, 6%, and 4%, describing a natural fall-off in enhancement with increasing PTV size. However, three PTVs bounded to the lung wall showed no significant increase, due to lack of dose enhancement in the denser PTV volume. In general, at 0.5 T, the GTV mean dose enhancement is around 60% lower than that at 1 T, while at 3 T, it is 5%–60% higher than 1 T. Conclusions: Monte Carlo methods have described significant and predictable dose enhancement effects in small lung tumor plans for 6 MV radiotherapy when an external inline magnetic field is included. Results of this study indicate that future clinical inline MRI-guided radiotherapy systems will be able to deliver a dosimetrically superior treatment to small (PTV < 15 cm{sup 3}), isolated lung tumors over non-MRI-Linac systems. This increased efficacy coincides with the reimbursement in the United States of lung CT screening and the likely rapid growth in the number of patients with small lung tumors to be treated with radiotherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-T-165: Daily QA Analysis for Spot Scanning Beamline

Purpose: The dosimetric results of our daily quality assurance over the last 8 years for discrete pencil beam scanning proton therapy will be presented. Methods: To perform the dosimetric checks, a multi-ion chamber detector is used, which consists of an array of 5 single parallel plate ion chambers that are aligned as a cross separated by 10cm each. The Tracker is snapped into a jig, which is placed on the tabletop. Different amounts of Solid Water buildup are added to shift the dose distribution. The dosimetric checks consist of 3 parts: position check, range check and volume dose check. Results: The average deviation of all position-check data were 0.2±1.3%. For the range check, the average deviation was 0.1%±1.2%, which also corresponds to a range stability of better than 1 mm over all measurements. The volumetric dose output readings were all within ±1% with the exception of 2 occasions when the cable to the dose monitor was being repaired. Conclusion: Morning QA using the Tracker device gives very stable dosimetric readings but is also sensitive to mechanical and output changes in the proton therapy delivery system.

07 ISOTOPE AND RADIATION SOURCES↗

Design of ultrasonic attenuation sensor with focused transmitter for density measurements of a slurry in a large steel pipeline

To design an ultrasonic sensor to measure the attenuation and density of a slurry carried by a large steel pipeline (diameter up to 70 cm) is the goal of this research. The pitch-catch attenuation sensor, placed in a small section of the pipeline, contains a send unit with a focused transducer that focuses the ultrasound to a small region of the receive unit on the opposite wall. The focused transducer consists of a section of a sphere (base ~12 cm) on the outer side of the send unit and a 500 kHz piezoelectric shell of PZT5A epoxied to it. The Rayleigh surface integral is used to calculate the pressure in steel and in water (slurry). An incremental method to plot the paths of ultrasonic rays shows that the rays focus where expected. Further, there is a region where the parallel rays are perpendicular to the wall of the receive unit. Designs for pipeline diameters of 25 cm and 71 cm show that the pressure in water at the receive transducer is about 17 times that for a pitch-catch system using 5 cm diameter disk transducers. The enhanced signal increases the sensitivity of the measurements and improves the signal-to-noise ratio.

Greenwood, Margaret Stautberg↗

Specific features of magnetic phase diagram of an MnSi helimagnet

Magnetization curves of MnSi single crystals have been measured in a range of temperatures T = 5.5-35 K and magnetic field strengths H {<=} 11 kOe for H Parallel-To [1 1 1], [0 0 1], and [1 1 0]. Special attention has been paid to the temperature interval near T{sub N} = 28.8 K, where MnSi exhibits a transition to the state with a long-period helical magnetic structure. Some new features in the magnetic behavior of MnSi have been found. In particular, in an intermediate temperature region above the transition (28.8 K = T{sub N} {<=} T < 31.5 K), the dM(H)/dH curves exhibit anomalies that are not characteristic of the typical paramagnetic state. It is established that the line of the characteristic field H*(T) of this anomaly is a natural extrapolation of the temperature dependence of the field of the transition from a conical phase to an induced ferromagnetic phase observed at T < T{sub N}. It is concluded that the properties of MnSi in the indicated intermediate region are related to and governed by those of the conical phase (rather than of the A phase). Based on these data, magnetic phase diagrams of MnSi for H Parallel-To [1 1 1], [0 0 1], and [1 1 0] are plotted and compared to diagrams obtained earlier by other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements and analytics of wide-area file transfers over dedicated connections

Distributed scientific and big-data computations are becoming increasingly dependent on access to remote files. Wide-area file transfers are supported by two basic schemes: (i) application-level tools, such as GridFTP, that provide transport services between file systems housed at geographically separated sites, and (ii) file systems mounted over wide-area networks, using mechanisms such as LNet routers that make them transparently available. In both cases, the file transfer performance critically depends on the configuration consisting of host, file, IO, and disk subsystems, which are complex by themselves, as well as on their complex compositions implemented using buffers and IO-network data transitions. We present extensive file transfer rate measurements collected over dedicated 10 Gbps connections with 0-366 ms round-trip times, using GridFTP and XDD file transfer tools, and Lustre file system extended over wide-area networks using LNet routers. Our test configurations are composed of: three types of host systems; XFS, Lustre, and ext3 file systems; and Ethernet and SONET wide-area connections. We present analytics based on the convexity-concavity of throughput profiles which provide insights into throughput and its superior or inferior trend compared to linear interpolations. We propose the utilization-concavity coefficient, a scalar metric that characterizes the overall performance of any file transfer method consisting of specific configuration and scheme. Our results enable performance optimizations by highlighting the significant roles of (i) buffer sizes and parallelism in GridFTP and XDD, and (ii) buffer utilization and credit mechanism in LNet routers.

Rao, Nageswara↗

DASSA v0.0.1

DASSA (Parallel DAS Data Storage and Analysis) provides a data storage engine and analysis engine for data from distributed acoustic sensing and other methods. It supports various data analysis operations, from FFT , signal filter, cross-correlation, compression, stacking etc.

Dong, Bin↗

BM3DORNL

BM3DORNL is a high-performance, open-source library for removing streak and ring artifacts from computed-tomography (CT) data, developed for neutron imaging at Oak Ridge National Laboratory's Spallation Neutron Source (VENUS beamline) and applicable to X-ray CT as well. Ring artifacts — concentric rings in reconstructed slices caused by detector pixel-to-pixel response non-uniformities — appear as vertical streaks in the sinogram and degrade both image quality and quantitative analysis. BM3DORNL operates in the sinogram domain using an adaptation of the BM3D (block-matching and 3D collaborative filtering) algorithm (Dabov et al., 2007). It provides a dedicated streak-removal mode, a true multi-scale BM3D variant (after Mäkinen et al., 2021) that suppresses wide streaks single-scale methods miss, and an alternative Fourier–SVD method (~2.6× faster) combining FFT-based energy detection with rank-1 SVD. The computationally intensive core is implemented in Rust with parallel (Rayon) block matching, integral-image pre-screening, and optimized transforms, and is exposed through a simple Python API (with an optional GUI) so it integrates directly into existing tomography reconstruction pipelines. It processes both 2D sinograms and 3D sinogram stacks, is pip-installable for Linux and macOS, and is documented at https://bm3dornl.readthedocs.io.

Zhang, Chen [Oak Ridge National Laboratory (ORNL),↗