Search NASA⌕ Search

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.

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

At least 271 records · Page 15

Correction factory techniques for improving aerodynamic prediction methods

A method for correcting discrete element lifting surface theory to reflect given experimental data is presented. Theoretical pressures are modified such that imposed constraints are satisfied while minimizing the changes to the pressures. Several types of correction procedures are presented and correlated; (1) scaling of pressures; (2) scaling of downwash values; and (3) addition of an increment to the downwash that is proportioned to pressure. Some special features are included in these methods and they include: (1) consideration of experimental data from multiple deflection modes, (2) limitation of the amplitudes of the correction factors, and (3) the use of correction factor mode shapes. These methods are correlated for cases involving all three Mach Number ranges using a FORTRAN IV computer program. Subsonically, a wing with an oscillating partial span control surface and a wing with a leading edge droop are presented. Transonically a two-dimensional airfoil with an oscillating flap is considered. Supersonically an arrow wing with and without camber is analyzed. In addition to correction factor methods an investigation is presented dealing with a new simplified transonic modification of the two-dimensional subsonic lifting surface theory. Correlations are presented for an airfoil with an oscillating flap.

Giesing, J. P.↗

HSCT Ref-H Transonic Flap Data Base: Wind-Tunnel Test and Comparison with Theory

In cooperation with personnel from the Boeing ANP Laboratory and NASA Langley, a performance test was conducted using the Reference-H 1.675% model ("NASA Modular Model") without nacelles at the NASA Langley 16-Ft Transonic Tunnel. The main objective of the test was to determine the drag reduction achievable with leading-edge and trailing-edge flaps deflected along the outboard wing span at transonic Mach numbers (M = 0.9 to 1.2) for purpose of preliminary design and for comparison with computational predictions. The obtained drag data with flap deflections for Mach numbers of 1.07 to 1.20 are unique for the Reference H wing. Four leading-edge and two trailing-edge flap deflection angles were tested at a mean-wing chord-Reynolds number of about 5.7 million. An outboard-wing leading-edge flap deflection of 81 provides a 4.5 percent drag reduction at M = 1.2 A = 0.2), and much larger values at lower Mach numbers with larger flap deflections. The present results for the baseline (no flaps deflected) compare reasonably well with previous Boeing and NASA Ref-H tunnel tests, including high-Reynolds number NTF results. Viscous CFD simulations using the OVERFLOW thin-layer N.S. method properly predict the observed trend in drag reduction at M = 1.2 as function of leading-edge flap deflection. Modified linear theory properly predicts the flap effects on drag at subsonic conditions (Aero2S code), and properly predicts the absolute drag for the 40 and 80 leading-edge deflection at M = 1.2 (A389 code).

Vijgen, Paul M.↗

A simple formula for estimating Stark widths of neutral lines

A simple formula for the prediction of Stark widths of neutral lines similar to the semiempirical method of Griem (1968) for ion lines is presented. This formula is a simplification of the quantum-mechanical classical path impact theory and can be used for complicated atoms for which detailed calculations are not readily available, provided that the effective position of the closest interacting level is known. The expression does not require the use of a computer. The formula has been applied to a limited number of neutral lines of interest, and the width obtained is compared with the much more complete calculations of Bennett and Griem (1971). The agreement generally is well within 50% of the published value for the lines investigated. Comparisons with other formulas are also made. In addition, a simple estimate for the ion-broadening parameter is given.

Freudenstein, S. A.↗

Satellite formation in capillary jet breakup

Computations of finite-amplitude, spatially periodic wave growth on a cylindrical jet have been carried out using a boundary integral method. The initial wave growth is in agreement with Rayleigh's linear theory. When followed to completion these waves pinch off large drops separated by smaller satellite drops (spherules) that decrease in size with decreasing wavelength. The computed sizes of both drops and satellites agree with experiment. It is found that satellites will form for all unstable wave numbers. The small satellites that are computed at wave numbers near the critical wave number were not predicted by near-linear analysis but are observed in experimental photographs of jet breakup. Computation of the collapse of elongated satellites shows short waves propagating on their surfaces.

Mansour, Nagi N.↗

Waves produced from a harmonic point source in a supersonic boundary layer

The disturbance wave pattern produced by a harmonic point source in a compressible flat-plate boundary layer is computed using linear stability theory and direct numerical integration approach. Receptivity coefficients are computed for the spectrum of spanwise modes generated at the source. The effect of boundary layer growth on the development of linear waves is determined by using the method of multiple scales. Results are presented for Mach numbers of 0, 2, 4.5, and 7. It is found that disturbances spread in wedge-shaped regions behind the source and the wedge angle decreases with Mach number. The lateral spreading angle for the instability waves turns out to be quite close to the angle found experimentally for turbulence lateral contamination.

Balakumar, P.↗

Irradiance from stars based on blackbody theory

