Search NASA⌕ Search

SEARCH · Search NASA

Results for “data ordering”

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 937 records · Page 52

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer’s cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗

Pairing symmetry of twisted bilayer graphene: A phenomenological synthesis

One of the outstanding questions in the study of twisted bilayer graphene—from both experimental and theoretical points of view—is the nature of its superconducting phase. In this paper we perform a comprehensive synthesis of existing experiments, and argue that experimental constraints are strong enough to allow the structure of the superconducting order parameter to be nearly uniquely determined. In particular, we argue that the order parameter is nodal, and is formed from an admixture of spin-singlet and spin-triplet Cooper pairs. This argument is made on phenomenological grounds, without committing to any particular microscopic model of the superconductor. Existing data is insufficient to determine the orbital parity of the order parameter, which could be either p wave or d wave. We propose a way in which the measurement of Andreev edge states can be used to distinguish between the two.

36 MATERIALS SCIENCE↗

PLANC: Parallel Low-rank Approximation with Nonnegativity Constraints

In this work, we consider the problem of low-rank approximation of massive dense nonnegative tensor data, for example, to discover latent patterns in video and imaging applications. As the size of data sets grows, single workstations are hitting bottlenecks in both computation time and available memory. We propose a distributed-memory parallel computing solution to handle massive data sets, loading the input data across the memories of multiple nodes, and performing efficient and scalable parallel algorithms to compute the low-rank approximation. We present a software package called Parallel Low-rank Approximation with Nonnegativity Constraints, which implements our solution and allows for extension in terms of data (dense or sparse, matrices or tensors of any order), algorithm (e.g., from multiplicative updating techniques to alternating direction method of multipliers), and architecture (we exploit GPUs to accelerate the computation in this work). We describe our parallel distributions and algorithms, which are careful to avoid unnecessary communication and computation, show how to extend the software to include new algorithms and/or constraints, and report efficiency and scalability results for both synthetic and real-world data sets.

97 MATHEMATICS AND COMPUTING↗

Order-Disorder Transitions in ( Ca x Sr 1 - x ) 3 Rh 4 Sn 13

The classification of structural phase transitions as displacive or order-disorder in character is usually based on spectroscopic data above the transition. Here, we use single crystal x-ray diffraction to investigate structural correlations in the quasiskutterudites, (Ca x Sr 1-x ) 3 Rh 4 Sn 13 , which have a quantum phase transition at x ~ 0.9. Three-dimensional pair distribution functions show that the amplitudes of local atomic displacements are temperature independent below the transition and persist to well above the transition, a signature of order-disorder behavior. The implications for the associated electronic transitions are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Field-induced intermediate ordered phase and anisotropic interlayer interactions in α-RuCl 3

In α-RuCl 3 , an external magnetic field applied within the honeycomb plane can induce a transition from a magnetically ordered state to a disordered state that is potentially related to the Kitaev quantum spin liquid. In zero field, single crystals with minimal stacking faults display a low-temperature state with in-plane zigzag antiferromagnetic order and a three-layer periodicity in the direction perpendicular to the honeycomb planes. In this work, we present angle-dependent magnetization, ac susceptibility, and thermal transport data that demonstrate the presence of an additional intermediate-field ordered state at fields below the transition to the disordered phase. Neutron-diffraction results show that the magnetic structure in this phase is characterized by a six-layer periodicity in the direction perpendicular to the honeycomb planes. Theoretically, the intermediate ordered phase can be accounted for by including spin-anisotropic couplings between the layers in a three-dimensional spin model. Together, this demonstrates the importance of interlayer exchange interactions in α-RuCl 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Using Bayesian Methodology to Estimate Liquefied Natural Gas Leak Frequencies

This analysis provides estimates on the leak frequencies of nine components found in liquefied natural gas (LNG) facilities. Data was taken from a variety of sources, with 25 different data sets included in the analysis. A hierarchical Bayesian model was used that assumes that the log leak frequency follows a normal distribution and the logarithm of the mean of this normal distribution is a linear function of the logarithm of the fractional leak area. This type of model uses uninformed prior distributions that are updated with applicable data. Separate models are fit for each component listed. Five order-of-magnitude fractional leak areas are considered, based on the flow area of the component. Three types of supporting analyses were performed: sensitivity of the model to the data set used, sensitivity of the leak frequency estimates to differences in the model structure or prior distributions, and sufficiency of sample sized used for convergence. Recommended leak frequency distributions for all component types and leak sizes are given. These leak frequency predictions can be used for quantitative risk assessments in the future.

03 NATURAL GAS↗

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗

Filament-Free Bulk Resistive Memory Enables Deterministic Analogue Switching

Digital computing is nearing its physical limits as computing needs and energy consumption rapidly increase. Analogue-memory-based neuromorphic computing can be orders of magnitude more energy efficient at data-intensive tasks like deep neural networks, but has been limited by the inaccurate and unpredictable switching of analogue resistive memory. Filamentary resistive random access memory (RRAM) suffers from stochastic switching due to the random kinetic motion of discrete defects in the nanometer-sized filament. Here, this stochasticity is overcome by incorporating a solid electrolyte interlayer, in this case, yttria-stabilized zirconia (YSZ), toward eliminating filaments. Filament-free, bulk-RRAM cells instead store analogue states using the bulk point defect concentration, yielding predictable switching because the statistical ensemble behavior of oxygen vacancy defects is deterministic even when individual defects are stochastic. Both experiments and modeling show bulk-RRAM devices using TiO2-X switching layers and YSZ electrolytes yield deterministic and linear analogue switching for efficient inference and training. Bulk-RRAM solves many outstanding issues with memristor unpredictability that have inhibited commercialization, and can, therefore, enable unprecedented new applications for energy-efficient neuromorphic computing. Beyond RRAM, this work shows how harnessing bulk point defects in ionic materials can be used to engineer deterministic nanoelectronic materials and devices.

36 MATERIALS SCIENCE↗

Graph identification of proteins in tomograms ( GRIP‐Tomo )

Abstract In this study, we present a method of pattern mining based on network theory that enables the identification of protein structures or complexes from synthetic volume densities, without the knowledge of predefined templates or human biases for refinement. We hypothesized that the topological connectivity of protein structures is invariant, and they are distinctive for the purpose of protein identification from distorted data presented in volume densities. Three‐dimensional densities of a protein or a complex from simulated tomographic volumes were transformed into mathematical graphs as observables. We systematically introduced data distortion or defects such as missing fullness of data, the tumbling effect, and the missing wedge effect into the simulated volumes, and varied the distance cutoffs in pixels to capture the varying connectivity between the density cluster centroids in the presence of defects. A similarity score between the graphs from the simulated volumes and the graphs transformed from the physical protein structures in point data was calculated by comparing their network theory order parameters including node degrees, betweenness centrality, and graph densities. By capturing the essential topological features defining the heterogeneous morphologies of a network, we were able to accurately identify proteins and homo‐multimeric complexes from 10 topologically distinctive samples without realistic noise added. Our approach empowers future developments of tomogram processing by providing pattern mining with interpretability, to enable the classification of single‐domain protein native topologies as well as distinct single‐domain proteins from multimeric complexes within noisy volumes.

59 BASIC BIOLOGICAL SCIENCES↗

Solar neutrinos and the strongest oscillation constraints on scalar NSI

Scalar non-standard neutrino interactions (sNSI) is a scenario where neutrinos can develop a medium dependent contribution to their mass due to a new scalar mediator. This scenario differs from the commonly discussed vector mediator case in that the oscillation effect scales with density rather than density and neutrino energy. Thus the strongest oscillation constraint comes from solar neutrinos which experience the largest density in a neutrino oscillation experiment. We derive constraints on all the sNSI parameters as well as the absolute neutrino mass scale by combining solar and reactor data and find solar neutrinos to be > 1 order of magnitude more sensitive to sNSI than terrestrial probes such as long-baseline experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Infrared triggered dwell and active cooling thermal control effects on microstructural uniformity in DED

Directed Energy Deposition (DED) offers rapid large scale fabrication, but difficulty in delivering consistent microstructures and properties hinders the use of DED fabricated components in safety or performance critical applications. Variability stems from the complex thermal cycles generated by the toolpath used to print the required geometry. Several practical methods have become established in DED to regulate overheating, such as active cooling of the baseplate structure or the use of an infrared camera to inject interlayer pauses to ensure the top layer of the component cools to a set temperature, which have been shown to affect microstructure. However, no critical assessment has been performed as to how effective these controls are in promoting microstructural uniformity in the context of complex layer timing commonly generated by non-prismatic geometries. Here we show how controls influence the thermal field, phase transformations, and dynamic annealing of a low-temperature transformation steel using infrared imaging and operando neutron diffraction. Counterintuitively, common thermal homogenization process controls can reduce microstructural uniformity because these approaches stabilize peak temperature while overlooking temperatures near the solid-state phase transformation fronts. Instead, the cyclic reheating induces spatially-variant dynamically annealed regions which can be modulated via control parameters. We show that these controls have spatially linked effects centimeters away from the active weld, which implies that microstructure control must co-optimize thermal input across many subsequent layers. In conclusion, our results demonstrate the pressing need for higher order controls that integrate predictive elements of simulation data to stabilize printed properties for future qualification of DED components.

