SEARCH · Search NASA
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.
Using the well-known “algebra of multifractality,” we derive the functional equation for anomalous dimensions Δ{sub q}, whose solution Δ = χq(q–1) corresponds to strict parabolicity of the multifractal spectrum. This result demonstrates clearly that a correspondence of the nonlinear σ-models with the initial disordered systems is not exact.
A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.
TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.
The local geometric properties of a nonlinear control system defined by vector fields E(sub 1),..., E(sub M) are determined by the algebraic properties of the iterated Lie brackets of the E(sub j)'s. The problem of rewriting expressions is examined involving the E(sub j)'s in terms of the D(sub mu)'s in such a way as to handle effectively any cancellation occurring due to commuting of the D(sub mu)'s. A data structure is introduced which allows us to organize the computation to take advantage of the symmetries in the expression and reduce the operation count.
The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.
The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.
In this work, we develop our idea on the construction of a system of combinatorial generators in a T-ideal of a free associative algebra, which is a full analogy of a Gröbner–Shirshov basis in a polynomial ideal. We prove a theorem on multilinear monomials that enables us to establish the existence of a finite set of combinatorial generators in a metabelian T-ideal.
The Jordan-Wigner transformation establishes a duality between $su(2)$ and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Here, our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicate that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams, and provides a better starting point for addressing challenging spin problems.
The scalable iterative solution of strongly coupled three-dimensional incompressible resistive magnetohydrodynamics (MHD) equations is quite challenging because disparate time scales arise from the electromagnetics, the hydrodynamics, as well as the coupling between these systems. This study considers a mixed finite element discretization of a dual saddle point formulation of the incompressible resistive MHD equations using a stable nodal (Q2/Q1) discretization for the hydrodynamics and a stable edge-node discretization of a reduced form of the Maxwell equations. This paper introduces new approximate block factorization preconditioners for this system which reduce the system to approximate Schur complement systems that can be solved using algebraic multilevel methods. These preconditioners include a new augmentation-based approximation for the magnetic induction saddle point system as well as efficient approximations of the Schur complements that arise from the complex coupling between the Navier--Stokes equations and the Maxwell equations.
In this report, the performance of the algebraic multigrid solver used as a preconditioner for conjugate gradient is investigated on 2 nodes of Spock, an early access system at ORNL with 4 MI-100 GPUs and a 64-core Rome CPU per node. We compare GPU and CPU performance for three different diffusion problems using increasing problem sizes.
The General Interpolants Method (GIM) code which solves the multidimensional Navier-Stokes equations for arbitrary geometric domains is described. The geometry module in the GIM code generates two and three dimensional grids over specified flow regimes, establishes boundary condition information and computes finite difference analogs for use in the GIM code numerical solution module. The technique can be classified as an algebraic equation approach. The geometry package uses multivariate blending function interpolation of vector-values functions which define the shapes of the edges and surfaces bounding the flow domain. By employing blending functions which conform to the cardinality conditions the flow domain may be mapped onto a unit square (2-D) or unit cube (3-D), thus producing an intrinsic coordinate system for the region of interest. The intrinsic coordinate system facilitates grid spacing control to allow for optimum distribution of nodes in the flow domain.
This paper presents a novel technique to validate and predict the rover slips on Martian surface for NASA's Mars Exploration Rovers. Different from the traditional approach, the proposed method uses the actual velocity profile of the wheels and the digital elevation map (DEM) from the stereo images of the terrain to formulate simplified equations of motion of the rovers. A weighted factor to the wheel-ground speed from the empirical data comprises the velocity equations of the simplified differential-algebraic system of the rover motion. Applying the discretization operator to these equations, the full kinematics state of the rover is then resolved by the configuration kinematics solution in the Robot Sequencing and Visualization Program (RSVP). This method produced accurate simulation of the rover movements compared with these of the earth testing vehicle. Using the telemetry from the onboard Visual Odometry, the simulated rover path also compares well with the actual track of the vehicle. Preliminary results indicated that the proposed computational method is very effective in planning the path of the rovers on the high-slope areas.
This work aims at bringing advances in discrete optimization algorithms to solving practical engineering problems at scale. Often times, in many engineering design problems, there is a need to select a small set of influential or representative elements from a large ground set of entities in an optimal fashion. Submodular optimization provides for a formal way to solve such problems. Common examples with infrastructure systems involve sensor placement and identification of key entities with certain objectives. However, scaling these approaches to large infrastructure systems can be challenging because of the high computational complexity of the overall framework that include the optimization algorithms as well as high-complexity compute-oracles that provide the necessary objective function values. In this work, we explore a well-studied and widely-applicable paradigm, namely leader-selection in a multi-agent networked setting in the context of scalable methodologies. We demonstrate novel frameworks that utilize variations of accelerated submodular optimization algorithms along with linear-algebraic methods that can help accelerate the oracle computations. We further explore this combination in conjunction with graph partitioning paradigms to take advantage of the accelerated algorithms in a distributed setting. Finally we demonstrate the key findings on a practical problem in an operational setting. For this, we leverage an example road network with approximately 18k nodes and 27k edges in a traffic control application, where we seek a limited number of k=200 key intersections. This problem can be solved in a serial setting in just under 5 hours providing more than 2 orders of magnitude speed-up over methods that do not consider acceleration techniques.
The NASA Glenn Research Center program in megawatt level electric propulsion is centered on electromagnetic acceleration of quasi-neutral plasmas. Specific concepts currently being examined are the Magnetoplasmadynamic (MPD) thruster and the Pulsed Inductive Thruster (PIT). In the case of the MPD thruster, a multifaceted approach of experiments, computational modeling, and systems-level models of self field MPD thrusters is underway. The MPD thruster experimental research consists of a 1-10 MWe, 2 ms pulse-forming-network, a vacuum chamber with two 32 diffusion pumps, and voltage, current, mass flow rate, and thrust stand diagnostics. Current focus is on obtaining repeatable thrust measurements of a Princeton Benchmark type self field thruster operating at 0.5-1 gls of argon. Operation with hydrogen is the ultimate goal to realize the increased efficiency anticipated using the lighter gas. Computational modeling is done using the MACH2 MHD code, which can include real gas effects for propellants of interest to MPD operation. The MACH2 code has been benchmarked against other MPD thruster data, and has been used to create a point design for a 3000 second specific impulse (Isp) MPD thruster. This design is awaiting testing in the experimental facility. For the PIT, a computational investigation using MACH2 has been initiated, with experiments awaiting further funding. Although the calculated results have been found to be sensitive to the initial ionization assumptions, recent results have agreed well with experimental data. Finally, a systems level self-field MPD thruster model has been developed that allows for a mission planner or system designer to input Isp and power level into the model equations and obtain values for efficiency, mass flow rate, and input current and voltage. This model emphasizes algebraic simplicity to allow its incorporation into larger trajectory or system optimization codes. The systems level approach will be extended to the pulsed inductive thruster and other electrodeless thrusters at a future date.
We study linear systems of ordinary differential equations of higher order with an identically singular matrix at the higher derivative of the required vector-valued function in the domain of definition. We give definitions of the index and singular points for these systems, formulate the conditions of solvability, and deduce a formula for the general solution. The algorithms for finding the index and singular points are presented.
Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.
The spiking activity of principal cells in mammalian hippocampus encodes an internalized neuronal representation of the ambient space-a cognitive map. Once learned, such a map enables the animal to navigate a given environment for a long period. However, the neuronal substrate that produces this map is transient: the synaptic connections in the hippocampus and in the downstream neuronal networks never cease to form and to deteriorate at a rapid rate. How can the brain maintain a robust, reliable representation of space using a network that constantly changes its architecture? We address this question using a computational framework that allows evaluating the effect produced by the decaying connections between simulated hippocampal neurons on the properties of the cognitive map. Using novel Algebraic Topology techniques, we demonstrate that emergence of stable cognitive maps produced by networks with transient architectures is a generic phenomenon. The model also points out that deterioration of the cognitive map caused by weakening or lost connections between neurons may be compensated by simulating the neuronal activity. Lastly, the model explicates the importance of the complementary learning systems for processing spatial information at different levels of spatiotemporal granularity.
We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data.We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of Oiii and Fe ii and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe ii]. Key words: atomic data - atomic processes - line: formation - methods: data analysis - molecular data - molecular processes - techniques: spectroscopic