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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Charging in the ac Conductance of a Double Barrier Resonant Tunneling Structure

There have been many studies of the linear response ac conductance of a double barrier resonant tunneling structure (DBRTS), both at zero and finite dc biases. While these studies are important, they fail to self consistently include the effect of the time dependent charge density in the well. In this paper, we calculate the ac conductance at both zero and finite do biases by including the effect of the time dependent charge density in the well in a self consistent manner. The charge density in the well contributes to both the flow of displacement currents in the contacts and the time dependent potential in the well. We find that including these effects can make a significant difference to the ac conductance and the total ac current is not equal to the simple average of the non-selfconsistently calculated conduction currents in the two contacts. This is illustrated by comparing the results obtained with and without the effect of the time dependent charge density included correctly. Some possible experimental scenarios to observe these effects are suggested.

Anantram, M. P.↗
DOE OSTI2025

Pauli potential formalism at finite temperature

At zero temperature, the Pauli potential—the functional derivative of the Pauli kinetic energy density functional—is the key to the accuracy of the orbital-free density functional theory (OFDFT) as it is supposed to capture all the effects associated with the Pauli exclusion principle. Here, we extend this concept to finite temperature by defining Pauli free energy and the modified Pauli free energy, both representing the natural generalizations of the Pauli term from zero- T to finite- T . We discuss their physical interpretation, the mathematical nuances, and the applicability, arguing that the modified Pauli potential should be used as an extension of the zero- T counterpart within the OF-DFT framework. Through analytical and numerical methods, we then analyze some of the exact properties concerning the modified Pauli free- and kinetic-energy terms and examine the temperature dependence of the modified Pauli potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modeling inter- and intra-granular dislocation transport using crystal plasticity

Here, this work presents the development of a crystal plasticity material model that incorporates both dislocation transport within grains and dislocation transfer across grain boundaries. This model has been implemented in the open-source finite element code MOOSE. In addition, a novel geometry-based criterion is developed to determine the direction of dislocation transfer across grain boundaries. The transfer criterion incorporates the geometric features of the grain boundary, such as the grain boundary plane normal, and its misorientation, which is accounted for through the orientation of the incoming and outgoing slip systems. The model is tested with several cases, including a copper single crystal, bi-crystal, and polycrystal. The development of the transfer criterion, implementation of the model, and its application to these test cases are discussed in detail.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Geometrically Nonlinear Finite Element Analysis of a Composite Space Reflector

Lightweight aerospace structures, such as low areal density composite space reflectors, are highly flexible and may undergo large deflection under applied loading, especially during the launch phase. Accordingly, geometrically nonlinear analysis that takes into account the effect of finite rotation may be needed to determine the deformed shape for a clearance check and the stress and strain state to ensure structural integrity. In this study, deformation of the space reflector is determined under static conditions using a geometrically nonlinear solid shell finite element model. For the solid shell element formulation, the kinematics of deformation is described by six variables that are purely vector components. Because rotational angles are not used, this approach is free of the limitations of small angle increments. This also allows easy connections between substructures and large load increments with respect to the conventional shell formulation using rotational parameters. Geometrically nonlinear analyses were carried out for three cases of static point loads applied at selected points. A chart shows results for a case when the load is applied at the center point of the reflector dish. The computed results capture the nonlinear behavior of the composite reflector as the applied load increases. Also, they are in good agreement with the data obtained by experiments.

Lee, Kee-Joo↗
DOE OSTI2023

FINITE ELEMENT MODEL MESH REFINEMENT EFFECTS ON QUALIFICATION OF NUCLEAR GRADE GRAPHITE COMPONENTS

The American Society of Mechanical Engineers (ASME) provides the full and simplified design-by-analysis probabilistic assessments for determining acceptance of nuclear grade graphite core components. The assessments can be characterized by three parts: (1) a component stress distribution, often determined by a finite element (FE) model; (2) a Weibull probability density function (pdf) that characterizes the experimental tensile strength distribution; and (3) the post-processor, which combines the FE model and the Weibull strength distribution in accordance with the full and simplified assessments to determine component acceptance. It is known that the level of mesh refinement in FE models can affect the modeled component’s calculated stress distribution. Depending on the component geometry, the stress distribution may converge with sufficient refinement. It was previously unknown whether the acceptance decision resulting from the full and simplified assessments might change even with sufficient mesh refinement. This study explores that question using experimental strength results for a dog-bone geometry for two graphite grades, IG-110 and PCEA. The simplified assessment has two criteria that must be met, the first limits the combined membrane stress by the allowable stress and the second limits the peak equivalent stress by the allowable stress scaled by the ratio of flexural to tensile strength. In the application of the simplified assessment, convergence of the peak equivalent stress required extreme mesh refinement, however, the acceptance decision was not affected. It is hypothesized that more complex geometries with stress concentrations may present mesh refinement effects on the simplified assessment acceptance decision. Mesh refinement did affect the acceptance decision in the full assessment for the applied pressure loadings in this study. This work suggests component stress distribution convergence is not a sufficient criteria for POF convergence in the full assessment and that mesh refinement should continue until the POF has converged, especially where the resulting POF is bordering the SRC acceptable POF limit.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Zero-pairing and zero-temperature limits of finite-temperature Hartree-Fock-Bogoliubov theory

Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is non-trivial for open-shell systems such that HFB theory does {\it not} reduce to the Hartree-Fock (HF) formalism in all cases. The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits. The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction. While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e. for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the non-trivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e. the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases. In conclusion, while in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a non trivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.

