SEARCH · Search NASA

Results for “Computational Number Theory”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Shear and bulk viscosity for a pure glue theory using an effective matrix model

At nonzero temperatures, the deconfining phase transition can be analyzed using an effective matrix model to characterize the change in holonomy. The model includes gluons and two-dimensional ghost fields in the adjoint representation, or “teens.” As ghosts, the teen fields are responsible for the decrease of the pressure as 𝑇 →𝑇 𝑑 , with 𝑇 𝑑 the transition temperature for deconfinement. Using the solution of this matrix model for a large number of colors, the parameters of the teen fields are adjusted so that the expectation value of the Polyakov loop is close to the values from the lattice. The shear, 𝜂, and bulk, 𝜁, viscosities are computed at nonzero holonomy to leading logarithmic order in weak coupling. In the pure glue theory, the value of the Polyakov loop is relatively large in the deconfined phase, ≈1/2 at 𝑇 𝑑 . Consequently, if 𝑠 is the entropy density, while 𝜂/𝑠 decreases as 𝑇 →𝑇 𝑑 , it is still well above the conformal bound. In contrast, 𝜁/𝑠 is largest at 𝑇 𝑑 , comparable to 𝜂/𝑠, then falls off rapidly with increasing temperature and is negligible by ∼2⁢𝑇 𝑑 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1972

An innovative approach to compensator design

The primary goal is to present for a control system a computer-aided-compensator design technique from a frequency domain point of view. The thesis for developing this technique is to describe the open loop frequency response by n discrete frequency points which result in n functions of the compensator coefficients. Several of these functions are chosen so that the system specifications are properly portrayed; then mathematical programming is used to improve all of these functions which have values below minimum standards. In order to do this several definitions in regard to measuring the performance of a system in the frequency domain are given. Next, theorems which govern the number of compensator coefficients necessary to make improvements in a certain number of functions are proved. After this a mathematical programming tool for aiding in the solution of the problem is developed. Then for applying the constraint improvement algorithm generalized gradients for the constraints are derived. Finally, the necessary theory is incorporated in a computer program called CIP (compensator improvement program).

Mitchell, J. R.↗

DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗

NASA NTRS1985

A faster 'transition' to laminar flow

A discussion is given of the ongoing research related to laminar flow airfoils, nacelles, and wings where the laminar flow is maintained by a favorable pressure gradient, surface suction or a combination of the two. Design methologies for natural laminar flow airfoil sections and wings for both low and high speed applications are outlined. Tests of a 7-foot chord, 23-deg sweep laminar flow-control airfoil at high subsonic Mach numbers are described, along with the associated stability theory used to design the suction system. The state-of-the-art of stability theory is simply stated and a typical calculation illustrated. In addition, recent computer simulations of transition using the time dependent Navier-Stokes equations are briefly described. Advances in wind tunnel capabilities and instrumentation will be reviewed, followed by the presentation of a few results from both wind tunnels and flight. Finally, some suggestions for future work will complete the paper.

Bobbitt, P. J.↗

NASA NTRS1988

Experimental Test Of Aerodynamic Computer Program

Linear potential-flow theory applied to advanced airplane design. Report discusses use of PANAIR computer program to predict airflow about advanced fighter aircraft. Predictions compared with measurements about scale model of aircraft in wind tunnel at mach numbers of 0.6, 0.9, and 1.2 and at angles-of-attack from 0 to 10 degree.

Madson, M. D.↗

NASA NTRS1988

The pros and cons of code validation

Computational and wind tunnel error sources are examined and quantified using specific calculations or experimental data, and a substantial comparison of theoretical and experimental results, or a code validation, is discussed. Wind tunnel error sources considered include wall interference, sting effects, Reynolds number effects, flow quality and transition, and instrumentation such as strain gage balances, electronically scanned pressure systems, hot film gages, hot wire anemometers, and laser velocimeters. Computational error sources include math model equation sets, the solution algorithm, artificial viscosity/dissipation, boundary conditions, the uniqueness of solutions, grid resolution, turbulence modeling, and Reynolds number effects. It is concluded that although improvements in theory are being made more quickly than in experiments, wind tunnel research has the advantage of the more realistic transition process of a right turbulence model in a free-transition test.

Bobbitt, Percy J.↗

NASA NTRS1988

The pros and cons of code validation

Computational and wind tunnel error sources are examined and quantified using specific calculations of experimental data, and a substantial comparison of theoretical and experimental results, or a code validation, is discussed. Wind tunnel error sources considered include wall interference, sting effects, Reynolds number effects, flow quality and transition, and instrumentation such as strain gage balances, electronically scanned pressure systems, hot film gages, hot wire anemometers, and laser velocimeters. Computational error sources include math model equation sets, the solution algorithm, artificial viscosity/dissipation, boundary conditions, the uniqueness of solutions, grid resolution, turbulence modeling, and Reynolds number effects. It is concluded that, although improvements in theory are being made more quickly than in experiments, wind tunnel research has the advantage of the more realistic transition process of a right turbulence model in a free-transition test.

