Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three-dimensional reconstruction”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Curvature Induced Modifications of Chirality and Magnetic Configuration in Perpendicular Films

Designing curvature in three-dimensional (3D) magnetic nanostructures enables controlled manipulation of local energy landscapes, allowing for the modification of noncollinear spin textures relevant for next-generation spintronic devices. In this study, we experimentally investigate 3D magnetization textures in a Co/Pd multilayer film, exhibiting strong perpendicular magnetic anisotropy (PMA), deposited onto curved Cu nanowire meshes with diameters as small as 50 nm and lengths of several microns. Utilizing magnetic soft X-ray nanotomography, we achieve reconstructions of 3D magnetic domain patterns at approximately 30 nm spatial resolution. This approach provides detailed information on both the orientation and magnitude of magnetization within the film. Our results reveal that interfacial anisotropy in the Co/Pd multilayers drives the magnetization toward the local surface normal. In contrast to typical labyrinth domains observed in planar films, the presence of curved nanowires significantly alters the domain structure, with domains preferentially aligning along the nanowire axis in close proximity, while adopting random orientations farther away. We report direct experimental observation of a curvature-induced Dzyaloshinskii-Moriya interaction (DMI), which is quantified to be approximately one-third of the intrinsic DMI in Co/Pd stacks. The curvature induced DMI enhances stability of Néel-type domain walls. These experimental observations are further supported by micromagnetic simulations. Altogether, our findings demonstrate that introducing curvature into magnetic nanostructures provides a powerful strategy for tailoring complex magnetic behaviors, paving the way for the design of advanced 3D racetrack memory and neuromorphic computing devices.

Raftrey, David↗
DOE OSTI2024

Laboratory-Based Micro-X-ray Computed Tomography of Energy Materials at Idaho National Laboratory

Abstract The Idaho National Laboratory (INL) has implemented laboratory-based micro-X-ray computed tomography in a laboratory equipped for the examination of highly radioactive samples. This capability provides nondestructive three-dimensional volumetric information on samples to inform subsequent traditional destructive examinations as well as real-world inputs for high-fidelity scientific modeling. Samples can be imaged with spatial resolutions ranging from several hundred nm/voxel up to ~ 100 µm/voxel. The best usable spatial resolution achieved to date is 384 nm/voxel with this instrument, while the highest radiological dose rate of a sample imaged is ~ 60 R/h β/γ on contact. Advanced data analysis, including custom tomographic reconstruction and segmentation methods, have also been developed to support this capability. In addition to traditional digital X-ray radiography and tomography, this instrument is also able to visualize in situ tensile and compression testing as well as perform diffraction contrast tomography. This work describes the X-ray computed tomography post-irradiation examination capabilities at INL, as well as detailing a variety of applications this instrument has examined.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultra-High-Speed X-ray Tomography: Bridging the Gaps to See the Unknown

Rapid, time-varying, three-dimensional physics underpin numerous engineering challenges. Often, these physics occur within opaque environments, internal to a component, severely limiting applicable diagnostics. Development of novel diagnostics is necessary to understand and predict transient three-dimensional (3D) phenomena within opaque environments. This report highlights progress in four key areas leading to advancements in high-speed X-ray radiography and tomography. The first area is enabling MHz-rate imaging of energetics at the Advanced Photon Source at Argonne National Laboratory. The second is modeling a high-flux, rotating-anode X-ray source to understand the heat loads on the anode. The third effort was to develop a novel reconstruction algorithm that is validated by ground experimental tomography data and synthetic tomography data. The fourth is the development of a novel approach to two-color X-ray imaging.

42 ENGINEERING↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental study of multiple-orientation muon tomography with image optimization in sparse data environments

Due to the high penetrating power of cosmic-ray muons, they can be used to probe very thick and dense objects. As muons are charged particles, they can be tracked by ionization detectors, determining the position and direction of the muons. With detectors on either side of an object to measure particle direction change, scattering information within the object can be found. This can be used to produce a scattering-intensity image within the object related to density and atomic number. Such imaging is typically performed with a single detector-object orientation, taking advantage of the more intense downward flux of muons, producing planar imaging with some depth-of-field information in the third dimension. Several simulation studies were published with multiorientation tomography, which can form a three-dimensional representation faster than a single-orientation view. In this study, experimental muon-scatter-based tomography was performed using a concrete filled steel drum with several different metal wedges inside, with the drum between detector planes. Data were collected from different detector-object orientations by rotating the steel drum. The data collected from each orientation were combined using two different tomographic methods. A traditional inverse Radon transform approach used for computed tomography and a combination of multiple depth-of-field reconstructions were applied to the data. As cosmic-ray muon flux imaging is rate limited, the imaging techniques were compared for sparse data. Using the combined depth-of-field reconstruction technique, fewer detector-object orientations were needed to reconstruct images that could be used to differentiate the metal wedges.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗
DOE OSTI2025

