Search NASA⌕ Search

SEARCH · Search NASA

Results for “Image coding”

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,585 records · Page 88

Overview on Shielding Analyses for the VENUS Instrument at SNS

VENUS, a world‐class versatile neutron imaging instrument, is under construction and is expected to be completed and ready to start commissioning in 2023. The range of cold to epithermal neutrons at SNS will give users of VENUS access to novel imaging methods, as well as to significantly improved existing methods. The instrument is being built on beam line 10 at Spallation Neutron Source (SNS) First Target Station (FTS) facing a decoupled poisoned hydrogen moderator. Instrument design, which includes optics, Front-End components and instrument enclosure components and shape started over 10 years ago and during this time had significant changes. All changes were supported by neutronics analyses to provide adequate shielding for both Front-End and instrument enclosure.Final optics design will provide the field of view (FOV) at the detector position (where the image is formed) to be as high as 0.20 by 0.20 m. Both the beam, which contains a large fraction of high-energy neutrons, and the desire for a large footprint on the detector are challenges for shielding design and budget, because the driving cost for the instrument is the beam line and enclosure shielding. For cost reason VENUS baseline design is optimized with respect to both the instrument cave footprint and its wall thickness, and Front-End shielding is tailored along the beam. Analyses are performed with the Monte Carlo particle transport code MCNPX version 2.7.0 to make choices on materials, thicknesses and configurations of the shielding. Numerous calculations are performed to meet space and coast constrain and to satisfy instrument physics needs and to comply with radiation protection requirements.

Gallmeier, Franz X.↗

Using extended MHD to explore lasers as a trigger for x-pinches

X-pinches have been shown to be a source of extremely intense x-ray emissions useful for diagnosing plasma dynamics and imaging biological objects. The most striking feature of an x-pinch is the hotspot, the point source from where all the x rays come from. Unfortunately, the exact timing and location of the hotspot are still unpredictable. Since an x-pinch hotspot forms from instabilities (like an m = 0 mode), we will computationally explore whether we can use a high-power laser to control the timing and location of these instabilities. Our goal is to reduce the temporal and spatial jitter associated with the x-ray burst. Using an extended MHD (XMHD) code, we explore the non-relativistic instability generation using a current profile of a 250 kA linear transformer driver system and laser characteristics of SLAC's Matter in Extreme Conditions laboratory. Our results include both laser-penetration results using a boundary-defined electromagnetic (EM) wave and instability results from a power-deposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

TOPOS: A new topometric patient positioning and tracking system for radiation therapy based on structured white light

Purpose: A patient positioning system for radiation therapy based on structured white light and using off-the-shelf hardware components for flexibility and cost-effectiveness has been developed in house. Increased accuracy, patient comfort, abandonment of any skin marks, accelerated workflow, objective reading/recording, better usability and robust sensor design, compared to other positioning approaches, were the main goals of this work. Another aim was the application of a 6 degrees of freedom tracking system working without dose deposition. Methods: Two optical sensors are the main parts of the TOPOS Registered-Sign system (Topometrical Positioning, cyberTECHNOLOGIES, Germany). The components: cameras, projectors, and computers are commercial off-the-shelf products, allowing for low production costs. The black/white cameras of the prototype are capable of taking up to 240 frames per second (resolution: 640 Multiplication-Sign 488 pixels). The projector has a resolution of 1024 Multiplication-Sign 768 and a refresh rate of 120 Hz. The patient's body surface is measured continuously and registered to a reference surface, providing a transformation to superimpose the patient's surface to the reference (planning CT) surface as best as possible. The execution of the calculated transformation provides the correct patient position before the treatment starts. Due to the high-speed acquisition of the surfaces, a surveillance of the patient's (respiration) motion during treatment is also accomplished. The accuracy of the system was determined using a male mannequin. Two treatment sites were evaluated: one simulating a head and neck treatment and the other simulating a thoracic wall treatment. The mannequin was moved to predefined positions, and shift vectors given by the surface registration were evaluated. Additionally manual positioning using a color-coding system was evaluated. Results: Two prototypes have been developed, each allowing a continuous high density scan of a 500 Multiplication-Sign 500 Multiplication-Sign 400 mm{sup 3} (L Multiplication-Sign W Multiplication-Sign D) large volume with a refresh rate of 10 Hz (extendible to 20 Hz for a single sensor system). Surface and position correction display, as well as respiratory motion, is shown in real-time (delay < 200 ms) using present graphical hardware acceleration. For an intuitive view of the patient's misalignment, a fast surface registration algorithm has been developed and tested and a real-time color-coding technique is proposed and verified that allows the user to easily verify the position of the patient. Using first the surface registration and then the color coding the best results were obtained: for the head and neck case, the mean difference between the actual zero position and the final match was 0.1 {+-} 0.4, -0.2 {+-} 0.7, and -0.1 {+-} 0.3 mm in vertical, longitudinal, and lateral direction. For the thoracic case, the mean differences were 0.3 {+-} 0.5, -0.6 {+-} 1.9, 0.0 {+-} 0.4 mm. Conclusions: The presented system copes with the increasing demand for more accurate patient positioning due to more precise irradiation technologies and minimizes the preparation times for correct patient alignment, therefore optimizing the treatment workflow. Moreover, TOPOS is a versatile and cost effective image guided radiation therapy device. It allows an objective rating of the patient's position before and during the irradiation and could also be used for respiratory gating or tracking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with GEANT 4

