Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Algebraic structures”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Viscous three-dimensional calculations of transonic fan performance

A 3-D flow analysis code was used to compute the design speed operating line of a transonic fan rotor, and the results were compared with experimental data. The code is an explicit finite difference code with an algebraic turbulence model. The transonic fan, called rotor 67, was tested experimentally at NASA-Lewis with conventional aerodynamic probes and with user anemometry and was included as one of the AGARD test cases for the computation of internal flows. The experimental data are described. Maps of total pressure ratio and adiabatic efficiency versus mass flow were computed and are compared with the experimental maps, with good agreement. Detailed comparisons between calculations and experiment are made at two operating points, one near peak efficiency and the other near stall. Blade-to-blade contour plots are used to show the shock structure. Comparisons of circumferentially integrated flow quantities downstream of the rotor show spanwise distributions of several aerodynamic parameters. Calculated Mach number distributions are compared with laser anemometer data within the blade row and the wake to quantify the accuracy of the calculations. Particle traces are used to show the nature of secondary flow.

Chima, Rodrick V.↗
NASA NTRS1997

Sensitivity Analysis and Optimization of Aerodynamic Configurations with Blend Surfaces

A novel (geometrical) parametrization procedure using solutions to a suitably chosen fourth order partial differential equation is used to define a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. The general airplane configuration has wing, fuselage, vertical tail and horizontal tail. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point Form method. A graphic interface software is developed which dynamically changes the surface of the airplane configuration with the change in input design variable. The software is made user friendly and is targeted towards the initial conceptual development of any aerodynamic configurations. Grid sensitivity with respect to surface design parameters and aerodynamic sensitivity coefficients based on potential flow is obtained using an Automatic Differentiation precompiler software tool ADIFOR. Aerodynamic shape optimization of the complete aircraft with twenty four design variables is performed. Unstructured and structured volume grids and Euler solutions are obtained with standard software to demonstrate the feasibility of the new surface definition.

Thomas, A. M.↗
NASA NTRS1999

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗
NASA NTRS1995

Pressure correlations at a fluid/structure interface

The structure of pressure-pressure correlations at the interface of an incompressible steady-state turbulent flow with a rigid boundary was investigated. For the sake of completeness, the absolute value of the correlation between two random varying functions is herein defined as a number greater than or equal to zero and less than or equal to unity which is a measure of that fraction of one of the functions that 'follows' the second function (or vice versa). It was found that the soughtafter correlations can be determined by consideration of the high Re Navier-Stokes equation, but that the complexity of boundary layer turbulence, in particular the inhomogeneity perpendicular to the boundary and the anisotropy due to convective flow gradients, makes the structure of said correlations extremely difficult to assess. One of the earlier researchers in this field described the quantity under present consideration as 'a quantity which is beyond assessment.' Nonetheless, it was found that under some rather simplifying assumptions the determination of the required structure necessitates the formulation of the related structure of second order two-point correlations of turbulent velocity gradients, as well as third order two-point correlations of velocity gradients. The presence of these latter gradients is due to the nonlinearity in the turbulence ('turbulence self-interaction'). Both of these correlations are scaled, although not similarly, by factors dependent upon the magnitude of the convective flow, which can be modeled using a log law approximation. Fourth order correlations, although present, can be ignored, since they constitute 'higher order terms.' In a slightly more complex situation, it was found that convective flow gradients also have to be incorporated. At the moment, no definitive algebraic information peculiar to pressure-pressure correlations is available in the most highly idealized cases.

Trevino, George↗
NASA NTRS1990

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗
DOE OSTI2026

Tracking discontinuities in parameter space

We develop a geometric framework in Feynman-parameter space to determine constraints on the sequential discontinuities of Feynman integrals. Our method is based on tracking the deformation of the integration contour as external kinematics are analytically continued. This procedure imposes powerful constraints on the analytic structure of Feynman integrals, providing crucial inputs for their bootstrap. We demonstrate the usefulness of this framework by applying it to integrals in dimensional regularization, with higher propagator powers, and to examples with non-uniform transcendental weight. The method is illustrated with several one- and two-loop calculations.

Differential and Algebraic Geometry↗
DOE OSTI2024

Holographic tensor networks with bulk gauge symmetries

