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 109 records · Page 6

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Watson, W. R.↗

Recovery of sounding rocket payloads by center-of-gravity position control.

A computerized procedure has been developed to predict dynamic pressure at parachute deployment altitude. Inputs to the program are payload geometry, center-of-gravity position, velocity, and altitude at start of reentry. Payload is assumed to be trimmed, and potential and cross-flow drag theory is used to compute trim angle-of-attack and trim drag coefficient as a function of free-stream Mach number. The computations involve only simple geometric quantities and closed-form equations. Payload flight path is assumed to be vertical, and equations of motion of the payload are integrated from start of reentry to chute deployment altitude.

Mcgarvey, J. F.↗

Theory and computation of optimal low- and medium-thrust transfers

This report presents two numerical methods considered for the computation of fuel-optimal, low-thrust orbit transfers in large numbers of burns. The origins of these methods are observations made with the extremal solutions of transfers in small numbers of burns; there seems to exist a trend such that the longer the time allowed to perform an optimal transfer the less fuel that is used. These longer transfers are obviously of interest since they require a motor of low thrust; however, we also find a trend that the longer the time allowed to perform the optimal transfer the more burns are required to satisfy optimality. Unfortunately, this usually increases the difficulty of computation. Both of the methods described use small-numbered burn solutions to determine solutions in large numbers of burns. One method is a homotopy method that corrects for problems that arise when a solution requires a new burn or coast arc for optimality. The other method is to simply patch together long transfers from smaller ones. An orbit correction problem is solved to develop this method. This method may also lead to a good guidance law for transfer orbits with long transfer times.

Chuang, C.-H.↗

Generalized Hall conductivities in local commuting projector models: Generalized symmetries and protected surface modes

Hall conductivities are important characterizations of phases of matter. It is known that nonzero Hall conductivities are difficult to realize in local commuting projector lattice models due to no-go theorems in (2+1)⁢D. In this work we construct local commuting projector models in (2+1)⁢D and (3+1)⁢D with nonzero generalized Hall conductivities for ordinary and higher-form continuous symmetries on tensor product Hilbert space of finite local dimension. The model is given by a standard ℤ 𝑁 toric code, but the symmetries do not admit expression in terms of on-site charge operators. The symmetry do not have local charges or currents on the lattice in the absence of boundaries, but there is still a notion of Hall conductivities that coincide with the continuum field theories. We construct protected gapless boundaries of the lattice models using modified Villain formalism. The generalized Hall conductivities are computed by surface currents as well as bulk flux insertion and many-body Chern number.

Anomalies↗

Flame-Vortex Interactions Imaged in Microgravity - To Assess the Theory Flame Stretch

The goals of this research are to: 1) Assess the Theory of Flame Stretch by operating a unique flame-vortex experiment under microgravity conditions in the NASA Glenn 2.2 Second Drop Tower (drops to identify operating conditions have been completed); 2) Obtain high speed shadowgraph images (500-1000 frames/s) using the drop rig (images were obtained at one-g, and the NASA Kodak RO camera is being mounted on the drop rig); 3) Obtain shadowgraph and PIV images at 1-g while varying the effects of buoyancy by controlling the Froude number (completed); 4) Numerically model the inwardly-propagating spherical flame that is observed in the experiment using full chemistry and the RUN 1DL code (completed); 5) Send images of the flame shape to Dr. G. Patniak at NRL who is numerically simulating the entire flame-vortex interaction of the present experiment (data transfer completed); and 6) Assess the feasibility of obtaining PIV velocity field images in the drop rig, which would be useful (but not required) for our assessment of the Theory of Flame Stretch (PIV images were obtained at one-g using same low laser power that is available from fiber optic cable in drop tower). The motivation for the work is to obtain novel measurement needed to develop a physically accurate model of turbulent combustion that can help in the control of engine pollutants. The unique experiment allows, for the first time, the detailed study of a negatively-curved (negatively stretched) flame, which is one of the five fundamental types of premixed flames. While there have been studies of flat flames, positively-curved (outwardly-propagating) cases and positively-strained (counterflow) cases, this is the first detailed study of a negatively-curved (inwardly-propagating) flame. The first set of drops in the 2.2 Second Drop Tower showed that microgravity provides more favorable conditions for achieving inwardly-propagating flames (IPFs) than 1-g. A vortex interacts with a flame and creates a spherical pocket, which burns inwardly. Shadowgraphs at 1000 frames/sec quantify the Markstein number and flame speed. A Low-Laser Power PIV System was developed and is being added to the drop package. Numerical computations were required to explain why the Markstein numbers measured for the inwardly-propagating flames differ from those of outward propagating flames; this is an important research issue in the assessment of the Theory of Flame Stretch. The RUN-1DL code (developed by Prof. B. Rogg) was run for IPF and OPFs with complex methane and propane chemistry. Results confirmed that Ma for the IPFs are larger than for OPFs as was observed experimentally. Physical reasons for these new findings about the Theory of Flame Stretch are being determined from the experiments and the computations. Several journal papers have been published; the drop package is described in the AIAA Journal, while the one-g results appear in three other journal papers.