Purpose: This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and GEANT4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field. Methods: A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm{sup 2}. Dose scoring was performed using a 3D grid comprising 1 mm{sup 3} voxels. The production thresholds were equivalent for both codes. Results were analyzed based upon a voxel by voxel dose difference between the two codes and also using a volumetric gamma analysis. Results: Comparisons were drawn from central axis depth doses, cross beam profiles, and isodose contours. Both in the presence and absence of a 1.5 T static magnetic field the relative differences in doses scored along the beam central axis were less than 1% for the homogeneous water phantom and all results matched within a maximum of ±2% for heterogeneous phantoms. Volumetric gamma analysis indicated that more than 99% of the examined volume passed gamma criteria of 2%—2 mm (dose difference and distance to agreement, respectively). These criteria were chosen because the minimum primary statistical uncertainty in dose scoring voxels was 0.5%. The presence of the magnetic field affects the dose at the interface depending upon the density of the material on either sides of the interface. This effect varies with the field size. For example, at the water-lung interface a 33.94% increase in dose was observed (relative to the D{sub max}), by both GPUMCD and GEANT4 for the field size of 2 × 2 cm{sup 2} (compared to no B-field case), which increased to 47.83% for the field size of 5 × 5 cm{sup 2} in the presence of the magnetic field. Similarly, at the lung-water interface, the dose decreased by 19.21% (relative to D{sub max}) for a field size of 2 × 2 cm{sup 2} and by 30.01% for 5 × 5 cm{sup 2} field size. For more complex combinations of materials the dose deposition also becomes more complex. Conclusions: The GPUMCD algorithm showed good agreement against GEANT4 both in the presence and absence of a 1.5 T external magnetic field. The application of 1.5 T magnetic field significantly alters the dose at the interfaces by either increasing or decreasing the dose depending upon the density of the material on either side of the interfaces.

61 RADIATION PROTECTION AND DOSIMETRY↗

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

Numerical Simulation of Isothermal Flow Across Slant Five-Tube Bundle with Spectral Element Method Code Nek5000

The helical coil steam generator (HCSG) is a specific type of shell-and-tube heat exchanger known for having a higher heat transfer coefficient than many other designs. For this reason, they are considered in small modular reactor and high temperature reactor designs. Investigation of flow behavior in HCSGs is important for better design. In this paper we present our study for modeling the coolant flow in the primary side of the HCSG. We used Nek5000, an open source, high-order spectral element computational fluid dynamics code developed in Argonne National Laboratory. Nek5000 accepts only hexahedral mesh, which makes the meshing process for the complicated HCSG geometry very challenging. A tetrahedral-to-hexahedral meshing strategy was applied to bypass the geometric complexities. In this study large eddy simulation (LES) was performed at the Reynolds number of 9000 based on the inlet velocity and the tube diameter. The employed subgrid-scale model for LES relies on explicit filtering. First- and second-order statistics were compared to available experimental data. We report overall velocity and turbulent kinetic energy showed good agreement with particle image velocimetry data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The resolved outflow from 3C 48