Bobbitt, Percy J.↗

NASA NTRS1953

An analysis of laminar free-convection flow and heat transfer about a flat plate paralled to the direction of the generating body force

The free-convection flow and heat transfer (generated by a body force) about a flat plate parallel to the direction of the body force are formally analyzed and the type of flow is found to be dependent on the Grashof number alone. For large Grashof numbers (which are of interest in aeronautics), the flow is of the boundary-layer type and the problem is reduced in a formal manner, which is analogous to Prandtl's forced-flow boundary-layer theory, to the simultaneous solution of two ordinary differential equations subject to the proper boundary conditions. Velocity and temperature distributions for Prandtl numbers of 0.01, 0.72, 0.733, 1, 1, 10, 100, and 1000 are computed, and it is shown that velocities and Nusselt numbers of the order of magnitude of those encountered in forced-convection flows may be obtained in free-convection flows. The theoretical and experimental velocity and temperature distributions are in good agreement. A flow and a heat-transfer parameter, from which the important physical quantities such as shear stress and heat-transfer rate can be computed, are derived as functions of Prandtl number alone.

Ostrach, Simon↗

NASA NTRS1996

Effect of Surface Waviness on Transition in Three-Dimensional Boundary-Layer Flow

The effect of a surface wave on transition in three-dimensional boundary-layer flow over an infinite swept wing was studied. The mean flow computed using interacting boundary-layer theory, and transition was predicted using linear stability theory coupled with the empirical eN method. It was found that decreasing the wave height, sweep angle, or freestream unit Reynolds number, and increasing the freestream Mach number or suction level all stabilized the flow and moved transition onset to downstream locations.

Masad, Jamal A.↗

NASA NTRS1982

A Numerical Method for Solving the Equations of Compressible Viscous Flow

Although much progress has already been made In solving problems in aerodynamic design, many new developments are still needed before the equations for unsteady compressible viscous flow can be solved routinely. This paper describes one such development. A new method for solving these equations has been devised that 1) is second-order accurate in space and time, 2) is unconditionally stable, 3) preserves conservation form, 4) requires no block or scalar tridiagonal inversions, 5) is simple and straightforward to program (estimated 10% modification for the update of many existing programs), 6) is more efficient than present methods, and 7) should easily adapt to current and future computer architectures. Computational results for laminar and turbulent flows at Reynolds numbers from 3 x 10(exp 5) to 3 x 10(exp 7) and at CFL numbers as high as 10(exp 3) are compared with theory and experiment.

MacCormack, R. W.↗

NASA NTRS2002

Efficient Reformulation of HOTFGM: Heat Conduction with Variable Thermal Conductivity

Functionally graded materials (FGMs) have become one of the major research topics in the mechanics of materials community during the past fifteen years. FGMs are heterogeneous materials, characterized by spatially variable microstructure, and thus spatially variable macroscopic properties, introduced to enhance material or structural performance. The spatially variable material properties make FGMs challenging to analyze. The review of the various techniques employed to analyze the thermodynamical response of FGMs reveals two distinct and fundamentally different computational strategies, called uncoupled macromechanical and coupled micromechanical approaches by some investigators. The uncoupled macromechanical approaches ignore the effect of microstructural gradation by employing specific spatial variations of material properties, which are either assumed or obtained by local homogenization, thereby resulting in erroneous results under certain circumstances. In contrast, the coupled approaches explicitly account for the micro-macrostructural interaction, albeit at a significantly higher computational cost. The higher-order theory for functionally graded materials (HOTFGM) developed by Aboudi et al. is representative of the coupled approach. However, despite its demonstrated utility in applications where micro-macrostructural coupling effects are important, the theory's full potential is yet to be realized because the original formulation of HOTFGM is computationally intensive. This, in turn, limits the size of problems that can be solved due to the large number of equations required to mimic realistic material microstructures. Therefore, a basis for an efficient reformulation of HOTFGM, referred to as user-friendly formulation, is developed herein, and subsequently employed in the construction of the efficient reformulation using the local/global conductivity matrix approach. In order to extend HOTFGM's range of applicability, spatially variable thermal conductivity capability at the local level is incorporated into the efficient reformulation. Analytical solutions to validate both the user-friendly and efficient reformulations am also developed. Volume discretization sensitivity and validation studies, as well as a practical application of the developed efficient reformulation are subsequently carried out. The presented results illustrate the accuracy and implementability of both the user-friendly formulation and the efficient reformulation of HOTFGM.

Zhong, Yi↗

NASA NTRS1979

NASF transposition network: A computing network for unscrambling p-ordered vectors

The viewpoints of design, programming, and application of the transportation network (TN) is presented. The TN is a programmable combinational logic network that connects 521 memory modules to 512 processors. The unscrambling of p-ordered vectors to 1-ordered vectors in one cycle is described. The TN design is based upon the concept of cyclic groups from abstract algebra and primitive roots and indices from number theory. The programming of the TN is very simple, requiring only 20 bits: 10 bits for offset control and 10 bits for barrel switch shift control. This simple control is executed by the control unit (CU), not the processors. Any memory access by a processor must be coordinated with the CU and wait for all other processors to come to a synchronization point. These wait and synchronization events can be a degradation in performance to a computation. The TN application is for multidimensional data manipulation, matrix processing, and data sorting, and can also perform a perfect shuffle. Unlike other more complicated and powerful permutation networks, the TN cannot, if possible at all, unscramble non-p-ordered vectors in one cycle.

