SEARCH · Engineering Papers
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.
In this paper we present results from experiments with a system consisting of two immiscible fluid layers in rectangular and annular geometries, driven by a vertical temperature gradient. Time-dependent variations in the type of coupling observed between the two layers are described and characterized.
Lunar laser ranging (LLR) data and Apollo seismic data analyses, revealed independent evidence for the presence of a fluid lunar core. However, the size of the lunar fluid core remained uncertain by _55 km (encompassing two contrasting 2011 Apollo seismic data analyses). Here we show that a new description of the lunar interior's dynamical model provides a determination of the radius and geometry of the lunar core-mantle boundary (CMB) from the LLR observations. We compare the present day lunar core oblateness obtained from LLR analysis with the expected hydrostatic model values, over a range of previously expected CMB radii. The findings suggest a core oblateness (f(c) = (2.2 ± 0.6) x 10(exp 14)) that satisfies the assumption of hydrostatic equilibrium over a tight range of lunar CMB radii (R(CMB) = 381 ± 12 km). Our estimates of a presently-relaxed lunar CMB translates to a core mass fraction in the range of 1.59 - 1.77% with a present-day Free Core Nutation (FCN) within (367±100) years.
Previous studies found that it is easier for observers to spatially update displays during imagined self-rotation versus array rotation. The present study examined whether either the physics of gravity or the geometric relationship between the viewer and array guided this self-rotation advantage. Experiments 1-3 preserved a real or imagined orthogonal relationship between the viewer and the array, requiring a rotation in the observer's transverse plane. Despite imagined self-rotations that defied gravity, a viewer advantage remained. Without this orthogonal relationship (Experiment 4), the viewer advantage was lost. We suggest that efficient transformation of the egocentric reference frame relies on the representation of body-environment relations that allow rotation around the observer's principal axis. This efficiency persists across different and conflicting physical and imagined postures.
The exact contraction of a generic two-dimensional (2D) tensor network state (TNS) is known to be exponentially hard, making simulation of 2D systems difficult. The recently introduced class of isometric TNS (isoTNS) represents a subset of TNS that allows for efficient simulation of such systems on finite square lattices. The isoTNS ansatz requires the identification of an “orthogonality column” of tensors, within which one-dimensional matrix product state (MPS) methods can be used for calculation of observables and optimization of tensors. Here we extend isoTNS to infinitely long strip geometries and introduce an infinite version of the Moses Move algorithm for moving the orthogonality column around the network. Using this algorithm, we iteratively transform an infinite MPS representation of a 2D quantum state into a strip isoTNS and investigate the entanglement properties of the resulting state. In addition, we demonstrate that the local observables can be evaluated efficiently. Lastly, we introduce an infinite time-evolving block decimation algorithm (iTEBD 2 ) and use it to approximate the ground state of the 2D transverse field Ising model on lattices of infinite strip geometry.
Long-term measurements of the hard X-ray spectrum from 3 keV to 8 MeV of the black-hole candidate Cygnus X-1 in its low state are reported. Observations were made from October 26 to November 18, 1977 with the A2 (Cosmic X-ray) and A4 (Hard X-ray and Low-Energy Gamma-Ray) experiments on board HEAO 1 in the spacecraft's scanning mode. The measured spectrum below 200 keV is found to agree well with previous spectra which have been fit by a model of the Compton scattering of optical or UV photons in a very hot plasma of electron temperature 32.4 keV and optical depth 3.9 or 1.6 for spherical or disk geometry, respectively. At energies above 300 keV, however, flux excess is observed which may be accounted for by a distribution of electron temperatures from 15 to about 100 keV.
Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.
We explore the connections between various coordinate systems associated with observersmoving inwardly along radial geodesics in the Schwarzschild geometry. Painleve-Gullstrand (PG) time is adapted to freely falling observers dropped from rest from infinity; Lake-Martel-Poisson (LMP) time coordinates are adapted to observers who start at infinity with non-zero initial inward velocity; Gautreau-Hoffmann time coordinates are adapted to observers dropped from rest from a finite distance from the black hole horizon.We construct from these an LMP family and a proper-time family of time coordinates, the intersection of which is PG time. We demonstrate that these coordinate families are distinct, but related, one-parameter generalizations of PG time, and show linkage to Lemaître coordinates as well.
Bidirectional reflectance measurements are the only type of reflectance data available to the remote observer. For compositional interpretations, data are desired not only for identification of possible mineral components but also for modal abundance. The latter requires detailed information about the strength of absorption features. Using a new laboratory facility, the RELAB, laboratory data in the near infrared are presented that document effects of particle size, mineral mixtures, and viewing geometry for selected materials with well-developed absorption bands. The commonly observed increase in reflectance with decrease in particle size is also observed for absorption bands as well as a related decrease in absorption strength. For small particles in parts of the spectrum of maximum reflectance, however, a minor decrease in reflectance with a decrese in particle size is sometimes observed. Small particles dominate the observed characteristics of particulate surfaces, which contain a range of particle sizes. The mean optical path length (transmission through particles) of reflected radiation measured for a variety of particle sizes has an apparent upper limit of about 2 mm for particles of less than 250 microns. The typical number of particles involved in the optical path is less than 50.
Precision measurements involve the accurate determination of parameters through repeated measurements of identically prepared experimental setups. For many parameters there is a 'natural' choice for the quantum observable which is expected to give optimal information; and from this observable one can construct an Heinsenberg uncertainty principle (HUP) bound on the precision attainable for the parameter. However, the classical statistics of multiple sampling directly gives us tools to construct bounds for the precision available for the parameters of interest (even when no obvious natural quantum observable exists, such as for phase, or time); it is found that these direct bounds are more restrictive than those of the HUP. The implication is that the natural quantum observables typically do not encode the optimal information (even for observables such as position, and momentum); we show how this can be understood simply in terms of the Hilbert space geometry. Another striking feature of these bounds to parameter uncertainty is that for a large enough number of repetitions of the measurements all V quantum states are 'minimum uncertainty' states - not just Gaussian wave-packets. Thus, these bounds tell us what precision is achievable as well as merely what is allowed.
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of approx 240 km/s. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of approx 750 +/- 50 km/s, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
Emission gradient curves for extreme-ultraviolet resonance lines of lithiumlike ions have been constructed from spectroheliograms of a northern polar coronal hole observed on August 14, 1973, with the Harvard experiment on Skylab. An emission-measure analysis indicates both reduced density and coronal temperature in the coronal hole. The boundary geometry of the coronal hole is determined, and a temperature-density model that is consistent with the observed intensities is constructed. The model gives a conductive flux of 60,000 erg/sq cm per sec at 1.03 solar radii. The boundary geometry and density distribution are combined with typical solar-wind parameters at the earth to determine an outflow velocity of 15 km/s at 1.08 solar radii.
In this work, we employ thermoreflectance thermal imaging to directly measure the steady-state two-dimensional (2D) temperature field generated by nanostructured heat sources deposited on silicon substrate with different geometrical configurations and characteristic sizes down to 400nm. The analysis of the results using Fourier’s law not only breaks down as size scales down, but it also fails to capture the impact of the geometry of the heat source. The substrate effective Fourier thermal conductivities fitted to wire-shaped and circular-shaped structures with identical characteristic lengths are found to display up to 40% mismatch. Remarkably, a hydrodynamic heat transport model reproduces the observed temperature fields for all device sizes and shapes using just intrinsic Si parameters, i.e., a geometry and size-independent thermal conductivity and nonlocal length scale. The hydrodynamic model provides insight into the observed thermal response and of the contradictory Fourier predictions. We discuss the substantial Silicon hydrodynamic behavior at room temperature and contrast it to InGaAs, which shows less hydrodynamic effects due to dominant phonon-impurity scattering.
Asteroid discoveries are essential for planetary-defense efforts aiming to prevent impacts with Earth, including the more frequent megaton explosions from decameter impactors. While large asteroids (≥100 km) have remained in the main belt since their formation, small asteroids are commonly transported to the near-Earth object (NEO) population. However, due to the lack of direct observational constraints, their size-frequency distribution — which informs our understanding of the NEOs and the delivery of meteorite samples to Earth —varies significantly among models. Here, we report 138 detections of the smallest asteroids (⪆10 m) ever observed in the main belt, which were enabled by JWST’s infrared capabilities covering the asteroids’ emission peaks and synthetic tracking techniques. Despite small orbital arcs, we constrain the objects’ distances and phase angles using known asteroids as proxies, allowing us to derive sizes via radiometric techniques. Their size-frequency distribution exhibits a break at∼100 m (debiased cumulative slopes of q=−2.66±0.60 and−0.97±0.14 for diameters smaller and larger than∼100 m, respectively), suggestive of a population driven by collisional cascade. These asteroids were sampled from multiple asteroid families —most likely Nysa, Polana and Massalia — according to the geometry of pointings considered here. Through additional long-stare infrared observations, JWST is poised to serendipitously detect thousands of decameter-scale asteroids across the sky, probing individual asteroid families and the source regions of meteorites “in-situ”.
The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.
Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.
This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.
Neutral and charged excitons (trions) in atomically thin materials offer important capabilities for photonics, from ultrafast photodetectors to highly efficient light-emitting diodes and lasers. Recent studies of van der Waals (vdW) heterostructures comprised of dissimilar monolayer materials have uncovered a wealth of optical phenomena that are predominantly governed by interlayer interactions. Here, we examine the optical properties in NbSe 2 -MoSe 2 vdW hetero- structures, which provide an important model system to study metal–semiconductor interfaces, a common element in optoelectronics. Through low-temperature photoluminescence (PL) microscopy, we discover a sharp emission feature, L1, that is localized at the NbSe 2 -capped regions of MoSe 2 . L1 is observed at energies below the commonly studied MoSe 2 excitons and trions and exhibits temperature- and power-dependent PL consistent with exciton localization in a confining potential. This PL feature is robust, observed in a variety of samples fabricated with different stacking geometries and cleaning procedures. Using first-principles calculations, we reveal that the confinement potential required for exciton localization naturally arises from the in-plane band bending due to the changes in the electron affinity between pristine MoSe 2 and NbSe 2 -MoSe 2 heterostructure. Here, we discuss the implications of our studies for atomically thin optoelectronics devices with atomically sharp interfaces and tunable electronic structures.
Using ab initio density functional theory, we study the electronic and magnetic properties of the van der Waals chain material OsCl 4 . In the nonmagnetic state, a strongly anisotropic band structure was observed, in agreement with its anticipated one-dimensional crystal geometry. Based on Wannier functions, we found that the four electrons of the 5d Os atom form a low-spin S = 1 state, with a large crystal field between the d xz/yz and d xy orbitals, corresponding to a strong Jahn–Teller distortion (Q 3 < 0). As a consequence, the magnetic properties are mainly contributed by the d xz/yz states. Furthermore, when a Mott gap develops after the introduction of the Hubbard U and Hund coupling J, we found that the staggered spin order is the most likely magnetic state, namely, spins arranged as (↑-↓-↑-↓) with π wavevector along the chain. In addition, the energy differences between various spin states are small, suggesting a weak magnetic exchange coupling along the chain. Overall, our results provide guidance to experimentalists and theorists working on quasi-one-dimensional osmium halides chain materials.