Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “positive definite functions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Decoherence of multimode thermal squeezed coherent states

It is well known that any multimode positive definite quadratic Hamiltonian can be transformed into a Hamiltonian of uncoupled harmonic oscillators. Based on this theorem, the multimode thermal squeezed coherent states are constructed in terms of density operators. Decoherence of multimode thermal squeezed coherent states is investigated via the characteristic function and it is shown that the decohered (reduced) states are still thermal squeezed coherent states in general.

Yeh, Leehwa↗
NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗
DOE OSTI2022

Deconvolving the components of the sign problem

Auxiliary field quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field h introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices $M_σ(h)$ where σ labels the species, and hence is typically not positive definite. Indeed, the average sign $\langle \mathscr {L} \rangle$ of the determinants goes to zero exponentially with increasing spatial size and decreasing temperature for most Hamiltonians of interest. This statement, however, does not explicitly separate two possible origins for the vanishing of $\langle \mathscr {L} \rangle$. Does $\langle \mathscr {L} \rangle$ → 0 because randomly chosen field configurations have det[M(h)] < 0, or does the sign problem arise because the specific subset of configurations chosen by the weighting function have a greater preponderance of negative values? In the latter case, the process of weighting the configurations with |det[M(h)]| might steer the simulation to a region of configuration space of h where positive and negative determinants are equally likely, even though randomly chosen h would preferentially have determinants with a single dominant sign. Here in this paper, we address the relative importance of these two mechanisms for the vanishing of $\langle \mathscr {L} \rangle$ in quantum simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stringy canonical forms

Canonical forms of positive geometries play an important role in revealing hidden structures of scattering amplitudes, from amplituhedra to associahedra. In this paper, we introduce “stringy canonical forms”, which provide a natural definition and extension of canonical forms for general polytopes, deformed by a parameter α'. They are defined by real or complex integrals regulated with polynomials with exponents, and are meromorphic functions of the exponents, sharing various properties of string amplitudes. As α'→ 0, they reduce to the usual canonical form of a polytope given by the Minkowski sum of the Newton polytopes of the regulating polynomials, or equivalently the volume of the dual of this polytope, naturally determined by tropical functions. At finite α', they have simple poles corresponding to the facets of the polytope, with the residue on the pole given by the stringy canonical form of the facet. There is the remarkable connection between the α'→ 0 limit of tree-level string amplitudes, and scattering equations that appear when studying the α'→ ∞ limit. We show that there is a simple conceptual understanding of this phenomenon for any stringy canonical form: the saddle-point equations provide a diffeomorphism from the integration domain to the interior of the polytope, and thus the canonical form can be obtained as a pushforward via summing over saddle points. When the stringy canonical form is applied to the ABHY associahedron in kinematic space, it produces the usual Koba-Nielsen string integral, giving a direct path from particle to string amplitudes without an a priori reference to the string worldsheet. We also discuss a number of other examples, including stringy canonical forms for finite-type cluster algebras (with type A corresponding to usual string amplitudes), and other natural integrals over the positive Grassmannian.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Common 5S rRNA variants are likely to be accepted in many sequence contexts

Over evolutionary time RNA sequences which are successfully fixed in a population are selected from among those that satisfy the structural and chemical requirements imposed by the function of the RNA. These sequences together comprise the structure space of the RNA. In principle, a comprehensive understanding of RNA structure and function would make it possible to enumerate which specific RNA sequences belong to a particular structure space and which do not. We are using bacterial 5S rRNA as a model system to attempt to identify principles that can be used to predict which sequences do or do not belong to the 5S rRNA structure space. One promising idea is the very intuitive notion that frequently seen sequence changes in an aligned data set of naturally occurring 5S rRNAs would be widely accepted in many other 5S rRNA sequence contexts. To test this hypothesis, we first developed well-defined operational definitions for a Vibrio region of the 5S rRNA structure space and what is meant by a highly variable position. Fourteen sequence variants (10 point changes and 4 base-pair changes) were identified in this way, which, by the hypothesis, would be expected to incorporate successfully in any of the known sequences in the Vibrio region. All 14 of these changes were constructed and separately introduced into the Vibrio proteolyticus 5S rRNA sequence where they are not normally found. Each variant was evaluated for its ability to function as a valid 5S rRNA in an E. coli cellular context. It was found that 93% (13/14) of the variants tested are likely valid 5S rRNAs in this context. In addition, seven variants were constructed that, although present in the Vibrio region, did not meet the stringent criteria for a highly variable position. In this case, 86% (6/7) are likely valid. As a control we also examined seven variants that are seldom or never seen in the Vibrio region of 5S rRNA sequence space. In this case only two of seven were found to be potentially valid. The results demonstrate that changes that occur multiple times in a local region of RNA sequence space in fact usually will be accepted in any sequence context in that same local region.