Lim, R. S.↗

NASA NTRS1999

Development of a Three-Dimensional PSE Code for Compressible Flows: Stability of Three-Dimensional Compressible Boundary Layers

A program is developed to investigate the linear stability of three-dimensional compressible boundary layer flows over bodies of revolutions. The problem is formulated as a two dimensional (2D) eigenvalue problem incorporating the meanflow variations in the normal and azimuthal directions. Normal mode solutions are sought in the whole plane rather than in a line normal to the wall as is done in the classical one dimensional (1D) stability theory. The stability characteristics of a supersonic boundary layer over a sharp cone with 50 half-angle at 2 degrees angle of attack is investigated. The 1D eigenvalue computations showed that the most amplified disturbances occur around x(sub 2) = 90 degrees and the azimuthal mode number for the most amplified disturbances range between m = -30 to -40. The frequencies of the most amplified waves are smaller in the middle region where the crossflow dominates the instability than the most amplified frequencies near the windward and leeward planes. The 2D eigenvalue computations showed that due to the variations in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the most amplified disturbances are shifted to 120 degrees compared to 90 degrees for the parallel theory. It is also observed that the nonparallel amplification rates are smaller than that is obtained from the parallel theory.

Balakumar, P.↗

NASA NTRS1984

Fine-cut bluntness effects on the hypersonic static stability of a 10 degree cone

A systematic wind-tunnel investigation has been performed on a series of spherically blunted 10-degree cones to determine Mach number and angle-of-attack effects on the hypersonic static stability. The cones, which range in nose-to-base radius ratios from 0 to 0.5, have been tested from -4 to 20 degrees angle of attack at Mach numbers of 6 and 10 in air and 20 in helium. Relatively large excursions in the center-of-pressure location were observed for small changes in cone nose bluntness. The movement of the center of pressure was also noted to become more pronounced as Mach number increased; increases in angle of attack tended to lessen the magnitude of these excursions. Computational predictions of basic aerodynamic coefficients as well as the center-of-pressure locations using both viscous and inviscid theory compare very well with the experimental data.

Bryan, M. D.↗

NASA NTRS1984

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗

DOE OSTI2025

String breaking in the heavy quark limit with scalable circuits

Quantum simulations of non-Abelian gauge theories require efficient mappings onto quantum computers and practical state preparation and measurement procedures. A truncation of the Hilbert space of non-Abelian lattice gauge theories with matter in the heavy quark limit is developed. This truncation is applied to SU(2) lattice gauge theory in 1+1D to map the theory efficiently onto a quantum computer. Scalable variational circuits are found to prepare the vacuum and single meson states. It is also shown how these state preparation circuits can be used to perform measurements of the number of mesons produced during the system’s time evolution. A state with a single qq¯ pair is prepared on quantum hardware and the inelastic production of qq¯ pairs is observed using 104 qubits on IBM’s Heron quantum computer ibm_torino.

Ciavarella, Anthony N↗

NASA NTRS1974

Space Shuttle flutter as affected by wing-body aerodynamic interaction

In the NASA Langley Research Center 26-inch transonic blowdown wind-tunnel, flutter speeds were measured on 1/80-th scale semispan models of the orbiter wing, the complete Space Shuttle, and intermediate component combinations. Using the doublet lattice method combined with slender body theory to calculate unsteady aerodynamic forces, subsonic flutter speeds were computed for comparison. Aerodynamic interaction was found by test and analysis to raise the flutter speed in some configurations while lowering it in others. Although at Mach number less than 0.7, predicted speeds correlated to within 6% of those measured, rapid deterioration of the agreement occurred at higher subsonic Mach numbers, especially on the more complicated configurations. Additional analysis showed that aerodynamic forces arising from body flexibility potentially can have a large effect on flutter speed, but that the current shuttle design is not so affected.

Chipman, R. R.↗

NASA NTRS1975

Analytical and experimental study of the effects of wing-body aerodynamic interaction on space shuttle subsonic flutter

The effects on flutter of the aerodynamic interaction between the space shuttle bodies and wing, 1/80th-scale semispan models of the orbiter wing, the complete shuttle and intermediate component combinations were tested in the NASA Langley Research Center 26-inch Transonic Blowdown Wind Tunnel. Using the double lattice method combined with slender body theory to calculate unsteady aerodynamic forces, subsonic flutter speeds were computed for comparison. Using calculated complete vehicle modes, flutter speed trends were computed for the full scale vehicle at an altitude of 15,200 meters and a Mach number of 0.6. Consistent with findings of the model studies, analysis shows the shuttle to have the same flutter speed as an isolated cantilevered wing.

Chipman, R. R.↗