Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Transition form factors: γ * + p → Δ ( 1232 ) , Δ ( 1600 )

We report that electroproduction form factors describing the γ * p → Δ + ( 1232 ) , Δ + ( 1600 ) transitions are computed using a fully dynamical diquark-quark approximation to the Poincaré-covariant three-body bound-state problem in relativistic quantum field theory. In this approach, the Δ ( 1600 ) is an analogue of the Roper resonance in the nucleon sector, appearing as the simplest radial excitation of the Δ ( 1232 ) . Precise measurements of the γ * p → Δ + ( 1232 ) transition already exist on 0 ≤ Q 2 ≲ 8 GeV 2 , and the calculated results compare favorably with the data outside the meson-cloud domain. The predictions for the γ * p → Δ + ( 1600 ) magnetic dipole and electric quadrupole transition form factors are consistent with the empirical values at the real photon point, and extend to Q 2 ≈ 6 m p 2 , enabling a meaningful direct comparison with experiment once analysis of existing data is completed. In both cases, the electric quadrupole form factor is particularly sensitive to deformation of the Δ -baryons. Interestingly, while the γ * p → Δ + ( 1232 ) transition form factors are larger in magnitude than those for γ * p → Δ + ( 1600 ) in some neighborhood of the real photon point, this ordering is reversed on Q 2 ≳ 2 m p 2 , suggesting that the γ * p → Δ + ( 1600 ) transition is more localized in configuration space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Charge localization and ordering in A 2 Mn 8 O 16 hollandite group oxides: Impact of density functional theory approaches

The phases of A 2 Mn 8 O 16 hollandite group oxides emerge from the competition between ionic interactions, Jahn-Teller effects, charge ordering, and magnetic interactions. Their balanced treatment with feasible computational approaches can be challenging for commonly used approximations in density functional theory. Three examples (A = Ag, Li, and K) are studied with a sequence of different approximate exchange-correlation functionals. Starting from a generalized gradient approximation (GGA), an extension to include van der Waals interactions and a recently proposed meta-GGA are considered. Then local Coulomb interactions for the Mn 3d electrons are more explicitly considered with the DFT + U approach. Finally, selected results from a hybrid functional approach provide a reference. Results for the binding energy of the A species in the parent oxide highlight the role of van der Waals interactions. Relatively accurate results for insertion energies can be achieved with a low-U and a high-U approach. In the low-U case, the materials are described as band metals with a high-symmetry, tetragonal crystal structure. In the high-U case, the electrons donated by A result in formation of local Mn 3+ centers and corresponding Jahn-Teller distortions characterized by a local order parameter. The resulting degree of monoclinic distortion depends on charge ordering and magnetic interactions in the phase formed. The reference hybrid functional results show charge localization and ordering. Comparison to low-temperature experiments of related compounds suggests that charge localization is the physically correct result for the hollandite group oxides studied here. Lastly, while competing effects in the local magnetic coupling are subtle, the fully anisotropic implementation of DFT + U gives the best overall agreement with results from the hybrid functional.

meta-GGA↗
DOE OSTI2023

SuperNeuro: A Fast and Scalable Simulator for Neuromorphic Computing

In many neuromorphic workflows, simulators play a vital role for important tasks such as training spiking neural networks, running neuroscience simulations, and designing, implementing, and testing neuromorphic algorithms. Currently available simulators cater to either neuroscience workflows (e.g., NEST and Brian2) or deep learning workflows (e.g., BindsNET). Problematically, the neuroscience-based simulators are slow and not very scalable, and the deep learning-based simulators do not support certain functionalities that are typical of neuromorphic workloads (e.g., synaptic delay). In this paper, we address this gap in the literature and present SuperNeuro, which is a fast and scalable simulator for neuromorphic computing capable of both homogeneous and heterogeneous simulations as well as GPU acceleration. We also present preliminary results that compare SuperNeuro to widely used neuromorphic simulators such as NEST, Brian2, and BindsNET in terms of computation times. We demonstrate that SuperNeuro can be approximately 10×--300× faster than some of the other simulators for small sparse networks. On large sparse and large dense networks, SuperNeuro can be approximately 2.2×--3.4× faster than the other simulators, respectively.

Date, Prasanna↗
DOE OSTI2018

Viscosity Control Experiment Feasibility Study