We investigate the properties of the high-velocity outflow driven by the young radio jet of 3C 48, a compact-steep-spectrum source. We use the Space Telescope Imaging Spectrograph on board the Hubble Space Telecope to obtain (1) low-resolution UV and optical spectra and (2) multi-slit medium-resolution spectra of the ionized outflow. With supporting data from ground-based spectrographs, we are able to accurately measure the ratios of diagnostic emission lines such as [O III] λ5007, [O III] λ3727, [N II] λ6548, Hα, Hβ, [Ne V] λ3425, and [Ne III] λ3869. We fit the observed emission-line ratios using a range of ionization models, powered by active galactic nucleus (AGN) radiation and shocks, produced by the MAPPINGS code. We have determined that AGN radiation is likely the dominant ionization source. The outflow's density is estimated to be in the range n = 10{sup 3}-10{sup 4} cm{sup –3}, the mass is ∼6 × 10{sup 6} M {sub ☉}, and the metallicity is likely equal to or higher than solar. Compared with the typical outflows associated with more evolved radio jets, this young outflow is denser, less massive, and more metal rich. Multi-slit observations allow us to construct a two-dimensional velocity map of the outflow that shows a wide range of velocities with distinct velocity components, suggesting a wide-angle clumpy outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

Analyzing line-of-sight selection biases in galaxy-scale strong lensing with external convergence and shear

The upcoming Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will dramatically increase the number of strong gravitational lensing systems, requiring precise modeling of line-of-sight (LOS) effects to mitigate biases in lensing observations and cosmological inferences. We develop a method to construct joint distributions of external convergence (κ ext ) and shear (γ ext ) for strong lensing LOS by aggregating large-scale structure simulations with high-resolution halo renderings and non-linear correction. Our approach captures both smooth background matter and perturbations from halos, enabling accurate modeling of LOS effects. Here, we apply non-linear LOS corrections to κ ext and γ ext that address the non-additive lensing effects caused by objects along the LOS in strong lensing. We find that, with a minimum image separation of 1.0'', non-linear LOS correction due to the presence of a dominant deflector slightly increases the ratio of quadruple to double lenses; this non-linear LOS correction also introduces systematic biases of ∼ 0.1% for galaxy-AGN lens in the inferred Hubble constant (H 0 ) if not accounted for. We also observe a 0.66% bias for galaxy-galaxy lenses on H 0 , and even larger biases 1.02% for galaxy-AGN systems if LOS effects are not accounted for. These results highlight the importance of LOS for precision cosmology. The publicly available code and datasets provide tools for incorporating LOS effects in future analyses.

Hubble constant↗

Determining the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5

