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,243 records · Page 69

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗

Progress in computational unsteady aerodynamics

After vigorous development for over twenty years, Computational Fluid Dynamics (CFD) in the field of aerospace engineering has arrived at a turning point toward maturity. This paper discusses issues related to algorithm development for the Euler/Navier Stokes equations, code validation and recent applications of CFD for unsteady aerodynamics. Algorithm development is a fundamental element for a good CFD program. Code validation tries to bridge the reliability gap between CFD and experiment. Many of the recent applications also take a multidisciplinary approach, which is a future trend for CFD applications. As computers become more affordable, CFD is expected to be a better scientific and engineering tool.

Obayashi, Shigeru↗

Horizon radiation reaction forces

Using Effective Field Theory (EFT) methods, we compute the effects of horizon dissipation on the gravitational interactions of relativistic binary black hole systems. We assume that the dynamics is perturbative, i.e it admits an expansion in powers of Newton’s constant (post-Minkowskian, or PM, approximation). As applications, we compute corrections to the scattering angle in a black hole collision due to dissipative effects to leading PM order, as well as the post-Newtonian (PN) corrections to the equations of motion of binary black holes in non-relativistic orbits, which represents the leading order finite size effect in the equations of motion. The methods developed here are also applicable to the case of more general compact objects, eg. neutron stars, where the magnitude of the dissipative effects depends on non-gravitational physics (e.g, the equation of state for nuclear matter).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A computational study of yttria-stabilized zirconia: I. Using crystal chemistry to search for the ground state on a glassy energy landscape

Yttria-stabilized zirconia (YSZ), a ZrO 2 -Y 2 O 3 solid solution that contains a large population of oxygen vacancies, is widely used in energy and industrial applications. Past computational studies correctly predicted the anion diffusivity but not the cation diffusivity, which is important for material processing and stability. One of the challenges lies in identifying a plausible configuration akin to the ground state in a glassy landscape. This is unlikely to come from random sampling of even a very large sample space, but the odds are much improved by incorporating packing preferences revealed by a modest sized configurational library established from empirical potential calculations. Ab initio calculations corroborated these preferences, which prove remarkably robust extending to the fifth cation-oxygen shell about 8 Å away. Yet because of frustration there are still rampant violations of packing preferences and charge neutrality in the ground state, and the approach toward it bears a close analogy to glass relaxations. Fast relaxations proceed by fast oxygen movement around cations, while slow relaxations require slow cation diffusion. The latter is necessarily cooperative because of strong coupling imposed by the long-range packing preferences.

36 MATERIALS SCIENCE↗

Global Transformations for Legacy Parallel Applications via Structural Analysis and Rewriting

Performance and scalability optimization of large HPC applications is currently a labor-intensive, manual process with very low productivity. Major difficulties come from the disaggregated environment for HPC application development: the compiler is only involved in local decisions (core or multithreaded domain), while a library-based, communication-oriented programming model realizes whole-machine parallelism. Realizing any major global change in such a disaggregated environment is very difficult and involves changing large portions of the source code. We present semi-automated techniques, based on structural analysis and rewriting, for performing global transformations on an HPC application source code. We present two case studies using the Self-Consistent Field (SCF) standalone benchmark as well as the Coupled Cluster (CCSD) module (2.9 million lines of Fortran code), a key module of the NWChem computational chemistry application. We demonstrate how structural rewriting techniques can be used to automate transformations that affect multiple sections of the application's source code. We show that the transformations can be applied in a systematic fashion across the source code bases with minimal manual effort. These transformations improve the scalability of the SCF benchmark by more than two orders of magnitude and the performance of the full CCSD module by a factor of 4.0.

legacy parallel applications↗

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗

SparseLU, A Novel Algorithm and Math Library for Sparse LU Factorization

Decomposing sparse matrices into lower and upper triangular matrices (sparse LU factorization) is a key operation in many computational scientific applications. We developed SparseLU, a sparse linear algebra library that implements a new algorithm for LU factorization on general sparse matrices. The new algorithm divides the input matrix into tiles to which OpenMP tasks are created for factorization computation, where only tiles that contain nonzero elements are computed. For comparative performance analysis, we used the reference library SuperLU. Testing was performed on synthetically generated matrices which replicate the conditions of the real-world matrices. SparseLU is able to reach a mean speedup of ~29× compared to SuperLU.

Valero Lara, Pedro↗

Kinetic description of electron beams in the solar chromosphere