Driscoll, James F.↗

Transonic aerodynamics - History and statement of the problem

An historical account is given of the state of transonic aerodynamics research and theorization before 1940, and at selected points between that year and 1975, in order to provide a basis for the understanding and evaluation of the most recent developments in the field. Intensive development began during the Second World War, when fighter aircraft diving at speeds in excess of Mach 0.8 underwent buffeting and longitudinal trim shift effects. By 1960, transonic small disturbance theory had been established as a mathematical model for the prediction of the aerodynamic properties of airfoils, wings, and bodies throughout the Mach number range from zero to 1.5. In the years leading up to 1975, which saw the formulation of the transonic equivalence rule, improvements in the computational speed and memory capacity of computers became the most significant factor in the development of theory.

Spreiter, J. R.↗

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TLD measurements of gamma heating in heavy elements.

Measurements and calculations of gamma heating in polyethylene and lead containers were done and compared. The objective was to provide a workable method of getting good values for gamma heating in in-pile experiments containing materials of high atomic numbers. It was inferred that a combination of thermoluminescent dosimeter measurements, using Bragg-Gray theory, with photon transport calculations using the ANISN computer program, would meet this objective.

Reilly, H. J.↗

Prediction of static aerodynamic characteristics for slender bodies alone and with lifting surfaces to very high angles of attack

An engineering-type method is presented for computing normal-force and pitching-moment coefficients for slender bodies of circular and noncircular cross section alone and with lifting surfaces. In this method, a semi-empirical term representing viscous-separation crossflow is added to a term representing potential-theory crossflow. For many bodies of revolution, computed aerodynamic characteristics are shown to agree with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. The angles of attack extend from 0 deg to 180 deg for M = 2.9 from 0 deg to 60 deg for M = 0.6 to 2.0. For several bodies of elliptic cross section, measured results are also predicted reasonably well over the investigated Mach number range from 0.6 to 2.0 and at angles of attack from 0 deg to 60 deg. As for the bodies of revolution, the predictions are best for supersonic Mach numbers. For body-wing and body-wing-tail configurations with wings of aspect ratios 3 and 4, measured normal-force coefficients and centers are predicted reasonably well at the upper test Mach number of 2.0. Vapor-screen and oil-flow pictures are shown for many body, body-wing and body-wing-tail configurations. When spearation and vortex patterns are asymmetric, undesirable side forces are measured for the models even at zero sideslip angle. Generally, the side-force coefficients decrease or vanish with the following: increase in Mach number, decrease in nose fineness ratio, change from sharp to blunt nose, and flattening of body cross section (particularly the body nose).

Jorgensen, L. H.↗

Effect of angle of attack and Mach number on slender wing aerodynamics

An analytic theory is presented in which the classical slender wing theory is modified to account for the combined effects of large angle of attack and nonsonic Mach number on the unsteady aerodynamics. The computed results agree well with available static and dynamic experimental data for slender delta wings in the freestream Mach number range between 0 and 2.8. The method was extended to compute the unsteady aerodynamics of the space shuttle orbiter by defining an equivalent slender wing using static experimental data. The results obtained in this manner are in good agreement with dynamic experimental results for the freestream Mach number range between 0.3 and 1.2.

Ericsson, L. E.↗

Navy propeller section characteristics as used in propeller design

This report contains artificial aerodynamic characteristics of a set of propeller sections to be used in designing propellers by means of the blade element theory. Characteristics computed from model propeller tests for a single section are extended to cover sections of Navy propeller sections at high Reynolds Number in the variable density tunnel of the NACA.

Weick, Fred E↗

Design of centrifugal impeller blades

