SEARCH · Engineering Papers

Results for “numeric invariants”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Shielding, the bulk chemical potential, and cohesion in alloys

It is shown that the bulk chemical potential in alloys is intimately related to the spatial dependence of the shielding cloud that results when the electronic charge rearranges itself as one atom type is replaced by another at a given site. Such a relationship fixes the relative energy scale between the alloy and its pure constituents, important in determining the stability of alloys. A correct treatment of shielding is thus essential to quantitative calculations of alloy stability. A model calculation of the bulk chemical potential and cohesion of alloys in the tight-binding approximation is presented as a numerical example. In the course of this investigation a general invariant of an integral over the shielding cloud is derived.

Stern, E. A.↗

NASA NTRS1978

Microwave brightness of polar firn as measured by Nimbus 5 and 6 ESMR

The microwave emission from a half-space medium characterized by coordinate dependent scattering and absorbing centers was calculated by numerically solving the radiative transfer equation by the method of invariant imbedding. A Mie scattering phase function and surface polarization was included in the calculation. Also included are the physical temperature profile and the temperature variation of the index of refraction for ice. Using published values of grain size and temperature profile data of polar firn, the brightness temperature was calculated for the 1.55 cm and 0.8 cm wavelengths. For selected regions in Greenland and Antarctica, the results are in reasonable agreement with the observed Nimbus-5 and Nimbus-6 ESMR data.

Chang, A. T. C.↗

NASA NTRS1983

Manipulator interactive design with interconnected flexible elements

This paper describes the development of an analysis tool for the interactive design of control systems for manipulators and similar electro-mechanical systems amenable to representation as structures in a topological chain. The chain consists of a series of elastic bodies subject to small deformations and arbitrary displacements. The bodies are connected by hinges which permit kinematic constraints, control, or relative motion with six degrees of freedom. The equations of motion for the chain configuration are derived via Kane's method, extended for application to interconnected flexible bodies with time-varying boundary conditions. A corresponding set of modal coordinates has been selected. The motion equations are imbedded within a simulation that transforms the vector-dyadic equations into scalar form for numerical integration. The simulation also includes a linear, time-invariant controler specified in transfer function format and a set of sensors and actuators that interface between the structure and controller. The simulation is driven by an interactive set-up program resulting in an easy-to-use analysis tool.

Singh, R. P.↗

DOE OSTI2015

Wigner phase space distribution via classical adiabatic switching

Evaluation of the Wigner phase space density for systems of many degrees of freedom presents an extremely demanding task because of the oscillatory nature of the Fourier-type integral. We propose a simple and efficient, approximate procedure for generating the Wigner distribution that avoids the computational difficulties associated with the Wigner transform. Starting from a suitable zeroth-order Hamiltonian, for which the Wigner density is available (either analytically or numerically), the phase space distribution is propagated in time via classical trajectories, while the perturbation is gradually switched on. According to the classical adiabatic theorem, each trajectory maintains a constant action if the perturbation is switched on infinitely slowly. We show that the adiabatic switching procedure produces the exact Wigner density for harmonic oscillator eigenstates and also for eigenstates of anharmonic Hamiltonians within the Wentzel-Kramers-Brillouin (WKB) approximation. We generalize the approach to finite temperature by introducing a density rescaling factor that depends on the energy of each trajectory. Time-dependent properties are obtained simply by continuing the integration of each trajectory under the full target Hamiltonian. Further, by construction, the generated approximate Wigner distribution is invariant under classical propagation, and thus, thermodynamic properties are strictly preserved. Numerical tests on one-dimensional and dissipative systems indicate that the method produces results in very good agreement with those obtained by full quantum mechanical methods over a wide temperature range. The method is simple and efficient, as it requires no input besides the force fields required for classical trajectory integration, and is ideal for use in quasiclassical trajectory calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1998

Numerical Simulation of the Generation of Axisymmetric Mode Jet Screech Tones

