Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed space systems”

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 919 records · Page 51

Experimental evaluation of electrical conductivity imaging of anisotropic brain tissues using a combination of diffusion tensor imaging and magnetic resonance electrical impedance tomography

Anisotropy of biological tissues is a low-frequency phenomenon that is associated with the function and structure of cell membranes. Imaging of anisotropic conductivity has potential for the analysis of interactions between electromagnetic fields and biological systems, such as the prediction of current pathways in electrical stimulation therapy. To improve application to the clinical environment, precise approaches are required to understand the exact responses inside the human body subjected to the stimulated currents. In this study, we experimentally evaluate the anisotropic conductivity tensor distribution of canine brain tissues, using a recently developed diffusion tensor-magnetic resonance electrical impedance tomography method. At low frequency, electrical conductivity of the biological tissues can be expressed as a product of the mobility and concentration of ions in the extracellular space. From diffusion tensor images of the brain, we can obtain directional information on diffusive movements of water molecules, which correspond to the mobility of ions. The position dependent scale factor, which provides information on ion concentration, was successfully calculated from the magnetic flux density, to obtain the equivalent conductivity tensor. By combining the information from both techniques, we can finally reconstruct the anisotropic conductivity tensor images of brain tissues. The reconstructed conductivity images better demonstrate the enhanced signal intensity in strongly anisotropic brain regions, compared with those resulting from previous methods using a global scale factor.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

In-Memory Graph Databases for Web-Scale Data

RDF databases have emerged as one of the most relevant way for organizing, integrating, and managing expo- nentially growing, often heterogeneous, and not rigidly structured data for a variety of scientific and commercial fields. In this paper we discuss the solutions integrated in GEMS (Graph database Engine for Multithreaded Systems), a software framework for implementing RDF databases on commodity, distributed-memory high-performance clusters. Unlike the majority of current RDF databases, GEMS has been designed from the ground up to primarily employ graph-based methods. This is reflected in all the layers of its stack. The GEMS framework is composed of: a SPARQL-to-C++ compiler, a library of data structures and related methods to access and modify them, and a custom runtime providing lightweight software multithreading, network messages aggregation and a partitioned global address space. We provide an overview of the framework, detailing its component and how they have been closely designed and customized to address issues of graph methods applied to large-scale datasets on clusters. We discuss in details the principles that enable automatic translation of the queries (expressed in SPARQL, the query language of choice for RDF databases) to graph methods, and identify differences with respect to other RDF databases.

GEMS↗

Machine learning applied to classifying neutron resonances

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this project, we apply a machine learning technique, namely decision trees, to automate the quantum number assignments. The tree is trained from simulated data generated to mimic the errors found in real data. We explore the use of several physics-motivated features for training our trees, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results using random matrix theory motivated fits which demonstrated that we can determine resonance spin groups somewhat reliably. If we use these fits as features in our trees, we can train them to spot outliers corresponding to misassigned resonances. We found that with the large number of features used in this project that the decision tree tended to over t training data resulting in poor performance with respect to the test data. By reducing the number of features, we can achieve nearly perfect assignment of quantum numbers with our training data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗

THE MUSCLES TREASURY SURVEY. I. MOTIVATION AND OVERVIEW

Ground- and space-based planet searches employing radial velocity techniques and transit photometry have detected thousands of planet-hosting stars in the Milky Way. With so many planets discovered, the next step toward identifying potentially habitable planets is atmospheric characterization. While the Sun–Earth system provides a good framework for understanding the atmospheric chemistry of Earth-like planets around solar-type stars, the observational and theoretical constraints on the atmospheres of rocky planets in the habitable zones (HZs) around low-mass stars (K and M dwarfs) are relatively few. The chemistry of these atmospheres is controlled by the shape and absolute flux of the stellar spectral energy distribution (SED), however, flux distributions of relatively inactive low-mass stars are poorly understood at present. To address this issue, we have executed a panchromatic (X-ray to mid-IR) study of the SEDs of 11 nearby planet-hosting stars, the Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanetary Systems (MUSCLES) Treasury Survey. The MUSCLES program consists visible observations from Hubble and ground-based observatories. Infrared and astrophysically inaccessible wavelengths (EUV and Lyα) are reconstructed using stellar model spectra to fill in gaps in the observational data. In this overview and the companion papers describing the MUSCLES survey, we show that energetic radiation (X-ray and ultraviolet) is present from magnetically active stellar atmospheres at all times for stars as late as M6. The emission line luminosities of C iv and Mg ii are strongly correlated with band-integrated luminosities and we present empirical relations that can be used to estimate broadband FUV and XUV (≡X-ray + EUV) fluxes from individual stellar emission line measurements. We find that while the slope of the SED, FUV/NUV, increases by approximately two orders of magnitude form early K to late M dwarfs (≈0.01–1), the absolute FUV and XUV flux levels at their corresponding HZ distances are constant to within factors of a few, spanning the range 10–70 erg cm{sup −2} s{sup −1} in the HZ. Despite the lack of strong stellar activity indicators in their optical spectra, several of the M dwarfs in our sample show spectacular UV flare emission in their light curves. We present an example with flare/quiescent ultraviolet flux ratios of the order of 100:1 where the transition region energy output during the flare is comparable to the total quiescent luminosity of the star E{sub flare}(UV) ∼ 0.3 L{sub *}Δt (Δt = 1 s). Finally, we interpret enhanced L(line)/L{sub Bol} ratios for C iv and N v as tentative observational evidence for the interaction of planets with large planetary mass-to-orbital distance ratios (M{sub plan}/a{sub plan}) with the transition regions of their host stars.

