Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Algebra”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

PyMTensor 1.0.0

SAND2021-1654 O The Python Material Tensor code (PyMTensor) applies crystal symmetries to arbitrary material tensors. It uses symbolic algebra to exactly determine the zero tensor components and the relationships between the nonzero components for all crystallographic point groups. It is capable of computing components of higher-order material tensors needed when treating nonlinear effects. It has a simple and flexible interface for inputting tensor rank information and any known symmetries between tensor indices.Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Jensen, DanielS↗
DOE OSTI2013

Trilinos Developers SQE Guide: ASC Software Quality Engineering Practices Version 3.0

The Trilinos Project is an effort to develop algorithms and enabling technologies within an object-oriented software framework for the solution of large-scale, complex multi-physics engineering and scientific problems. A new software capability is introduced into Trilinos as a package. A Trilinos package is an integral unit and, although there are exceptions such as utility packages, each package is typically developed by a small team of experts in a particular algorithms area such as algebraic preconditioners, nonlinear solvers, etc. The Trilinos Developers SQE Guide is a resource for Trilinos package developers who are working under Advanced Simulation and Computing (ASC) and are therefore subject to the ASC Software Quality Engineering Practices as described in the Sandia National Laboratories Advanced Simulation and Computing (ASC) Software Quality Plan: ASC Software Quality Engineering Practices Version 3.0 document. The Trilinos Developer Policies webpage contains a lot of detailed information that is essential for all Trilinos developers. The Trilinos Software Lifecycle Model defines the default lifecycle model for Trilinos packages and provides a context for many of the practices listed in this document.

Willenbring, James Michael↗
DOE OSTI2023

Efficient Integration of Algebraic Constraint for Exponential Time Integration

This work is a continuation of a previous project where the high efficiency of exponential integration methods for DRMHD systems was demonstrated. Often algebraic constraints must also be enforced for DRMHD systems of interest. Straightforward application of exponential methods to DRMHD equations with constraints leads to prohibitively computationally expensive methods. In this work, we propose new exponential schemes that allow the constraints to be removed from evaluation of exponential matrix functions which drastically reduces computational complexity by eliminating the need to perform computations enforcing constraints in exponential calculations while still preserving a high order of accuracy and allowing for a large time step even when the problem is stiff. This idea is similar to the W-methods and is achieved by carefully designing a method with the desired order of accuracy, even with an incomplete Jacobian matrix used as an argument of exponential-like functions. The constraints are accounted for by including them in the evaluation of the right-hand-side forcing function of the spatially discretized system. We study performance of the new methods on test problems and outline future research directions that this work opens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Non-commutative disintegrations: Existence and uniqueness in finite dimensions

Motivated by advances in categorical probability, we introduce non-commutative almost everywhere (a.e.) equivalence and disintegrations in the setting of C*-algebras. We show that C*-algebras (resp. W*-algebras) and a.e. equivalence classes of 2-positive (resp. positive) unital maps form a category. We prove non-commutative disintegrations are a.e. unique whenever they exist. We provide an explicit characterization for when disintegrations exist in the setting of finite-dimensional C*-algebras, and we give formulas for the associated disintegrations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

Study-simulation of space station dynamics

Matrix algebra translator and executor /MATE/ takes equations describing structural control system environmental interaction problem for flexible spacecraft components and loads them into self programming computer.

Gaitens, M. J.↗
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. I - Efficiency comparison

A comparison of the efficiency of several recently developed numerical techniques for solving chemical kinetic rate equations is presented. The solution procedures examined include two general-purpose codes, EPISODE and LSODE, developed as multipurpose differential equation solvers, and three specialzed codes, CHEMEQ, CREK1D, and GCKP84, developed specifically for chemical kinetics. The efficiency comparison is made by applying these codes to two practical combustion kinetics problems. Both problems describe adiabatic, constant-pressure, gas-phase chemical reactions and include all three combustion regimes: induction, heat release, and equilibration. The comparison shows that LSODE is the fastest routine currently available for solving chemical kinetic rate equations. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for temperature can be significantly faster than evaluation of the temperature by integration of its time derivative. Significant increases in computational speed are realized by updating the reaction rate constants only when the temperature change exceeds an amount Delta-T that is problem dependent. An approximate expression for the automatic evaluation of Delta-T is presented and is shown to result in increased computational speed.

Radhakrishnan, K.↗
NASA NTRS1990

On the receptivity and non-parallel stability of travelling disturbances in rotating disk flow

