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 325 records · Page 18

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing the Neutron Skin with Ultrarelativistic Isobaric Collisions

Particle production in ultrarelativistic heavy ion collisions depends on the details of the nucleon density distributions in the colliding nuclei. We demonstrate that the charged hadron multiplicity distributions in isobaric collisions at ultrarelativistic energies provide a novel approach to determine the poorly known neutron density distributions and thus the neutron skin thickness in finite nuclei, which can in turn put stringent constraints on the nuclear symmetry energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stability, metallicity, and magnetism in niobium silicide nanofilms

Modern superconducting qubits based on two-dimensional (2D) transmons typically involve the growth of Nb thin films on high-resistivity Si substrates. Since imperfections at the Nb-Si heterointerface have been implicated as a source of two-level systems that limit quantum coherence times, detailed characterization and understanding of niobium silicide interfacial layers are critical to improving superconducting qubit technology. While bulk binary intermetallic niobium silicide phases are well understood, the thermodynamic phase stability and properties of ultrathin niobium silicides, such as those found at the Nb-Si heterointerface in 2D transmons, have not yet been explored. Here, we report finite-sized effects for ultrathin niobium silicide films using density functional theory calculations and predict nanoscale stabilization of Nb 6 Si 5 over the bulk α-Nb 5 Si 3 phase. This result is consistent with our experimental observations of a niobium silicide interfacial layer between a sputtered Nb thin film and the underlying Si substrate. Furthermore, our calculations show that Nb 6 Si 5 nanofilms are nonmagnetic, making them superior to nanofilms of α-Nb 5 Si 3 that exhibit antiferromagnetic correlations detrimental to long coherence times in superconducting qubits. Furthermore, by providing atomic-scale insight into niobium silicide nanofilms, this paper can help guide ongoing efforts to optimize Nb-Si heterointerfaces for long coherence times in superconducting qubits.

36 MATERIALS SCIENCE↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗

Vacuum polarization in a strong external field

Vacuum polarization in the field of a high-Z finite-size nucleus is examined, and the polarization charge density in coordinate space of order alpha(Z alpha)exp n not less than 3 is calculated. Energy-level shifts of K- and L-shell electrons in hydrogenlike systems are given.

Soff, Gerhard↗

A simple algebraic grid adaptation scheme with applications to two- and three-dimensional flow problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented which provides high flowfield resolution when applied to two-dimensional and three-dimensional flow problems. The present scheme locally adjusts the grid density using gradients of selected flow variables from either finite difference or finite volume calculations, and it can specify user-prescribed grid stretching such that control of the grid spacing can be maintained in areas of known flowfield behavior. A robust and efficient grid smoothing technique has been incorporated into the adaptive grid routine.

Hsu, Andrew T.↗

Recent advances in PDF modeling of turbulent reacting flows

This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Leonard, Andrew D.↗

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↗

Comparison of Fatigue Life Estimation Using Equivalent Linearization and Time Domain Simulation Methods

The Monte Carlo simulation method in conjunction with the finite element large deflection modal formulation are used to estimate fatigue life of aircraft panels subjected to stationary Gaussian band-limited white-noise excitations. Ten loading cases varying from 106 dB to 160 dB OASPL with bandwidth 1024 Hz are considered. For each load case, response statistics are obtained from an ensemble of 10 response time histories. The finite element nonlinear modal procedure yields time histories, probability density functions (PDF), power spectral densities and higher statistical moments of the maximum deflection and stress/strain. The method of moments of PSD with Dirlik's approach is employed to estimate the panel fatigue life.

Mei, Chuh↗

Current sheet Formation in a Conical Theta Pinch Faraday Accelerator with Radio-Frequency Assisted Discharge

The inductive formation of current sheets in a conical theta pinch FARAD (Faraday Accelerator with Radio-frequency Assisted Discharge) thruster is investigated experimentally with time-integrated photography. The goal is to help in understanding the mechanisms and conditions controlling the strength and extent of the current sheet, which are two indices important for FARAD as a propulsion concept. The profiles of these two indices along the inside walls of the conical acceleration coil are assumed to be related to the profiles of the strength and extent of the luminosity pattern derived from photographs of the discharge. The variations of these profiles as a function of uniform back-fill neutral pressure (with no background magnetic field and all parameters held constant) provided the first clues on the nature and qualitative dependencies of current sheet formation. It was found that there is an optimal pressure for which both indices reach a maximum and that the rate of change in these indices with pressure differs on either side of this optimal pressure. This allowed the inference that current sheet formation follows a Townsend-like breakdown mechanism modified by the existence of a finite pressure-dependent radio-frequency-generated electron density background. The observation that the effective location of the luminosity pattern favors the exit-half of the conical coil is explained as the result of the tendency of the inductive discharge circuit to operate near its minimal self-inductance. Movement of the peak in the luminosity pattern towards the upstream side of the cone with increasing pressure is believed to result from the need of the circuit to compensate for the increase in background plasma resistivity due to increasing pressure.

Hallock, Ashley K.↗

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comprehensive Studies on Hot Compaction and Vibration-Assisted Compaction Tests of Aluminum Powder

Abstract Aluminum powder compaction was studied using both test and simulation. Cold compaction, hot compaction, and vibration-assisted (cold) compaction tests were conducted to achieve different density ratios. First, the hot compaction test (at 300 °C, compression pressure 140 MPa) improved about 6% compared with cold compaction under the same compression pressure. Second, although the relative density ratio does not obviously improve at a vibration-assisted (cold) compaction, the strength of the specimens made under vibration loading is much better than those of cold compaction. Additionally, finite element models with well-calibrated Drucker–Prager Cap (DPC) material constitutive model were built in abaqus/standard to simulate the powder compaction process. The results of the finite element model have very good correlations with test results up to the tested range, and this finite element model further predicts the loading conditions needed to achieve the higher density ratios. Two exponential equations of the predicted density ratio were obtained by combining the test data and the simulation results. A new analytical solution was developed to predict the axial pressure versus the density ratio for the powder compaction according to DPC material model. The results between the analytical solution and the simulation model have a very good match.

Engineering↗

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗

Nonequilibrium spintronic transport through Kondo impurities

In this work we analyze the nonequilibrium transport through a quantum impurity (quantum dot or molecule) attached to ferromagnetic leads by using a hybrid numerical renormalization group–time-dependent density matrix renormalization group thermofield quench approach. For this, we study the bias dependence of the differential conductance through the system, which shows a finite zero-bias peak, characteristic of the Kondo resonance and reminiscent of the equilibrium local density of states. In the nonequilibrium settings, the resonance in the differential conductance is also found to decrease with increasing the lead spin polarization. The latter induces an effective exchange field that lifts the spin degeneracy of the dot level. Therefore, as we demonstrate, the Kondo resonance can be restored by counteracting the exchange field with a finite external magnetic field applied to the system. Lastly, we investigate the influence of temperature on the nonequilibrium conductance, focusing on the split Kondo resonance. Our work thus provides an accurate quantitative description of the spin-resolved transport properties relevant for quantum dots and molecules embedded in magnetic tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kinetically constrained freezing transition in a dipole-conserving system

Here, we study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fractonlike conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase, where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