Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗

Nuclear Matrix Elements for Neutrinoless Double-Beta Decay

Neutrinoless double-beta decay ($0\nu\beta\beta$) is a rare hypothesised process that, if discovered, would establish that the neutrino is Majorana, that is, it is its own antiparticle. Interpretation of experimental results relies on knowledge of nuclear matrix elements, whose large model uncertainty is the limiting factor in comparing measured (bounds on) half-lives to the neutrino mass. Nuclear effective field theory and lattice QCD have the potential to compute these matrix elements with better control over uncertainties, enhancing the discovery potential of next-generation $0\nu\beta\beta$ experiments. This work will survey various lattice QCD double-beta decay calculations and discuss their implications.

Grebe, Anthony V. [Fermilab] (ORCID:00000003103201↗

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗

Application of transport-based metric for continuous interpolation between cryo-EM density maps

Cryogenic electron microscopy (cryo-EM) has become widely used for the past few years in structural biology, to collect single images of macromolecules "frozen in time". As this technique facilitates the identification of multiple conformational states adopted by the same molecule, a direct product of it is a set of 3D volumes, also called EM maps. To gain more insights on the possible mechanisms that govern transitions between different states, and hence the mode of action of a molecule, we recently introduced a bioinformatic tool that interpolates and generates morphing trajectories joining two given EM maps. This tool is based on recent advances made in optimal transport, that allow efficient evaluation of Wasserstein barycenters of 3D shapes. As the overall performance of the method depends on various key parameters, including the sensitivity of the regularization parameter, we performed various numerical experiments to demonstrate how MorphOT can be applied in different contexts and settings. Finally, we discuss current limitations and further potential connections between other optimal transport theories and the conformational heterogeneity problem inherent with cryo-EM data.

3D shapes↗

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

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

Mcgarvey, J. F.↗

Formation of ion-acoustic shocks

The formation of an ion-acoustic shock in a numerically modeled two-chamber double plasma device is investigated for a plasma of cold ions and isothermal Boltzmann electrons. An initial potential ramp applied to the driver chamber launches an ion-acoustic pulse into the target chamber which steepens into a shock. The quasi-steady shock structure agrees with observed double plasma shocks. An upper limit of Mach 1.6 is observed independent of the potential ramp magnitude, in agreement with theory.

White, R. B.↗

The collisional interaction of a beam of charged particles with a hydrogen target of arbitrary ionization level

The classical theory of scattering under the Coulomb potential of both charged and neutral particles is used to derive formulae for the energy deposition rate and mean scattering of a beam of charged particles interacting with a cold hydrogen target of arbitrary ionization level as a function of the column density traversed by the beam. These general results hold for any form of stable injection energy spectrum, and their relevance to the existing literature on chromospheric heating during solar flares is discussed.

Emslie, A. G.↗

The analysis of convolutional codes via the extended Smith algorithm

Convolutional codes have been the central part of most error-control systems in deep-space communication for many years. Almost all such applications, however, have used the restricted class of (n,1), also known as 'rate 1/n,' convolutional codes. The more general class of (n,k) convolutional codes contains many potentially useful codes, but their algebraic theory is difficult and has proved to be a stumbling block in the evolution of convolutional coding systems. In this article, the situation is improved by describing a set of practical algorithms for computing certain basic things about a convolutional code (among them the degree, the Forney indices, a minimal generator matrix, and a parity-check matrix), which are usually needed before a system using the code can be built. The approach is based on the classic Forney theory for convolutional codes, together with the extended Smith algorithm for polynomial matrices, which is introduced in this article.

Mceliece, R. J.↗

Unsteady aerodynamics and flutter of propfans using a three-dimensional Full-Potential Solver

A full-potential solver coupled with a linear structural dynamics model is used to calculate the unsteady aerodynamics and aeroelasticity of propfans. The solver allows calculations for arbitrary interblade phase angles. Results are presented for two propfan configurations. Good agreement is seen between the full-potential results and results from linear theory since the flow is subsonic and the thickness of the propfan blades is small. Some difficulty is encountered due to wave reflections from outer computational boundaries; however, this does not affect the results in the range of frequencies of aeroelastic interest.

Bakhle, Milind A.↗

Unsteady aerodynamics and flutter based on the potential equation

A time-domain three-dimensional full-potential solver is coupled with a linear structural dynamics model to investigate the unsteady aerodynamics and aeroelasticity of propfans. The solver allows calculations in multiple blade passages with independent blade motions. Aeroelastic calculations are performed in both frequency and time domains. Results are presented for two propfan configurations. Good agreement is seen between the full-potential results and results from linear theory since the flow is subsonic and the thickness of the propfan blades is small; the agreement is not as good for the cases in which the angle of attack is high. Some difficulty is encountered due to wave reflections from outer computational boundaries; however, this does not affect the results in the range of frequencies of interest.

Bakhle, Milind A.↗

Stability derivatives of triangular wings at supersonic speeds

The analysis of the stability derivatives of low-aspect-ratio triangular wings at subsonic and supersonic speeds, given in NACA TN no. 1423, is extended to apply to triangular wings having large vertex angles and traveling at supersonic speeds. The lift, rolling moment due to sideslip, and damping in roll and pitch for this more general case have been treated elsewhere on the basis of the theory of small disturbances. The surface potentials for angle of attack and rolling taken therefrom are used to obtain the several side-force and yawing-moment derivatives that depend on leading-edge suction, and a tentative value for the rolling moment due to yawing. The lift and moment due to downward acceleration are obtained on the basis of an unpublished unsteady-flow solution. All the known stability derivatives of the triangular wing at supersonic speeds, regardless of source, are summarized for convenience and presented with respect to both body axes and stability axes. The results are limited to Mach numbers for which the triangular wing is contained within the Mach cone for its vertex. The spanwise variation of Mach number in the case of yawing is neglected, although the effect must be of importance.

Ribner, Herbert S↗

Stability derivatives at supersonic speeds of thin rectangular wings with diagonals ahead of tip Mach lines

The investigation includes steady and accelerated vertical and longitudinal motions and steady rolling, yawing, sideslipping, and pitching for Mach numbers and aspect ratios greater than those for which the Mach line from the leading edge of the tip section intersects the trailing edge of the opposite tip section. The stability derivatives are derived with respect to principal body axes and then transformed to a system of stability axes. Theoretical results are obtained, by means of the linearized theory, for the surface-velocity-potential functions, surface-pressure distributions, and stability derivatives for various motions at supersonic speeds of thin flat rectangular wings without dihedral. In the case of yawing, a treatment for the infinitely long wing which takes account of the spanwise variation in the stream Mach number is extended to the finite wing, and a plausible, although not rigorous, solution is obtained for the wing tip effects.

Harmon, Sidney M↗

Gravity waves in Titan's atmosphere

Scintillations (high frequency variations) observed in the radio signal during the occultation of Voyager 1 by Titan (Hinson and Tyler, 1983) provide information concerning neutral atmospheric density fluctuations on scales on hundreds of meters to a few kilometers. Those seen at altitudes higher than 25 km above the surface were interpreted by Hinson and Tyler as being caused by linear, freely propagating (energy-conserving) gravity waves, but this interpretation was found to be inconsistent with the scintillation data below the 25-km altitude level. Here an attempt is made to interpret the entire scintillation profile between the surface and the 90-km altitude level in terms of gravity waves generated at the surface. Numerical calculations of the density fluctuations caused by two-dimensional, nonhydrostatic, finite-amplitude gravity waves propagating vertically through Titan's atmosphere are performed to produce synthetic scintillation profiles for comparison with the observations. The numerical model accurately treats the effects of wave transience, nonlinearity, and breakdown due to convective instability in the overturned part of the wave. The high-altitude scintillation data were accurately recovered with a freely propagating wave solution, confirming the analytic model of Hinson and Tyler. It is found that the low-altitude scintillation data can be fit by a model where a component of the gravity waves becomes convectively unstable and breaks near the 15 km level. The large-scale structure of the observed scintillation profile in the entire altitude range between 5 and 85 km can be simulated by a model where the freely propagating and breaking waves are forced at the surface simultaneously. Further analysis of the Voyager 1 Titan low-altitude scintillation data, using inversion theory appropriate for strong scattering, could potentially remove some of the ambiguities remaining in this analysis and allow a better determination of the strength and source of the waves.

Friedson, A. James↗