The generation and evolution of small amplitude wavelength traveling disturbances in rotating disk flow is discussed. The steady rotational speed of the disk is perturbed so as to introduce high frequency oscillations in the flow field. Secondly, surface imperfections are introduced on the disk such as roughness elements. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations that are solved numerically. For the class of disturbances considered (wavelength on the order of Reynolds number), it is found that eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are excited the most. The width of these roughness elements in the radial direction is about .1 r(sub 0) where r(sub 0) is the radial location of the roughness. It is also found that the receptivity coefficient is larger for a negative spanwise wavenumber than for a positive one. Typical wave angles found for these disturbances are about -26 degrees.

Balakumar, P.↗
NASA NTRS1991

On the receptivity and nonparallel stability of traveling disturbances in rotating-disk flow

The generation and evolution of small amplitude long wavelength traveling disturbances in rotating-disk flow are discussed. The steady rotational speed of the disk is perturbed so as to introduce high-frequency oscillations in the flow field. Secondly, surface imperfections on the disk such as roughness elements are introduced. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations which are solved numerically. It is found that, for the class of disturbances considered here (wavelength on the order of Reynolds number), eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are preferentially excited.

Balakumar, P.↗
NASA NTRS1994

Modification of the MML turbulence model for adverse pressure gradient flows

Computational fluid dynamics is being used increasingly to predict flows for aerospace propulsion applications, yet there is still a need for an easy to use, computationally inexpensive turbulence model capable of accurately predicting a wide range of turbulent flows. The Baldwin-Lomax model is the most widely used algebraic model, even though it has known difficulties calculating flows with strong adverse pressure gradients and large regions of separation. The modified mixing length model (MML) was developed specifically to handle the separation which occurs on airfoils and has given significantly better results than the Baldwin-Lomax model. The success of these calculations warrants further evaluation and development of MML. The objective of this work was to evaluate the performance of MML for zero and adverse pressure gradient flows, and modify it as needed. The Proteus Navier-Stokes code was used for this study and all results were compared with experimental data and with calculations made using the Baldwin-Lomax algebraic model, which is currently available in Proteus. The MML model was first evaluated for zero pressure gradient flow over a flat plate, then modified to produce the proper boundary layer growth. Additional modifications, based on experimental data for three adverse pressure gradient flows, were also implemented. The adapted model, called MMLPG (modified mixing length model for pressure gradient flows), was then evaluated for a typical propulsion flow problem, flow through a transonic diffuser. Three cases were examined: flow with no shock, a weak shock and a strong shock. The results of these calculations indicate that the objectives of this study have been met. Overall, MMLPG is capable of accurately predicting the adverse pressure gradient flows examined in this study, giving generally better agreement with experimental data than the Baldwin-Lomax model.

Conley, Julianne M.↗
NASA NTRS1994

Influence of cross-flow on nonlinear Tollmien-Schlichting/vortex interaction

The transition of an incompressible three-dimensional boundary layer with strong cross-flow is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two lower-branch Tollmien-Schlichting waves which mutually interact nonlinearly to induce a longitudinal vortex flow. The vortex motion in turn gives rise to significant wave modulation via wall-shear forcing. The characteristic Reynolds number is large and, as a consequence, the waves' and the vortex motion are governed primarily by triple deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Following analysis and computation over a wide range of parameters, three distinct responses are found to emerge in the nonlinear behavior of the flow solution downstream: an algebraic finite-distance singularity, far-downstream saturation or far-downstream wave decay leaving pure vortex flow. These depend on the input conditions, the wave angles and the size of the cross flow.

Davis, D. A. R.↗
NASA NTRS1995

The analysis of control trajectories using symbolic and database computing

This final report comprises the formal semi-annual status reports for this grant for the periods June 30-December 31, 1993, January 1-June 30, 1994, and June 1-December 31, 1994. The research supported by this grant is broadly concerned with the symbolic computation, mixed numeric-symbolic computation, and database computation of trajectories of dynamical systems, especially control systems. A review of work during the report period covers: trajectories and approximating series, the Cayley algebra of trees, actions of differential operators, geometrically stable integration algorithms, hybrid systems, trajectory stores, PTool, and other activities. A list of publications written during the report period is attached.

Grossman, Robert↗
NASA NTRS2016

A Nonlinear Modal Aeroelastic Solver for FUN3D

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Goldman, Benjamin D.↗
NASA NTRS1978

Finite element concepts in computational aerodynamics