20 ≤ A ≤ 38↗
DOE OSTI2023

Microscopic Origin of Polarization Charges at Ga N / ( Al , Ga ) N Interfaces

Ga⁢N/(Al,Ga)⁢N heterojunctions are at the heart of high-electron-mobility transistors that are being adopted for high-power and high-frequency applications. The strong polarization fields present at this interface significantly enhance the density of the two-dimensional electron gas (2DEG) that is confined on the Ga⁢ N side of the junction. The microscopic origin of these electrons has been debated over the years: after excluding that they would be contributed by bulk donors, a model that identifies surface states on the (Al, Ga)⁢N surface as the source of electrons has become widely adopted. Recently, it has become clear, however, that the measured density of surface states is insufficient to account for the high electron density in the 2DEG. Here we demonstrate, based on state-of-the-art first-principles calculations, that the source of electrons is intrinsic to the overall structure and that the negative charge in the 2DEG is balanced by fixed charge on the surface. We perform a rigorous study of polarization, using our recently developed methodology for quantifying polarization fields within the finite-sized systems that can be addressed with density-functional calculations. The results show that the electrons that appear in the 2DEG originate locally at the interface, and that the net charge at the interface is predominantly compensated by fixed charge on the surface, rather than surface states. We elucidate the source of this fixed charge and associate it with surface reconstructions or the presence of heterovalent impurities (such as oxygen). Our results force a reassessment of the impact of surface states on the density of the 2DEG: rather than serving as a supply of electrons, the surface states mainly act to pin the Fermi level. Our conclusions allow a fresh interpretation of experimental observations and allow devising guidelines for optimizing carrier densities in the 2DEG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2021

Multiple fermion scattering in the weakly coupled spin-chain compound YbAlO 3

The Heisenberg antiferromagnetic spin-1/2 chain, originally introduced almost a century ago, is one of the best studied models in quantum mechanics due to its exact solution, but nevertheless it continues to present new discoveries. Its low-energy physics is described by the Tomonaga-Luttinger liquid of spinless fermions, similar to the conduction electrons in one-dimensional metals. In this work we investigate the Heisenberg spin-chain compound YbAlO3 and show that the weak interchain coupling causes Umklapp scattering between the left- and right-moving fermions and stabilizes an incommensurate spin-density wave order at q = 2k F under finite magnetic fields. These Umklapp processes open a route to multiple coherent scattering of fermions, which results in the formation of satellites at integer multiples of the incommensurate fundamental wavevector Q = nq. Our work provides surprising and profound insight into bandstructure control for emergent fermions in quantum materials, and shows how neutron diffraction can be applied to investigate the phenomenon of coherent multiple scattering in metals through the proxy of quantum magnetic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

A phenomenological pulsar model

Particle injection energies and rates previously calculated for the stellar wind generation by rotating magnetized neutron stars are adopted. It is assumed that the ambient space-charge density being emitted to form this wind is bunched. These considerations immediately place the coherent radio frequency luminosity from such bunches near 10 to the 28th erg/s for typical pulsar parameters. A comparable amount of incoherent radiation is emitted for typical (1 second) pulsars. For very rapid pulsars, however, the latter component grows more rapidly than the available energy sources. The comparatively low radio luminosity of the Crab and Vela pulsars is attributed to both components being limited in the same ratio. The incoherent radiation essentially has a synchotron spectrum and extends to gamma-ray energies; consequently the small part of the total luminosity that is at optical wavelengths is unobservable. Assuming full coherence at all wavelengths short of a critical length gives a spectral index for the flux density of -8/3 at higher frequencies. The finite energy available from the injected particles would force the spectrum to roll over below about 100 MHz, although intrinsic morphological factors probably enter for any specific pulsar as well.

Michel, F. C.↗
NASA NTRS1987

Ground-state properties of the periodic Anderson model

The ground-state energy, hybridization matrix element, local moment, and spin-density correlations of a one-dimensional, finite-chain, periodic, symmetric Anderson model are obtained by numerical simulations and compared with perturbation theory and strong-coupling results. It is found that the local f-electron spins are compensated by correlation with other f-electrons as well as band electrons leading to a nonmagnetic ground state.

Blankenbecler, R.↗
NASA NTRS1989

MHD surface waves in high- and low-beta plasmas. I - Normal-mode solutions