Turbulent mix has been invoked to explain many results in Inertial Confinement Fusion (ICF) and High Energy Density (HED) physics, such as reduced yield in capsule implosions. Many ICF capsule implosions exhibit interfacial instabilities seeded by the drive shock, but it is not clear that fully developed turbulence results from this. Many simulations use turbulent mix models to help match simulation results to data, but this is not appropriate if turbulence is not present. It would be useful to have an experiment where turbulent mixing could be turned on or off by design. The use of high-Z dopants to modify viscosity and the resulting influence on turbulence is considered here. A complicating factor is that the plasma in some implosions can become strongly coupled, which makes the Spitzer expression for viscosity invalid. We first consider equations that cover a broad parameter space in temperature and density to address regimes for various experimental applications. Next, a previous shock-tube and other ICF experiments that investigate viscosity or use doping to examine the effects on yield are reviewed. How viscosity and dopants play a role in capsule yield depends on the region and process under consideration. Experiments and simulations have been performed to study the effects of viscosity on both the hot spot and the fuel/ablator mix. Increases in yield have been seen for some designs, but not all. We then discuss the effect of adding krypton dopant to the gas region of a typical OMEGA and a 2-shock NIF implosion to determine approximately the effect of adding dopant on the computed Reynolds number. Recommendations for a path forward for possible experiments using high-Z dopants to affect viscosity and turbulence are made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bolt: A Fast Solver for Kinetic Theories Using a High-Resolution Constrained Transport Scheme [Slides]

Our understanding of collisionless and semi-collisional plasmas in the nonlinear regime is limited by the expense of computing solutions numerically. Bolt is a fast, GPU-accelerated code for rapidly computing such solutions with accurate transport and an approximate collision operator. Such calculations are relevant to both problems in astrophysics, such as heat conduction and magnetic reconnection in accretion disks around black holes, and also to programmatic interests at LANL. The bolt code paper, demonstrating accuracy via a suite of test problems calculated on kodiak, is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of the uncertainty quantification toolkit's python interface and surrogate construction tutorial

The uncertainty quantification toolkit (UQTk) is a collection of c++ libraries that assess the confidence of numerical models. Surrogate approximations, often polynomial chaos expansions (PCEs), lessen the computational cost of these assessments. I developed a Python interface in UQTk for regression and Bayesian compressive sensing to add to the existing Galerkin projection method. These methods receive an object containing the polynomial basis information and NumPy arrays of sample points, call c++ methods, and return the PCE coefficients in a NumPy array. To demonstrate these methods, I wrote a tutorial in which I use them to construct surrogates for Genz functions and calculate the resulting error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Quantifying How Density Gradients and Front Curvature Affect Carbon Detonation Strength during SNe Ia

Accurately reproducing the physics behind the detonations of SNe Ia and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full-star simulations in the present and near future. In the lower density regions of the star, the curvature of the detonation front will slow the detonation, affecting the production of intermediate mass elements. We find that shock strengthening due to the density gradient present in the outer layers of the progenitor is essential for understanding the nucleosynthesis there, with burning extending well below the density at which a steady-state detonation is extinct. We show that a complete reaction network is not sufficient to obtain physical detonations at high densities and modest resolution due to numerical mixing at the unresolved reaction front. At low densities, below 6 × 10 5 g cm -3 , it is possible to achieve high enough resolution to separate the shock and the reaction region, and the abundance structure predicted by fully resolved quasi-steady-state calculations is obtained. For our best current benchmark yields, we utilize a method in which the unresolved portion of Lagrangian histories are reconstructed based on fully resolved quasi-steady-state detonation calculations. These computations demonstrate that under-resolved simulations agree approximately, ~10% in post-shock values of temperature, pressure, density, and abundances, with expected detonation structures sufficiently far from the under-resolved region, but that there is still room for some improvement in the treatment of subgrid reactions in the hydrodynamics to before better than 1% can be achieved at all densities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Power series evaluation of transition and covariance matrices.

Reexamination power series solutions to the matrix covariance differential equation and the transition differential equation. Truncation error bounds are derived which are computationally attractive and which extend previous results. Polynomial approximations are obtained by exploiting the functional equations satisfied by the transition and covariance matrices. The series-functional equation propagation technique represents a fast and accurate alternative to the numerical integration of the time-invariant transition and covariance equations.

Bierman, G. J.↗
NASA NTRS1974

Phase-bunching and other non-linear processes occurring in gyroresonant wave-particle interactions

Consideration of the movement of individual electrons in a magnetized plasma in which a monochromatic wave is propagating in the whistler mode. Simple expressions are derived which give the displacement of the electrons as a function of time, the phase angle that their velocity vector makes with the magnetic component of the wave, their pitch angle, and energy changes. A useful formula is obtained which gives the velocity range over which particles remain trapped inside the wave, as a function of the wave intensity and of the initial phase angle of the particle. It is shown that even strictly resonant particles can escape from the wave when their initial phase angle is very small. From the derived expressions, it is possible to compute the phase-bunching effect which occurs approximately at one trapping wavelength behind the leading edge of the interaction region.

Gendrin, R.↗
NASA NTRS1974

A Monte Carlo investigation of thrust imbalance of solid rocket motor pairs