Purpose: Recent advances in the imaging of {sup 90}Y using positron emission tomography (PET) and improved uncertainty in the branching ratio for the internal pair production component of {sup 90}Y decay allow for a more accurate determination of the activity distribution of {sup 90}Y microspheres within a patient. This improved activity distribution can be convolved with the dose kernel of {sup 90}Y to calculate the dose distribution within a patient. This work investigates the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5 and compares the results of these two transport codes. Methods: Monte Carlo simulations were performed with egsnrc and mcnp5 to calculate the dose rate at multiple radial distances around various {sup 90}Y sources. Point source simulations were completed with mcnp5 to determine the optimal electron transport settings for this work. After determining the optimal settings, point source simulations were completed using egsnrc (user code edknrc) and mcnp5 in water and liver [as defined by the International Commission on Radiation Units and Measurements (ICRU) Report 44]. The results were compared to ICRU Report 72 reference data. Point source simulations were also completed in water with a density of 1.06 g{center_dot}cm{sup -3} to evaluate the effect of the density of the surrounding material. Glass and resin microsphere simulations were performed with average and maximum diameter and density values (based on values given in the literature) in water and in liver. The results were compared to point source simulation results using the same transport code and in the same surrounding material. All simulations had statistical uncertainties less than 1%. Results: The optimal transport settings in mcnp5 for this work included using the energy-and step-specific algorithm (DBCN 17J 2) and ESTEP set to 10. These settings were used for all subsequent simulations with mcnp5. The point source simulations in water for both egsnrc and mcnp5 were found to agree within 2% of the ICRU 72 reference data over the investigated range. Point source simulations in liver had large differences relative to ICRU 72, approaching -60% near the maximum range of {sup 90}Y. These differences are mostly attributed to the difference in density between water (1.0 g{center_dot}cm{sup -3}) and liver (1.06 g{center_dot}cm{sup -3}). Glass and resin microsphere simulations showed a slight decrease in the dose rate near the maximum range of {sup 90}Y relative to the point source simulations. The largest relative differences were approximately -4.2% and -2.8% for the glass and resin microspheres, respectively. Agreement between the egsnrc and mcnp5 simulations results was generally good. Conclusions: The presence of the microsphere material causes slight differences in the {sup 90}Y dose kernel compared to those calculated with point sources. Large differences were seen between simulations in water and those in liver. For the most accurate calculation of the dose distribution, the density of the patient's liver should be accounted for in the calculation of the dose kernel. Lastly, due to the need to determine the optimal transport settings with mcnp5, electron transport with this code should be used with caution.

61 RADIATION PROTECTION AND DOSIMETRY↗

High-resolution time-resolved experiments on mixing and entrainment of buoyant jets in stratified environments

Fluid jets interacting with a stratified layer play an important role in the safety of several reactor designs. In the containment of nuclear power plants, fluid jets dominate the transport and mixing of gaseous species and consequent hydrogen distribution in case of a severe accident. The mixing phenomena in the containment are driven by buoyant high-momentum injections (jets) and low momentum injection plumes. Mixing near the postulated break is initially dominated by high flow velocities. Plumes with moderate flow velocities are instead relevant in the break compartment during the long-term pressurization phase, or in any of the apertures between two connected compartments if the mass flows are sufficiently high and the density differences between efflux and ambient are sufficiently low. Phenomena of interest include free plumes (as produced by the efflux from the break compartment in a larger room or directly from a break flow), wall plumes (such those produced by low mass flows through inter-compartment apertures), and propagating stratification fronts in the ambient (for any stably stratified conditions). These phenomena have been highly ranked about nuclear reactor design, especially regarding of safety protocols. During a Pressurized Thermal Shock (PTS) scenario, the interaction between the cold ECCS injection plume and the stratified fluid present in the cold (or hot) leg is important in order to determine the temperature at the time-dependent temperature at the inlet of the reactor pressure vessel (RPV) and the potential to cause a thermal shock on the RPV wall. In sodium-cooled fast reactors (SFRs), core channels are typically hydro-dynamically isolated so that there exists a considerable temperature variation at the exit of adjacent fuel assemblies. All the above phenomena are characterized by the interaction of buoyant jets with the stratified flow. In stratified layers baroclinic forces create significant redistribution of turbulent kinetic energy and scales, which leads to anisotropy. This important physical phenomenon is highly three dimensional and is challenging to capture even with high-fidelity CFD simulations, due in part to lack of sufficiently resolved validation data. Furthermore, the experimental data available in the open literature do not feature the level of fidelity needed for an extensive validation of turbulence models in lower order CFD. To shed new lights into the crucial phenomena object of the present research project, it was proposed to conduct coordinated experiments and simulations at the University of Michigan and the George Washington University. The project has resulted in an experimental database of high-resolution time-resolved measurements of jets in uniform and stratified environments. The novel experimental data will be used to validate computational fluid dynamic (CFD) codes, including both Large Eddy Simulations (LES) and unsteady Reynolds-averaged Navier-Stokes equations (URANS) methodologies. In the Experimental and Multiphase flow (ECMF) laboratory at Univerisity of Michigan, we built two experimental facilities to investigate also the effect of scaling. The first facility, DESTROJER (DEnsity Stratified Turbulent ROund free Jet ExpeRiment), featuring a contoured jet nozzle with a diameter of D=12.7mm and a 1m×1m×1m cubic tank, which is made of acrylic glass for optical access. The ratio between the tank width and the nozzle diameter is equal to 78, which ensures that there is no direct interaction between the jet and the side walls. A second, modular experimental facility, features three different tank sizes of size 10×10×30, 20×20×30, 30×30×30 cm3 respectively (all tanks have the same height of 30 cm), and a jet diameter of 2mm. For the modular facility, tank-to-nozzle ratio of 50, 100, and 150 are obtained respectively. Experiments with different-density jet impactions and sharp interface with a density difference of 3.16% In the Laboratory at George Washington University, complementary experiments have been carried out at a facility featuring the same dimensions as DESTROYER. Experiments with different jet diameters have been carried out as well. High-resolution time-resolved measurements have been performed at all facilities using high-speed particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF). In this way, simultaneous measurements of velocity and density fields have been carried out. At the University of Michigan, wire-mesh sensors have been designed and built to gather high-resolution, time-dependent data of the stratified front in the tank. These measurements are important to provide accurate initial conditions for the CFD simulations. An important outcome of the project has been the development of novel techniques to achieve refractive index matching (RIM) of acqueous solutions with high density differences. RIM is crucial if optical techniques such as PIV are used to measure velocity fields. Previous to this project, experiments up to 3% density difference had been reported in the literature. With the methodology developed in this project, we have been able to push the envelop to density differences up to 9%. The experiments have been used to validate Reynolds Averaged Navier-Stokes (RANS) turbulence models and Large Eddy Simulations (LES) models. Both NEK5000 and the commercial CFD code STAR-CCM+ have been used. The experimental data will be made available to the community through the NE-KAMS databasehave been successfully conducted in the scaled facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array