79 ASTRONOMY AND ASTROPHYSICS↗

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗

TARGET SELECTION FOR THE LBTI EXOZODI KEY SCIENCE PROGRAM

The Hunt for Observable Signatures of Terrestrial planetary Systems (HOSTS) on the Large Binocular Telescope Interferometer will survey nearby stars for faint emission arising from ∼300 K dust (exozodiacal dust), and aims to determine the exozodiacal dust luminosity function. HOSTS results will enable planning for future space telescopes aimed at direct spectroscopy of habitable zone terrestrial planets, as well as greater understanding of the evolution of exozodiacal disks and planetary systems. We lay out here the considerations that lead to the final HOSTS target list. Our target selection strategy maximizes the ability of the survey to constrain the exozodi luminosity function by selecting a combination of stars selected for suitability as targets of future missions and as sensitive exozodi probes. With a survey of approximately 50 stars, we show that HOSTS can enable an understanding of the statistical distribution of warm dust around various types of stars and is robust to the effects of varying levels of survey sensitivity induced by weather conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron-Phonon Systems on a Universal Quantum Computer

We present an algorithm that extends existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed matter physics to include bosons in general and phonons in particular. We introduce a qubit representation for the low-energy subspace of phonons which allows an efficient simulation of the evolution operator of the electron-phonon systems. As a consequence of the Nyquist-Shannon sampling theorem, the phonons are represented with exponential accuracy on a discretized Hilbert space with a size that increases linearly with the cutoff of the maximum phonon number. The additional number of qubits required by the presence of phonons scales linearly with the size of the system. The additional circuit depth is constant for systems with finite-range electron-phonon and phonon-phonon interactions and linear for long-range electron-phonon interactions. Our algorithm for a Holstein polaron problem was implemented on an Atos quantum learning machine quantum simulator employing the quantum phase estimation method. The energy and the phonon number distribution of the polaron state agree with exact diagonalization results for weak, intermediate, and strong electron-phonon coupling regimes.

97 MATHEMATICS AND COMPUTING↗

Accelerating the Global Arrays ComEx Runtime using Multiple Progress Ranks

Partitioned Global Address Space (PGAS) models are part of system software that are being designed to support communication runtime for exascale applications. Suitability of using MPI as a scalable PGAS communication subsystem has been shown as a viable option due to its standardization and higher performance. We used MPI two-sided semantics with a combination of automatic and user defined splitting of MPI communicators to achieve asynchronous progress. Our implementation can make use of multiple asynchronous progress ranks (PR) per node that can be mapped to the computing architecture of a node in a distributed cluster. We are able to show significant speed up of over 2.0X and scaling of communication bound computational chemistry application over 1024 nodes of state- of-the-art HPC clusters. Our results show that while running a communication bound application workload on a certain number of cluster nodes, a use of optimum number of ranks dedicated for communication may be used to achieve asynchronous communication progress and obtain highest performance.

PGAS, global-arrays, NWCHEM, NWChemEx, HPC, GPGPU↗

Flexibility from Energy Systems Integration: Supporting Synergies Among Sectors

Energy systems integration, or sector coupling, has several drivers that span climate impact mitigation and economics to social and regulatory considerations. A key question is what is sector coupling, and how does it impact the flexibility of the energy system? Here, the energy system includes several sectors - electricity, gas, heat, and transportation - that have been independent for decades in most countries except for their coupling via combined heat and power (CHP) units. In energy systems integration, some sectors may provide flexibility to other sectors, while other sectors will require flexibility when interlinking. To support these synergies among sectors, it is important to explore and quantify mutual interactions as well as seek examples of how these integrations can provide flexibility and other benefits. From the perspective of the electricity sector, it is important to ensure that there is enough flexibility in the interconnected systems to support decarbonization goals, such as those set in the Paris Agreement, while ensuring operational reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗

A tri-level optimization model for interdependent infrastructure network resilience against compound hazard events

