SEARCH · Search NASA

Results for “3D 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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2021

On the use of neural networks for energy reconstruction in high-granularity calorimeters

We contrasted the performance of deep neural networks — Convolutional Neural Network (CNN) and Graph Neural Network (GNN) — to current state of the art energy regression methods in a finely 3D-segmented calorimeter simulated by GEANT4. This comparative benchmark gives us some insight to assess the particular latent signals neural network methods exploit to achieve superior resolution. A CNN trained solely on a pure sample of pions achieved substantial improvement in the energy resolution for both single pions and jets over the conventional approaches. It maintained good performance for electron and photon reconstruction. We also used the Graph Neural Network (GNN) with edge convolution to assess the importance of timing information in the shower development for improved energy reconstruction. We implement a simple simulation based correction to the energy sum derived from the fraction of energy deposited in the electromagnetic shower component. This serves as an approximate dual-readout analogue for our benchmark comparison. Although this study does not include the simulation of detector effects, such as electronic noise, the margin of improvement seems robust enough to suggest these benefits will endure in real-world application. We also find reason to infer that the CNN/GNN methods leverage latent features that concur with our current understanding of the physics of calorimeter measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Multi-dimensional incoherent Thomson scattering system in PHAse Space MApping (PHASMA) facility

A multi-dimensional incoherent Thomson scattering diagnostic system capable of measuring electron temperature anisotropies at the level of the electron velocity distribution function (EVDF) is implemented on the PHAse Space MApping facility to investigate electron energization mechanisms during magnetic reconnection. This system incorporates two injection paths (perpendicular and parallel to the axial magnetic field) and two collection paths, providing four independent EVDF measurements along four velocity space directions. For strongly magnetized electrons, a 3D EVDF comprised of two characteristic electron temperatures perpendicular and parallel to the local magnetic field line is reconstructed from the four measured EVDFs. As a result, validation of isotropic electrons in a single magnetic flux rope and a steady-state helicon plasma is presented.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Three dimensions, two microscopes, one code: Automatic differentiation for x-ray nanotomography beyond the depth of focus limit

Conventional tomographic reconstruction algorithms assume that one has obtained pure projection images, involving no within-specimen diffraction effects nor multiple scattering. Advances in x-ray nanotomography are leading toward the violation of these assumptions, by combining the high penetration power of x-rays, which enables thick specimens to be imaged, with improved spatial resolution that decreases the depth of focus of the imaging system. We describe a reconstruction method where multiple scattering and diffraction effects in thick samples are modeled by multislice propagation and the 3D object function is retrieved through iterative optimization. We show that the same proposed method works for both full-field microscopy and for coherent scanning techniques like ptychography. Our implementation uses the optimization toolbox and the automatic differentiation capability of the open-source deep learning package TensorFlow, demonstrating a straightforward way to solve optimization problems in computational imaging with flexibility and portability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2024

Simulated Feasibility of 3D Lightning Mapping from Space

In addition to the awe it inspires, lightning can illuminate the microphysical processes hidden away within deep convection. The current generation of space-based lightning mapping uses mostly 2D optical imaging to connect overall flash characteristics to their parent storm dynamics, but are missing a dimension’s worth of information. With lightning now classified as an essential climate variable, future spaceborne mappers will need improved capabilities to take advantage of the 3D structure of lightning flashes to support meteorological and climate modeling. We report here on a study of the feasibility of implementing a radio frequency (RF)–based network of satellites for 3D lightning mapping in high resolution from low Earth orbit (LEO). Lightning sources are simulated using existing tools of the Lightning Mapping Array (LMA), modified for orbital detection, and spatially reconstructed using a Levenberg-Marquardt geolocation algorithm to assess sources of uncertainty in these solutions. We analyze the benefits and limitations of this approach compared to existing orbital and ground-based methods. Results of this study show that lightning can be mapped in 3D with a vertical location accuracy better than 2 km using as few as five satellites in LEO capable of measuring the time-of-arrival of impulsive RF signals in the very high frequency (VHF) band. The consequence of this study is that high-resolution, spaceborne 3D mapping of lightning is achievable across most of the globe, having crucial implications for our understanding of not only lightning, but also severe weather development, climate science, and more.

CubeSpark↗

DOE OSTI2024

Subject-specific multi-scale modeling of the fate of inhaled aerosols