We formulate the relativistic Fokker-Plank equation for a beam of accelerated electrons interacting with a partially ionized plasma. In our derivation we conserved those terms contributing to velocity diffusion and found that this effect cannot be neglected a priori. We compute the terms accounting for elastic and inelastic collisions with neutral hydrogen and helium. Collisions with neutral hydrogen are found to be dominant throughout the chromosphere, except at the uppermost layers close to the transition region. As an application, we compute the loss of energy and momentum for a power-law beam impinging on the solar chromosphere, for a particular case in which the Fokker-Planck equation can be integrated analytically. We find that most of the beam energy is deposited in a relatively thin region of the chromosphere, a result which is largely insensitive to the theoretical method employed to compute the energy deposition rate.

Gomez, Daniel O.↗

A Computational Study for the Utilization of Jet Pulsations in Gas Turbine Film Cooling and Flow Control

This report is the second part of a three-part final report of research performed under an NRA cooperative Agreement contract. The first part is NASA/CR-2012-217415. The third part is NASA/CR-2012-217417. Jets have been utilized in various turbomachinery applications in order to improve gas turbines performance. Jet pulsation is a promising technique because of the reduction in the amount of air removed from compressor. In this work two areas of pulsed jets applications were computationally investigated using the commercial code Fluent (ANSYS, Inc.); the first one is film cooling of High Pressure Turbine (HPT) blades and second one is flow separation control over Low Pressure Turbine (LPT) airfoil using Vortex Generator Jets (VGJ). Using pulsed jets for film cooling purposes can help to improve the effectiveness and thus allow higher turbine inlet temperature. Effects of the film hole geometry, blowing ratio and density ratio of the jet, pulsation frequency and duty cycle of blowing on the film cooling effectiveness were investigated. As for the low-pressure turbine (LPT) stages, the boundary layer separation on the suction side of airfoils can occur due to strong adverse pressure gradients. The problem is exacerbated as airfoil loading is increased. Active flow control could provide a means for minimizing separation under conditions where it is most severe (low Reynolds number), without causing additional losses under other conditions (high Reynolds number). The effects of the jet geometry, blowing ratio, density ratio, pulsation frequency and duty cycle on the size of the separated region were examined in this work. The results from Reynolds Averaged Navier-Stokes and Large Eddy Simulation computational approaches were compared with the experimental data.

Kartuzova, Olga V.↗

An operational matrix method for solving variable-order fractional biharmonic equation

In this paper, we focus on the fourth-order biharmonic equation which also has been appeared with fractional derivative order. Then, we generalize the fractional biharmonic equation using the concept of variable-order (V-O) fractional derivative in the type of Caputo namely V-O fractional biharmonic equation (V-OFBE). To establish a method for solving V-OFBE, we firstly derive the operational matrix of V-O fractional derivatives for the shifted second kind Chebyshev polynomials. For obtaining this operational matrix with a general procedure, we implement the analytical form of these polynomials and some properties of the V-O fractional derivatives. This approach reduces the V-OFBE to a system of algebraic equations which considerably decreases the following computations. The applicability and reliability of the proposed method are investigated by solving some test problems. The experimental results show the spectral rate of convergence.

97 MATHEMATICS AND COMPUTING↗

Validation of transfer map calculation for electrostatic deflectors in the code COSY INFINITY

The code COSY INFINITY uses a beamline coordinate system with a Frenet–Serret frame relative to the reference particle, and calculates differential algebra-valued transfer maps by integrating the ODEs of motion in the respective vector space over a differential algebra (DA). In our work, we described and performed computation of the DA transfer map of an electrostatic spherical deflector in a laboratory coordinate system using two conventional methods: (1) by integrating the ODEs of motion using a numerical integrator and (2) by computing analytically and in closed form the properties of the respective elliptical orbits from Kepler theory. We compared the resulting transfer maps with (3) the DA transfer map of COSY INFINITY’s built-in electrostatic spherical deflector element ESP and (4) the transfer map of the electrostatic spherical deflector computed using the program GIOS, which uses analytic formulas from a paper1 by Hermann Wollnik regarding second-order aberrations. In addition to the electrostatic spherical deflector, we studied an electrostatic cylindrical deflector, where the Kepler theory is not applicable. We computed the DA transfer map by the ODE integration method (1), and we compared it with the transfer maps by (3) COSY INFINITY’s built-in electrostatic cylindrical deflector element ECL and (4) GIOS. The transfer maps of electrostatic spherical and cylindrical deflectors obtained using the direct calculation methods (1) and (2) are in excellent agreement with those computed using (3) COSY INFINITY. On the other hand, we found a significant discrepancy with (4) the program GIOS.

43 PARTICLE ACCELERATORS↗

On Holder-Brascamp-Lieb inequalities for torsion-free discrete Abelian groups