An imperfectly expanded supersonic jet, invariably, radiates both broadband noise and discrete frequency sound called screech tones. Screech tones are known to be generated by a feedback loop driven by the large scale instability waves of the jet flow. Inside the jet plume is a quasi-periodic shock cell structure. The interaction of the instability waves and the shock cell structure, as the former propagates through the latter, is responsible for the generation of the tones. Presently, there are formulas that can predict the tone frequency fairly accurately. However, there is no known way to predict the screech tone intensity. In this work, the screech phenomenon of an axisymmetric jet at low supersonic Mach number is reproduced by numerical simulation. The computed mean velocity profiles and the shock cell pressure distribution of the jet are found to be in good agreement with experimental measurements. The same is true with the simulated screech frequency. Calculated screech tone intensity and directivity at selected jet Mach number are reported in this paper. The present results demonstrate that numerical simulation using computational aeroacoustics methods offers not only a reliable way to determine the screech tone intensity and directivity but also an opportunity to study the physics and detailed mechanisms of the phenomenon by an entirely new approach.

Shen, Hao↗

DOE OSTI2014

Simulation of sharp interface multi-material flows involving an arbitrary number of components through an extended five-equation model

In this paper, we present an anti-diffusive method dedicated to the simulation of interface flows on Cartesian grids involving an arbitrary number m of compressible components. Our work is two-fold: first, we introduce a m-component flow model that generalizes a classic two material five-equation model. In that way, interfaces are localized using color function discontinuities and a pressure equilibrium closure law is used to complete this new model. The resulting model is demonstrated to be hyperbolic under simple assumptions and consistent. Second, we present a discretization strategy for this model relying on a Lagrange–Remap scheme. Here, the projection step involves an anti-dissipative mechanism allowing to prevent numerical diffusion of the material interfaces. The proposed solver is built ensuring consistency and stability properties but also that the sum of the color functions remains equal to one. The resulting scheme is first order accurate and conservative for the mass, momentum, energy and partial masses. Furthermore, the obtained discretization preserves Riemann invariants like pressure and velocity at the interfaces. Finally, validation computations of this numerical method are performed on several tests in one and two dimensions. The accuracy of the method is also compared to results obtained with the upwind Lagrange–Remap scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Adiabatic model of field reversal by fast ions in an axisymmetric open trap

A model of field reversal by fast ions has been developed under the assumption of preservation of fast-ion adiabatic invariants. Analytical solutions obtained in the approximation of a narrow fast-ion layer and numerical solutions to the evolutionary problem are presented. The solutions demonstrate the process of formation of a field reversed configuration with parameters close to those of the planned experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2024

DDFRG3: Double-Differential FRaGmentation model for Neutron production

A Double-Differential FRaGmentation (DDFRG) model for neutron production from nucleus - nucleus collisions is developed for space radiation applications. DDFRG1 was a previous model developed for production of protons and light ions, and DDFRG2 was a model developed for pion production. A new model (DDFRG3) for neutron production is developed in the present work, and is based upon thermal production of neutrons from the projectile, the target and central fireball sources. The Lorentz-invariant double-differential cross sections are calculated in the various source frames, and are Lorentz transformed to the laboratory frame, resulting in a closed-form analytic formula involving no numerical integration, and which can be run very efficiently in radiation transport codes. The Lorentz-invariant double-differential cross section is then integrated over angle to give the single-differential energy spectral distribution. The DDFRG3 neutron model compares very well to an extensive experimental data set.

neutrons↗

DOE OSTI2015

Spontaneous supersymmetry breaking in two dimensional lattice super QCD

We report on a non-perturbative study of two dimensional N=(2,2) super QCD. Our lattice formulation retains a single exact supersymmetry at non-zero lattice spacing, and contains N f fermions in the fundamental representation of a U(N c ) gauge group. The lattice action we employ contains an additional Fayet-Iliopoulos term which is also invariant under the exact lattice supersymmetry. This work constitutes the first numerical study of this theory which serves as a toy model for understanding some of the issues that are expected to arise in four dimensional super QCD. As a result, we present evidence that the exact supersymmetry breaks spontaneously when N f < N c in agreement with theoretical expectations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

thornado-transport: IMEX schemes for two-moment neutrino transport respecting Fermi-Dirac statistics