NASA Discipline Exobiology↗
NASA NTRS1999

Flight Test of an Adaptive Configuration Optimization System for Transport Aircraft

A NASA Dryden Flight Research Center program explores the practical application of real-time adaptive configuration optimization for enhanced transport performance on an L-1011 aircraft. This approach is based on calculation of incremental drag from forced-response, symmetric, outboard aileron maneuvers. In real-time operation, the symmetric outboard aileron deflection is directly optimized, and the horizontal stabilator and angle of attack are indirectly optimized. A flight experiment has been conducted from an onboard research engineering test station, and flight research results are presented herein. The optimization system has demonstrated the capability of determining the minimum drag configuration of the aircraft in real time. The drag-minimization algorithm is capable of identifying drag to approximately a one-drag-count level. Optimizing the symmetric outboard aileron position realizes a drag reduction of 2-3 drag counts (approximately 1 percent). Algorithm analysis of maneuvers indicate that two-sided raised-cosine maneuvers improve definition of the symmetric outboard aileron drag effect, thereby improving analysis results and consistency. Ramp maneuvers provide a more even distribution of data collection as a function of excitation deflection than raised-cosine maneuvers provide. A commercial operational system would require airdata calculations and normal output of current inertial navigation systems; engine pressure ratio measurements would be optional.

Gilyard, Glenn B.↗
NASA NTRS2008

Space Adaptation Back Pain: A Retrospective Study

Back pain is frequently reported by astronauts during the early phase of space flight as they adapt to the microgravity environment. However, the epidemiology of space adaptation back pain has not been well defined. The purpose of this retrospective study was to develop a case definition of space adaptation back pain, determine the incidence of space adaptation back pain, and determine the effectiveness of available treatments. Medical records from the Mercury, Apollo, Apollo-Soyuz Test Project (ASTP), Skylab, Mir, International Space Station (ISS), and Shuttle programs were reviewed. All episodes of in-flight back pain that met the criteria for space adaptation back pain were recorded. Pain characteristics, including intensity, location, and duration of the pain were noted. The effectiveness of specific treatments also was recorded. The incidence of space adaptation back pain among astronauts was determined to be 53% (384/722). Most of the affected astronauts reported mild pain (85%). Moderate pain was reported by 11% of the affected astronauts and severe pain was reported by only 4% of the affected astronauts. The most effective treatments were fetal positioning (91% effective) and the use of analgesic medications (85% effective). This retrospective study aids in the development of a case definition of space adaptation back pain and examines the epidemiology of space adaptation back pain. Space adaptation back pain is usually mild and self-limited. However, there is a risk of functional impairment and mission impact in cases of moderate or severe pain that do not respond to currently available treatments. Therefore, the development of preventive measures and more effective treatments should be pursued.

Kerstman, E. L.↗
NASA NTRS1973

Free vibration analysis of spinning structural systems.

This article presents an efficient digital computer procedure, along with the complete listing of the associated computer program, which may be conveniently utilized for the accurate solution of a wide range of practical eigenvalue problems. Important applications of the present work are envisaged in the natural frequency analysis of spinning structures discretized by the finite element technique, and in the determination of transfer functions associated with the dynamic blocks of control systems of spacecraft utilizing gas jets or reaction wheels for attitude control, as well as of spin-stabilized and dual-spin-stabilized satellites. The validity of the Sturm sequence property is first established for the related matrix formulation involving Hermitian and real symmetric, positive-definite matrices, both being usually of highly banded configuration. A numerically stable algorithm based on the Sturm sequence method is then developed which fully exploits the banded form of the associated matrices.

Gupta, K. K.↗
NASA NTRS1994

LQG optimal compensator transfer function for the NASA LaRC CSI Evolutionary Model

Following the general form for LQG optimal compensators for flexible structures with collocated rate sensors we develop an explicit compensator transfer function for the NASA LaRC CSI Evolutionary model in the form: psi(i omega) = g i omega B(sub u)(sup *)(-M(sub b)omega(exp 2) + T(i omega) + i gamma omega B(sub u)Bu(sub u)(sup *))(exp -1)B(sub u) where T(i omega) is a 48 x 48 positive definite matrix whose derivation is the main result of this report. The undamped mode frequencies can be expressed in terms of T(i omega) as the zeros of Det (-omega(exp 2)M(sub b) + T(i omega)) while 'clamped-clamped' modes of the structure (with all nodes clamped) are the poles.