Resilient operation of interdependent infrastructures against compound hazard events is essential for maintaining societal well-being. To address consequence assessment challenges in this problem space, we propose a novel policy-guided tri-level optimization model applied to a proof-of-concept case study with fuel distribution and transportation networks – encompassing one realistic network; one fictitious, yet realistic network; as well as networks drawn from three synthetic distributions. Mathematically, our approach takes the form of a defender-attacker-defender (DAD) model—a multi-agent tri-level optimization, comprised of a defender, attacker, and an operator acting in sequence. Here, in this study, our notional operator may choose proxy actions to operate an interdependent system comprised of fuel terminals and gas stations (functioning as supplies) and a transportation network with traffic flow (functioning as demand) to minimize unmet demand at gas stations. A notional attacker aims to hypothetically disrupt normal operations by reducing supply at the supply terminals, and the notional defender aims to identify best proxy defense policy options which include hardening supply terminals or allowing alternative distribution methods such as trucking reserve supplies. We solve our DAD formulation at a metropolitan scale and present practical defense policy insights against hypothetical compound hazards. We demonstrate the generalizability of our framework by presenting results for a realistic network; a fictitious, yet realistic network; as well as for three networks drawn from synthetic distributions. Additionally, we demonstrate the scalability of the framework by investigating runtime performance as a function of the network size. Steps for future research are also discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r,v r ,v θ ) phase space coordinates

In this paper, a fourth-order finite-volume Vlasov–Poisson algorithm is developed for simulating axisymmetric plasma configurations in ( r , v r , v θ ) phase space coordinates. The Vlasov equation for cylindrical phase space coordinates is cast into conservation-law form and is discretized on a structured grid. The conservative finite-volume discretization is based on fifth-order upwind reconstructions of the distribution function and a fourth-order quadrature rule that accounts for transverse variations of fluxes along control volume surfaces. High-order specular reflection boundary conditions enable high-fidelity treatment of plasma distribution functions at the axis and at wall boundaries. The numerical method is applied to simulate a confined uniform neutral gas to assess convergence properties for an equilibrium system in which effects of finite-temperature and acceleration due to centrifugal and Coriolis forces are present. The discretization is also applied to a Z-pinch configuration to study electrostatic ion confinement and the dynamics of the associated breathing mode. Simulations show that the ion distribution function exhibits non-Maxwellian features and that a temperature anisotropy develops and is sustained. Finally, the finite-volume implementation is demonstrated to converge at fourth-order for both applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Supramolecular Assembly of Multielement Ribbon-like Structures Derived from Halide Perovskites

Halide perovskites are crucial materials with broad applications owing to their exceptional optoelectronic properties. Vacancy-ordered double perovskites, featuring highly tunable transition metal sites, enable controllable optoelectronic properties through multielement compositional design. Here, in this study, we introduced 18-crown-6 into the vacancy-ordered double perovskites system and developed two-dimensional ribbon-like single crystals (18C6@K) 2 {PtSnTeIrRe} 1 Cl 6 via an antisolvent supramolecular assembly method, demonstrating morphology modulation through multielement composition design. The crystals crystallize in the centrosymmetric space group $P\bar{1}$. The dumbbell-shaped structural units (crown ether@A)2MX6 pack along a and b axes to form a 2-dimensional (2D) monolayer, and these monolayers further stack along the c axis to generate the ribbon-like single crystals. Energy-dispersive X-ray spectroscopy (EDX) qualitatively confirmed the uniform distribution of the five transition metals throughout the crystal, while inductively coupled plasma atomic emission spectroscopy (ICP–AES) quantitatively verified their atomic ratios. We further investigated the origin of the morphology, distinct from the previously reported cube-like single crystals with the $R\bar{3}$ space group. When acetonitrile was used as the solvent, three-dimensional crystals with $R\bar{3}$ symmetry were obtained, whereas dimethylformamide (DMF) was essential for forming two-dimensional ribbon-like single crystals. The essential role of DMF could be ascribed to its capability to maintain a higher concentration of the building blocks. Moreover, Ir 4+ and Pt 4+ cations also played critical roles in inducing the two-dimensional ribbon-like morphology. The three-element (18C6@K) 2 {PtSnTe}1Cl 6 single crystals exhibited bright yellow emission under 375 nm laser excitation, demonstrating the tunability of optoelectronic properties of this class of material.

Zhu, Heqing [University of California, Berkeley, C↗

The effects of phase decorrelation on the dynamics of the bump-on-tail instability

The dynamics of the bump-on-tail instability has been studied. The novelty of the work is the analysis of how the bump-on-tail dynamics is affected by an extrinsic stochastisation of the phase of the wave-particle interaction; here referred to as phase decorrelation. For this purpose, a nonlinear Monte Carlo model has been developed. When the characteristic time scale for macroscopic phase decorrelation becomes shorter than time scales of nonlinear wave-particle dynamics, the system may be described quasilinearly, with the phase decorrelation being replaced by a quasilinear diffusion coefficient in particle energy. A purely quasilinear Monte Carlo model, which is typically less computationally demanding than the fully nonlinear description due to the reduced dimensionality of phase space, has been developed for comparison. In this paper, parameter regimes, where the nonlinear and the quasilinear descriptions quantitatively agree on a macroscopic level, have been investigated, using combined theoretical and numerical analyses. Qualitative effects on the macroscopic dynamics by the presence of phase decorrelation and/or by structures of the energy distribution function in the proximity of the wave-particle resonance are also studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