We develop implicit-explicit (IMEX) schemes for neutrino transport in a background material in the context of a two-moment model that evolves the angular moments of a neutrino phase-space distribution function. Considering the upper and lower bounds that are introduced by Pauli's exclusion principle on the moments, an algebraic moment closure based on Fermi-Dirac statistics and a convex-invariant time integrator both are demanded. A finite-volume/first-order discontinuous Galerkin(DG) method is used to illustrate how an algebraic moment closure based on Fermi-Dirac statistics is needed to satisfy the bounds. Several algebraic closures are compared with these bounds in mind, and the Cernohorsky and Bludman closure, which satisfies the bounds, is chosen for our IMEX schemes. For the convex-invariant time integrator, two IMEX schemes named PD-ARS have been proposed. PD-ARS denotes a convex-invariant IMEX Runge-Kutta scheme that is high-order accurate in the streaming limit, and works well in the diffusion limit. Our two PD-ARS schemes use second-and third-order, explicit, strong-stability-preserving Runge-Kutta methods as their explicit part, respectively, and therefore are second-and third-order accurate in the streaming limit, respectively. The accuracy and convex-invariance of our PD-ARS schemes are demonstrated in the numerical tests with a third-order DG method for spatial discretization and a simple Lax-Friedrichs flux. The method has been implemented in our high-order neutrino-radiation hydrodynamics (thornado) toolkit. We show preliminary results employing tabulated neutrino opacities.

Chu, Ran↗

DOE OSTI2017

Novel phases in strongly coupled four-fermion theories