A Gigaparsec-scale Hydrodynamic Volume Reconstructed with Deep Learning

The next generation of spectroscopic surveys will map the large-scale structure of the Universe at high redshifts (2 ≤ z ≤ 5) using millions of quasar spectra, enabling major advances in constraining both the standard cosmological model and its extensions. Robust cosmological analyses of these data sets require numerical simulations that both cover gigaparsec volumes and resolve features on ∼10 kpc scales and smaller. However, running such large-volume, high-resolution hydrodynamic simulations is computationally prohibitive. We present a generative deep learning model that enhances a low-resolution, gigaparsec-scale (960 h −1 Mpc) hydrodynamic simulation using a smaller (80 h −1 Mpc) high-resolution input hydrodynamic simulation as training data. The resulting enhanced simulation reproduces the line-of-sight power spectrum to within ∼10% and the three-dimensional power spectrum at the ∼20% level at intermediate to small scales (k ≲ 2 h Mpc −1 ). Our method shows strong promise for producing realistic simulations for cosmological analyses with current surveys such as the Dark Energy Spectroscopic Instrument and upcoming next-generation experiments, but further improvements are needed to accurately recover the large-scale modes. We publicly release the enhanced hydrodynamic simulation, along with a halo catalog from a companion N-body dark matter simulation to support the calibration of data analysis pipelines for these large-scale surveys.

Convolutional neural networks↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Laser-driven high-resolution MeV x-ray tomography

The need for high-resolution MeV x-ray tomography to observe the three-dimensional structure of dense, large-sized objects is rapidly increasing for the non-destructive evaluation of critical additively manufactured parts, national security, and other applications. We report a demonstration of high-resolution MeV computed tomography of a dense, large object with a laser-driven x-ray source. A record detector-limited MeV radiograph resolution of < 200 µm as determined with the Bennett approximation of the point spread function was achieved by irradiating millimeter-thick tungsten targets with 300 TW femtosecond laser pulses at a 0.5 Hz repetition rate. A tungsten alloy step wedge spectrometer indicates that the peak of the x-ray emission is between 1 and 2 MeV, with an endpoint energy of 19 MeV. To illustrate the radiographic imaging capability of the system, a tomographic reconstruction of a nickel superalloy turbine blade ( maximum $ρr$ = 139 g/cm 2 ) with sub-millimeter resolution was performed using 2160 individual radiographs. The small x-ray source size opens the prospect of extremely high-resolution tomographs of large, dense objects. This laser-driven approach has major advantages for non-destructive evaluation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Using multimodal X-ray computed tomography to advance 3D petrography: A non-destructive investigation of olivine inside a carbonaceous chondrite

Rocks form in three dimensions through time and studying them provides information from inside dynamic systems we cannot otherwise observe. Yet how we typically access the interior of the rocks themselves to gain that information may limit our understanding and influence how we reconstruct the processes that formed them. Here, we demonstrate combined non-destructive 3D X-ray imaging techniques that produce quantitative densitometric and crystallographic maps of entire individual grains inside a rock. Olivine grains throughout a sample of the carbonaceous chondrite Northwest Africa (NWA) 11346 were each characterized by size, shape, composition, zoning intensity, and crystallographic orientation. The addition of 3D crystallographic mapping to calibrated 3D densitometric analysis—used to calculate chemical composition—demonstrates a fully non-destructive petrographic method and provides unique insight. For instance, in our case, using crystallographic data to delineate individual grains and then measuring the 3D size, shape, and composition of each distinguishes variably reset relict grains from those later crystallized after a melting event. Intersection in a 2D slice could not have led to this interpretation because the integration of three-dimensional size, rounding, composition, location, and crystallographic orientation measured from each grain forms the key patterns. Multimodal laboratory X-ray imaging has strong potential to advance 3D petrography.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nano-laminography with a transmission X-ray microscope

Nano-laminography combines the penetrating power of hard X-rays with a tilted rotational geometry to deliver high-resolution, three-dimensional images of laterally extended, flat specimens that are otherwise incompatible with, or difficult to image using, conventional nano-tomography. In this work, we demonstrate a full-field, X-ray nano-laminography system implemented with the transmission X-ray microscope at beamline 32-ID of the upgraded Advanced Photon Source at Argonne National Laboratory, USA. By rotating the sample around an axis inclined by 20° to the incident beam, the technique minimizes the long optical path lengths that would otherwise generate excessive artifacts when planar samples are imaged edge-on. The efficiency of the technique is demonstrated with 50 nm spatial resolution and minute-scale temporal resolution 3D imaging of a planar integrated circuit sample and targeted imaging of an individual particle within a powder sample, where mounting procedures are typically challenging in regular nano-tomography. The sample mounting strategy, data acquisition, and reconstruction method will also be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Comparison of automated chemical-guided segmentation and human annotation of soil organic matter in X-ray microcomputed tomography imaging in contrasted soil types