Purpose: Feasibility of targeted and volumetric hyperthermia (40–45 °C) delivery to the prostate with a commercial MR-guided endorectal ultrasound phased array system, designed specifically for thermal ablation and approved for ablation trials (ExAblate 2100, Insightec Ltd.), was assessed through computer simulations and tissue-equivalent phantom experiments with the intention of fast clinical translation for targeted hyperthermia in conjunction with radiotherapy and chemotherapy. Methods: The simulations included a 3D finite element method based biothermal model, and acoustic field calculations for the ExAblate ERUS phased array (2.3 MHz, 2.3 × 4.0 cm{sup 2}, ∼1000 channels) using the rectangular radiator method. Array beamforming strategies were investigated to deliver protracted, continuous-wave hyperthermia to focal prostate cancer targets identified from representative patient cases. Constraints on power densities, sonication durations and switching speeds imposed by ExAblate hardware and software were incorporated in the models. Preliminary experiments included beamformed sonications in tissue mimicking phantoms under MR temperature monitoring at 3 T (GE Discovery MR750W). Results: Acoustic intensities considered during simulation were limited to ensure mild hyperthermia (T{sub max} < 45 °C) and fail-safe operation of the ExAblate array (spatial and time averaged acoustic intensity I{sub SATA} < 3.4 W/cm{sup 2}). Tissue volumes with therapeutic temperature levels (T > 41 °C) were estimated. Numerical simulations indicated that T > 41 °C was calculated in 13–23 cm{sup 3} volumes for sonications with planar or diverging beam patterns at 0.9–1.2 W/cm{sup 2}, in 4.5–5.8 cm{sup 3} volumes for simultaneous multipoint focus beam patterns at ∼0.7 W/cm{sup 2}, and in ∼6.0 cm{sup 3} for curvilinear (cylindrical) beam patterns at 0.75 W/cm{sup 2}. Focused heating patterns may be practical for treating focal disease in a single posterior quadrant of the prostate and diffused heating patterns may be useful for heating quadrants, hemigland volumes or even bilateral targets. Treatable volumes may be limited by pubic bone heating. Therapeutic temperatures were estimated for a range of physiological parameters, sonication duty cycles and rectal cooling. Hyperthermia specific phasing patterns were implemented on the ExAblate prostate array and continuous-wave sonications (∼0.88 W/cm{sup 2}, 15 min) were performed in tissue-mimicking material with real-time MR-based temperature imaging (PRFS imaging at 3.0 T). Shapes of heating patterns observed during experiments were consistent with simulations. Conclusions: The ExAblate 2100, designed specifically for thermal ablation, can be controlled for delivering continuous hyperthermia in prostate while working within operational constraints.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Advancing Artificial Intelligence with Liquid Argon Neutrino Experiments (Technical Report)