Finite element theory was employed to establish an implicit numerical solution algorithm for the time averaged unsteady Navier-Stokes equations. Both the multidimensional and a time-split form of the algorithm were considered, the latter of particular interest for problem specification on a regular mesh. A Newton matrix iteration procedure is outlined for solving the resultant nonlinear algebraic equation systems. Multidimensional discretization procedures are discussed with emphasis on automated generation of specific nonuniform solution grids and accounting of curved surfaces. The time-split algorithm was evaluated with regards to accuracy and convergence properties for hyperbolic equations on rectangular coordinates. An overall assessment of the viability of the finite element concept for computational aerodynamics is made.

Baker, A. J.↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

ExaTN

A challenging scientific problem that will benefit from exascale computing is predictive simulation of two- and three-dimensional quantum many-body Hamiltonians in condensed matter physics, together with the related strong entanglement problems in quantum chemistry and quantum computing. Efficient numerical simulations will require new insights into wavefunction compression algorithms that scale on extremely heterogeneous HPC architectures. To address this need, we propose to develop a scalable math library of parallel multi-linear algebraic numerical primitives for hierarchical tensor compression. Our library will subsequently be used for building an efficient, performance-portable HPC framework for simulating strongly correlated quantum systems. In particular, our will will be the first to deliver a massively parallel implementation of the multiscale entanglement renormalization ansatz (MERA), a hierarchical tensor decomposition quantum scheme capable of expressing local expectation values of strongly entangled higher-dimensional quantum systems efficiently. To address the extreme complexity and high computational cost of such schemes, we will automate their parallel implementation by delivering an integrated software framework, ExaTN (Exascale Tensor Networks), that will process hierarchical tensor representations on exascale HPC systems via an asymchronous task-based programming model. As a specific application, our framework will provide scientists with a capability of simulating strongly entangled quantum systems, thus positioning ORNL as a strategic leader in modeling of strongly correlated materials.

McCaskey, AlexanderJ↗
NASA NTRS2021

Sulfur Molecules in Space by X-rays: A Computational Study

X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur Kedge X-ray absorption spectra of H2S, SO2, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core−valence separation (CVS)approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.

Goranka Bilalbegovic↗
DOE OSTI2019

Higher spin de Sitter Hilbert space

We propose a complete microscopic definition of the Hilbert space of minimal higher spin de Sitter quantum gravity and its Hartle-Hawking vacuum state. The fundamental degrees of freedom are 2N bosonic fields living on the future conformal boundary, where N is proportional to the de Sitter horizon entropy. The vacuum state is normalizable. The model agrees in perturbation theory with expectations from a previously proposed dS-CFT description in terms of a fermionic Sp(N) model, but it goes beyond this, both in its conceptual scope and in its computational power. In particular it resolves the apparent pathologies affecting the Sp(N) model, and it provides an exact formula for late time vacuum correlation functions. We illustrate this by computing probabilities for arbitrarily large field excursions, and by giving fully explicit examples of vacuum 3- and 4-point functions. We discuss bulk reconstruction and show the perturbative bulk QFT canonical commutations relations can be reproduced from the fundamental operator algebra, but only up to a minimal error term ~e -O(N) , and only if the operators are coarse grained in such a way that the number of accessible “pixels” is less than O(N). Independent of this, we show that upon gauging the higher spin symmetry group, one is left with 2N physical degrees of freedom, and that all gauge invariant quantities can be computed by a 2N × 2N matrix model. This suggests a concrete realization of the idea of cosmological complementarity.

1/N expansion↗
DOE OSTI2023

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this article, we present algorithms and implementations for the end-to-end GPU acceleration of matrix-free low-order-refined preconditioning of high-order finite element problems. The methods described here allow for the construction of effective preconditioners for high-order problems with optimal memory usage and computational complexity. The preconditioners are based on the construction of a spectrally equivalent low-order discretization on a refined mesh, which is then amenable to, for example, algebraic multigrid preconditioning. The constants of equivalence are independent of mesh size and polynomial degree. For vector finite element problems in H(curl) and H(div) (e.g., for electromagnetic or radiation diffusion problems), a specially constructed interpolation–histopolation basis is used to ensure fast convergence. Detailed performance studies are carried out to analyze the efficiency of the GPU algorithms. The kernel throughput of each of the main algorithmic components is measured, and the strong and weak parallel scalability of the methods is demonstrated. The different relative weighting and significance of the algorithmic components on GPUs and CPUs is discussed. Results on problems involving adaptively refined nonconforming meshes are shown, and the use of the preconditioners on a large-scale magnetic diffusion problem using all spaces of the finite element de Rham complex is illustrated.

97 MATHEMATICS AND COMPUTING↗