Determining the fate of inhaled aerosols in the respiratory system is essential in assessing the potential toxicity of inhaled airborne materials, responses to airborne pathogens, or in improving inhaled drug delivery. The availability of high-resolution clinical lung imaging and advances in the reconstruction of lung airways from CT images have led to the development of subject-specific in-silico 3D models of aerosol dosimetry, often referred to as computational fluid-particle-dynamics (CFPD) models. As CFPD models require extensive computing resources, they are typically confined to the upper and large airways. These models can be combined with lower-dimensional models to form multiscale models that predict the transport and deposition of inhaled aerosols in the entire respiratory tract. Understanding where aerosols deposit is only the first of potentially several key events necessary to predict an outcome, being a detrimental health effect or a therapeutic response. To that end, multiscale approaches that combine CFPD with physiologically-based pharmacokinetics (PBPK) models have been developed to evaluate the absorption, distribution, metabolism, and excretion (ADME) of toxic or medicinal chemicals in one or more compartments of the human body. CFPD models can also be combined with host cell dynamics (HCD) models to assess regional immune system responses. Here, this paper reviews the state of the art of these different multiscale approaches and discusses the potential role of personalized or subject-specific modeling in respiratory health.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

DOE OSTI2020

Three-dimensional super line-localization in low signal-to-noise microscope images via prior-apprised unsupervised learning (PAUL)

Biological processes such as processive enzyme turnover and intracellular cargo tracking involve the dynamic motion of a small "article" along a curvilinear biopolymer track. To understand these processes that occur across multiple length and time scales, one must acquire both the trajectory of the particle and the position of the track along which it moves, possibly by combining high-resolution single-particle tracking with conventional microscopy. Yet, usually there is a significant resolution mismatch between these modalities: while the tracked particle is localized with a precision of 10 nm, the image of the surroundings is limited by optical difraction, with 200 nm lateral and 500 nm axial resolutions. Compared to the particle's trajectory, the surrounding curvilinear structure appears as a blurred and noisy image. This disparity in the spatial resolutions of the particle trajectory and the surrounding curvilinear structure image makes data reconstruction, as well as interpretation, particularly challenging. Analysis is further complicated when the curvilinear structures are oriented arbitrarily in 3D space. Here, we present a prior-apprised unsupervised learning (PAUL) approach to extract information from 3D images where the underlying features resemble a curved line such as a filament or microtubule. This three-stage framework starts with a Hessian-based feature enhancement, which is followed by feature registration, where local line segments are detected on repetitively sampled subimage tiles. In the final stage, statistical learning, segments are clustered based on their geometric relationships. Principal curves are then approximated from each segment group via statistical tools including principal component analysis, bootstrap and kernel transformation. This procedure is characterized on simulated images, where sub-voxel medium deviations from true curves have been achieved. The 3D PAUL approach has also been implemented for successful line localization in experimental 3D images of gold nanowires obtained using a multifocal microscope. Lastly, this work not only bridges the resolution gap between two microscopy modalities, but also allows us to conduct 3D super line-localization imaging experiments, without using super-resolution techniques.

42 ENGINEERING↗

DOE OSTI2023

Fast 2D Bicephalous Convolutional Autoencoder for Compressing 3D Time Projection Chamber Data

High-energy large-scale particle colliders produce data at high speed in the order of 1 terabytes per second in nuclear physics and petabytes per second in high energy physics. Developing real-time data compression algorithms to reduce such data at high throughput to fit permanent storage has drawn increasing attention. Specifically, at the newly constructed sPHENIX experiment at the Relativistic Heavy Ion Collider (RHIC), a time projection chamber is used as the main tracking detector, which records particle trajectories in a volume of three-dimensional (3D) cylinder. The resulting data are usually very sparse with occupancy around 10.8%. Such sparsity presents a challenge to conventional learning-free lossy compression algorithms, such as SZ, ZFP, and MGARD. The 3D convolutional neural network (CNN)-based approach, Bicephalous Convolutional Autoencoder (BCAE), outperforms traditional methods both in compression rate and reconstruction accuracy. BCAE can also utilize the computation power of graphical processing units suitable for deployment in a modern heterogeneous highperformance computing environment. This work introduces two BCAE variants: BCAE++ and BCAE-2D. BCAE++ achieves a 15% better compression ratio and a 77% better reconstruction accuracy measured in mean absolute error compared with BCAE. BCAE-2D treats the radial direction as the channel dimension of an image, resulting in a 3× speedup in compression throughput. In addition, we demonstrate an unbalanced autoencoder with a larger decoder can improve reconstruction accuracy without significantly sacrificing throughput. Lastly, we observe both the BCAE++ and BCAE-2D can benefit more from using half-precision mode in throughput (76 - 79% increase) without loss in reconstruction accuracy. The source code and links to data and pretrained models can be found at https://github.com/BNL-DAQ-LDRD/NeuralCompression_v2

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Analysis of Defects in Metal Additive Manufacturing with Augmented Data Generation