The grant allowed two main contributions: 1) The development of a first successful demonstration of the employment of Optimal Transport in liquid argon time projection chamber neutrino detectors. Optimal Transport, used in other contexts and specifically with LHC calorimetric data, was adapted to address a key particle identification challenge in LArTPCs: the separation of pi0 backgrounds from single-electrons produced in charged-current electron neutrino interactions. The work, leveraging ML methods such as k-nearest-neighbor (kNN) and support-vector-machine (SVM), showed an increase in background rejection of a factor of two or more. Work is now ongoing to incorporate this development in physics analyses for LArTPC experiments and more broadly expand the use of OT in LArTPC detectors including DUNE. This work was done in collaboration with the phenomenology group led by Nathaniel Craig at UCSB. 2) The deployment of NuGraph2, a graph neural network developed for LArTPC reconstruction, in the MicroBooNE experiment. NuGraph2 uses novel graph-neural-network methods on the rather simple LArTPC inputs of reconstructed hits, greatly simplifying the workflow compared to the use of waveform or signal-deconvolved wire ROIs. The network performed particle classification and was shown to address many challenging problems in LArTPC imaging including track-shower separation and the identification of protons and charged pions from primary muons. Our group collaborated with Giuseppe Cerati (FNAL scientist) who is one of the core developers of NuGraph2 to integrate this tool in MicroBooNE’s analysis framework. This consisted in tow key contributions: a) Studying performance on real data, which came with several months of iterations because the MC-trained version of the network was found to show significant bias that our group investigated and addressed. b) Integrating the output hit labeling of NuGraph2 into the existing particle tracking and shower reconstruction code. As a result of this work led by our team NuGraph2 is now enabling a suite of new analyses which benefit from enhanced capabilities and thus broader physics reach. The grant supported primarily the salary of UCSB graduate student Chuyue “Michaelia” Fang as well as partial summer salary support for PI Caratelli. Some funds were used for travel by Michaelia to ML related schools and conferences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the late-time growth of the two-dimensional Richtmyer–Meshkov instability in shock tube experiments [On the late-time growth of the 2D Richtmyer–Meshkov instability in shock tube experiments]

In the present study, shock tube experiments are used to study the very late-time development of the Richtmyer–Meshkov instability from a diffuse, nearly sinusoidal, initial perturbation into a fully turbulent flow. The interface is generated by two opposing gas flows and a perturbation is formed on the interface by transversely oscillating the shock tube to create a standing wave. The puncturing of a diaphragm generates a Mach $1. 2$ shock wave that then impacts a density gradient composed of air and SF 6 , causing the Richtmyer–Meshkov instability to develop in the 2.0 m long test section. The instability is visualized with planar Mie scattering in which smoke particles in the air are illuminated by a Nd:YLF laser sheet, and images are recorded using four high-speed video cameras operating at 6 kHz that allow the recording of the time history of the instability. In addition, particle image velocimetry (PIV) is implemented using a double-pulsed Nd:YAG laser with images recorded using a single CCD camera. Initial modal content, amplitude, and growth rates are reported from the Mie scattering experiments while vorticity and circulation measurements are made using PIV. Amplitude measurements show good early-time agreement but relatively poor late-time agreement with existing nonlinear models. The model of Goncharov agrees with growth rate measurements at intermediate times but fails at late experimental times. Measured background acceleration present in the experiment suggests that the late-time growth rate may be influenced by Rayleigh–Taylor instability induced by the interfacial acceleration. Numerical simulations conducted using the LLNL codes Ares and Miranda show that this acceleration may be caused by the growth of boundary layers, and must be accounted for to produce good agreement with models and simulations. Adding acceleration to the Richtmyer–Meshkov buoyancy–drag model produces improved agreement. It is found that the growth rate and amplitude trends are also modelled well by the Likhachev–Jacobs vortex model. Here, circulation measurements also show good agreement with the circulation value extracted by fitting the vortex model to the experimental data.