This paper restricts itself to radial impellers with cylindrical blades since, as Prasil has shown, the flow about an arbitrarily curved surface of revolution may be reduced to this normal form we have chosen by a relatively simple conformal transformation. This method starts from the simple hypotheses of the older centrifugal impeller theory by first assuming an impeller with an infinite number of blades. How the flow is then modified is then investigated. For the computation of flow for a finite number of blades, the approximation method as developed by Munk, Prandtl and Birnbaum, or Glauert is found suitable. The essential idea of this method is to replace the wing by a vortex sheet and compute the flow as the field of these vortices. The shape of the blades is then obtained from the condition that the flow must be along the surface of the blade.

Betz, A↗

Thermal analysis of a functionally graded material subject to a thermal gradient using the boundary element method

The boundary element method is utilized in this study to conduct thermal analysis of functionally graded composites, materials in which the internal microstructure or properties are explicitly tailored in order to obtain an optimal response, on the micromechanical (constituent) scale. A unique feature of the boundary element formulations used here is the use of circular shape functions to convert the two-dimensional integrations of the composite fibers to one dimensional integrations. Using the computer code BEST-CMS, the through the thickness temperature profiles are computed for a representative material with varying numbers of fibers and fiber spacing in the thickness direction. The computed temperature profiles are compared to those obtained using an alternate analytical theory which explicitly couples the heterogeneous microstructure to the global analysis. The boundary element results compared favorably to the analytical calculations, with discrepancies that are explainable based on the boundary element formulation. The results serve both to demonstrate the ability of the boundary element method to analyze these types of materials, and to verify the accuracy of the analytical theory.

Goldberg, Robert K.↗

Corner Wrinkling at a Square Membrane Due to Symmetric Mechanical Loads

Thin-film membrane structures are under consideration for use in many future gossamer spacecraft systems. Examples include sunshields for large aperture telescopes, solar sails, and membrane optics. The development of capabilities for testing and analyzing pre-tensioned, thin film membrane structures is an important and challenging aspect of gossamer spacecraft technology development. This paper presents results from experimental and computational studies performed to characterize the wrinkling behavior of thin-fi[m membranes under mechanical loading. The test article is a 500 mm square membrane subjected to symmetric comer loads. Data is presented for loads ranging from 0.49 N to 4.91 N. The experimental results show that as the load increases the number of wrinkles increases, while the wrinkle amplitude decreases. The computational model uses a finite element implementation of Stein-Hedgepeth membrane wrinkling theory to predict the behavior of the membrane. Comparisons were made with experimental results for the wrinkle angle and wrinkled region. There was reasonably good agreement between the measured wrinkle angle and the predicted directions of the major principle stresses. The shape of the wrinkle region predicted by the finite element model matches that observed in the experiments; however, the size of the predicted region is smaller than that determined in the experiments.

Blandino, Joseph R.↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

Knowledge Driven Image Mining with Mixture Density Mercer Kernels

This paper presents a new methodology for automatic knowledge driven image mining based on the theory of Mercer Kernels; which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly infinite dimensional feature space. In that high dimensional feature space, linear clustering, prediction, and classification algorithms can be applied and the results can be mapped back down to the original image space. Thus, highly nonlinear structure in the image can be recovered through the use of well-known linear mathematics in the feature space. This process has a number of advantages over traditional methods in that it allows for nonlinear interactions to be modelled with only a marginal increase in computational costs. In this paper, we present the theory of Mercer Kernels, describe its use in image mining, discuss a new method to generate Mercer Kernels directly from data, and compare the results with existing algorithms on data from the MODIS (Moderate Resolution Spectral Radiometer) instrument taken over the Arctic region. We also discuss the potential application of these methods on the Intelligent Archive, a NASA initiative for developing a tagged image data warehouse for the Earth Sciences.

Srivastava, Ashok N.↗

Knowledge Driven Image Mining with Mixture Density Mercer Kernals

This paper presents a new methodology for automatic knowledge driven image mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly infinite dimensional feature space. In that high dimensional feature space, linear clustering, prediction, and classification algorithms can be applied and the results can be mapped back down to the original image space. Thus, highly nonlinear structure in the image can be recovered through the use of well-known linear mathematics in the feature space. This process has a number of advantages over traditional methods in that it allows for nonlinear interactions to be modelled with only a marginal increase in computational costs. In this paper we present the theory of Mercer Kernels; describe its use in image mining, discuss a new method to generate Mercer Kernels directly from data, and compare the results with existing algorithms on data from the MODIS (Moderate Resolution Spectral Radiometer) instrument taken over the Arctic region. We also discuss the potential application of these methods on the Intelligent Archive, a NASA initiative for developing a tagged image data warehouse for the Earth Sciences.

Srivastava, Ashok N.↗

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