Abstract Tensor networks are useful toy models for understanding the structure of entanglement in holographic states and reconstruction of bulk operators within the entanglement wedge. They are, however, constrained to only prepare so-called “fixed-area states” with flat entanglement spectra, limiting their utility in understanding general features of holographic entanglement. Here, we overcome this limitation by constructing a variant of random tensor networks that enjoys bulk gauge symmetries. Our model includes a gauge theory on a general graph, whose gauge-invariant states are fed into a random tensor network. We show that the model satisfies the quantum-corrected Ryu-Takayanagi formula with a nontrivial area operator living in the center of a gauge-invariant algebra. We also demonstrate nontrivial,n-dependent contributions to the Rényi entropy and Rényi mutual information from this area operator, a feature shared by general holographic states.

Physics↗
NASA NTRS1993

The harmonic oscillator and nuclear physics

The three-dimensional harmonic oscillator plays a central role in nuclear physics. It provides the underlying structure of the independent-particle shell model and gives rise to the dynamical group structures on which models of nuclear collective motion are based. It is shown that the three-dimensional harmonic oscillator features a rich variety of coherent states, including vibrations of the monopole, dipole, and quadrupole types, and rotations of the rigid flow, vortex flow, and irrotational flow types. Nuclear collective states exhibit all of these flows. It is also shown that the coherent state representations, which have their origins in applications to the dynamical groups of the simple harmonic oscillator, can be extended to vector coherent state representations with a much wider range of applicability. As a result, coherent state theory and vector coherent state theory become powerful tools in the application of algebraic methods in physics.

Rowe, D. J.↗
DOE OSTI2023

Decorated TQFTs and their Hilbert spaces

We discuss topological quantum field theories that compute topological invariants which depend on additional structures (or decorations) on three-manifolds. The q-series invariant $\hat{Z}$(q) proposed by Gukov, Pei, Putrov, and Vafa is an example of such an invariant. We describe how to obtain these decorated invariants by cutting and gluing and make a proposal for Hilbert spaces that are assigned to two-dimensional surfaces in the $\hat{Z}$-TQFT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fibration structure in toric hypersurface Calabi-Yau threefolds

We find through a systematic analysis that all but 29,223 of the 473.8 million 4D reflexive polytopes found by Kreuzer and Skarke have a 2D reflexive subpolytope. Such a subpolytope is generally associated with the presence of an elliptic or genus one fibration in the corresponding birational equivalence class of Calabi-Yau threefolds. This extends the growing body of evidence that most Calabi-Yau threefolds have an elliptically fibered phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

The symmetric turbulent plane wake downstream of a sharp trailing edge

The analysis and modeling of the symmetric turbulent plane wake downstream of a sharp trailing edge is addressed. A compact description of the flow near the trailing edge is formulated using the results of a previous asymptotic analysis. The new description retains the two-layered structure identified in the previous work, and it clarifies the principal dynamics of the flow in the near-wake outer layer, away from the wake centerline. For zero-pressure-gradient flow, the near-wake outer layer is shown to be represented to leading order by the similarity solution that governs the outer region of the surface boundary layer. The leading perturbation in the outer layer due to the developing near-wake inner-layer flow is identified, and this is shown to be asymptotically smaller than undetermined higher-order terms associated with the surface boundary-layer flow. Results of the new near-wake analysis are used to formulate an algebraic eddy viscosity model for wake flow predictions at arbitrary distances from the trailing edge. The model is used in a numerical solution of the boundary layer equations, and computed velocity and Reynolds stress profiles are shown to compare well with experimental data.

Bogucz, E. A.↗
DOE OSTI2021

Generalising G 2 geometry: involutivity, moment maps and moduli

We analyse the geometry of generic Minkowski N = 1, D = 4 flux compactifications in string theory, the default backgrounds for string model building. In M-theory they are the natural string theoretic extensions of G2 holonomy manifolds. In type II theories, they extend the notion of Calabi-Yau geometry and include the class of flux backgrounds based on generalised complex structures first considered by Graña et al. (GMPT). Using E 7(7) × $\mathbb{R}$ + generalised geometry we show that these compactifications are characterised by an SU(7) ⊂ E 7(7) structure defining an involutive subbundle of the generalised tangent space, and with a vanishing moment map, corresponding to the action of the diffeomorphism and gauge symmetries of the theory. The Kähler potential on the space of structures defines a natural extension of Hitchin’s G 2 functional. Using this framework we are able to count, for the first time, the massless scalar moduli of GMPT solutions in terms of generalised geometry cohomology groups. It also provides an intriguing new perspective on the existence of G 2 manifolds, suggesting possible connections to Geometrical Invariant Theory and stability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗
DOE OSTI2020

Performance and Feature Improvements in Parareal-based Power System Dynamic Simulation