Soil organic matter (OM) formation and persistence is strongly influenced by the spatial distribution of organic substrates and microscale soil heterogeneity by dictating OM accessibility to microorganisms. However, traditional size and/or density fractionation techniques disrupt aggregate architecture, eliminating spatial information needed to fully understand intra-aggregate OM distribution. To quantify three-dimensional OM spatial distribution and automate segmentation in X-ray microcomputed tomography (µCT) imaging without human annotation bias, we developed an iodine gas vapor (I2) based staining workflow that eliminates labor-intensive manual annotation while maintaining segmentation accuracy, using aggregates from four taxonomically diverse soils (Xerofluvent, Haploxeroll Sphagnofibrist, Palehumult) with an 8-fold range of soil organic carbon. Human annotation of 10 µCT slices by the experienced and inexperienced annotators resulted in variations up to 3% in the Dice similarity coefficient (DSC), reflecting a degree of inherent subjectivity of manual labeling. Such inconsistencies are expected to compound as the number of manually annotated slices increases. Dual-energy µCT imaging at 33.1 keV (below the iodine (I) K-edge) and 33.2 keV (above the I K-edge) was used to resolve aggregate microstructure following I2 staining. The automated image subtraction pipeline identified OM regions by the I Kedge induced brightness increases, achieving DSC values of 0.58–0.83 relative to an experienced annotator. Sensitivity analyses revealed that the reconstruction alpha value—optimized via the open-source tool TomocuPy—and the 3D registration slice count were the primary determinants of accuracy, providing a novel benchmark for dual-energy soil imaging. The pipeline without GPU acceleration achieved 9.6 to 43.2 times faster than manual annotation. Using GPU-accelerated image post-processing and affine transformation matrices, the pipeline successfully segmented OM elements for large-scale datasets (3232×3232 pixel, 2048 slices) within ~5200 s from raw file acquisition to segmented output. The high-throughput approach enables the quantification of OM spatial distribution across diverse and heterogeneous soil.

Soil microbial biomass↗
DOE OSTI2025

3D pore structure reconstruction of fluff region in top of EBR-II irradiated fuel by X-ray micro-computed tomography

Metallic fuels comprised of U-Zr binary and U-Pu-Zr ternary alloys are a promising fuel candidate for sodium-cooled fast reactors and possess many advantages over ceramic fuels such as higher uranium density and improved thermal conductivity. Fluff is a highly porous structure forming at the top of metallic fuel slugs during irradiation and was extensively observed in past sodium-cooled fast reactor (SFR) experiments such as Idaho National Laboratory’s Experimental Breeder Reactor II (EBR-II). Previous investigations have shown that fluff is rich in fissile atoms and hence could have neutronics implications, but its formation mechanisms are not understood. In this work, two cube-shaped lift outs were taken from the fluff and bulk fuel regions of a high burnup U-19Pu-10Zr composition EBR-II fuel slug and characterized by micro computed tomography (Micro-CT) and energy dispersive x-ray spectroscopy (EDS). This was used to characterize differences in microscale porosity and composition between fluff and bulk fuel, and to support or oppose postulated fluff formation mechanisms. Surface EDS showed minimal differences between the elemental compositions of the fluff and fuel cubes. However, three-dimensional (3D) pore analysis from the tomography scans showed significant differences between the fluff and fuel cube pore size distributions at the microscale level. The high porosity and irregular pore distribution seen in macroscale fluff was also present within the micro-scale fluff cube, having a much higher porosity than the fuel cube which was concentrated in a few large pore networks. In conclusion, these results indicate support for a creep driven mechanism of fluff formation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Generalized parton distributions from the pseudodistribution approach on the lattice

Generalized parton distributions (GPDs) are key quantities for the description of a hadron’s three-dimensional structure. They are the current focus of all areas of hadronic physics—phenomenological, experimental and theoretical, including lattice QCD. Synergies between these areas are desirable and essential to achieve precise quantification and understanding of the structure of, particularly, nucleons, as the basic ingredients of matter. In this paper, we investigate, for the first time, the numerical implementation of the pseudodistribution approach for the extraction of zero-skewness GPDs for unpolarized quarks. Pseudodistributions are Euclidean parton correlators computable in lattice QCD that can be perturbatively matched to the light-cone parton distributions of interest. Although they are closely related to the quasidistributions and come from the same lattice-extracted matrix elements, they are, however, subject to different systematic effects. We use the data previously utilized for quasi-GPDs and extend it with other momentum transfers and nucleon boosts, in particular a higher one ( P 3 = 1.67 GeV ) with eightfold larger statistics than the largest one used for quasidistributions ( P 3 = 1.25 GeV ). We renormalize the matrix elements with a ratio scheme and match the resulting Ioffe time distributions to the light cone in coordinate space. The matched distributions are then used to reconstruct the x dependence with a fitting . We investigate some systematic effects related to this procedure, and we also compare the results with the ones obtained in the framework of quasi-GPDs. Our final results involve the invariant four-momentum transfer squared ( − t ) dependence of the flavor nonsinglet ( u − d ) H and E GPDs. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗
DOE OSTI2024

