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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗
DOE OSTI2023

Enhancement of charge density wave correlations in a Holstein model with an anharmonic phonon potential

The Holstein Hamiltonian describes itinerant electrons whose site density couples to local phonon degrees of freedom. In the single-site limit, at half filling, the electron-phonon coupling results in a double-well structure for the lattice displacement, favoring empty or doubly occupied sites. In two dimensions and on a bipartite lattice in d ≥ 2, intersite hopping causes these doubly occupied and empty sites to alternate in a charge density wave (CDW) pattern when the temperature is lowered. Because a discrete symmetry is broken, this occurs in a conventional second-order transition at finite T cdw . In this paper, we investigate the effect of changing the phonon potential energy to one with an intrinsic double-well structure even in the absence of electron-phonon coupling. While this aids in the initial process of pair formation, the implications for subsequent CDW order are nontrivial. One expects that, when the electron-phonon coupling is too strong, the double wells become deep and the polaron mass gets large, an effect which reduces T cdw . Here, we show here the existence of regions of parameter space where the double-well potential, while aiding local pair formation, does so in a way which also substantially enhances long-range CDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

On the mechanical behavior of a Nb 3 Sn superconducting coil during a quench: Two-dimensional finite element analysis of a quench heater protected magnet

New high-field accelerator magnets based on $Nb_3Sn$ superconductors are pushing the boundaries of magnet design and quench protection towards new limits. While their large stored energies and current densities result in a very challenging scenario for magnet protection, their great electromagnetic forces create also new requirements in terms of magnet design and stress management techniques. Furthermore, the strain sensitivity of $Nb_3Sn$ cables turns the electro-mechanical limits of the conductor into a parameter of the highest importance, where conductor degradation becomes a critical aspect in magnet operation. The coupling of all the above-mentioned considerations during quench is a case of special interest that adds further complexity to the design of $Nb_3Sn$ magnets. The objective of this paper is to provide a complete two-dimensional investigation of the coil and magnet structure mechanics during a quench event. Here, the analysis is performed using a combination of finite element codes that provide the necessary input for the mechanical study. The core of the modelling strategy relies on a thermal-electric model, whose results are directly used as loads in the mechanical simulation. In doing so, the stress evolution during and after a quench is obtained. We focus for this time in the analysis of a quench heater protected magnet, where dynamic effects are of less importance in contrast to other protection systems like CLIQ.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene

We employ a functional renormalization group approach to ascertain the pairing mechanism and symmetry of the superconducting phase observed in rhombohedral trilayer graphene. Superconductivity in this system occurs in a regime of carrier density and displacement field with a weakly distorted annular Fermi sea. We find that repulsive Coulomb interactions can induce electron pairing on the Fermi surface by taking advantage of momentum-space structure associated with the finite width of the Fermi sea annulus. Furthermore, the degeneracy between spin-singlet and spin-triplet pairing is lifted by valley-exchange interactions that strengthen under the RG flow and develop nontrivial momentum-space structure. We find that the leading pairing instability is d-wave-like and spin-singlet, and that the theoretical phase diagram versus carrier density and displacement field agrees qualitatively with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling of plasma beta effects on the island divertor transport in the standard configuration of W7-X

The influence of plasma beta effects on the edge plasma transport in the Wendelstein 7-X standard configuration is studied systematically by using EMC3-EIRENE combined with a 3D equilibrium code named HINT. The magnetic topology changes induced by plasma beta effects are significantly reflected in plasma transport behaviors and heat flux patterns on divertor targets. After validating the modeling strategy by comparisons with experimental data, the extended simulations for high performance plasmas show that the threshold separatrix density for accessing the power detachment is reduced in higher beta plasmas. Compared with the vacuum field case, the impurity radiation distributions with finite beta effects are modified in the magnetic island region. The divertor heat flux is distributed more evenly along the toroidal direction on the strike line at the vertical target. The strike line on the horizontal target moves towards the pumping gap with an increase in the plasma beta. In addition, the different pressure profiles with the same central beta also result in a modified heat flux pattern on the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hall viscosity of the composite-fermion Fermi seas for fermions and bosons