42 ENGINEERING↗

Detection and characterization of detached tidal dwarf galaxies

Tidal interactions between galaxies often give rise to tidal tails, which can harbor concentrations of stars and interstellar gas resembling dwarf galaxies. Some of these tidal dwarf galaxies (TDGs) have the potential to detach from their parent galaxies and become independent entities, but their long-term survival is uncertain. In this study, we conducted a search for detached TDGs associated with a sample of 39 interacting galaxy pairs in the local Universe using infrared, ultraviolet, and optical images. We employed IR colors and UV/optical/IR spectral energy distributions to identify potential interlopers, such as foreground stars or background quasars. Through spectroscopic observations using the Boller and Chivens spectrograph at San Pedro Mártir Observatory, we confirmed that six candidate TDGs are at the same redshift as their putative parent galaxy pairs. We identified and measured emission lines in the optical spectra and calculated nebular oxygen abundances, which range from log(O/H) = 8.10 ± 0.01 to 8.51 ± 0.02. We have serendipitously discovered an additional detached TDG candidate in Arp72 using available spectra from SDSS. Utilizing the photometric data and the CIGALE code for stellar population and dust emission fitting, we derived the stellar masses, stellar population ages, and stellar metallicities for these detached TDGs. Compared to standard mass-metallicity relations for dwarf galaxies, five of the seven candidates have higher than expected metallicities, confirming their tidal origins. One of the seven candidates remains unclear due to large uncertainties in metallicity, and another has stellar and nebular metallicities compatible with those of a preexisting dwarf galaxy. The latter object is relatively compact in the optical relative to its stellar mass, in contrast to the other candidate TDGs, which have large diameters for their stellar masses compared to most dwarf galaxies. The derived stellar population ages range from 100 Myr to 900 Myr, while the inferred stellar masses are between 2 × 10 6 M ⊙ and 8 × 10 7 M ⊙ . Four of the six TDGs are associated with the gas-rich M51-like pair Arp 72, one TDG is associated with a second M51-like pair Arp 86, and another is associated with Arp 65, an approximately equal mass pair. In spite of the relatively low stellar masses of these TDGs, they have survived for at least 100–900 Myrs, suggesting that they are stable and in dynamical equilibrium. We conclude that encounters with a relatively low-mass companion (1/10th–1/4th of the mass of the primary) can also produce long-lasting TDGs.

Astronomy & Astrophysics↗

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

Highly Extensible X-ray Diffraction Toolkit

This code has been previous released as open source under LLNL-CODE-529294. The HEXRD program comprises a library of X-ray diffraction analysis tools and a stand-alone GUI application. While there are many software packages that provide similar capabilities, HEXRD was designed to provide an extremely flexible and extensible basis for describing a broad array of X-ray diffraction instruments; this includes both mono- and poly-chromatic modalities as well as poly- and single-crystal samples. The high-level objects include Material: an abstraction for defining crystalline materials, including unit cell parameters and anisotropic elastic moduli Imageseries: an interface for series of diffraction images containing one or more frames, including provisions for metadata. Instrument: an abstraction of an X-ray diffraction instrument including support for the definition of multiple detector elements as well as powder, Laue, and mono-chromatic single/multi-grain diffraction analysis modalities (including attendant calibration routines). The hexrdgui package is a standalone PySide2/Qt5 GUI built on top of hexrd. The 0.8.1 release includes various bug fixes as well as some new functionality in the GUI, including: o Enhancements to the import tools for LLNL diffraction instruemtns (PXRDIP and TARDIS) o The ability to save 1-d LeBail fits of integrated diffraction spectra.

