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.
Banner Highlight for Structure-Dependent Doppler Broadening Using a Generalized Thermal Scattering Law [DOI 10.3390/jne2020013].
Polar vortex domains have recently become an emergent research field due to the abundant physical phenomena and potential applications in high-density memories. Here, we explore the mechanisms of creating polar vortex domains in the BiFeO3 (BFO) membranes subjected to different boundary conditions using phase-field simulations. A major difference is that the vortex in membrane can be stabilized even under short-circuit electrical boundary conditions compared to vortex in other systems, such as thin film or superlattice. We found that (a) the formation of polar vortex domains at the membrane interior under bending is mainly driven by the reduction of elastic energy under short-circuit boundary condition, and the vortex chirality (namely, clockwise and counterclockwise) could be identified by n-shape and u-shape bending; (b) in the unbent open-circuit BFO membrane case, exotic trapezoid-shaped vortex nanodomains form at the terminations of 109 degrees domain walls (DWs) and partially charged 71 degrees DWs, which is driven by the local depolarization field and the interplay among electrostatic, elastic, and gradient and Landau energies. We also examine Kittel's law by establishing the dependence of vortex periods on the membrane thickness. Furthermore, these results give further understanding of the effect of boundary conditions on the formation of polar vortex domains, guiding experimental designs of vortex-based high-density memories.
Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.
ABSTRACT The Milky Way halo is one of the few galactic haloes that provides a unique insight into galaxy formation by resolved stellar populations. Here, we present a catalogue of ∼47 million halo stars selected independent of parallax and line-of-sight velocities, using a combination of Gaia DR3 proper motion and photometry by means of their reduced proper motion. We select high tangential velocity (halo) main sequence stars and fit distances to them using their simple colour-absolute-magnitude relation. This sample reaches out to ∼21 kpc with a median distance of 6.6 kpc thereby probing much further out than would be possible using reliable Gaia parallaxes. The typical uncertainty in their distances is $0.57_{-0.26}^{+0.56}$ kpc. Using the colour range 0.45 < (G0 − GRP, 0) < 0.715, where the main sequence is narrower, gives an even better accuracy down to $0.39_{-0.12}^{+0.18}$ kpc in distance. The median velocity uncertainty for stars within this colour range is 15.5 km s−1. The distribution of these sources in the sky, together with their tangential component velocities, are very well-suited to study retrograde substructures. We explore the selection of two complex retrograde streams: GD-1 and Jhelum. For these streams, we resolve the gaps, wiggles and density breaks reported in the literature more clearly. We also illustrate the effect of the kinematic selection bias towards high proper motion stars and incompleteness at larger distances due to Gaia’s scanning law. These examples showcase how the full RPM catalogue made available here can help us paint a more detailed picture of the build-up of the Milky Way halo.
Quantum spin liquids are highly entangled ground states of insulating spin systems, in which magnetic ordering is prevented down to the lowest temperatures due to quantum fluctuations. One of the most extraordinary characteristics of quantum spin liquid phases is their ability to support fractionalized, low-energy quasiparticles known as spinons, which carry spin-1/2 but bear no charge. Relaxometry based on color centers in crystalline materials—of which nitrogen-vacancy (NV) centers in diamond are a well-explored example—provides an exciting new platform to probe the spin spectral functions of magnetic materials with both energy and momentum resolution and to search for signatures of these elusive, fractionalized excitations. In this work, we theoretically investigate the color-center relaxometry of two archetypal quantum spin liquids: the two-dimensional U(1) quantum spin liquid with a spinon Fermi surface and the spin-1/2 antiferromagnetic spin chain. The former is characterized by a metallic, spin-split ground state of mobile, interacting spinons, which closely resembles a spin-polarized Fermi liquid ground state but with neutral quasiparticles. We show that the observation of the Stoner continuum and the collective spin wave mode in the spin spectral function would provide a strong evidence for the existence of spinons and fractionalization. In one dimension, mobile spinons form a Luttinger liquid ground state. We show that the spin spectral function exhibits strong features representing the collective density and spin-wave modes, which are broadened in an algebraic fashion with an exponent characterized by the Luttinger parameter. The possibilities of measuring these collective modes and detecting the power-law decay of the spectral weight using NV relaxometry are discussed. We also examine how the transition rates are modified by marginally irrelevant operators in the Heisenberg limit. Published by the American Physical Society 2024
While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.
The commercially available Outotec HSC 9.9.0{sup R} provides excellent predictions for the solubilities of many species for Hanford waste tank concentrates. The HSC solubility predictions with a Harvie activity model agreed well with five boildowns that showed solubility limits for NaNO{sub 3}, NaNO{sub 2}, Na{sub 2}CO{sub 3}, NaF, Na{sub 2}SO{sub 4}, and Na{sub 3}PO{sub 4}. However, the HSC did not predict density and therefore water content very well and therefore we used the measured densities and derived water content from the difference between measured density and assays. Furthermore, we adjusted the HSC enthalpy for the minor species Na{sub 3}PO{sub 4}*NaF double salt in order to fit the five boildowns as well. Previous work has shown that gibbsite becomes and remains supersaturated in all of these boildowns and so there was no alumina precipitation despite supersaturation. Furthermore, an equilibrium model like HSC always shows complete oxidation of any TOC in the presence of nitrate and nitrite. We therefore suppressed the HSC oxidation of acetate and oxalate in order to fit measured acetate and oxalate solubilities as TOC. Space in Hanford waste tanks is very limited and therefore evaporator campaigns continue to recover needed tank space from dilute tank liquids. Prior to each evaporator campaign, a boildown of a sample of tank liquid for each campaign provides a basis for planning by measuring density, water activity, soluble species, and solids precipitation as a function of concentration at 18 deg. C. These boildowns then provide validations and calibrations for the HSC solubility prediction for Hanford wastes. Previous work has shown the solubilities of minor electrolytes in these five boildowns follow rather simple power law solubility products. Hanford's 222-S laboratory has reported each of five boildowns for concentrated liquids from five waste tanks and tank blends. All boildowns eventually showed the precipitation of the dominant electrolyte, NaNO{sub 3}, and overall solution activity largely drives all solubility. (authors)
Graphite is a classic neutronic material that has been used as both a reactor reflector and moderator in various nuclear reactor systems. The ability to accurately predict the slowing down and thermalization of neutrons in graphite can have significant implications on the safety and operation of such reactor systems. In reactors, the neutron thermalization process is quantified using the thermal scattering law (TSL) and related cross sections for a given moderator. An ideal approach to assess the validity of TSL data is using benchmark measurement based on the pulsed Slowing-Down-Time technique and its comparison with the appropriate graphite nuclear library. In this work, experimental measurement and computational Monte Carlo simulations were performed to benchmark the slowing down characteristics and thermalization of neutrons in nuclear (reactor-grade) graphite. Given the density of graphite, various graphite libraries were selected for the benchmark analysis.
Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.
The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.
A newly uncovered class of plasmons in the strongly excited limit opens access to unprecedented Petavolts per meter electromagnetic fields with wide-ranging, transformative impact. Unlike conventional plasmons, such plasmons are constituted by non-perturbative, large-amplitude oscillations of the ultradense, delocalized free electron Fermi gas inherent in conductive media. Here structured semiconductors doped to have an appropriate conduction electron density are introduced to tune the properties of the Fermi gas for matched excitation of large-amplitude plasmons using readily available electron beams which enables immediate experimental validation. Specifically, an electrostatic, surface “crunch-in” plasmon is collisionlessly excited by the beam launched inside a tube. Strong excitation due to matching results in relativistic oscillations of the electron gas and unravels unique phenomena. Relativistically induced ballistic electron transport comes about due to relativistic multifold increase in the mean free path and also leads to unconventional heat deposition beyond Ohm’s law. This explains the absence of observed damage or solid-plasma formation in past experiments on conductive samples interacting with electron bunches shorter than 10-13 seconds. Furthermore, relativistic momentum leads to copious tunneling of electron gas across the surface, which then crunches inside the tube. Relativistic effects along with large, localized electron density variations underlying these modes necessitate kinetic approach to theoretical and computational modeling. Kinetic model presented here demonstrates experimental viability of observing tens of gigavolts per meter plasmonic fields excited by matching readily available electron beams to plasmons in semiconductors with 1018cm-3 free electron density, and paves the way for Petavolts per meter plasmonics.
Hydrogen dense polymers, specifically polyethylene, polystyrene, and Plexiglas, have served as neutron moderator and reflector materials in hundreds of separate critical benchmark experiments because of their low cost and abundance of hydrogen. In order to accurately model and simulate these critical benchmarks, the thermal scattering law (TSL) evaluation that governs how neutrons will thermalize must be well understood and rigorously validated. To support this validation, researchers at Rensselaer Polytechnic Institute performed total neutron cross section measurements for high-density polyethylene & polystyrene over the energy range 0.0005–20 eV and for Plexiglas G & Plexiglas G-UVT over the energy range 0.0005–3 eV. Comparisons were made between the measured cross section and that predicted by the ENDF/B-VIII.0 and Oak Ridge National Laboratory/European Spallation Source/Rensselaer Polytechnic Institute TSL evaluations for polyethylene, Plexiglas, and polystyrene, as available. These experiments represent the first total neutron cross section measurements for polystyrene.
We report a recurrent neural network (RNN) based model is developed as a surrogate to predict nonlinear plastic response under multiaxial loading. The RNN-based model is trained and tested on stress versus strain curves generated using a numerical solution based on the classical radial return method. Besides simply learning the basic constitutive relationship, a novel approach is taken to enforce certain physical conditions. Specifically, regularization is employed to maintain non-negative plastic power density throughout the loading history thereby ensuring monotonically increasing plastic work and thermodynamic consistency. Enforcing physics in this manner permits coupling of the data-driven RNN approach with physics-based knowledge and laws. This has the effect of reducing the necessary amount of data and ensuring known physical laws are not violated. Since, once trained, the model need not perform the expensive task of solving nonlinear equations, its efficiency is orders of magnitude greater than its numerical counterpart. The RNN-based model has been trained on varied sets of data and the accuracy on test datasets validated. The developed model is general and robust and has widespread application such as in the simulation of metal forming, large scale plasticity, and part life prediction.
Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.
Motivated by recent cold-atom experiments, here we study the relaxation of spin helices in quantum XXZ spin chains. The experimentally observed relaxation of spin helices follows scaling laws that are qualitatively different from linear-response transport. We construct a theory of the relaxation of such spin helices, combining generalized hydrodynamics with diffusive corrections and a generalized form of the local density approximation. Although helices are far from local equilibrium, our hydrodynamic approach reproduces the experimentally observed relaxational dynamics and also predicts the late-time relaxation, which is outside the experimentally accessible timescales. In particular, our theory explains the existence of temporal regimes with apparent anomalous diffusion, as well as the asymmetry between positive and negative anisotropy regimes at short and intermediate times.
We explore the evolution of the X-ray luminosity function of quasars and the intrinsic correlation between the X-ray and 2500 Å ultraviolet luminosities, utilizing techniques verified in previous works and a sample of over 4000 quasars detected with Chandra and XMM-Newton in the range 0 < z < 5. We find that quasars have undergone significantly less evolution with redshift in their total X-ray luminosity than in other wave bands. We then determine that the best-fit intrinsic power-law correlation between the X-ray and ultraviolet luminosities, of the form ${L}_{{\rm{X}}}^{{\prime} }\propto {({L}_{\mathrm{UV}}^{{\prime} })}^{\gamma }$, is γ = 0.28 ± 0.03, and we derive the luminosity function and density evolution in the X-ray band. We discuss the implications of these results for models of quasar systems.
Bi-component protection-systems need to withstand a variety of extreme loading conditions such as high-velocity projectile impact to shock and blast mitigation. The large number of parameters at each stage of bi-component protection system design process results in a near-infinite number of possible variables for optimization. The aim of this project is to identify the scaling laws linking material and geometric scales of the problem to the final impact performance, further including the effects of wave propagation/interaction and inter-component bonding. Specifically, the impact dynamics for different optimal projectile/target density ratios and relative component mass fractions are quantified for bonded ceramic/titanium alloy systems. The impact behavior such as wave mechanics and damage propagation within the target system is characterized and quantified via digital image correlation or via x-ray phase-contrast imaging.
In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity $\overline{k_{L}}$, where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts $\overline{k_{L}}$ = α o Sc -1/2 ($v\bar{ϵ}$) 1/4 , where Sc is the molecular Schmidt number, $v$ is the water kinematic viscosity, and $\bar{ϵ}$ is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While α o = 0.39 - 0.46 has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous ϵ around $\bar{ϵ}$, a mechanism that was not previously considered. As the coefficient of variation (CV e ) in ϵ increases, α o must be adjusted by a multiplier (1 = CV e 2 ) -3/32 that was derived from a log‐normal model for the probability density function of ϵ. Reported variations in α o with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in ϵ. Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of ϵ. That α o varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade.