Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer applications”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,153 records · Page 64

A computational fluid dynamics simulation framework for ventricular catheter design optimization

OBJECTIVE Cerebrospinal fluid (CSF) shunts are the primary treatment for patients suffering from hydrocephalus. While proven effective in symptom relief, these shunt systems are plagued by high failure rates and often require repeated revision surgeries to replace malfunctioning components. One of the leading causes of CSF shunt failure is obstruction of the ventricular catheter by aggregations of cells, proteins, blood clots, or fronds of choroid plexus that occlude the catheter’s small inlet holes or even the full internal catheter lumen. Such obstructions can disrupt CSF diversion out of the ventricular system or impede it entirely. Previous studies have suggested that altering the catheter’s fluid dynamics may help to reduce the likelihood of complete ventricular catheter failure caused by obstruction. However, systematic correlation between a ventricular catheter’s design parameters and its performance, specifically its likelihood to become occluded, still remains unknown. Therefore, an automated, open-source computational fluid dynamics (CFD) simulation framework was developed for use in the medical community to determine optimized ventricular catheter designs and to rapidly explore parameter influence for a given flow objective. METHODS The computational framework was developed by coupling a 3D CFD solver and an iterative optimization algorithm and was implemented in a high-performance computing environment. The capabilities of the framework were demonstrated by computing an optimized ventricular catheter design that provides uniform flow rates through the catheter’s inlet holes, a common design objective in the literature. The baseline computational model was validated using 3D nuclear imaging to provide flow velocities at the inlet holes and through the catheter. RESULTS The optimized catheter design achieved through use of the automated simulation framework improved significantly on previous attempts to reach a uniform inlet flow rate distribution using the standard catheter hole configuration as a baseline. While the standard ventricular catheter design featuring uniform inlet hole diameters and hole spacing has a standard deviation of 14.27% for the inlet flow rates, the optimized design has a standard deviation of 0.30%. CONCLUSIONS This customizable framework, paired with high-performance computing, provides a rapid method of design testing to solve complex flow problems. While a relatively simplified ventricular catheter model was used to demonstrate the framework, the computational approach is applicable to any baseline catheter model, and it is easily adapted to optimize catheters for the unique needs of different patients as well as for other fluid-based medical devices.

Neurosciences & Neurology↗

A Comparison of Three Ways to Measure Time-Dependent Densities With Quantum Simulators

Quantum algorithms are touted as a way around some classically intractable problems such as the simulation of quantum mechanics. At the end of all quantum algorithms is a quantum measurement whereby classical data is extracted and utilized. In fact, many of the modern hybrid-classical approaches are essentially quantum measurements of states with short quantum circuit descriptions. Here, we compare and examine three methods of extracting the time-dependent one-particle probability density from a quantum simulation: direct Z -measurement, Bayesian phase estimation, and harmonic inversion. We have tested these methods in the context of the potential inversion problem of time-dependent density functional theory. Our test results suggest that direct measurement is the preferable method. We also highlight areas where the other two methods may be useful and report on tests using Rigetti's quantum virtual device. This study provides a starting point for imminent applications of quantum computing.

Yang, Jun↗

Kullback–Leibler Divergence of an Open-Queuing Network of a Cell-Signal-Transduction Cascade

Queuing networks (QNs) are essential models in operations research, with applications in cloud computing and healthcare systems. However, few studies have analyzed the cell’s biological signal transduction using QN theory. This study entailed the modeling of signal transduction as an open Jackson’s QN (JQN) to theoretically determine cell signal transduction, under the assumption that the signal mediator queues in the cytoplasm, and the mediator is exchanged from one signaling molecule to another through interactions between the signaling molecules. Each signaling molecule was regarded as a network node in the JQN. The JQN Kullback–Leibler divergence (KLD) was defined using the ratio of the queuing time (λ) to the exchange time (μ), λ/μ. The mitogen-activated protein kinase (MAPK) signal-cascade model was applied, and the KLD rate per signal-transduction-period was shown to be conserved when the KLD was maximized. Our experimental study on MAPK cascade supported this conclusion. This result is similar to the entropy-rate conservation of chemical kinetics and entropy coding reported in our previous studies. Thus, JQN can be used as a novel framework to analyze signal transduction.

97 MATHEMATICS AND COMPUTING↗

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗

The [N ii] 205 μm emission in local luminous infrared galaxies

In this paper, we present the measurements of the [N ii] 205 μm line for a flux-limited sample of 122 (ultra-)luminous infrared galaxies [(U)LIRGs] and 20 additional normal galaxies, obtained with the Herschel Space Observatory (Herschel). We explore the far-infrared (FIR) color dependence of the [N ii] 205 μm (L{sub [N ii]205 μm}) to the total infrared (L{sub IR}) luminosity ratio, and find that L{sub [N ii]205 μm}/L{sub IR} only depends modestly on the 70–160 μm flux density ratio (f{sub 70}/f{sub 160}) when f{sub 70}/f{sub 160}≲0.6, whereas such dependence becomes much steeper for f{sub 70}/f{sub 160}>0.6. We also investigate the relation between L{sub [N ii]205 μm} and star formation rate (SFR), and show that L{sub [N ii]205 μm} has a nearly linear correlation with SFR, albeit the intercept of such a relation varies somewhat with f{sub 60}/f{sub 100}, consistent with our previous conclusion that [N ii] 205 μm emission can serve as an SFR indicator with an accuracy of ∼0.4 dex, or ∼0.2 dex if f{sub 60}/f{sub 100} is known independently. Furthermore, together with the Infrared Space Observatory measurements of [N ii], we use a total of ∼200 galaxies to derive the local [N ii] 205 μm luminosity function (LF) by tying it to the known IR LF with a bivariate method. As a practical application, we also compute the local SFR volume density ( ρ-dot {sub SFR}) using the newly derived SFR calibrator and LF. The resulting log ρ-dot {sub SFR}=−1.96±0.11 M{sub ⊙} yr{sup −1} Mpc{sup −3} agrees well with previous studies. Finally, we determine the electron densities (n{sub e}) of the ionized medium for a subsample of 12 (U)LIRGs with both [N ii] 205 μm and [N ii] 122 μm data, and find that n{sub e} is in the range of ∼1–100 cm{sup −3}, with a median value of 22 cm{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal spin-qubit hallmarks of sulfur-vacancy defects in 4H-SiC: Design from first principles

We propose neutral defects in 4H-SiC comprising a silicon vacancy and a sulfur atom dopant substituting a carbon atom as an optically controllable spin qubit with application in quantum computing and quantum telecommunications. According to our state-of-the-art first-principles calculation, the proposed defect possess very promising qubit functionalities. This system once again confirms our hypothesis for the rational design of spin qubits and single-photon emitters. Importantly, the atoms of this system have high-abundance isotopes with zero nuclear spin ensuring high spin-coherence time of the qubit.

36 MATERIALS SCIENCE↗

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗

Using Whisper-1.1 to Guide Improvements to Nuclear Data Evaluations [Slides]

This presentation covers the use of Whisper-1.1 to guide improvements to nuclear data evaluations. This presentation provides a brief overview on what Whisper-1.1 is. Additionally covered is the benchmark selection: using benchmark and application correlations (ck) computed using nuclear data sensitivity profiles and covariance data. Next the presentation touches on compute bias and bias uncertainty based on selected benchmarks using Extreme Value Theory Finally the presentation shows the additional margin for nuclear data uncertainty estimated by GLLS.

97 MATHEMATICS AND COMPUTING↗

Toward Improved Modeling of Spectral Solar Irradiance for Solar Energy Applications: Preprint

This study introduces the National Renewable Energy Laboratory's (NREL's) recent efforts to extend the capability of the Fast All-sky Radiation Model for Solar applications (FARMS) by computing spectral solar irradiances over both horizontal and inclined surfaces. A new model is developed by computing the optical thickness of the atmosphere using a spectral irradiance model for clear-sky conditions, SMARTS2. A comprehensive lookup table (LUT) of cloud bidirectional transmittance distribution functions (BTDFs) is precomputed for 2002 wavelength bands using an atmospheric radiative transfer model, libRadtran. The solar radiation transmitted through the atmosphere is given by considering all possible paths of photon transmission and the relevent scattering and absorption attenuation. Our results indicate that this new model has an accuracy that is similar to that of state-of-the-art radiative transfer models, but it is significantly more efficient.

14 SOLAR ENERGY↗

Parallel Asynchronous Distributed-Memory Maximal Independent Set Algorithm with Work Ordering

The maximal independent set (MIS) graph problem arises in many applications such as computer vision, information theory, molecular biology, and process scheduling. The growing scale of graph data suggests the use of distributed memory hardware as a cost-effective approach to providing necessary compute and memory resources. Existing distributed memory parallel MIS algorithms rely on synchronous communication and use techniques such as subgraph computations. In this paper, we present an asynchronous distributed-memory parallel graph algorithm that relies on a virtual directed acyclic graph (DAG) that is created during the algorithm execution. We introduce two additional algorithms that save computations by ordering generated work. The first algorithm applies ordering globally to reduce computations, and the second algorithm applies ordering locally at the level of threads to minimize the synchronization overhead. We use two different implementations of Luby’s algorithm variants as baseline to compare the performance of the presented algorithms: (1) vertex-centric Luby A and Luby B implementations, and (2) the CombBLAS linear-algebra Luby A implementation. Results show that proposed algorithms outperform both implementations of Luby algorithms, especially in distributed execution. Furthermore, we show that for low-diameter graphs the algorithm that applies global ordering scales better than other algorithms and for high diameter graphs the original asynchronous algorithm and thread-level ordering algorithm show better performance.

maximal independent set (MIS), distributed memory ↗

Emergence of excitonic superfluid at topological-insulator surfaces

Excitons are spin integer particles that are predicted to condense into a coherent quantum state at sufficiently low temperature, and exciton condensates can be realized at much higher temperature than condensates of atoms because of strong Coulomb binding and small mass. Signatures of exciton condensation have been reported in double quantum wells, microcavities, graphene, and transition metal dichalcogenides. Nonetheless, transport of exciton condensates is not yet understood and it is unclear whether an exciton condensate is a superfluid or an insulating electronic crystal. Topological insulators (TIs) with massless particles and unique spin textures have been theoretically predicted as a promising platform for achieving exciton condensation. Here we report experimental evidence of excitonic superfluid phase on the surface of three-dimensional (3D) TIs. We unambiguously confirmed that electrons and holes are paired into charge neutral bound states by the electric field independent photocurrent distributions. And we observed a millimetre-long transport distance of these excitons up to 40 K, which strongly suggests dissipationless propagation. The robust macroscopic quantum states achieved with simple device architecture and broadband photoexcitation at relatively high temperature are expected to find novel applications in quantum computations and spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Classical Optimizers for Noisy Intermediate-Scale Quantum Devices

We present a collection of optimizers tuned for usage on Noisy Inter-mediate-Scale Quantum (NISQ) devices. Optimizers have a range of applications in quantum computing, including the Variational Quantum Eigensolver (VQE) and Quantum Approximate Optimization (QAOA) algorithms. They have further uses in calibration, hyperparameter tuning, machine learning, etc. We employ the VQE algorithm as a case study. VQE is a hybrid algorithm, with a classical minimizer step driving the next evaluation on the quantum processor. While most results to date concentrated on tuning the quantum VQE circuit, our study indicates that in the presence of quantum noise the classical minimizer step is a weak link and a careful choice combined with tuning is required for correct results. We explore state-of-the-art gradient-free optimizers capable of handling noisy, black-box, cost functions and stress-test them using a quantum circuit simulation environment with noise injection capabilities on individual gates. Our results indicate that specifically tuned optimizers are crucial to obtaining valid science results on NISQ hardware, as well as projecting forward on fault tolerant circuits.

Lavrijsen, Wim↗

Controlling AM spatter and conduction

An intelligent feed forward model to control additive manufacturing (AM) laser powder bed fusion process and reduce spattering whereby defects are eliminated by controlling the laser power and reducing spattering through a computer model. This application describes using a proportional integral derivative (PID) controller to create a power map that reduces spattering.

Khairallah, Saad A.↗

PersGNN: Applying Topological Data Analysis and Geometric Deep Learning to Structure-Based Protein Function Prediction

Understanding protein structure-function relationships is a key challenge in computational biology, with applications across the biotechnology and pharmaceutical industries. While it is known that protein structure directly impacts protein function, many functional prediction tasks use only protein sequence. In this work, we isolate protein structure to make functional annotations for proteins in the Protein Data Bank in order to study the expressiveness of different structure-based prediction schemes. We present PersGNN - an end-to-end trainable deep learning model that combines graph representation learning with topological data analysis to capture a complex set of both local and global structural features. While variations of these techniques have been successfully applied to proteins before, we demonstrate that our hybridized approach, PersGNN, outperforms either method on its own as well as a baseline neural network that learns from the same information. PersGNN achieves a 9.3% boost in area under the precision recall curve (AUPR) compared to the best individual model, as well as high F1 scores across different gene ontology categories, indicating the transferability of this approach.

Swenson, Nicolas↗

Enabling Combustion Science Simulations for Future Exascale Machines

Reacting flow simulations for combustion applications require extensive computing capabilities. Leveraging the AMReX library, the Pele suite of combustion simulation tools targets the largest supercomputers available and future exascale machines. We introduce PeleC, the compressible solver in the Pele suite, and detail its capabilities, including complex geometry representation, chemistry integration, and discretization. We present a comparison of development efforts using both OpenACC and AMReX's C++ performance portability framework for execution on multiple GPU architectures. We discuss relevant details that have allowed PeleC to achieve high performance and scalability. PeleC's performance characteristics are measured through relevant simulations on multiple supercomputers. The success of PeleC's design for exascale is exhibited through demonstration of a 160 billion cell simulation and weak scaling onto 100\% of Summit, an NVIDIA-based GPU supercomputer at Oak Ridge National Laboratory. Our results provide confidence that PeleC will enable future combustion science simulations with unprecedented fidelity.

combustion↗

Reducing the cost of energy estimation in the variational quantum eigensolver algorithm with robust amplitude estimation

Quantum chemistry and materials is one of the most promising applications of quantum computing. Yet much work is still to be done in matching industry-relevant problems in these areas with quantum algorithms that can solve them. Most previous efforts have carried out resource estimations for quantum algorithms run on large-scale fault-tolerant architectures, which include the quantum phase estimation algorithm. In contrast, few have assessed the performance of near-term quantum algorithms, which include the variational quantum eigensolver (VQE) algorithm. Recently, a large-scale benchmark study [Gonthier et al. 2020] found evidence that the performance of the variational quantum eigensolver for a set of industry-relevant molecules may be too inefficient to be of practical use. This motivates the need for developing and assessing methods that improve the efficiency of VQE. In this work, we predict the runtime of the energy estimation subroutine of VQE when using robust amplitude estimation (RAE) to estimate Pauli expectation values. Under conservative assumptions, our resource estimation predicts that RAE can reduce the runtime over the standard estimation method in VQE by one to two orders of magnitude. Despite this improvement, we find that the runtimes are still too large to be practical. These findings motivate two complementary efforts towards quantum advantage: 1) the investigation of more efficient near-term methods for ground state energy estimation and 2) the development of problem instances that are of industrial value and classically challenging, but better suited to quantum computation.

Johnson, Peter D.↗

Task 2.1: Adsorption-Based ISPR for BETO-Relevant Bioproducts

This task focuses on the development of adsorption-based in situ product recovery (ISPR) integrated with simulated moving bed chromatography for the recovery and purification of carboxylate products that are relevant to BETO. ISPR has been pursued previously in the Separations Consortium to recover carboxylic acids near or below their pKa values with liquid-liquid extraction coupled to downstream distillation. However, there are many acid products in the BETO portfolio that require neutralization well above their pKa values wherein ISPR could still be a major benefit to the bioprocess performance, including muconic acid, beta-ketoadipic acid, 3-hydroxypropionic acid, itaconic acid, butyric acid, and others. In this task, we are combining dynamic filtration with a rotating ceramic disk, resin capacity measurements, tailored resin synthesis, and simulated moving bed chromatography into an ISPR system that can be used to recover BETO-relevant carboxylates from bioreactor cultivations. We are working across process scales and using computational modeling where applicable alongside techno-economic analysis and life cycle assessment to understand major cost, energy, and GHG emissions drivers. The impact of this project will be a bench-scale integrated approach to recover carboxylate products in situ, which will reduce the waste generation from biological carboxylate production processes and improve the productivities of biological systems.

bio-based acid↗

Residual loss due to magnetic flux trapping in niobium superconducting microwave cavities

One of the significant source of radio frequency loss in superconducting microwave cavities is the trapped residual magnetic field during the cooldown due to the incomplete Messier effect. The sensitivity of flux trapping depends on the distribution of impurities, network of dislocations and defects which pin vortices. Recent advancement in cavity processing techniques resulted in very high quality factors with impurity doping, however the rf loss related to the flux trapping increases. High quality factor in SRF cavities is not only limited to particle accelerators but also emerging as an application to quantum computing and quantum information science. In this presentation, we present the results of measurements of flux trapping sensitivity on several 1.3 GHz elliptical accelerating cavities with different grain size, impurity content and surface preparations.

Dhakal, P.↗