Laser powder bed fusion (LPBF) is a method of additive manufacturing (AM) that selectively melts and fuses together microscopic metallic powder. LPBF offers the benefit of producing custom structures out of high strength metals that can be difficult to fabricate with conventional methods. The challenge of LPBF is that 3D printed structures often have internal pores due to process flaws. Pulsed thermal tomography (PTT) is a method for reconstructing the depth profile of materials, allowing the visualization internal voids in solids. In prior work, we developed a convolutional neural network (CNN) which, having been trained on simulated 2D PTT images of subsurface elliptical defects, was able to classify the semi-major radii, semi-minor radii, and angular orientation of the best-fit ellipses in previously unseen PTT images. The unseen PTT images contained subsurface irregular defect shapes imported from scanning electron microscopy (SEM) images of metallic LPBF-printed specimens. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Self-consistent models for triaxial galaxies with flat rotation curves - The disk case

We examine the possibility of constructing scale-free triaxial logarithmic potentials self-consistently, using Schwarzschild's linear programing method. In particular, we explore the limit of nonaxisymmetric disks. In this case it is possible to reduce the problem to the self-consistent reconstruction of the disk surface density on the unit circle, a considerably simpler problem than the usual 2D or 3D one. Models with surface densities of the form Sigma = (x exp n + (y/q) exp n) exp - 1/n with n = 2 or 4 are investigated. We show that the complicated shapes of the 'boxlet' orbit families (which replace the box orbit family found in potentials with smooth cores) limit the possibility of building self-consistent models, though elliptical disks of axis ratio above 0.7 and a restricted range of boxier models can be constructed. This result relies on using sufficiently fine bins, smaller than the 10 deg bins commonly used in 2D or 3D investigations. It also indicates the need for caution in interpreting N-body models of triaxial halos in which the core of the potential is numerically smoothed.

Kuijken, Konrad↗

NASA NTRS2023

Solar Activity Relations in Energetic Electron Events Measured By the MESSENGER Mission

Aims. We perform a statistical study of the relations between the properties of solar energetic electron (SEE) events measured by the MESSENGER mission from 2010 to 2015 and the parameters of the respective parent solar activity phenomena in order to identify the potential correlations between them. During the time of analysis, the MESSENGER heliocentric distance varied between 0.31 and 0.47 au. Methods. We used a published list of 61 SEE events measured by MESSENGER, which includes information on the near-relativistic electron peak intensities, the peak-intensity energy spectral indices, and the measured X-ray peak intensity of the flares related to the SEE events. Taking advantage of multi-viewpoint remote-sensing observations, we reconstructed, whenever possible, the associated coronal mass ejections (CMEs) and shock waves; and we determined the three-dimensional (3D) properties (location, speed, and width) of the CMEs and the maximum speed of the 3D CME-driven shocks in the corona. We used different methods (Spearman, Pearson, and a Bayesian approach, namely the Kelly method to linear regression) to estimate the correlation coefficients between the flare intensity, maximum speed at the apex of the CME-driven shock, CME speed at the apex, and CME width with the electron peak intensities and with the energy spectral indices. In this statistical study, we considered and addressed the limitations of the particle instrument on board MESSENGER (elevated background intensity level, anti-Sun pointing). Results. There is an asymmetry to the east in the range of connection angles (CAs) for which the SEE events present the highest peak intensities, where the CA is the longitudinal separation between the footpoint of the magnetic field connecting to the spacecraft and the flare location. Based on this asymmetry, we define a subsample of well-connected events as when −65° ≤ CA ≤ +33°. For the well-connected sample, we find moderate to strong correlations between the near-relativistic electron peak intensity and the 3D CME-driven shock maximum speed at the apex (Spearman: cc = 0.53 ± 0.05; Pearson: cc = 0.65 ± 0.04; Kelly: cc = 0.87 ± 0.20), the flare peak intensity (Spearman: cc = 0.63 ± 0.03; Pearson: cc = 0.59 ± 0.03; Kelly: cc = 0.74 ± 0.30), and the 3D CME speed at the apex (Spearman: cc = 0.50 ± 0.04; Pearson: cc = 0.46 ± 0.03; Kelly: cc = 0.60 ± 0.39). When including poorly connected events (full sample), the relations between the peak intensities and the solar-activity phenomena are blurred, showing lower correlation coefficients. Conclusions. Based on the comparison of the correlation coefficients presented in this study using near 0.4 au data, (1) both flare and shock-related processes may contribute to the acceleration of near relativistic electrons in large SEE events, in agreement with previous studies based on near 1 au data; and (2) the maximum speed of the CME-driven shock is a better parameter to investigate particle-acceleration-related mechanisms than the average CME speed, as suggested by the stronger correlation with the SEE peak intensities.