Balakrishnan, A. V.↗
DOE OSTI2026

Framework for X-ray mirror surface shape fitting

For accurate characterization of grazing-incidence X-ray mirrors, we present a comprehensive framework to fit measured surface shapes (either slope or height) of X-ray mirrors used in synchrotron radiation and free-electron laser facilities. We summarize the closed-form expressions of some typical surface shapes of X-ray mirrors including elliptic cylinders, hyperbolic cylinders, ellipsoids, hyperboloids, and diaboloids. This framework is composed of four layers: definition of standard shapes with closed-form expressions, generation of theoretical surface with pose parameters (six degrees of freedom defining an object's position and orientation relative to a coordinate system), parameter optimization with the ability to select which parameters are fit and which are held constant, and the development of user-friendly fitting function wrappers for particular fitting tasks. A few practical fitting examples are demonstrated to verify the effectiveness of the proposed fitting framework. We discuss the physical meanings of the fitting parameters, and provide several examples using the elliptic cylinder and ellipsoid shapes to highlight some features of the framework. Moreover, we provide the presented framework as open-source codes (MATLAB and Python codes available at https://github.com/nsls2omf/xmf) to the community to encourage academic collaboration and further improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Cooperative Design of Systems of Systems Against Attack on One Subsystem

Future large-scale control systems will be cyberphysical System of Systems (SoS)s each consisting of many interacting subsystems operating harmoniously and making decisions that benefit the entire enterprise. In terms of accomplishing their own objectives, the subsystems in an SoS will be designed so that their decision-making process achieves a Nash equilibrium for the entire SoS. This equilibrium ensures that no one subsystem can improve its position by deviating from the agreed decision. One of the weaknesses of this equilibrium is that it is vulnerable to cyber-attacks directed towards the objectives of one subsystem. In this paper, we consider an attack on a linear quadratic SoS that consists of altering some parameters in the objective function of one subsystem by injecting negative definite matrices into the attacked subsystem's objective function. Such an attack would be designed to disrupt the operation of the entire SoS by pitting one subsystem against the others and causing a change in the behavior and trajectory of the entire SoS. We consider two possible outcomes of the attack. The first is when the attack is not detected and all subsystems continue to operate as designed but now reaching a different Nash equilibrium that may or may not yield the outcome desired by the attacker. In this case, existence results related to the new Nash equilibrium will be derived. The second is when the attack is detected causing the remaining subsystems to react by forming a cooperative team to counter the actions of the attacked subsystem. The team will then establish a new Nash equilibrium with the attacked subsystem. Both of these possible outcomes of the cyber-attack will be analyzed. A 3-subsystems SoS example illustrating the results will also be presented. This example illustrates the ability of an attacker to drive the SoS's state trajectory away from the origin if the attack on one subsystem is not detected. It also illustrates the advantage that the remaining subsystems have in restoring the trajectory to one that is close to the pre-attack equilibrium trajectory by cooperating as a team if the attack is detected.

Talebi, Shahriar↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

The terminal area automated path generation problem

The automated terminal area path generation problem in the advanced Air Traffic Control System (ATC), has been studied. Definitions, input, output and the interrelationships with other ATC functions have been discussed. Alternatives in modeling the problem have been identified. Problem formulations and solution techniques are presented. In particular, the solution of a minimum effort path stretching problem (path generation on a given schedule) has been carried out using the Newton-Raphson trajectory optimization method. Discussions are presented on the effect of different delivery time, aircraft entry position, initial guess on the boundary conditions, etc. Recommendations are made on real-world implementations.

Hsin, C.-C.↗
DOE OSTI2020

Nucleon localization function in rotating nuclei

An electron localization function was originally introduced to visualize in positional space bond structures in molecules. It became a useful tool to describe electron configurations in atoms, molecules, and solids. In nuclear physics, a nucleon localization function (NLF) has been used to characterize cluster structures in light nuclei, formation of fragments in fission, and pasta phases appearing in the inner crust of neutron stars. In the case of a deformed rotating nucleus, several NLFs can be introduced, depending on the definition of the spin-quantization axis, direction of the total angular momentum, and self-consistent symmetries of the system. Contributions to the NLF from the current density, spin-current tensor density, and density gradient terms are negligible in the nuclear interior. The oscillating pattern of the simplified NLF can be explained in terms of a constructive interference between kinetic-energy and particle densities. The characteristic nodal pattern seen in the NLF in the direction of major axis of a rotating nucleus comes from single-particle orbits carrying large aligned angular momentum. The variation of the NLF along the minor axis of the nucleus can be traced back to deformation-aligned orbits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