Directed energy deposition↗

Integration of renewable energy generation and storage systems for emissions reduction in an islanded campus microgrid

Microgrid connected building communities are projected to play an integral part in the clean energy transition. These types of systems, when integrated with distributed energy resources (DERs) such as combined heat and power (CHP), district heating and cooling, renewable generation, and energy storage, can provide clean, reliable power to critical facilities and vulnerable communities. The intermittent nature of renewable generation is a challenge when integrating renewables into any grid system, but particularly in islanded microgrids. The University of Texas at Austin (UT) operates an islanded microgrid powered by a CHP plant, while also utilizing district heating and cooling systems and thermal energy storage (TES). High fidelity operating data was used to develop a validated reduced order model of UT’s integrated campus energy systems to serve as a testbed for use in a case study. Hypothetical renewable energy installations on land owned by UT in west Texas were modeled and integrated into the validated campus models along with battery energy storage (BES). Simulations showed that a combination of renewable energy from wind, and optimally controlled 24-hour thermal and battery storage systems could reduce carbon dioxide emissions on campus by 45.4%. The additional retrofit of burner systems to utilize hydrogen natural gas blends resulted in an overall annual emissions reduction of 54.7%. Carbon capture and storage eliminated the majority of the remaining emissions with increased plant energy expenditure. The presented simulations display the practical limitations of a CHP system complemented by renewable generation and short-term storage in eliminating emissions. Furthermore, results highlight the need for further research and development in long duration storage technologies and hydrogen fueled turbines to increase penetration of renewable energy and reduce emissions.

CHP↗

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗

Typological representation of the offshore oceanographic environment along the Alaskan North Slope

Erosion and flooding impacts to Arctic coastal environments are intensifying with nearshore oceanographic conditions acting as a key environmental driver. Robust and comprehensive assessment of the nearshore oceanographic conditions require knowledge of the following boundary conditions: incident wave energy, water level, incident wind energy, ocean temperature and salinity, bathymetry, and shoreline orientation. The number of offshore oceanographic boundary conditions can be large, requiring a significant computational investment to reproduce nearshore conditions. This present study develops location-independent typologies to reduce the number of boundary conditions needed to assess nearshore oceanographic environments in both a Historical (2007–2019) and Future (2020–2040) timespan along the Alaskan North Slope. We used WAVEWATCH III® and Delft3D Flexible Mesh model output from six oceanographic sites located along a constant ~50 m bathymetric line spanning the Chukchi to Beaufort Seas. K-means clustering was applied to the energy-weighted joint-probability distribution of significant wave height (H s ) and peak period (T p ). Distributions of wave and wind direction, wind speed, and water level associated with location-independent centroids were assigned single values to describe a reduced order, typological rendition of offshore oceanographic conditions. Reanalysis data (e.g., ASRv2, ERA5, and GOFS) grounded the historical simulations while projected conditions were obtained from downscaled GFDL-CM3 forced under RCP8.5 conditions. Location-dependence for each site is established through the occurrence joint-probability distribution in the form of unique scaling factors representing the fraction of time that the typology would occupy over a representative year. As anticipated, these typologies show increasingly energetic ocean conditions in the future. They also enable computationally efficient simulation of the nearshore oceanographic environment along the North Slope of Alaska for better characterization of coastal processes (e.g., erosion, flooding, or sediment transport).

54 ENVIRONMENTAL SCIENCES↗

Accelerating high-order mesh optimization using finite element partial assembly on GPUs

In this paper we present a new GPU-oriented mesh optimization method based on high order finite elements. Our approach relies on node movement with fixed topology, through the Target-Matrix Optimization Paradigm (TMOP) and uses a global nonlinear solve over the whole computational mesh, i.e., all mesh nodes are moved together. A key property of the method is that the mesh optimization process is recast in terms of finite element operations, which allows us to utilize recent advances in the field of GPU-accelerated high order finite element algorithms. For example, we reduce data motion by using tensor factorization and matrix-free methods, which have superior performance characteristics compared to traditional full finite element matrix assembly and offer advantages for GPU based HPC hardware. Furthermore, we describe the major mathematical components of the method along with their efficient GPU-oriented implementation. In addition, we propose an easily reproducible mesh optimization test that can serve as a performance benchmark for the mesh optimization community.

97 MATHEMATICS AND COMPUTING↗