A computer program which calculates a theoretical backbody curve for a star when given only the star's temperature and visual magnitude is presented. The program was written to produce approximate stellar spectra to be used in the design of the Space Telescope Fine Guidance System (FGS). In addition to the blackbody curve, the program calculates the number of photons tht would be received by the Space Telescope per 25 MSEC interval in wavelength increments of 0.2 microns. The equations used for the computations work best for stars that emit primarily in the range of near infrared to ultraviolet.

Jacobs, W. A.↗

Computer modeling of solar wind interaction with Venus and Mars

The fundamental ideas, successes, and shortcomings of fluid models based on MHD theory for the interaction of the solar wind and a planetary magnetic field and/or atmosphere are reviewed, with special attention to applications to Venus and Mars. An account, including new quantitative results, is given of a number of enhancements of a basic gasdynamic convected field approximation to the magnetohydrodynamic model that incorporate features of solar wind flow past a planet not treated previously. These permit solutions to be calculated for 3D, as opposed to axisymmetric, obstacle shapes, the effects of a magnetic barrier region just exterior to the ionopause, the complete MHD equations for the special case in which the interplanetary magnetic field is aligned with the flow, the effects of mass loading by pickup of ions newly created by photoionization of atmospheric neutrals, and for the unsteady effects associated with the passage of an interplanetary shock wave. Particular emphasis is placed on the information furnished by these models that can be compared with observations, and conversely on the observational data that provide critical tests of the models.

Spreiter, John R.↗

Upper-surface-blowing jet wing interaction

A linear, inviscid, subsonic compressible flow theory is formulated for predicting the aerodynamic characteristics of upper-surface-blowing configurations. The effect of the thick jet is represented by a two-vortex-sheet model in order to account for the Mach number nonuniformity. The wing loading with the jet interaction effects is computed by satisfying boundary conditions on the wing and the jet surfaces. The vortex model is discussed in detail.

Lan, C. E.↗

Static and dynamic stability analysis of the space shuttle vehicle-orbiter

The longitudinal static and dynamic stability of a Space Shuttle Vehicle-Orbiter (SSV Orbiter) model is analyzed using the FLEXSTAB computer program. Nonlinear effects are accounted for by application of a correction technique in the FLEXSTAB system; the technique incorporates experimental force and pressure data into the linear aerodynamic theory. A flexible Orbiter model is treated in the static stability analysis for the flight conditions of Mach number 0.9 for rectilinear flight (1 g) and for a pull-up maneuver (2.5 g) at an altitude of 15.24 km. Static stability parameters and structural deformations of the Orbiter are calculated at trim conditions for the dynamic stability analysis, and the characteristics of damping in pitch are investigated for a Mach number range of 0.3 to 1.2. The calculated results for both the static and dynamic stabilities are compared with the available experimental data.

Chyu, W. J.↗

Prediction and control of transition in hypersonic boundary layers

In this paper, the role of compressible linear stability theory in prediction of boundary layer transition at supersonic and hypersonic speeds is investigated. Computations for sharp cones, using the e exp N method with N = 10, show that the first oblique Tollmien-Schlichting mode is responsible for transition at adiabatic wall conditions for freestream Mach numbers up to 7. For cold walls, the two-dimensional second mode dominates the transition process at lower hypersonic Mach numbers due to the well-known destabilizing effect of cooling on the second mode. It is shown that pressure gradient and suction may be used to stabilize this mode. Some results on the real gas effects on hypersonic boundary-layer stability are presented.

Malik, Mujeeb R.↗

Interactive orbital proximity operations planning system instruction and training guide

This guide instructs users in the operation of a Proximity Operations Planning System. This system uses an interactive graphical method for planning fuel-efficient rendezvous trajectories in the multi-spacecraft environment of the space station and allows the operator to compose a multi-burn transfer trajectory between orbit initial chaser and target trajectories. The available task time (window) of the mission is predetermined and the maneuver is subject to various operational constraints, such as departure, arrival, spatial, plume impingement, and en route passage constraints. The maneuvers are described in terms of the relative motion experienced in a space station centered coordinate system. Both in-orbital plane as well as out-of-orbital plane maneuvering is considered. A number of visual optimization aids are used for assisting the operator in reaching fuel-efficient solutions. These optimization aids are based on the Primer Vector theory. The visual feedback of trajectory shapes, operational constraints, and optimization functions, provided by user-transparent and continuously active background computations, allows the operator to make fast, iterative design changes that rapidly converge to fuel-efficient solutions. The planning tool is an example of operator-assisted optimization of nonlinear cost functions.

Grunwald, Arthur J.↗

An Efficient Scheme for Updating Sparse Cholesky Factors