Novel Instrumentation for Advanced Quantum Chromodynamics Studies with Flavor Sensitivity

The true nature of ordinary matter has always fascinated humankind, whose imagination has been pointed to the right direction since the 5th century BC with the first known atomic theory by Leucippus and Democritus. They described the matter as formed by small, invisible, indivisible, and eternal particles: the atoms. A long investigation followed them, including contributions from several of history?s greatest philosophers and scientists. The progress led to the understanding that atoms have an internal structure consisting of a central nucleus composed of protons and neutrons, surrounded by a cloud of electrons. About 50 years ago, scientists discovered that nucleons also have an internal structure. Through the Deep Inelastic Scattering (DIS) experiments, on which a lepton is scattered off a nucleon, it was possible to access the distribution of partons inside the nucleon along the longitudinal direction defined by the hard-scale probe. It was a crucial step in developing Quantum Chromodynamics (QCD), currently the best theory describing the subatomic matter. Ultimately, DIS was recognized as a limited tool to investigate the nucleon structure because it is mainly sensitive to observables at the probe?s energy scale. In recent decades, theoretical developments have allowed defining soft-scale observables by measuring more complex processes. Reactions such as Semi-Inclusive Deep Inelastic Scattering (SIDIS) with single or di-hadron production, can provide information about the momentum of parton inside the nucleon. In particular, polarized SIDIS allows access to the Transverse-Momentum-Dependent (TMD) parton distributions describing the correlations of the quark and gluon degrees of freedom in transverse momentum and spin. This study can provide three-dimensional imaging of nucleon structure and inner dynamics, as long as the hadronic component in the final state can be measured. Identifying the species of the produced particles provides information on the quark flavor and assumes a crucial role in this context. The thesis describes the work carried out by the author in preparing novel instrumentation for particle identification in the final state of DIS experiments, to support modern analysis of the parton dynamics within the basic confined object (nucleon) with flavor sensitivity. The work has been focused on the CLAS12 spectrometer in operation at Jefferson Laboratory (JLab), Newport News, VA, USA, and the ePIC experiment in preparation for the future Electron-Ion Collider (EIC) at Brookhaven National Laboratory (BNL), Long Island, NY, USA. CLAS12 is a fixed-target experiment using a 12 GeV beam of polarized electrons scattering off polarized nuclear targets. In 2018 and 2022, two Ring Imaging Cherenkov (RICH) detector modules were installed to improve the ?±/K± separation in the high-momentum region (3÷8 GeV/c) of the experiment. The author contributed to the assembly, installation, and com- missioning of the second RICH module, to the efficiency studies of the first module, and to the first SIDIS analysis with high-momentum kaons identified by the RICH. These studies led to the observation of the first spin asymmetry with high-momentum kaons from the CLAS12 experiment, and an initial assessment of the systematic error associated with the hadron identification. ePIC will be the first experiment at the future EIC, the new collider designed to expand the frontiers of QCD. The author was involved in the development of the dual-radiator Ring Imaging Cherenkov (dRICH). This detector will interpolate the measurements of the Cherenkov angles of photons produced by relativistic particles crossing two different radiators to identify charged hadrons. The author contributed to the studies conducted with the dRICH prototype, developing the reconstruction and analysis software and simulations, characterizing the aerogel radiator samples, and being responsible for the tracking system used during the beam tests. The performance obtained by the prototype has been progressively improved and the results are now comparable with the expectation derived from the simulation and satisfy the requirements of the experiment. In Chapter 1, the SIDIS theory is briefly introduced, outlying the connection with the experimental measurement. Chapter 2 includes the description of the CLAS12 RICH, its assembly and installation, and the efficiency study. In Chapter 3, the analysis of the Beam-Spin Asymmetry associated with SIDIS kaons is described. Chapter 4 describes the contribution to the dRICH for EIC, starting from a description of the ePIC experiment and focusing on the studies performed for the dRICH prototype. The conclusions of this work are summarized at the end.

Vallarino, Simone↗
DOE OSTI2024

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