Here, we study a lattice model comprising four massless reduced staggered fermions in four dimensions coupled through an SU(4)-invariant four-fermion interaction. We present both theoretical arguments and numerical evidence that no bilinear fermion condensates are present for any value of the four-fermi coupling, in contrast to earlier studies of Higgs{Yukawa models with different exact lattice symmetries. At strong coupling we observe the formation of a four-fermion condensate and a mass gap in spite of the absence of bilinear condensates. Unlike those previously studied systems we do not find a ferromagnetic phase separating this strong-coupling phase from the massless weak-coupling phase. Instead we observe long-range correlations in a narrow region of the coupling, still with vanishing bilinear condensates. While our numerical results come from relatively small lattice volumes that call for caution in drawing conclusions, if this novel phase structure is verified by future investigations employing larger volumes it may offer the possibility for new continuum limits for strongly interacting fermions in four dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Reinforcement Learning of Structured Stabilizing Control for Linear Systems With Unknown State Matrix

This paper delves into designing feedback control gains for a continuous-time linear quadratic regulator (LQR) problem that is constrained to certain predefined structure with unknown state matrix. We bring forth the ideas from reinforcement learning (RL) in conjunction with sufficient stability and performance guarantees in order to design these structured gains using the trajectory measurements of states and controls. Here we first formulate a model-based framework using dynamic programming (DP) to embed the structural constraint to the LQR gain computation in the continuous-time setting, and then subsequently, formulate a policy iteration RL algorithm that can alleviate the requirement of known state matrix in conjunction with maintaining the feedback gain structure. The design enables a distributed learning control design which is necessary for many large-scale cyber-physical systems. Theoretical guarantees are provided for stability and convergence of the structured reinforcement learning (SRL) algorithm. We validate our theoretical results with numerical simulations on a multi-agent networked linear time-invariant (LTI) dynamic system.

42 ENGINEERING↗

NASA NTRS1989

Application of stochastic robustness to aircraft control systems

Stochastic robustness, a simple numerical procedure for estimating the stability robustness of linear, time-invariant systems, is applied to a forward-swept-wing aircraft control system. Based on Monte Carlo evaluation of the system's closed-loop eignevalues, this analysis approach introduces the probability of instability as a scalar stability robustness measure. The related stochastic root locus provides insight into robustness characteristics of the closed-loop system. Three Linear Quadratic controllers of decreasing robustness are chosen to demonstrate the use of stochastic robustness to analyze and compare control designs. Examples are presented illustrating the use of stochastic robustness analysis to address the effects of actuator dynamics and unmodeled dynamics on the stability robustness of the forward-swept-wing aircraft.

Stengel, Robert F.↗

DOE OSTI2025

More about the lattice Hamiltonian for Adjoint QCD 2

In our earlier work [1], we introduced a lattice Hamiltonian for Adjoint QCD 2 using staggered Majorana fermions. We found the gauge invariant space of states explicitly for the gauge group SU(2) and used them for numerical calculations of observables, such as the spectrum and the expectation value of the fermion bilinear. In this paper, we carry out a more in-depth study of our lattice model, extending it to any compact and simply-connected gauge group G. We show how to find the gauge invariant space of states and use it to study various observables. We also use the lattice model to calculate the mixed ’t Hooft anomalies of Adjoint QCD 2 for arbitrary G. We show that the matrix elements of the lattice Hamiltonian can be expressed in terms of the Wigner 6j-symbols of G. For G = SU(3), we perform exact diagonalization for lattices of up to six sites and study the low-lying spectrum, the fermion bilinear condensate, and the string tension. We also show how to write the lattice strong coupling expansion for ground state energies and operator expectation values in terms of the Wigner 6j-symbols. For SU(3) we carry this out explicitly and find good agreement with the exact diagonalizations, and for SU(4) we give expansions that can be compared with future numerical studies.

confinement↗

DOE OSTI2023

Spectral localizer for line-gapped non-Hermitian systems

Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

In many recent applications, particularly in the field of atom-centered descriptors for interatomic potentials, tensor products of spherical harmonics have been used to characterize complex atomic environments. When coupled with a radial basis, the atomic cluster expansion (ACE) basis is obtained. However, symmetrization with respect to both rotation and permutation results in an overcomplete set of ACE descriptors with linear dependencies occurring within blocks of functions corresponding to particular generalized Wigner symbols. All practical applications of ACE employ semi-numerical constructions to generate a complete, fully independent basis. While computationally tractable, the resultant basis cannot be expressed analytically, is susceptible to numerical instability, and thus has limited reproducibility. Here we present a procedure for generating explicit analytic expressions for a complete and independent set of ACE descriptors. The procedure uses a coupling scheme that is maximally symmetric w.r.t. permutation of the atoms, exposing the permutational symmetries of the generalized Wigner symbols, and yields a permutation-adapted rotationally and permutationally invariant basis (PA-RPI ACE). Theoretical support for the approach is presented, as well as numerical evidence of completeness and independence. A summary of explicit enumeration of PA-RPI functions up to rank 6 and polynomial degree 32 is provided. The PA-RPI blocks corresponding to particular generalized Wigner symbols may be either larger or smaller than the corresponding blocks in the simpler rotationally invariant basis. Finally, we demonstrate that basis functions of high polynomial degree persist under strong regularization, indicating the importance of not restricting the maximum degree of basis functions in ACE models a priori.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Thermal Hall conductance and a relative topological invariant of gapped two-dimensional systems

We derive a Kubo-like formula for the thermal Hall conductance of a 2d lattice systems which is free from ambiguities associated with the definition of energy magnetization. We use it to define a relative topological invariant of gapped 2d lattice systems at zero temperature. Up to a numerical factor, it can be identified with the difference of chiral central charges for the corresponding edge modes. This establishes the bulk-boundary correspondence for the chiral central charge. Here, we also show that for any local commuting projector Hamiltonian, the relative chiral central charge vanishes, while for free fermionic systems, it is related to the zero-temperature electric Hall conductance via the Wiedemann-Franz law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Shadow of complex fixed point: Approximate conformality of Q > 4 Potts model

We study the famous example of weakly first-order phase transitions in the (1 + 1)-dimensional quantum Q-state Potts model with Q > 4. We numerically show that these weakly first-order transitions have approximately conformal invariance. Specifically, we find that entanglement entropy on considerably large system sizes fits perfectly with the universal scaling law of this quantity in conformal field theories (CFTs). This supports that the weakly first-order transitions are proximate to complex fixed points, which are described by recent conjectured complex CFTs. Moreover, the central charge extracted from this fitting is close to the real part of the complex central charge of these complex CFTs. Furthermore, we also study the conformal towers and the drifting behaviors of these conformal data (e.g., central charge and scaling dimensions).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