Search NASA⌕ Search

SEARCH · Search NASA

Results for “steep”

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 163 records · Page 9

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗

A multifrequency study of sub-parsec jets with the Event Horizon Telescope

Context. The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN.Aims. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion.Methods. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132–1654, J0006–0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2–86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole.Results. Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factor δ ∝ r≤0.5. This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radius B ∝ r−3, hinting at jet acceleration or efficient magnetic energy dissipation.Key words: techniques: interferometric / galaxies: active / galaxies: jets / galaxies: nuclei / quasars: general / quasars: supermassive black holes⋆⋆ Deceased.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Onset of High-Harmonic Fast Wave Propagation in Annulus Geometry

Previous experiments on NSTX have reported high-harmonic fast wave (HHFW) power losses in the scrape-off layer (SOL), with dependencies on the wave parallel wavenumber, plasma density, and magnetic field. This study presents an analysis of HHFW propagation in the SOL plasma, focusing on parasitic coupling to the waveguide mode bounded by the steep gradient pedestal and the outer wall boundary. The analytic and finite element wave models are discussed using a 2D annulus geometry, which serves as a simplified representation of the NSTX-U plasma. The model shows that a minimum radial gap width is required to support wave propagation, analogous to the higher-order modes in coaxial geometries. The impact of several factors, including gap width, wave parallel phasing, and magnetic field strength, on parasitic wave coupling is examined, highlighting the role of the geometric dimension in relation to the wave perpendicular wavelength. The findings suggest that controlling the edge/SOL density profile (gap control) is crucial for mitigating parasitic coupling when the density in front of the antenna exceeds the density cutoff.

