SEARCH · Engineering Papers
Results for “density distribution”
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.
High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations
The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.
Long-Range Fit: A Software Package for the Representation and Study of Long-Range Molecular Interactions
Describing intermolecular forces is fundamental to modeling and predicting the behavior of molecular systems. In particular, long-range molecular interactions─with electrostatic, induction, and dispersion as the main components─play a critical role, especially for low-temperature and low-density regimes. Long-range interactions are often described through perturbation theory, representing the electronic charge distribution via a multipolar series of the moments and polarizability tensors corresponding to each molecule. However, while the theory is well established, obtaining the resulting analytical expressions (and their practical implementation) constitutes a highly complex and system-dependent task. To address this challenge, we developed long-range-fit (LRF), an interactive and user-friendly software package designed to automate the generation and fitting of long-range interaction terms for arbitrary molecules in nondegenerate (ground or excited) electronic states. We have derived and implemented all terms up to 15th order, without approximations, via a spherical tensor representation, with symmetry adaptation to all molecular point-group symmetries. The resulting potential energy surface is compatible with most representations of the close interaction region.
Atomic Three-Dimensional Investigations of Pd Nanocatalysts for Acetylene Semi-hydrogenation
Deciphering the three-dimensional (3D) insight into nanocatalyst surfaces at the atomic level is crucial to understanding catalytic reaction mechanisms and developing high-performance catalysts. Nevertheless, better understanding the inherent insufficiency of a long-range ordered lattice in nanocatalysts is a big challenge. Here, in this work, we report the local structure of Pd nanocatalysts, which is beneficial for demonstrating the shape-structure-adsorption relationship in acetylene hydrogenation. The 5.27 nm spherical Pd catalyst (Pd sph ) shows an ethylene selectivity of 88% at complete acetylene conversion, which is much higher than those of the Pd octahedron and Pd cube and superior to other reported monometallic Pd nanocatalysts so far. By virtue of the local structure revelation combined with the atomic pair distribution function (PDF) and reverse Monte Carlo (RMC) simulation, the atomic surface distribution of the unique compressed strain of Pd-Pd pairs in Pd sph was revealed. Density functional theory calculations verified the obvious weakening of the ethylene adsorption energy on account of the surface strain of Pd sph . It is the main factor to avoid the over-hydrogenation of acetylene. The present work, entailing shape-induced surface strain manipulation and atomic 3D insight, opens a new path to understand and optimize chemical activity and selectivity in the heterogeneous catalysis process.
Can electric fields drive chemistry for an aqueous microdroplet?
Abstract Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O–H bonds at the surface are ~16 MV/cm larger on average than that found for O–H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.
Dynamic characteristics of multi-charged ions emitted from nanosecond laser produced molybdenum plasmas
The charge state distribution and temporal evolution of laser-produced molybdenum plasma were investigated at various laser power densities.
New interaction potentials for borate glasses with mixed network formers
Here, we adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination, and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reliably reproduce the structure, coordination, and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce the properties of various boroaluminate and borosilicate glasses.
Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method
In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.
Poisson_CCD: A dedicated simulator for modeling CCDs
A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.
Compact galaxies and the size–mass galaxy distribution from a colour-selected sample at 0.04 < z < 0.15 supplemented by ugrizYJHK photometric redshifts
ABSTRACT The size–mass galaxy distribution is a key diagnostic for galaxy evolution. Massive compact galaxies are potential surviving relics of a high-redshift phase of star formation. Some of these could be nearly unresolved in Sloan Digital Sky Survey (SDSS) imaging and thus not included in galaxy samples. To overcome this, a sample was selected from the combination of SDSS and UKIRT Infrared Deep Sky Survey (UKIDSS) photometry to r < 17.8. This was done using colour–colour selection, and then by obtaining accurate photometric redshifts (photo-z) using scaled flux matching (SFM). Compared to spectroscopic redshifts (spec-z), SFM obtained a 1σ scatter of 0.0125 with only 0.3 per cent outliers (|Δln (1 + z)| > 0.06). A sample of 163 186 galaxies was obtained with 0.04 < z < 0.15 over $2300\, {\rm deg}^2$ using a combination of spec-z and photo-z. Following Barro et al. log Σ1.5 = log M* − 1.5log r50, maj was used to define compactness. The spectroscopic completeness was 76 per cent for compact galaxies (log Σ1.5 > 10.5) compared to 92 per cent for normal-sized galaxies. This difference is primarily attributed to SDSS ‘fibre collisions’ and not the completeness of the main galaxy sample selection. Using environmental overdensities, this confirms that compact quiescent galaxies are significantly more likely to be found in high-density environments compared to normal-sized galaxies. By comparison with a high-redshift sample from 3D-HST, log Σ1.5 distribution functions show significant evolution, with this being a compelling way to compare with simulations such as EAGLE. The number density of compact quiescent galaxies drops by a factor of about 30 from z ∼ 2 to log (n/Mpc−3) = − 5.3 ± 0.4 in the SDSS–UKIDSS sample. The uncertainty is dominated by the steep cut off in log Σ1.5, which is demonstrated conclusively using this complete sample.
First simultaneous global QCD analysis of kaon and pion parton distributions with lattice QCD constraints
We perform the first simultaneous global QCD analysis of pion and kaon parton distribution functions (PDFs) constrained by pion- and kaon-induced Drell-Yan (DY) and leading neutron electroproduction data together with lattice QCD data on pion and kaon PDF moments. The analysis indicates a softer valence u ¯ distribution in the K − than in the π − and a significantly more peaked valence s -quark density in K − compared with the u ¯ . The effective exponent governing the high- x behavior of the PDF is found to be larger for u ¯ in the kaon, β u ¯ K − = 1.6 ( 2 ) , than in the pion, β u ¯ π − = 1.16 ( 4 ) , in the range 0.7 ≤ x ≤ 0.95 . From the gluon momentum fractions, we find the pion’s gluon content accounts for ≈ 1 / 3 of the mass budget of the pion at μ = 2 GeV , but only ≈ 1 / 4 for the kaon.
Hierarchy among the crystal lattice, charge density wave, and superconducting orders in transition metal dichalcogenides
Using high-energy x-ray scattering and large-scale three-dimensional (3D) structure modeling, in this study we investigate the relationship between the crystal lattice, charge density wave (CDW), and superconducting (SC) orders in transition metal dichalcogenides (TMDs). In particular, we systematically substitute Te for Se in Ta-Se-Te solid solutions, determine changes in their crystal lattice, and relate them to changes in the CDW transition temperature, $T_{CDW}$, and SC critical temperature, $T_c$ . We find that strong lattice distortions such as buckling of Ta layers are detrimental to the CDW and SC orders. The presence of a perfect lattice order in two dimensions is a prerequisite to the emergence of CDWs but insufficient to achieve a SC ordered state. For the SC order to emerge, the Ta sublattice should also appear periodic in 3D. Local chemical disorder may promote the SC order, and the perfectness of Ta coordination polyhedra is a factor contributing to its strength. A hierarchical relationship among the crystal lattice, CDW, and SC orders thus appears to exist in TMDs in a sense that different degrees of crystal lattice order-disorder promote and maintain the CDW and SC orders to a different extent, offering an opportunity to control the latter through modifying the former by rational design. Our findings are a step towards a better understanding of the correlation between lattice and electronic degrees of freedom in TMDs. We also demonstrate an efficient experimental approach to study them in fine detail.
Pauli energy contribution to the nucleus-nucleus interaction
Background: The Pauli exclusion principle plays a crucial role as a building block of many-body quantal systems comprised of fermions. It also induces a “Pauli repulsion” in the interaction between di-nuclear systems. It has been shown in [Phys. Rev. C 95, 031601(R) (2017)] that the Pauli repulsion widens the nucleus-nucleus potential barrier, thus hindering sub-barrier fusion. Purpose: To investigate the proton and neutron contributions to the Pauli repulsion, both in the bare potential neglecting shape polarization and transfer between the reactants, as well as in the dynamical potential obtained by accounting for such dynamical rearrangements. Methods: As the basis of our study we utilize the Pauli kinetic energy (PKE) obtained by studying the nuclear localization function (NLF). Recently this approach has been generalized to incorporate all of the dynamical and time-odd terms present in the nuclear energy density functional. This approach is employed in the density constrained frozen Hartree-Fock (DCFHF) and in the density constrained time-dependent Hartree-Fock (DC- TDHF) microscopic methods. Results: The PKE spatial distribution shows that a repulsion occurs in the neck between the nuclei when they first touch. Inside the barrier, neutrons can contribute significantly more to the Pauli repulsion in neutron-rich systems. Dynamical effects tend to lower the Pauli repulsion near the barrier. Proton and neutron dynamical contributions to the PKE significantly differ inside the barrier for asymmetric collisions, which is interpreted as an effect of multinucleon transfer. Conclusions: The PKE is shown to make a significant contribution to nuclear interaction potentials. Protons and neutrons can play very different roles in both the bare potential and in the dynamical rearrangement. Further microscopic studies are required to better understand the role of transfer and to investigate the effect of pairing and deformation.
Radius and equation of state constraints from massive neutron stars and GW190814
Motivated by the unknown nature of the 2.50–2.67M⊙ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. Here we demonstrate that a modification of the highly uncertain suprasaturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of 2.6M⊙ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of ≈ 2.0M⊙ neutron stars (but not necessarily the radii of 1.4 M⊙ neutron stars), which further motivates future NICER radius measurements of PSR J1614 – 2230 and PSR J0740+6620.
La 2 O 3 Mn 2 Se 2 : A correlated insulating layered d-wave altermagnet
Altermagnets represent a new class of magnetic phases without net magnetization, invariant under a combination of rotation and time reversal. Unlike conventional collinear antiferromagnets (AFM), altermagnets could lead to new correlated states and important material properties deriving from their nonrelativistic spin-split band structure. Indeed, they serve as the magnetic analogue of unconventional superconductors and can yield spin-polarized electrical currents in the absence of external magnetic fields, making them promising candidates for next-generation spintronics. Here, we report altermagnetism in the correlated insulator, magnetically ordered tetragonal oxychalcogenide, La 2 O 3 Mn 2 Se 2 . Symmetry analysis reveals a 𝑑 𝑥 2 −𝑦 2 -wave-like spin-momentum locking arising from the Mn 2 O Lieb lattice, supported by density functional theory (DFT) calculations. Magnetic measurements confirm the AFM transition below ∼166K while neutron pair distribution function analysis reveals a 2D short-range magnetic order that persists above the Néel temperature. Single crystals are grown and characterized using x-ray diffraction, optical and electron microscopy, and micro-Raman spectroscopy to confirm the crystal structure, stoichiometry, and uniformity. Furthermore, our findings establish La 2 O 3 Mn 2 Se 2 as a model altermagnetic system realized on a Lieb lattice.
Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water
The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.
Analysis of the Surface Morphology and Chemical Composition of Zr-Nb3Sn Alloys with Zr Different Concentrations
The inclusion of zirconium (Zr) in niobium-tin (Nb3Sn) significantly enhances the performance of Nb3Sn radiofrequency cavities in high magnetic fields. This research project is dedicated to characterizing the surface properties and chemical composition of Zr-doped Nb3Sn. Two different concentrations of Zr were used for doping: approximately 0.5% and 24%. Various spectroscopy techniques were employed to analyze how the surface morphology and chemical composition of Nb3Sn change with increasing Zr content. The findings of this study indicate that the grain size decreases as the Zr concentration increases. In samples with 24% Zr, nanometer-sized particles, likely oxides, were observed. X-ray photoelectron spectroscopy (XPS) data revealed that the thickness of niobium oxides (NbOx) decreases with increasing Zr, while the thickness of tin oxides (SnOx) increases. Additionally, grain size distribution analysis showed that the average grain size is around 5300 nm , with a grain area density of about 165 grains/μm .
Statistical Characterization of Hot Jupiter Atmospheres Using Spitzer’s Secondary Eclipses
We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6 and 4.5 μm using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for James Webb Space Telescope such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit (ecos ω=0.00614±0.00064), and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan & Agol (2011). Planets having maximum day-side temperatures exceeding ~2200 K are consistent with having zero albedo and a distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 ± 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800–2500 K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profiles of real planetary atmospheres as compared to models.