The Hall viscosity has been proposed as a topological property of incompressible fractional quantum Hall states and can be evaluated as Berry curvature. This paper reports on the Hall viscosities of composite-fermion Fermi seas at ν = 1/m, where m is even for fermions and odd for bosons. A well-defined value for the Hall viscosity is not obtained by viewing the 1/m composite-fermion Fermi seas as the n → ∞ limit of the Jain ν = n/(nm ± 1) states, whose Hall viscosities (±n + m)ℏρ/4 (ρ is the two-dimensional density) approach ±∞ in the limit n → ∞. A direct calculation shows that the Hall viscosities of the composite-fermion Fermi sea states are finite and also relatively stable with system size variation, although they are not topologically quantized in the entire τ space. Here, I find that the ν = 1/2 composite-fermion Fermi sea wave function for a square torus yields a Hall viscosity that is expected from particle-hole symmetry and is also consistent with the orbital spin of 1/2 for Dirac composite fermions. I compare my numerical results with some theoretical conjectures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Trap-Assisted Dopant Compensation Prevents Shunting in poly-Si Passivating Interdigitated Back Contact Silicon Solar Cells

Interdigitated back contact (IBC) solar cells achieve the highest efficiencies of single-junction architectures, but complicated patterning of the rear fingers and spreading of dopants during processing inhibit their mainstream adoption due to concerns of shunting between the IBC fingers. One method of simplifying patterning at the rear is by using contact masks combined with plasma-enhanced chemical vapor deposition (PECVD) or ion implantation. However, the intrinsic isolation region becomes contaminated during high-temperature annealing by lateral diffusion of dopants and during masked PECVD by spreading of dopant radicals through region between the mask and the substrate. Despite this contamination, we show through scanning spreading resistance microscopy and Kelvin probe force microscopy that a ~20 µm wide compensating region exists with high enough resistivity to prevent shunting. We model this p-i-n poly-Si system using two simulation models: a simple resistor model considering only the capture of charge carriers by trap defects in poly-Si to reduce the conductivity, and a more refined 1-dimensional finite element model using Poisson’s equation, drift-diffusion equations, and recombination of carriers. Using this model, we show that high defect density significantly decreases the current across the region between the p- and n-type fingers, preventing shunting.

atom probe↗
DOE OSTI2021

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗
DOE OSTI2024

Transducer Resolution Effect on Pressure Fluctuations Beneath Hypersonic Turbulent Boundary Layers

The size of a pressure transducer is known to affect the accuracy of measurements of wall-pressure fluctuations beneath a turbulent boundary layer because of spatial averaging over the sensing area of the transducer. In this paper, the effect of finite transducer size is investigated by applying spatial averaging or wavenumber filters to a database of hypersonic wall pressure generated from a direct numerical simulation (DNS) that simulates the turbulent portion of the boundary layer over a sharp 7° half-angle cone at nominally Mach 8. Here, a good comparison between the DNS and the experiment in the Sandia Hypersonic Wind Tunnel at Mach 8 is achieved after spatial averaging is applied to the DNS data over an area similar to the sensing area of the transducer. The study shows that a finite sensor size similar to that of the PCB132 transducer can cause significant attenuation in the root-mean-square and power spectral density (PSD) of wall-pressure fluctuations, and the attenuation effect is identical between cone and flat plate configurations at the same friction Reynolds number. The Corcos theory is found to successfully compensate for the attenuated high-frequency components of the wall-pressure PSD.

42 ENGINEERING↗
DOE OSTI2024

Simulated structure and thermodynamics of decagonal Al-Co-Cu quasicrystals