Baek, Seung-Gyou [Massachusetts Inst. of Technolog↗

Evaluation of the Passive Cooling Potential of Thermal Mass Inherent in Medium to Large Commercial Buildings

Typical commercial construction inherently contains large surface areas of relatively thin mass located in the floor slabs. This research asks how effective this mass could be for passive cooling, relative to other mass depths, if it were to be exposed and used as thermal mass. Although much of the rich literature on thermal mass has been conducted with far greater mass depths, a more limited amount of research suggests that because the surface heat transfer rate is limited, the focus should be on the surface area of exposed mass, not its depth. This presents an opportunity for more buildings to behave more thermally massively, if the mass inherent in typical floor slabs contains utility, particularly over diurnal cycles. In all climates analyzed, it was found that there is a pronounced shoulder where energy savings from passive cooling of increasing mass depths was steep until roughly 7.5–10 cm, and beyond this point the achieved energy savings diminished rapidly. When considering the embodied energy of concrete, the incremental benefit of added mass beyond a typical topping slab of 10 cm does not justify the incremental embodied energy cost from a total energy standpoint alone.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A New Modeling Approach for a Priori Uncertainties of Laser Tracker Angle Measurements

Methods for modeling the uncertainty in laser tracker angle measurements vary within the metrology industry, leading to confusion and questionable stochastic modeling for survey network adjustments and error propagation analysis. Interpreting the published laser tracker manufacturer performance specifications to determine an a priori sigma value for weighting azimuth and zenith angle measurements can be confusing and has led to differing implementations. Furthermore, this paper proposes a unique way to model survey network a priori laser tracker angular uncertainties based on laser tracker manufacturers’ published maximum permissible error (MPE) values referenced to current standards for weighting survey network measurements. This paper’s proposed model takes into account the disparate effects that pointing errors, target centering errors, and leveling errors have on azimuth and zenith angular uncertainties for measurements with steep sightings and at near ranges.

42 ENGINEERING↗

Mechanical Purcell filters for microwave quantum machines

In circuit quantum electrodynamics, measuring the state of a superconducting qubit introduces a loss channel, which can enhance spontaneous emission through the Purcell effect, thus decreasing the qubit lifetime. This decay can be mitigated by performing the measurement through a Purcell filter, which strongly suppresses signal propagation at the qubit transition frequency. If the filter is also well-matched at the readout cavity frequency, it will protect the qubit from decoherence channels without sacrificing measurement bandwidth. In this work, we propose and analyze design for a mechanical Purcell filter, which we also fabricate and characterize at room temperature. The filter is composed of an array of nanomechanical resonators in thin-film lithium niobate, connected in a ladder topology, with series and parallel resonances arranged to produce a bandpass response. Their modest footprint, steep band edges, and lack of cross talk make these filters an appealing alternative to analogous electromagnetic versions currently used in microwave quantum machines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational analysis of ion orbital loss in diverted positive- and negative-triangularity tokamaks

The mechanism of ion orbital loss in diverted negative-triangularity (NT) tokamaks is significantly different from that of positive-triangularity (PT) tokamaks. The orbital loss from the edge of NT tokamaks can be quite large, because the separatrix X-point is located on the outboard side, meaning that trapped ions crossing the separatrix are lost to the inboard divertor plate after a mirror reflection on the side of the strong magnetic field. NT tokamaks have difficulties in building up steep pressure gradients in the edge regions. However, the majority of trapped ions crossing the separatrix in PT tokamaks survive (through complete banana motions) without reaching the X-point region. For NT tokamaks, topological changes in the trapped ion trajectories can be achieved in the presence of the radial electric field, thus reducing the orbital loss, whereas for PT tokamaks, the orbital loss is relatively insensitive to the radial electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High quality epitaxial thin films and exchange bias of antiferromagnetic Dirac semimetal FeSn

FeSn is a topological semimetal (TSM) and kagome antiferromagnet (AFM) composed of alternating Fe 3 Sn kagome planes and honeycomb Sn planes. This unique structure gives rise to exotic features in the band structures such as the coexistence of Dirac cones and flatbands near the Fermi level, fully spin-polarized 2D surface Dirac fermions, and the ability to open a large gap in the Dirac cone by reorienting the Néel vector. In this paper, we report the synthesis of high-quality epitaxial (0001) FeSn films by magnetron sputtering. Using FeSn/Py heterostructures, we show a large exchange bias effect that reaches an exchange field of 220 Oe at 5 K, providing unambiguous evidence of antiferromagnetism and strong interlayer exchange coupling in our films. Field cycling studies show steep initial training effects, highlighting the complex magnetic interactions and anisotropy. Importantly, our work provides a simple, alternative means to fabricate FeSn films and heterostructures, making it easier to explore the topological physics of AFM TSMs and develop FeSn-based spintronics.

36 MATERIALS SCIENCE↗

Gunn threshold voltage characterization in GaAs devices with wedge-shaped tapering

Here, we fabricate gallium arsenide-based devices with a wedge-shaped tapering region connected to a rectangular-shaped region and measure the threshold voltage required to trigger the Gunn effect. The threshold voltage reduction is attributed to the focusing of the electric field toward the narrower end of the device and is effective when the device has a steep enough tapering. We also model the electric field profile for the tapered devices using an intuitive graphical approach and the finite element method and provide estimates for the threshold voltages of tapered devices. Finally, we compare the estimates to the measured values and provide possible reasons for the discrepancies. We believe the capability of threshold voltage reduction with the wedge-shaped tapering design could be useful in device applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electric field dynamics in an atmospheric pressure helium plasma jet impinging on a substrate

Time and spatially resolved electric field measurements by Stark polarization spectroscopy in a nanosecond pulsed atmospheric pressure helium jet operating in ambient air and impinging on a conductive indium tin oxide (ITO) coated glass slide are reported. An automatic fitting procedure of the Stark shifted spectra taking into consideration constraints regarding Stark components positions and intensities as well as molecular nitrogen emission subtraction was implemented. This allowed electric field vector components measurements both in the gas phase and at the interface when the jet impinges on the substrate and during the development of a surface ionization wave. The obtained results show an increase of the axial electric field in the jet effluent in the gas phase with a peak magnitude from 12 to 18 kVcm -1 before the ionization wave impinges on the substrate. A steep electric field enhancement to a peak value of about 24 kVcm -1 was observed when the ionization wave impinges on the surface. A peak radial electric field of about 27 kVcm -1 was measured off-axis in the surface ionization wave. These results are consistent with previously reported modelling predictions. While Stark polarization spectroscopy is limited to electric field measurement from regions with emission, we illustrate that the capability to measure near surface electric fields in helium makes it a valuable complementary technique for the electric field-induced second harmonic (E-FISH) technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temperature distribution in a laser-heated diamond anvil cell as described by finite element analysis

Finite element analysis (FEA) is a powerful tool for numerically solving partial differential equations over complex geometries and is thus useful for analyzing heat transport in laser-heated diamond anvil cell (LHDAC) experiments. Our models expand on previously published simulations by calculating the volume-averaged temperatures of both the sample and insulation/pressure media under steady-state heating to determine the thermal pressure of the hot sample. Our goal is to produce an accurate relationship between the measured surface temperature of the absorbing sample and the temperature of the transparent insulating media, which is used to determine thermal pressure but susceptible to steep temperature gradients. We find that in doing so, our FEA models of temperature within the pressure/insulation media can differ from simplified estimates of temperature gradients by more than a factor of 2. We also explore temperature-dependent and temperature-independent thermal conductivity models and find that the volume-averaged temperatures differ by up to a factor of 1.3, forcing the predicted thermal pressures determined to also differ by up to a factor of 1.5 at a temperature of 2000 K at 50 GPa for neon. Higher temperatures exacerbate this difference. We also find that unintentional asymmetric sample insertion and sample heating, which are common in LHDAC experiments, do not have a first-order effect on volume-averaged temperatures. The FEA models, available in both Python and FlexPDE, are versatile across different sample geometries, materials, and heat source laser shapes.

Farah, Frederick↗

Pedestal fluctuation measurements with charge exchange imaging at the DIII-D tokamak

A new high radial resolution 2D multichannel Charge eXchange Imaging (CXI) diagnostic is under development for deployment at DIII-D. Here, the diagnostic system will measure low-to-intermediate radial wavenumber carbon density fluctuations by observing the n = 8 – 7 (λ = 529.06 nm) C-VI emission line, resulting from charge exchange collisions between heating neutral beam atoms and the intrinsic carbon ion density. The new CXI diagnostic will provide measurements with ΔR ~ 0.4 cm to access higher k r instabilities (k r < 8 cm –1 ) predicted to arise in the steep-gradient region of the H-mode pedestal. The CXI system will feature 60 fiber bundles in a 12 × 5 arrangement, with each bundle consisting of four 1 mm fibers. A custom optical system has been designed to filter and image incoming signals onto an 8 × 8 avalanche photodiode array. Additionally, a novel electronics suite has been designed and commissioned to amplify and digitize the relatively low-intensity carbon signal at a 2 MHz bandwidth. Forward modeling results of the active C-VI emission suggest sufficient signal to noise ratios to resolve turbulent fluctuations. Prototype measurements demonstrate the ability to perform high frequency pedestal measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