Raghavan had earlier developed the software package DCSPACK which can be used for solving sparse linear systems where the coefficient matrix is symmetric and positive definite (this project was not funded by NASA but by agencies such as NSF). DSCPACK-S is the serial code and DSCPACK-P is a parallel implementation suitable for multiprocessors or networks-of-workstations with message passing using MCI. The main algorithm used is the Cholesky factorization of a sparse symmetric positive positive definite matrix A = LL(T). The code can also compute the factorization A = LDL(T). The complexity of the software arises from several factors relating to the sparsity of the matrix A. A sparse N x N matrix A has typically less that cN nonzeroes where c is a small constant. If the matrix were dense, it would have O(N2) nonzeroes. The most complicated part of such sparse Cholesky factorization relates to fill-in, i.e., zeroes in the original matrix that become nonzeroes in the factor L. An efficient implementation depends to a large extent on complex data structures and on techniques from graph theory to reduce, identify, and manage fill. DSCPACK is based on an efficient multifrontal implementation with fill-managing algorithms and implementation arising from earlier research by Raghavan and others. Sparse Cholesky factorization is typically a four step process: (1) ordering to compute a fill-reducing numbering, (2) symbolic factorization to determine the nonzero structure of L, (3) numeric factorization to compute L, and, (4) triangular solution to solve L(T)x = y and Ly = b. The first two steps are symbolic and are performed using the graph of the matrix. The numeric factorization step is of dominant cost and there are several schemes for improving performance by exploiting the nested and dense structure of groups of columns in the factor. The latter are aimed at better utilization of the cache-memory hierarchy on modem processors to prevent cache-misses and provide execution rates (operations/second) that are close to the peak rates for dense matrix computations. Currently, EPISCOPACY is being used in an application at NASA directed by J. Newman and M. James. We propose the implementation of efficient schemes for updating the LL(T) or LDL(T) factors computed in DSCPACK-S to meet the computational requirements of their project. A brief description is provided in the next section.

Raghavan, Padma↗

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗

Propeller design by numerical optimization

A computer program designed to optimize propeller characteristics was developed by combining two main programs: the first is the optimization program based on the gradient algorithm; the second is based on a propeller blade element theory and uses an aerodynamics subprogram to approximate the lift and drag characteristics of the NACA 16-series airfoil section. To evaluate the propeller program alone (with its aerodynamics subprogram), propeller characteristics were computed and compared to those from wind tunnel investigations conducted on three different NACA propellers. Although the thrust and power coefficients which were computed using the blade element theory were generally higher than the experimental results for two of the three propellers, the corresponding efficiencies showed good agreement for all three propellers. The propeller optimization program was then used to study the NACA 4-(5)(08)-03 propeller at various Mach numbers from 0.175 to 0.60. Improvements in propeller efficiency and thrust were obtained through the use of the propeller optimization program.

Mendoza, J. P.↗

Instability and transition in supersonic boundary layers

Recently obtained flight transition data for a 5-deg half-angle cone at Mach numbers of 1.2, 1.35, 1.6, and 1.92, together with 'quiet' tunnel data at Mach 3.5, are analyzed by linear compressible stability theory. The results obtained indicate that the exponent N ranges in value between 9 and 11, which is in keeping with values previously determined for subsonic flows. Computations show that the transition in both the flight and quiet wind tunnel experiments was due to oblique first-mode disturbances.

Malik, M. R.↗

Intercomparison among plasma wake models for plasmaspheric and ionospheric conditions

The angular distributions of ions in the wake of a body moving through a space plasma computed from three different models are compared in order to investigate wake current depletion ratios under conditions simulating the topside ionosphere and plasmasphere. Results demonstrate the importance of taking into account the thermal flux at low Mach numbers and the angular acceptance of ion detectors in making theory-experiment comparisons. For all models considered, gradients in the angular variations of the fluxes are shown to be steeper near the wake-ambient interface than closer to the maximum rarefaction region.

Samir, U.↗

Qudit Gate Decomposition Dependence for Lattice Gauge Theories

In this work, we investigate the effect of decomposition basis on primitive qudit gates on superconducting radio-frequency cavity-based quantum computers with applications to lattice gauge theory. Three approaches are tested: SNAP & Displacement gates, ECD & single-qubit rotations $R(\theta,\phi)$, and optimal pulse control. For all three decompositions, implementing the necessary sequence of rotations concurrently rather then sequentially can reduce the primitive gate run time. The number of blocks required for the faster ECD &$R_p(\theta)$ is found to scale $\mathcal{O}(d^2)$, while slower SNAP & Displacement set scales at worst $\mathcal{O}(d)$. For qudits with $d<10$, the resulting gate times for the decompositions is similar, but strongly-dependent on experimental design choices. Optimal control can outperforms both decompositions for small $d$ by a factor of 2-12 at the cost of higher classical resources. Lastly, we find that SNAP & Displacement are slightly more robust to a simplified noise model.

Kürkçüoglu, Doga Murat↗