Atomic structures of Al-Co-Cu decagonal quasicrystals (dQCs) are investigated using empirical oscillating pair potentials (EOPP) in molecular dynamic (MD) simulations that we enhance by Monte Carlo (MC) swapping of chemical species and replica exchange. Predicted structures exhibit planar decagonal tiling patterns and are periodic along the perpendicular direction. We then recalculate the energies of promising structures using first-principles density functional theory (DFT), along with energies of competing phases. We find that our τ -inflated sequence of QC approximants (QCAs) are energetically unstable at low temperature by at least 3 meV/atom. Extending our study to finite temperatures by calculating harmonic vibrational entropy, as well as anharmonic contributions that include chemical species swaps and tile flips, our results suggest that the quasicrystal phase is entropically stabilized at temperatures in the range 600-800 K and above. It decomposes into ordinary (though complex) crystal phases at low temperatures, including a partially disordered B2-type phase. We discuss the influence of density and composition on QC phase stability; we compare the structural differences between Co-rich and Cu-rich quasicrystals; and we analyze the role of entropy in stabilizing the quasicrystal, concluding with a discussion of the possible existence of “high entropy” quasicrystals. Published by the American Physical Society 2024

Huang, Yang (ORCID:0009000045917347)↗
DOE OSTI2021

Atomic-scale mechanisms of void strengthening in tungsten

Neutron and heavy-ion irradiation of tungsten produces nanometer-size vacancy voids, gas-filled bubbles and dislocation loops. These defect features can affect mechanical properties and the impact can be quite significant because of their high density. Understanding the basic mechanisms of mechanical properties degradation is necessary for predicting radiation effects. Predictions can be made using discrete dislocation dynamics or/and finite element approaches which, however, need local interaction mechanisms as inputs. Such knowledge can be provided only by atomic-scale modeling. This paper reports the results of an extensive atomic-scale modeling study of the interactions between moving edge dislocations and voids in tungsten. Furthermore, the main focus is on the effects of the void size and ambient temperature. Critical resolved shear stress was calculated for voids up to 9 nm in diameter. Atomistic results are compared with the theoretical approach and with those obtained earlier for voids in body centered cubic (bcc) iron. An important role of the void surface has been revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigation of O( 3 P) Initiated Oxidation Products of 2,3-Dimethylfuran Using Synchrotron Photoionization

Biofuels are an appealing alternative to augment conventional fossil fuels, given that traditional fossil fuels are finite in nature despite being the highest-consumed energy source. Furans have been studied as potential biofuels due to their renewability, high energy density, and ability to both blend with existing fossil fuels and serve as a stand-alone fuel. Although promising, more research efforts on the combustion of furan-based biofuels are needed. Here, in this study, the oxidation of 2,3-dimethylfuran (2,3-DMF) initiated by O( 3 P) was investigated using vacuum-ultraviolet synchrotron radiation from the Advanced Light Source at the Lawrence Berkeley National Laboratory. The reaction was studied at 550 K and 7 Torr. Major reaction products include 3-methyl-3-buten-2-one, 4,5-dimethyl-3(2H)-furanone, and 3-methyl-4-oxo-2(Z)-pentenal, with 3-methyl-3-buten-2-one being the most abundant.

CBS-QB3↗
DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Preparing quantum many-body scar states on quantum computers

Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to typical eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols that enable the use of quantum computers to study this question. We present protocols both for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we present both a system-size-linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we formulate an exact state preparation approach based on matrix product states that yields quasipolynomial-depth circuits, as well as a variational approach with a polynomial-depth ansatz circuit. We also provide proof of principle state-preparation demonstrations on superconducting quantum hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Searching for the QCD Critical Point Along the Pseudo-critical/Freeze-out Line Using Padé-resummed Taylor Expansions of Cumulants of Conserved Charge Fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, T pc ≃ 156.5 MeV, we find the estimator of the radius of convergence to be µ B /T ≳ 3 for strangeness-neutral matter. We also present results for the pole structure of the Pad´e approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Pad´e is consistent with not having a critical point at temperatures larger than 135 MeV and a baryon chemical potential smaller than µ B /T ~ 2.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Assessment of Numerical Diffusion in NRELAP5 Code for Density Wave Oscillations Applications

Density wave oscillation (DWO) in a boiling channel can be a delicate phenomenon, and the precision of numerical tools can influence predictions. In this paper, relevant numerical diffusion was analyzed for the finite difference schemes used in the NRELAP5 code. The numerical diffusion in NRELAP5 was investigated under both single-phase and two-phase conditions. A unique perturbation technique called the “V” ramping approach was introduced, and the NRELAP5-predicted DWO results were then compared to test data.

Numerical diffusion, NRELAP5, RELAP5-3D, density w↗