Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2025

Skyrmionium dynamics and stability on one dimensional anisotropy patterns

We examine a skyrmionium driven over a periodic anisotropy pattern, which consists of disorder free regions and disordered regions. For small defect densities, the skyrmionium flows for an extended range of currents, and there is a critical current above which it transforms into a skyrmion. For increased amounts of quenched disorder, the current needed for the skyrmionium to transform into a skyrmion decreases, and there is a critical disorder density above which a moving skyrmionium is not stable. In the moving state, the skyrmionium to skyrmion transformation leads to a drop in the velocity and the onset of a finite skyrmion Hall angle. We also find a reentrance effect in which the pinned skyrmionium transforms into a skyrmion just above depinning, restabilizes into skyrmionium at larger drives, and becomes unstable again at large currents. In conclusion, we also show that adding a transverse shaking drive can increase the lifetime of a moving skyrmionium by reducing the effect of the pinning in the direction of the drive.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tunable polar distortions and magnetism in Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb epitaxial films

Hexagonal 𝐴⁢𝐵⁢𝐶 intermetallics are predicted to have tunable ferroelectric, topological, and magnetic properties as a function of the polar buckling of 𝐵⁢𝐶 atomic planes. Here, we report the impact of isovalent lanthanide substitution on the buckling, structural phase transitions, and electronic and magnetic properties of Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films grown by molecular beam epitaxy (MBE) on 𝑐 plane sapphire substrates. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films form a solid solution from 𝑥=0 to 𝑥=1 and retain the polar hexagonal structure (𝑃⁢6 3 ⁢𝑚⁢𝑐) out to 𝑥 ≤ 0.95. With increasing 𝑥, the PtSb buckling increases and the out-of-plane lattice constant 𝑐 decreases due to the lanthanide contraction. While hexagonal LaPtSb is a highly conductive polar metal, the carrier density decreases with 𝑥 until an abrupt phase transition to a zero band overlap semimetal is found for cubic GdPtSb at 𝑥=1. The magnetic susceptibility peaks at small but finite 𝑥, which we attribute to Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling between localized 4⁢𝑓 moments, whose concentration increases with 𝑥, and free carriers that decrease with 𝑥. Samples with 𝑥≥0.3 show antiferromagnetic Curie-Weiss behavior and a Neel temperature that increases with 𝑥. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb system provides opportunities to dramatically alter the polar buckling and concentration of local 4⁢𝑓 moments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
NASA NTRS1975

High-beta turbulence in two-dimensional magnetohydrodynamics

Incompressible turbulent flows were investigated in the framework of ideal magnetohydrodynamics. Equilibrium canonical distributions are determined in a phase whose coordinates are the real and imaginary parts of the Fourier coefficients for the field variables. The magnetic field and fluid velocity have variable x and y components, and all field quantities are independent of z. Three constants of the motion are found which survive the truncation in Fourier space and permit the construction of canonical distributions with three independent temperatures. Spectral densities are calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long wavelength, self-generated, magnetic fields ("magnetic islands"). In the presence of finite dissipation, energy cascades to higher wave numbers can be accompanied by vector potential cascades to lower wave numbers, in much the same way that in the fluid dynamic case, energy cascades to lower wave numbers accompany entropy cascades to higher wave numbers.

Fyfe, D.↗
NASA NTRS1979

Numerical solution of transonic full stream function equations in conservation form

The stream function equation in conservation form is solved iteratively based on the artificial compressibility method. The density is not a unique function of the mass flux. In order to avoid the ambiguity near the sonic line, the density is updated in terms of the velocity, which is obtained through a simple integration of a first order equation step by step in the flow field. Iteration algorithms and finite difference approximations are discussed and numerical results of both conservative and nonconservative calculations are presented.

Hafez, M. M.↗
NASA NTRS1989

Transonic flutter calculations using the Euler equations

In transonic flutter problems where shock motion plays an important part, it is believed that accurate predictions of the flutter boundaries will require the use of codes based on the Euler equations. Only Euler codes can obtain the correct shock location and shock strength, and the crucially important shock excursion amplitude and phase lag. The present study is based on the finite volume scheme developed by Jameson and Venkatakrishnan for the 2-D unsteady Euler equations. The equations are solved in integral form on a moving grid. The variable are pressure, density, Cartesian velocity components, and total energy.

Bendiksen, Oddvar O.↗
NASA NTRS1990

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗
NASA NTRS1992

Energetic particle signatures of satellites and rings in Neptune's magnetosphere

The cosmic ray system on Voyager 2 found a trapped radiation environment in Neptune's inner magnetosphere which is controlled primarily by absorption at the rings and satellite surfaces. The intensity of electrons with kinetic energies approximately greater than 1 MeV shows particularly strong and narrow signatures associated with absorption by the satellite 1989N1 at an orbital radius of 4.75 Neptune radii. Closer to the planet are several signatures of the inner satellites and rings. Absorption limits the intensity of the inner radiation belt sufficiently for the maximum intensity to occur outside the orbit of 1989N1 at a magnetic L shell of about 7. Radial profiles of the electron phase space density show that electrons diffuse inward from a source in the outer magnetosphere. Many of the inward-diffusing electrons are absorbed upon reaching a satellite orbital radius, but the finite absorption efficiency allows some of the electrons to pass by unaffected. The locations of the satellite and ring signatures also provide constraints on the nondipolar components of the planetary magnetic field.

Selesnick, R. S.↗
NASA NTRS2012

Effects of Subscale Size and Shape on Global Energy Dissipation in a Multiscale Model of a Fiber-Reinforced Composite Exhibiting Post-Peak Strain Softening Using Abaqus and FEAMAC