A technique is described for theoretical, statistical evaluation of the thrust imbalance of pairs of solid-propellant rocket motors (SRMs) firing in parallel. Sets of the significant variables, determined as a part of the research, are selected using a random sampling technique and the imbalance calculated for a large number of motor pairs. The performance model is upgraded to include the effects of statistical variations in the ovality and alignment of the motor case and mandrel. Effects of cross-correlations of variables are minimized by selecting for the most part completely independent input variables, over forty in number. The imbalance is evaluated in terms of six time - varying parameters as well as eleven single valued ones which themselves are subject to statistical analysis. A sample study of the thrust imbalance of 50 pairs of 146 in. dia. SRMs of the type to be used on the space shuttle is presented. The FORTRAN IV computer program of the analysis and complete instructions for its use are included. Performance computation time for one pair of SRMs is approximately 35 seconds on the IBM 370/155 using the FORTRAN H compiler.

Sforzini, R. H.↗
NASA NTRS1975

Doppler measurements of the ionosphere on the occasion of the Apollo-Soyuz test project. Part 1: Computer simulation of ionospheric-induced Doppler shifts

A computer simulation of the ionospheric experiment of the Apollo-Soyuz Test Project (ASTP) was performed. ASTP is the first example of USA/USSR cooperation in space and is scheduled for summer 1975. The experiment consists of performing dual-frequency Doppler measurements (at 162 and 324 MHz) between the Apollo Command Service Module (CSM) and the ASTP Docking Module (DM), both orbiting at 221-km height and at a relative distance of 300 km. The computer simulation showed that, with the Doppler measurement resolution of approximately 3 mHz provided by the instrumentation (in 10-sec integration time), ionospheric-induced Doppler shifts will be measurable accurately at all times, with some rare exceptions occurring when the radio path crosses regions of minimum ionospheric density. The computer simulation evaluated the ability of the experiment to measure changes of columnar electron content between CSM and DM (from which horizontal gradients of electron density at 221-km height can be obtained) and to measure variations in DM-to-ground columnar content (from which an averaged columnar content and the electron density at the DM can be deduced, under some simplifying assumptions).

Grossi, M. D.↗
NASA NTRS1976

Combustion of hydrogen injected into a supersonic airstream (a guide to the HISS computer program)

A computer program based on a finite-difference, implicit numerical integration scheme is described for the prediction of hydrogen injected into a supersonic airstream at an angle ranging from normal to parallel to the airstream main flow direction. Results of calculations for flow and thermal property distributions were compared with 'cold flow data' taken by NASA/Langley and show excellent correlation. Typical results for equilibrium combustion are presented and exhibit qualitatively plausible behavior. Computer time required for a given case is approximately one minute on a CDC 7600. A discussion of the assumption of parabolic flow in the injection region is given which demonstrates that improvement in calculation in this region could be obtained by a partially-parabolic procedure which has been developed. It is concluded that the technique described provides an efficient and reliable means for analyzing hydrogen injection into supersonic airstreams and the subsequent combustion.

Dyer, D. F.↗
NASA NTRS1976

Proposed powered explicit guidance thrust integrals derivation/implementation. Mission planning, mission analysis and software formulation

A new exoatmospheric, powered explicit guidance (PEG) thrust integral formulation and a simple method of implementation are presented. The new thrust integral formulation is significantly simpler than that currently used in PEG. Preliminary estimates indicate a computer storage savings of 220 words, which is approximately 10 percent of the current PEG ascent program. Alternate methods of implementation that could produce even more savings are noted.

Jaggers, R. F.↗
NASA NTRS1976

Finite element representations for thin shell instability analysis

The development of finite element calculational procedures for thin-shell instability analysis has involved the definition of appropriate element representations (i.e. the geometric form of the element and the approximation of displacement and/or stress), constitutive expressions, and computational algorithms. Among these, the status of thin-shell finite element representations remains unsettled. This paper, therefore, emphasizes recent developments in the basic aspects of thin-shell finite element analysis and discusses one simplified approach in more detail. Formulative and computational procedures for elastic instability analysis are then summarized, and numerical results are shown for the simplified shell element formulation.

Gallagher, R. H.↗
NASA NTRS1977

Combustion of hydrogen injected into a supersonic airstream (the SHIP computer program)

The mathematical and physical basis of the SHIP computer program which embodies a finite-difference, implicit numerical procedure for the computation of hydrogen injected into a supersonic airstream at an angle ranging from normal to parallel to the airstream main flow direction is described. The physical hypotheses built into the program include: a two-equation turbulence model, and a chemical equilibrium model for the hydrogen-oxygen reaction. Typical results for equilibrium combustion are presented and exhibit qualitatively plausible behavior. The computer time required for a given case is approximately 1 minute on a CDC 7600 machine. A discussion of the assumption of parabolic flow in the injection region is given which suggests that improvement in calculation in this region could be obtained by use of the partially parabolic procedure of Pratap and Spalding. It is concluded that the technique described herein provides the basis for an efficient and reliable means for predicting the effects of hydrogen injection into supersonic airstreams and of its subsequent combustion.

Markatos, N. C.↗

Refine your search to explore more results.