H\"older-Brascamp-Lieb inequalities provide upper bounds for a class of multilinear expressions, in terms of $L^p$ norms of the functions involved. They have been extensively studied for functions defined on Euclidean spaces. Bennett-Carbery-Christ-Tao have initiated the study of these inequalities for discrete Abelian groups and, in terms of suitable data, have characterized the set of all tuples of exponents for which such an inequality holds for specified data, as the convex polyhedron defined by a particular finite set of affine inequalities. In this paper we advance the theory of such inequalities for torsion-free discrete Abelian groups in three respects. The optimal constant in any such inequality is shown to equal $1$ whenever it is finite. An algorithm that computes the admissible polyhedron of exponents is developed. It is shown that nonetheless, existence of an algorithm that computes the full list of inequalities in the Bennett-Carbery-Christ-Tao description of the admissible polyhedron for all data, is equivalent to an affirmative solution of Hilbert's Tenth Problem over the rationals. That problem remains open. Applications to computer science will be explored in a forthcoming companion paper.

Christ, Michael↗

Applications of CALPHAD modeling and databases in advanced lightweight metallic materials

Here this paper presents an overview on the application of CALPHAD (CALculation of Phase Diagrams) methodology in the design and development of advanced lightweight metallic materials including magnesium, aluminum, titanium, aluminum-based metal matrix composites, and high entropy alloys. In this work, CALPHAD methodology has been established and summarized from the construction of databases describing thermodynamics, atomic mobility, kinetics, thermo-physical properties (such as viscosity) to the application of computational design of lightweight materials. The examples in this paper have demonstrated the effectiveness and capability of CALPHAD methodology in accelerating the design of advanced lightweight materials by optimizing the compositions and various heat-treatment conditions, modifying the evolution of microstructures during processing, and finally predicting the mechanical properties (e.g., yield strength and hardness) of the lightweight components. Although the examples are given in lightweight alloys for structural applications, the fundamental methodology and modeling principles are applicable to all materials and engineering applications. Thus, the future of the advanced material design will be strongly based on development of CALPHAD methodology such as the construction of reliable databases coupled with CALPHAD-based models for various applications.

36 MATERIALS SCIENCE↗

Software systems for modeling articulated figures

Research in computer animation and simulation of human task performance requires sophisticated geometric modeling and user interface tools. The software for a research environment should present the programmer with a powerful but flexible substrate of facilities for displaying and manipulating geometric objects, yet insure that future tools have a consistent and friendly user interface. Jack is a system which provides a flexible and extensible programmer and user interface for displaying and manipulating complex geometric figures, particularly human figures in a 3D working environment. It is a basic software framework for high-performance Silicon Graphics IRIS workstations for modeling and manipulating geometric objects in a general but powerful way. It provides a consistent and user-friendly interface across various applications in computer animation and simulation of human task performance. Currently, Jack provides input and control for applications including lighting specification and image rendering, anthropometric modeling, figure positioning, inverse kinematics, dynamic simulation, and keyframe animation.

Phillips, Cary B.↗

Exploiting universal nonlocal dispersion in optically active materials for spectro-polarimetric computational imaging

Recent years have seen significant advancements in exploring novel light-matter interactions such as hyperbolic dispersion within natural crystals. However, current studies have predominantly concentrated on local optical response of materials characterized by a dielectric tensor without spatial dispersion. Here, we investigate the nonlocal response in optically-active crystals with screw symmetries, revealing their lossless, super-dispersive properties compared to traditional optical response functions. We leverage this universal nonlocal dispersion, i.e. the dispersion of optical rotatory power, to explore a novel spectral de-multiplexing scheme compared to conventional gratings, prisms and metasurfaces. We design and demonstrate an ‘Nonlocal-Cam’ - a camera that exploits nonlocal dispersion through sampling of polarized spectral states and the application of computational spectral reconstruction algorithms. The Nonlocal-Cam captures information in both laboratory and outdoor field experiments which is unavailable to traditional intensity cameras - the spectral texture of polarization. Merging the fields of nonlocal electrodynamics and computational imaging, our work paves the way for exploiting nonlocal optics of optically active materials in a variety of applications, from biological microscopy to physics-driven machine vision and remote sensing.

Wang, Xueji [Purdue Univ., West Lafayette, IN (Uni↗

Computer Simulation of Aircraft Aerodynamics

The role of Ames Research Center in conducting basic aerodynamics research through computer simulations is described. The computer facilities, including supercomputers and peripheral equipment that represent the state of the art, are described. The methodology of computational fluid dynamics is explained briefly. Fundamental studies of turbulence and transition are being pursued to understand these phenomena and to develop models that can be used in the solution of the Reynolds-averaged Navier-Stokes equations. Four applications of computer simulations for aerodynamics problems are described: subsonic flow around a fuselage at high angle of attack, subsonic flow through a turbine stator-rotor stage, transonic flow around a flexible swept wing, and transonic flow around a wing-body configuration that includes an inlet and a tail.

Inouye, Mamoru↗