Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parabolic”

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 217 records · Page 12

Thermodynamic properties and enhancement of diamagnetism in nitrogen doped lutetium hydride synthesized at high pressure

Nitrogen doped lutetium hydride has drawn global attention in the pursuit of room-temperature superconductivity near ambient pressure and temperature. However, variable synthesis techniques and uncertainty surrounding nitrogen concentration have contributed to extensive debate within the scientific community about this material and its properties. We used a solid-state approach to synthesize nitrogen doped lutetium hydride at high pressure and temperature (HPT) and analyzed the residual starting materials to determine its nitrogen content. High temperature oxide melt solution calorimetry determined the formation enthalpy of LuH 1.96 N 0.02 (LHN) from LuH 2 and LuN to be −28.4 ± 11.4 kJ/mol. Magnetic measurements indicated diamagnetism which increased with nitrogen content. Ambient pressure conductivity measurements observed metallic behavior from 5 to 350 K, and the constant and parabolic magnetoresistance changed with increasing temperature. High pressure conductivity measurements revealed that LHN does not exhibit superconductivity up to 26.6 GPa. We compressed LHN in a diamond anvil cell to 13.7 GPa and measured the Raman signal at each step, with no evidence of any phase transition. Despite the absence of superconductivity, a color change from blue to purple to red was observed with increasing pressure. Thus, our findings confirm the thermodynamic stability of LHN, do not support superconductivity, and provide insights into the origins of its diamagnetism.

Science & Technology - Other Topics↗

Exponential Time Differencing Schemes for Fuel Depletion and Transport in Molten Salt Reactors: Theory and Implementation

A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-temperature oxidation of transient-liquid phase bonded Ni-based alloys in 1 bar and 250 bar carbon dioxide

Compact microchannel heat exchanger designs for supercritical CO 2 power cycles are valued for equipment size reduction and increased heat transfer. They can be fabricated via transient-liquid phase (TLP) bonding, which utilizes a lower melting point interlayer for bonding facilitation of etched thin plates. Haynes 230 and 282 TLP bonded stacks were exposed to high-purity CO 2 at 1 bar/700 °C to 4000 h and 250 bar/720 °C (supercritical) to 1500 h. All samples obeyed parabolic kinetics after an oxide buildup stage of 500-1000 h. Mass gains were similar for the base and TLP bonded alloys, and oxide thicknesses and scales did not significantly differ around the TLP bond regions. Furthermore, there were some slight oxide compositional differences for the H282 TLP samples at 250 bar. There was a slight pressure effect on mass gain, as samples exposed at the higher pressure experienced higher initial mass gains, but then slower continued oxide growth.

36 MATERIALS SCIENCE↗

Sulfidation kinetics of titanium and Ti-6Al-4V with elemental sulfur

The sulfidation kinetics of titanium and Ti-6Al-4V (Ti64) using elemental sulfur were studied by measuring the amount of metal consumed instead of the change in mass or scale thickness. Here this metal consumption approach combines the parabolic oxidation law, Arrhenius equation, and Pilling-Bedworth ratio to determine the apparent activation energy and the pre-exponential factor of consumption of Ti64 during sulfidation. The estimations of apparent activation energy and pre-exponential factor were 110.4 kJ·mol -1 and 1.3×10 -9 m 2 s -1 for the titanium-sulfur, and 116.2 kJ·mol -1 and 1.2×10 -9 m2s -1 for the Ti64-sulfur system. Finally, a material removal predictor was developed to predict the amount of Ti64 consumed for a given sulfidation temperature and time.

36 MATERIALS SCIENCE↗

Restoring the discontinuous heat equation source using sparse boundary data and dynamic sensors

Abstract This study focuses on addressing the inverse source problem associated with the parabolic equation. We rely on sparse boundary flux data as our measurements, which are acquired from a restricted section of the boundary. While it has been established that utilizing sparse boundary flux data can enable source recovery, the presence of a limited number of observation sensors poses a challenge for accurately tracing the inverse quantity of interest. To overcome this limitation, we introduce a sampling algorithm grounded in Langevin dynamics that incorporates dynamic sensors to capture the flux information. Furthermore, we propose and discuss two distinct dynamic sensor migration strategies. Remarkably, our findings demonstrate that even with only two observation sensors at our disposal, it remains feasible to successfully reconstruct the high-dimensional unknown parameters.

Mathematics↗

The theory of transport in helical spin-structure crystals

Abstract We study helical structures in spin-spiral single crystals. In the continuum approach for the helicity potential energy the simple electronic band splits into two non-parabolic bands. For low exchange integrals, the lower band is described by a surface with a saddle shape in the direction of the helicity axis. Using the Boltzmann equation with the relaxation due to acoustic phonons, we discover the dependence of the current on the angle between the electric field and helicity axis leading to the both parallel and perpendicular to the electric field components in the electroconductivity. The latter can be interpreted as a planar Hall effect. In addition, we find that the transition rates depend on an electron spin allowing the transition between the bands. The electric conductivities exhibit nonlinear behaviors with respect to chemical potential µ . We explain this effect as the interference of the band anisotropy, spin conservation, and interband transitions. The proposed theory with the spherical model in the effective mass approximation for conduction electrons can elucidate nonlinear dependencies that can be identified in experiments. We find the excellent agreement between the theoretical and experimental data for parallel resistivity depending on temperature at the phase transition from helical to ferromagnetic state in a M n P single crystal. In addition, we predict that the perpendicular resistivity abruptly drops to zero in the ferromagnetic phase.