A mesh objective crack band model is implemented in the generalized method of cells (GMC) micromechanics model to predict failure of a composite repeating unit cell (RUC). The micromechanics calculations are achieved using the MAC/GMC core engine within the ImMAC suite of micromechanics codes, developed at the NASA Glenn Research Center. The microscale RUC is linked to a macroscale Abaqus/Standard finite element model using the FEAMAC multiscale framework (included in the ImMAC suite). The effects of the relationship between the characteristic length of the finite element and the size of the microscale RUC on the total energy dissipation of the multiscale model are investigated. A simple 2-D composite square subjected to uniaxial tension is used to demonstrate the effects of scaling the dimensions of the RUC such that the length of the sides of the RUC are equal to the characteristic length of the finite element. These results are compared to simulations where the size of the RUC is fixed, independent of the element size. Simulations are carried out for a variety of mesh densities and element shapes, including square and triangular. Results indicate that a consistent size and shape must be used to yield preserve energy dissipation across the scales.

Pineda, Evan, J.↗
NASA NTRS1970

Prediction of electron density in the boundary layer on entry vehicles with ablation

A rather general computer code has been developed for the numerical solution of the laminar boundary layer and thin shock layer (stagnation point only) equations for a multi-component gas mixture with finite reaction rates. The purpose of this paper is to indicate the capabilities of this computer program by presenting solutions for pure air flows and flows with carbon ablation. At the body surface the effects of oxidation and sublimation of carbon are taken into account with the following species included in the gas model: 0 2 , N 2 , O, N, NO, NO + , C 1 , C 2 , C 3 , CN, CO, and CO 2 . Results for the boundary layer along a 10° half-angle hyperboloid with ablation of carbon are compared to pure air results. Shock layer results for several wall temperatures and altitudes are presented and compared to pure air results. Other types of problems that can be solved with this computer code are indicated.

F. G. Blottner↗
DOE OSTI2021

Linear stability analysis of two fluid columns of different densities and viscosities in a gravity field

The linear stability of a vertical interface separating two miscible fluid columns of different densities and viscosities under the influence of gravity is investigated. This flow possesses a time-dependent reference state (each column accelerates at different rates owing to their different densities) and the interface thickness grows as the square root of time (by diffusion). Numerical integration of the linear initial-value problem is carried out and discussed in detail as a function of vertical and spanwise wavenumbers and the flow parameters. Adjoint-based optimization is performed in order to determine initial conditions that lead to maximum growth of disturbances in finite time. Results indicate that the rate of growth of the perturbation energy at small wavenumbers (less affected by viscosity initially) is dominated by two-dimensional modes (no spanwise variation). Substantial transient growth is observed at higher wave modes initially, followed by asymptotic decay of the perturbations at large time. Sensitivity of perturbation growth with respect to initial time, density and viscosity ratios is investigated. This work is complementary to previous inviscid analysis of this configuration, which showed that the interface was unconditionally unstable at all wave modes, even in the presence of surface tension, and that instability grew as the exponential of time squared.

Mechanics↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
NASA NTRS2021

East–West Proton Flux Anisotropy Observed with the PAMELA Mission

We present a study of the east–west anisotropy of trapped-proton fluxes in low-Earth orbit based on the measurements of the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics(PAMELA)experiment. The differential intensities of eastward- and westward-traveling protons detected in the South Atlantic Anomaly region were estimated as a function of equatorial pitch angle and drift shell, for six energy bins between80 MeV and 2 GeV. We found that, as a consequence of the strong atmospheric gradient coupled with the large gyro radius in this energy range, the intensities of eastward fluxes exceed those of westward fluxes by a factor of∼10–20. However, the reported directional asymmetry also depends on the sign of the local flux gradient, resulting in more intense westward fluxes beyond the radial distances where the inner belt peaks. PAMELA observations can be used to improve the description of the near-Earth radiation environment at lowest altitudes and highest trapping energies, where current theoretical and empirical models are affected by the largest uncertainties. Unified Astronomy Thesaurus concepts: Cosmic rays(329);Van Allen radiation belts(1758)1. Introduction Low-altitude inner-belt protons are strongly influenced by the density distribution of Earth’s atmosphere, mostly through interactions with its neutral constituents, which induce significant flux anisotropies. In prim is, the atmospheric loss cone results in a steep pitch-angle distribution, which becomes narrower for lower drift shells. A further, azimuthal anisotropy originates from finite gyro radius effects at proton energies in excess of a few tens of MeV(Haerendel1962; Lenchek & Singer1962). In fact, for a given spacecraft position, protons with the same pitch angle but different gyro phase—the azimuth angle associated with the gyration motion—have their guiding centers on different drift shells. In particular, protons from the west and from the east gyrate around magnetic field lines located at higher and lower altitudes, respectively. The guiding-center separationΔhincreaseswith increasing energy, so that protons moving eastward will encounter progressively lower drift-averaged densities, thus experiencing less atmospheric absorption; the opposite situation will occur for protons traveling westward, resulting in an east–west asymmetry of flux intensities whenΔhbecomes comparable to or larger than the flux scale height(Garmire1963; Heckman &Nakano1963).The trapped-flux anisotropy is a relevant aspect of the modeling of the low Earth orbit(LEO)radiation environment, given the significant engineering implications, especially for The Astrophysical Journal,919:114(6pp), 2021 October 1https://doi.org/10.3847/1538-4357/ac1677© 2021. The American Astronomical Society.

A. Bruno↗
DOE OSTI2023

Influence of the symmetry energy on the nuclear binding energies and the neutron drip line position

A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