In recent years, a novel Parareal-based approach has been developed for fast transient simulations of large power system interconnections. Parareal belongs to the class of Parallel-in-time algorithms for solution of systems of differential-algebraic equations in parallel over an interval of time. The selection of a reasonably fast and accurate coarse solution is crucial to improve the performance of Parareal algorithm. Semi-analytical solution methods are one promising approach to achieve this goal. They have been investigated, and some preliminary results are presented here. In addition, Parareal-based simulator has been expanded to enable co-simulation with OpenDSS, a widely used open-source distribution system simulator. Preserving the parallel nature of the Parareal approach and taking advantage of the parallel capabilities of the latest versions of OpenDSS, each distribution system can be solved in their entirety on different processors in parallel within the main Parareal simulator. This paper also presents the structure of the transmission and distribution co-simulation and some results with different dynamic models of inverter-based resources in the distribution systems.

Park, Byungkwon↗
NASA NTRS1998

Formal Modeling of Multi-Agent Systems using the Pi-Calculus and Epistemic Logic

Multi-agent systems have become important recently in computer science, especially in artificial intelligence (AI). We allow a broad sense of agent, but require at least that an agent has some measure of autonomy and interacts with other agents via some kind of agent communication language. We are concerned in this paper with formal modeling of multi-agent systems, with emphasis on communication. We propose for this purpose to use the pi-calculus, an extension of the process algebra CCS. Although the literature on the pi-calculus refers to agents, the term is used there in the sense of a process in general. It is our contention, however, that viewing agents in the AI sense as agents in the pi-calculus sense affords significant formal insight. One formalism that has been applied to agents in the AI sense is epistemic logic, the logic of knowledge. The success of epistemic logic in computer science in general has come in large part from its ability to handle concepts of knowledge that apply to groups. We maintain that the pi-calculus affords a natural yet rigorous means by which groups that are significant to epistemic logic may be identified, encapsulated, structured into hierarchies, and restructured in a principled way. This paper is organized as follows: Section 2 introduces the pi-calculus; Section 3 takes a scenario from the classical paper on agent-oriented programming [Sh93] and translates it into a very simple subset of the n-calculus; Section 4 then shows how more sophisticated features of the pi-calculus may bc brought into play; Section 5 discusses how the pi-calculus may be used to define groups for epistemic logic; and Section 6 is the conclusion.

Rorie, Toinette↗
NASA NTRS2011

AutoBayes Program Synthesis System System Internals

This lecture combines the theoretical background of schema based program synthesis with the hands-on study of a powerful, open-source program synthesis system (Auto-Bayes). Schema-based program synthesis is a popular approach toward program synthesis. The lecture will provide an introduction into this topic and discuss how this technology can be used to generate customized algorithms. The synthesis of advanced numerical algorithms requires the availability of a powerful symbolic (algebra) system. Its task is to symbolically solve equations, simplify expressions, or to symbolically calculate derivatives (among others) such that the synthesized algorithms become as efficient as possible. We will discuss the use and importance of the symbolic system for synthesis. Any synthesis system is a large and complex piece of code. In this lecture, we will study Autobayes in detail. AutoBayes has been developed at NASA Ames and has been made open source. It takes a compact statistical specification and generates a customized data analysis algorithm (in C/C++) from it. AutoBayes is written in SWI Prolog and many concepts from rewriting, logic, functional, and symbolic programming. We will discuss the system architecture, the schema libary and the extensive support infra-structure. Practical hands-on experiments and exercises will enable the student to get insight into a realistic program synthesis system and provides knowledge to use, modify, and extend Autobayes.

Statistical Algorithms↗
DOE OSTI2021

Decomposition of Feynman integrals by multivariate intersection numbers

We present a detailed description of the recent idea for a direct decomposition of Feynman integrals onto a basis of master integrals by projections, as well as a direct derivation of the differential equations satisfied by the master integrals, employing multivariate intersection numbers. We discuss a recursive algorithm for the computation of multivariate intersection numbers, and provide three different approaches for a direct decomposition of Feynman integrals, which we dub the straight decomposition, the bottom-up decomposition, and the top-down decomposition. These algorithms exploit the unitarity structure of Feynman integrals by computing intersection numbers supported on cuts, in various orders, thus showing the synthesis of the intersection-theory concepts with unitarity-based methods and integrand decomposition. We perform explicit computations to exemplify all of these approaches applied to Feynman integrals, paving a way towards potential applications to generic multi-loop integrals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