Physics↗

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗

Demonstration of TNSA proton radiography on the National Ignition Facility Advanced Radiographic Capability (NIF-ARC) Laser

We report that proton radiography using short-pulse laser drivers is an important tool in high-energy density (HED) science for dynamically diagnosing key characteristics in plasma interactions. Here we detail the first demonstration of target-normal sheath acceleration (TNSA)-based proton radiography the NIF-ARC laser system aided by the use of compound parabolic concentrators (CPCs). The multi-kJ energies available at the NIF-ARC laser allows for a high-brightness proton source for radiography and thus enabling a wide range of applications in HED science. In this demonstration, proton radiography of a physics package was performed and this work details the spectral properties of the TNSA proton probe as well as description of the resulting radiography quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing extreme laser intensities by ponderomotive acceleration of protons from rarified gas

Abstract A new method to diagnose extreme laser intensities through measurement of angular and spectral distributions of protons directly accelerated by the laser focused into a rarefied gas is proposed. We simulated a laser pulse focused by an off-axis parabolic mirror by Stratton–Chu integrals, that enables description of laser pulse with different spatial-temporal profiles focusing in a focal spot down to the diffraction limit, that makes our theoretical predictions be a basis for experimental realization. The relationship between characteristics of the proton distributions and parameters of the laser pulse have been analyzed. The analytical and numerical results obtained justify the new method of laser diagnostics. The proposed scheme should be valuable for the commissioning of new extreme intensity laser facilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-significance detection of correlation between the unresolved gamma-ray background and the large-scale cosmic structure

Our understanding of the γ-ray sky has improved dramatically in the past decade, however, the unresolved γ-ray background (UGRB) still has a potential wealth of information about the faintest γ-ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ-ray background. In this study, we analyze the angular correlation between the γ-ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi-Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10 14 M ⊙ ) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ-ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ-ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗

Dispersive interactions between standard and Dirac materials and the role of dimensionality

The van der Waals (vdW) interaction plays a prominent role between neutral objects at separations where short ranged chemical forces are negligible. This type of dispersive coupling is determined by the interplay between geometry and response properties of the materials making up the objects. Here, we investigate the vdW interaction between 1D, 2D, and 3D standard and Dirac materials within the Random Phase Approximation, which takes into account collective excitations originating from the electronic Coulomb potential. A comprehensive understanding of characteristic functionalities and scaling laws are obtained for systems with parabolic energy dispersion (standard materials) and crossing linear bands (Dirac materials). By comparing the quantum mechanical and thermal limits the onset of thermal fluctuations in the vdW interaction is discussed showing that thermal effects are significantly pronounced at smaller scales in reduced dimensions.

36 MATERIALS SCIENCE↗

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗

Giant Planet Scatterings and Collisions: Hydrodynamics, Merger-Ejection Branching Ratio, and Properties of the Remnants

Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple “sticky-sphere” prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and N-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with MJ and 2MJ, we find that quick mergers occur when the impact parameter rp (the pericenter separation between the planets) is less than 2RJ, and the merger conserved at least 97% of the initial mass. Strong tidal effects can affect the “binary-planet” orbit when rp is between 2RJ and 4RJ. We quantify these effects using a set of fitting formulae that can be implemented in N-body codes. We run a suite of N-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular and nearly coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the branching ratio change. We also find that the merger products have broad distributions of spin magnitudes and obliquities.

79 ASTRONOMY AND ASTROPHYSICS↗

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

79 ASTRONOMY AND ASTROPHYSICS↗

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗

Efficacy of the Method of Four Coefficients to Determine Charge-Carrier Scattering

The investigation of the electronic properties of semiconductors from transport measurements (i.e., resistivity, Hall, and Seebeck coefficient measurements) is challenging due to the averaging of charge-carrier dynamics inherent in such measurements. Here, we investigate the incorporation of a fourth measurement of electronic transport, the Nernst coefficient, into the analysis, termed the method of four coefficients. This approach yields the Fermi level, effective mass, scattering exponent, and relaxation time prefactor. We begin with a review of the underlying mathematics and investigate the mapping between the four-dimensional material property and transport coefficient spaces. We then investigate how the traditional single-parabolic band method yields a single, potentially incorrect point on the solution subspace. This uncertainty can be resolved through Nernst coefficient measurements and we map the span of the ensuing subspace. We conclude with an investigation of how sensitive the analysis of transport coefficients is to experimental error for different sample types.

charge carriers↗