Sun↗

DOE OSTI2026

Improving ICARUS track reconstruction algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino program at Fermilab. Its primary objective is to explore the possible existence of sterile neutrinos in the O(1 eV) mass range and to clarify the anomalies observed in the Liquid Scintillator Neutrino Detector and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber, capable of producing high-resolution 3D images and precise calorimetric measurements of ionizing particles. This technology allows for a detailed study of neutrino interactions across a broad energy range, from a few keV to several hundred GeV. The track reconstruction is achieved through a software framework that applies a series of pattern recognition algorithms, transforming raw detector signals into fully reconstructed event topologies. This process involves identifying interaction vertices, particle tracks, and electromagnetic showers within the TPC. However, in certain cases, these algorithms may mistakenly break a single particle track into several shorter segments, interpreting each as a distinct particle. Since track length is used to estimate the particle's energy, such fragmentation can result in an energy underestimation of several hundred MeV. Furthermore, when a track is split into multiple segments, the particle identification (which relies on analyzing the energy loss as a function of the residual range) may fail, potentially leading to the loss of the entire event. To mitigate this problem, we have developed a dedicated algorithm designed to identify and reconnect (“stitch”) the tracks that were erroneously divided into multiple segments.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗

DOE OSTI2020

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Assessing the feasibility of a spaceborne 3D lightning observing concept

The distribution of electrical charge in thunderclouds results from thermodynamic, microphysical, and kinematic processes, which also modulate thunderstorm evolution. It is no surprise that the connection between lightning and these physical processes is so strong that the increase and vertical growth of lightning activity closely follows the vertical growth of the thundercloud, but unraveling these connections is not trivial and requires observations of the three-dimensional (3D) structure of electrical activity in a cloud. Ground-based 3D lightning mapping networks give excellent 3D flash-level detail but are limited to regional coverage. Satellite-based optical lightning mappers give excellent global coverage but are largely limited to 2D summaries of flash rate and radiant intensity, albeit new flash products and stereographic techniques are chipping away this limitation. New observing strategies are needed to expand and diversify the corpus of 3D lightning datasets and motivate studies that unravel connections lightning has with these key physical processes and the surrounding environment. This study examines the feasibility of using a distributed network of orbing satellites with VHF-based lightning detectors to obtain global maps of 3D lightning activity and assess efficacy of this approach for use in a new, small satellite mission concept called CubeSpark. CubeSpark combines new VHF and high-resolution, bispectral optical instruments on a constellation of low-Earth orbiting (LEO) satellites to globally map the 3D electrical structure of thunderstorms and study how it relates to thunderstorm evolution, extreme weather, nitrogen oxide production and distribution, upper atmospheric electrical phenomena, and how 3D flash observations can complement existing satellite-based lightning mappers and improve decision support tools. To locate lightning discharges, CubeSpark seeks to use the VHF time-of-arrival technique, similar to ground-based total lightning mapping networks. The vertical location accuracy of these satellite retrievals will be of poorer quality compared to a ground-based network, which has non-trivial implications for lightning flash reconstruction and lightning-based interpretations of deep convection. We adapt a Lightning Mapping Array (LMA) simulation framework to an orbiting network and use it to address feasibility of 3D lightning detection from space with particular attention to location accuracy of VHF detections of lightning in the vertical. These simulations inform a constellation design study that defines a realistic orbital configuration and depicts the global coverage for CubeSpark. Results indicate that a 3D location accuracy of <1-2 km for each dimension can be achieved across 300-500 km wide swaths, which suggests that CubeSpark can resolve the charge structure of thunderclouds from the tropics to the mid- and high- latitudes.