Since the first paper by Barston (1964) on electrostatic oscillations in inhomogeneous cold plasmas, it has been commonly accepted that all finite layers with a continuous profile in pressure, density and magnetic field cannot support normal surface waves but instead the waves always decay through phase mixing (also called resonant absorption). The problem is reanalyzed by studying a compressible current sheet of a general structure with rotation of the magnetic field included. All inhomogeneous layers considered in the high-beta plasma limit do not support normal modes. However, in the limit of a low-beta plasma there are some cases when normal-mode solutions are recovered. The latter means that the process of resonant absorption is not common for all inhomogeneous layers.

Musielak, Z. E.↗
NASA NTRS1990

Pre-mixed flame simulations for non-unity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet region is to wrinkle the flame fronts. For non-unity Lewis numbers (Le), the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of the three dimensional, constant density, decaying isotropic turbulence with a single step, finite rate chemical reaction. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. The turbulent flame speed, S(sub T), increases as Le decreases. The correlation between S(sub T) and u prime found in previous Le = 1 simulations has a strong Lewis number dependency. The variance of the pdf of the flame curvature increases as Le decreases, indicating that the flames become more wrinkled. A strong correlation between local flame speed and curvature was found. For Le greater than 1, the flame speed increases in regions concave towards the products and decreases in convex regions. The opposite correlation was found for Le less than 1. The mean temperature of the products was also found to vary with Lewis number. For Le = 0.8, it is less than the adiabatic flame temperature and for Le = 1.2 it is greater.

Rutland, C. J.↗
NASA NTRS2014

Computational Design of Materials: Planetary Entry to Electric Aircraft and Beyond

NASA's projects and missions push the bounds of what is possible. To support the agency's work, materials development must stay on the cutting edge in order to keep pace. Today, researchers at NASA Ames Research Center perform multiscale modeling to aid the development of new materials and provide insight into existing ones. Multiscale modeling enables researchers to determine micro- and macroscale properties by connecting computational methods ranging from the atomic level (density functional theory, molecular dynamics) to the macroscale (finite element method). The output of one level is passed on as input to the next level, creating a powerful predictive model.

Materials Design↗
DOE OSTI2025

DROP DURABILITY ASSESSMENT OF ELECTRONIC ASSEMBLIES UNDER OFF-AXIS LOADING WITH SKEWED FIXTURES

This thesis studies drop durability of electronic assemblies when the acceleration vector is oriented at 45° to the out-of-plane direction of the circuit card. The off-axis drop tests are accomplished with a skewed fixture and are conducted as a proxy for multiaxial drop testing. Advanced shock testing and vibration test methods have been developed over the last few decades to better represent real-world field environments during ground-based laboratory testing. However, many of these test methods require expensive and specialized equipment not available in most laboratories. An alternative approach for approximating simultaneous loading along multiple axes on conventional equipment utilizes skewed fixtures which have seen use in off-axis random vibration and drop impact testing. These methods generally rely on the conversion of a uniaxial input load from the test equipment (using a uniaxial drop tower or shaker) into a multiaxial load when resolved in the reference frame of the test article (mounted on a skewed fixture). Skewed fixture design is presented and recommendations for conducting skewed angle drop testing are introduced based on local measurements along the skewed face of the fixture to accurately monitor the impact event. Characterization tests were performed with a skewed fixture, at simultaneous acceleration loads from 500 to 3,000 g in two (in-plane and out-of-plane) directions, while meeting standard time domain tolerances. Upon experimental characterization, drop shock durability tests were conducted on a printed circuit assembly (PCA). Mean drops-to-failure were measured and quantified with Weibull statistics. Dominant solder joint failure modes were identified via failure analysis. Prior work on inclined angle impact testing is limited, and the majority of solder joint interconnect level fatigue studies are conducted considering perpendicular loading normal the circuit card. Low-cycle fatigue curves are generated based on plastic strain and plastic work density within the solder joint. A multiscale nonlinear finite element model is used to relate board-level flexure to solder joint interconnect level plastic strain. A high strain rate solder constitutive model allows for accurate modeling of solder plasticity resulting from high-impact drop shock. Fatigue parameters are computed from the Coffin-Manson relation and Palmgren-Miner damage accumulation. This work serves to apply established low-cycle fatigue methods for conventional drop shock loading (impact normal to circuit card) to non-perpendicular loading with a skewed fixture.

Hower, Jonathan [Kansas City National Security Cam↗
NASA NTRS1985

Dynamic fracture mechanics

Dynamic fracture and crack propagation concepts for ductile materials are reviewed. The equations for calculating dynamic stress integrity and the dynamic energy release rate in order to study dynamic crack propagation are provided. The stress intensity factor versus crack velocity relation is investigated. The uses of optical experimental techniques and finite element methods for fracture analyses are described. The fracture criteria for a rapidly propagating crack under mixed mode conditions are discussed; crack extension and fracture criteria under combined tension and shear loading are based on maximum circumferential stress or energy criteria such as strain energy density. The development and use of a Dugdale model and finite element models to represent crack and fracture dynamics are examined.

Kobayashi, A. S.↗