Barton, NathanR.↗

Calculation of energetic dual-energy detector efficiency using MCNP and GEANT Monte Carlo codes

The use of dual-energy detectors is a powerful tool in homeland security to identify materials atomic number in objects submitted to an X-ray analysis. In this work, the Monte Carlo methodology has been applied to study the response of a dual-energy lineal detector made of two scintillator materials, irradiated with a 160 KV photon spectrum emitted by an X-ray tube. To that, two different Monte Carlo codes have been used in the simulations, the MCNP5 and GEANT codes. The study of the energetic efficiency of a dual-energy detector has been developed by the simulation of several monoenergetic beams (10, 20,..., 160 KeV) oriented towards the detector, and registering the absorbed energy at both scintillators in each energy bin analyzed. Different configurations of scintillators materials and thickness have been simulated and the results show that a dual-energy detector formed GOS-CdWO{sub 4} scintillators materials is a good combination in order to separate the low energy part of the spectrum from the high energy part of it, obtaining as a result, two different images very useful for the identification of the effective atomic number of an irradiated object. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Approaches to reducing photon dose calculation errors near metal implants

Purpose: Dose calculation errors near metal implants are caused by limitations of the dose calculation algorithm in modeling tissue/metal interface effects as well as density assignment errors caused by imaging artifacts. The purpose of this study was to investigate two strategies for reducing dose calculation errors near metal implants: implementation of metal-based energy deposition kernels in the convolution/superposition (C/S) dose calculation method and use of metal artifact reduction methods for computed tomography (CT) imaging. Methods: Both error reduction strategies were investigated using a simple geometric slab phantom with a rectangular metal insert (composed of titanium or Cerrobend), as well as two anthropomorphic phantoms (one with spinal hardware and one with dental fillings), designed to mimic relevant clinical scenarios. To assess the dosimetric impact of metal kernels, the authors implemented titanium and silver kernels in a commercial collapsed cone C/S algorithm. To assess the impact of CT metal artifact reduction methods, the authors performed dose calculations using baseline imaging techniques (uncorrected 120 kVp imaging) and three commercial metal artifact reduction methods: Philips Healthcare’s O-MAR, GE Healthcare’s monochromatic gemstone spectral imaging (GSI) using dual-energy CT, and GSI with metal artifact reduction software (MARS) applied. For the simple geometric phantom, radiochromic film was used to measure dose upstream and downstream of metal inserts. For the anthropomorphic phantoms, ion chambers and radiochromic film were used to quantify the benefit of the error reduction strategies. Results: Metal kernels did not universally improve accuracy but rather resulted in better accuracy upstream of metal implants and decreased accuracy directly downstream. For the clinical cases (spinal hardware and dental fillings), metal kernels had very little impact on the dose calculation accuracy (<1.0%). Of the commercial CT artifact reduction methods investigated, the authors found that O-MAR was the most consistent method, resulting in either improved dose calculation accuracy (dental case) or little impact on calculation accuracy (spine case). GSI was unsuccessful at reducing the severe artifacts caused by dental fillings and had very little impact on calculation accuracy. GSI with MARS on the other hand gave mixed results, sometimes introducing metal distortion and increasing calculation errors (titanium rectangular implant and titanium spinal hardware) but other times very successfully reducing artifacts (Cerrobend rectangular implant and dental fillings). Conclusions: Though successful at improving dose calculation accuracy upstream of metal implants, metal kernels were not found to substantially improve accuracy for clinical cases. Of the commercial artifact reduction methods investigated, O-MAR was found to be the most consistent candidate for all-purpose CT simulation imaging. The MARS algorithm for GSI should be used with caution for titanium implants, larger implants, and implants located near heterogeneities as it can distort the size and shape of implants and increase calculation errors.

60 APPLIED LIFE SCIENCES↗