Lightning↗

DOE OSTI2023

Hybrid achromatic microlenses with high numerical apertures and focusing efficiencies across the visible

Abstract Compact visible wavelength achromats are essential for miniaturized and lightweight optics. However, fabrication of such achromats has proved to be exceptionally challenging. Here, using subsurface 3D printing inside mesoporous hosts we densely integrate aligned refractive and diffractive elements, forming thin high performance hybrid achromatic imaging micro-optics. Focusing efficiencies of 51–70% are achieved for 15μm thick, 90μm diameter, 0.3 numerical aperture microlenses. Chromatic focal length errors of less than 3% allow these microlenses to form high-quality images under broadband illumination (400–700 nm). Numerical apertures upwards of 0.47 are also achieved at the cost of some focusing efficiency, demonstrating the flexibility of this approach. Furthermore, larger area images are reconstructed from an array of hybrid achromatic microlenses, laying the groundwork for achromatic light-field imagers and displays. The presented approach precisely combines optical components within 3D space to achieve thin lens systems with high focusing efficiencies, high numerical apertures, and low chromatic focusing errors, providing a pathway towards achromatic micro-optical systems.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

One-Step Synthesis and Operando Electrochemical Impedance Spectroscopic Characterization of Heterostructured MoP–Mo 2 N Electrocatalysts for Stable Hydrogen Evolution Reaction

Here this study presents a novel synthesis of self-standing MoP and Mo 2 N heterostructured electrocatalysts with enhanced stability and catalytic performance. Facilitated by the controlled phase and interfacial microstructure, the seamless structures of these catalysts minimize internal resistivity and prevent local corrosion, contributing to increased stability. The chemical synthesis proceeds with etching step to activate the surface, followed by phosphor-nitriding in a chemical vapor deposition chamber to produce MoP-Mo 2 N@Mo heterostructured electrocatalysts. X-ray diffraction analyses confirm the presence of MoP, Mo 2 N, and Mo phases in the electrocatalyst. Morphology studies using scanning electron microscopy characterizes the hierarchical growth of structures, indicating successful formation of the heterostructure. X-ray photoelectron spectroscopy (XPS) analyses of the as-synthesized and post-catalytic activity samples reveal the chemical shift in terms of binding energy (BE) of Mo 3d XPS peak, especially after catalytic activity. The XPS BE shifts attributed to changes in oxidation state, electron transfer, and surface reconstruction during catalysis. Electrochemical evaluation of the catalysts demonstrates the superior performance of the MoP-Mo 2 N@Mo heterostructured catalyst in hydrogen evolution reactions (HER), with lower overpotentials and enhanced Tafel slopes. Stability tests reveal changes in double layer capacitance over time, suggesting surface reconstruction and increased active surface area during catalysis. Operando electrochemical impedance spectroscopy (EIS) further elucidates the dynamic changes in resistance and charge transfer during HER. Overall, comprehensive understanding of the synthesis, characterization, and electrochemical behavior of the developed MoP-Mo 2 N@Mo heterostructured electrocatalyst as presented in this paper highlights their potential utilization in sustainable energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

First look at neutron emission shape characteristics of ignition hotspots at the National Ignition Facility (invited)

The nuclear imaging system has been capturing neutron images of inertial confinement fusion (ICF) driven implosions for over a decade at the National Ignition Facility. This imaging system has evolved from one to three nearly orthogonal lines-of-sight, allowing for the study of three-dimensional shape characteristics of ignition shots. Limited-view tomography algorithms help visualize the burning hotspot in 3D and assess neutron source geometry using Legendre mode parameters. With its neutron, gamma-ray, and x-ray image reconstruction capabilities, NIS has provided critical insight into mechanisms that have limited implosion performance, such as fill tube diameter for ignition-type targets. This comprehensive diagnostic suite opens a window into the shape characteristics of ignition shots and how symmetry affects ICF implosion performance. In more recent ignition shots, neutron yields have visibly increased. Analyzing the shape and size of the reconstructed neutron source has shown an expansion of the burn volume, which is indicative of more efficient alpha heating during the implosion process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