Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local trap”

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 829 records · Page 46

Average spatial distributions of energetic O(+), O(2+), O(6+), and C(6+) ions in the magnetosphere observed by AMPTE CCE

Measurements with the charge energy mass spectrometer of the AMPTE CCE spacecraft are used to determine the relative fluxes of these ions near the equatorial plane as a function of the drift shell parameter L, the magnetic activity index Kp, and the local time LT. The O(+) and O(2+) ions have radial profiles with maxima at about L = 5 and diurnal variations with a broad maximum on the dayside. The O(6+) and C(6+) ion fluxes increase with L between L = 5 and L = 7 and level off farther out. The diurnal variations of the relative O(6+) and C(6+) fluxes exhibit a pronounced minimum on the dayside. The observations can be interpreted in terms of ion convection from the tail onto quasi-trapped drift orbits and further radial diffusion onto closed drift shells. The maximum of the O(+) and O(2+) fluxes as well as the minimum of the O(6+) and C(6+) fluxes on the dayside can be explained by drift shell splitting. The inward transport is associated with ion losses, mainly by charge exchange. O(6+) and C(6+) ions are removed already at larger distances than O(+) and O(2+). As a result, the ionosphere contributes most of the oxygen ions in the magnetosphere. Important contributions from the solar wind are encountered in the outer magnetosphere.

Kremser, G.↗

Maintaining a Flat Li Surface during the Li Stripping Process via Interface Design

Electroplating has been the main focus in mitigating the dendrite growth on the Li-metal electrode; however, the stripping process is equally critical, since the nonsmooth Li surface during stripping will lead to nonuniform local current density, planting the seeds for dendrite growth. In this paper, density functional theory (DFT) and kinetic Monte Carlo (KMC) techniques were combined to investigate the vacancy evolution in Li interfaced with different solid–electrolyte interphase (SEI) materials. It was found that the lithiophilic interface, such as Li/Li 2 O, repels vacancies into the bulk Li, so Li atoms can quickly fill the Li vacancies near the Li/Li 2 O interface and maintain a smooth Li surface. In contrast, the lithiophobic interface, such as Li/LiF, traps Li vacancies toward the interface, and the accumulated Li vacancies form voids and roughen the surface. The predicted critical stripping current density, below which a smooth Li surface will be maintained, is therefore much faster at the lithiophilic interface than that at the lithiophobic interface. It was further revealed that the lithiophilicity at different SEI or coating materials can be ranked as Li/Li 2 O > Li/LiPON > Li/Li 2 CO 3 > Li/LiF based on the calculated interfacial adhesion and accumulation of electron density at the interface. Furthermore, this suggests that interface and coating design at nanoscale can be effective for maintaining a smooth Li surface during the stripping process, solving another challenge to achieving a dendrite-free Li-metal electrode in both liquid and solid electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗

Ray-tracing studies and path-integrated gains of ELF unducted whistler mode waves in the earth's magnetosphere

Gyroresonance and Landau resonance interactions between unducted low-frequency whistler waves and trapped electrons in the earth's plasmasphere have been studied. Ray paths for waves launched near the plasmapause have been traced. In agreement with recent findings by Thorne et al. (1979), waves have been found which return through the equatorial zone with field-aligned wave normal angles. However, when the growth along the ray path is calculated for such waves, assuming an electron distribution function of the form E exp -n sin exp m alpha, it is found that for all the waves considered, the local growth rate becomes negative before plasmapause reflection, limiting the total gain to small values. Most waves reach zero gain before reflection. This is the result of Landau damping at oblique propagation angles, which necessarily occurs before reflection can take place. It is concluded that the concept of cyclic ray paths does not provide an explanation for the generation of unguided plasmaspheric hiss.

Huang, C. Y.↗

Very Large Array-RATAN 600 observations of a solar active region

VLA-RATAN-600 observations of the sun around the partial solar eclipse in July 1990 are discussed. The relatively flat spectrum of the extended emission is attributed to the optically thin thermal bremsstrahlung of electrons trapped in a magnetic loop at coronal temperatures. Step-spectrum sunspot-associated emission is attributed to thermal gyroresonance radiation at different heights along the leg of a loop joining regions of positive magnetic polarity. Comparisons with predicted distributions of gyroresonance radio indicate that compact sunspot-associated sources lie at heights of 2500-17,500 km above the photosphere. Differences in the observed and predicted microwave emission suggest an inhomogeneous distribution of thermal plasma in the corona above the sunspot or the presence of localized flux tubes containing higher density plasma or nonpotential fields.

Bogod, V. M.↗

Two-Dimensional Distribution of Volatiles in the Lunar Regolith from Space Weathering Simulations

We present simulations of space weathering effects on ice deposits in regions of permanent shadow on the Moon. These Monte Carlo simulations follow the effects of space weathering processes on the distribution of the volatiles over time. The model output constrains the coherence of volatile deposits with depth, lateral separation, and time. The results suggest that ice sheets become broken and buried with time. As impacts begin to puncture an initially coherent surficial ice sheet, small areas with a deficit of ice compared to surrounding areas are formed first. As time progresses, holes become prevalent and the anomalous regions are local enhancements of ice concentration in a volume. The 3-D distribution is also heterogeneous because the ice is buried to varying depths in different locations. Analysis of the coherence of ice on 10 cm scales predicts that putative ice sheets in anomalous radar craters are < 100 Myr old. Surface frost becomes homogenized within 20 Myr. The simulations show the data from the LCROSS impact and surrounding region are consistent with the ice deposit in Cabeus being >1000 Myr old. For future in situ analysis of cold trap volatiles, a horizontal range of 10 m is sufficient to acquire surface-based measurements of heterogeneously distributed ice. These results also support previous analyses that Mercury's cold traps are young.

Hurley, Dana M.↗

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA'S ISRU Technology

Introduction: In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s) range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant. System Model: In FY24, the System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team developed a lunar water processing system model using the Mission Analysis and Integration Tool (MAIT) to estimate the M/V/P for ISRU subsystems operating under a wide range of Hydrogen (H2) and Oxygen (O2) production targets for the Space Technology Mission Directorate (STMD) [1]. Based on Japan Aerospace Exploration Agency’s (JAXA) surface operational requirements, this system architecture was modified to include the ability to excavate consolidated icy regolith (versus granular ice excavation using Kennedy Space Center’s (KSC) ISRU Pilot Excavator, IPEx) and explore the feasibility of processing the lunar water both inside and outside of the PSR. For the consolidated icy regolith case study, excavation was performed via a mobility transport chassis outfitted with The Regolith Ice Drill for Exploring New Terrain (TRIDENT) for drilling [2] and the Cold Operable Lunar Deployable Arm (COLDArm) [3] for regolith transfer. The system model determines the required rover and payload. M/V/P to handle the required regolith processing rates. The regolith is then sorted and heated to sublimate the ice (via an auger dryer). The exiting high temperature, low pressure vapor is cleaned of volatiles (via cold trap) and electrolyzed to produce H2 and O2. These products are then dried, liquified with 20 K and 90 K cryocoolers (for H2 and O2, respectively), and stored in cylindrical tanks. Study Goals: Due to the different ConOps options of regolith transport to the ridge for processing versus processing it directly inside the PSR, as well as the unknown regolith/water-ice composition, new excavation techniques and their power configurations are being evaluated within a ISRU system architecture for production targets less than NASA’s pilot plant (1 mT). This analysis investigates the feasibility of numerous excavation techniques, power architectures, and logistical operations and determines an optimal system configuration with regards to M/V/P. It aims to investigate which parameters, both locally and globally, have the greatest effect on each subsystem within the plant. This can be used to identify the most critical components of the plant, and guide future decisions on allocating funding for research and development. The results from this study may provide subsystem developers with appropriate interfaces with excavation subsystems and downstream processes, and assessing the overall feasibility of each excavation technique, power architecture, and logistical timeframe. References: [1] Carlson, A. et. al. (2024) ICES. [2] Zacny, K., et. al. (2024) “ASCE Earth and Space”. [3] McCormick, R., et. Al. (2024) IEEE Xplore.

ISRU↗

A case of poor substructure diagnostics

The NASTRAN Manuals in the substructuring area are all geared toward instant success, but the solution paths are fraught with many traps for human error. Thus, the probability of suffering a fatal abort is high. In such circumstances, the necessity for diagnostics that are user friendly is paramount. This paper is written in the spirit of improving the diagnostics as well as the documentation in one area where the author felt he was backed into a blind corner as a result of his having committed a data oversight. This topic is aired by referring to an analysis of a particular structure. The structure, under discussion, used a number of local coordinate systems that simplified the preparation of input data. The principle features of this problem are introduced by reference to a series of figures.

Butler, Thomas G.↗

Off-center nuclei in galaxies

The nucleus of a galaxy orbits around the mass centroid. Orbital motions appear overstable in numerical experiments started with a galaxy's nucleus at rest atop its mass centroid. The amplitude doubles in 6-10 orbital periods. Orbits precess, nutate, and change their amplitudes, but they keep fairly constant periods. Orbital periods are in resonance with local particle motions, and amplitudes reach a core radius. This resonance suggests that center motions are a local, rather than a global, phenomenon. The overstability implies that a galaxy cannot be formed in nature with its nucleus at rest atop its mass centroid, and that nuclei orbit the mass centroid in real galaxies. These center motions should show up observationally as a shift of the nucleus away from the center defined by nearby isophotes. Off-center nuclei have been reported in many galaxies (e.g., M33, M101, NGC 3379, NGC 3384). Other kinds of observations confirmed the picture of nonsteady galactic centers as well. Gas trapped in moving nuclear regions of a galaxy should show strange flow patterns with possible shocks. The nuclear regions of galaxies including Milky Way and of globular clusters are not likely to be in a static steady state.

Miller, R. H.↗

The cutoff frequency for fast-mode magnetohydrodynamic waves in an isothermal atmosphere with a uniform horizontal magnetic field

This study analytically examines conditions for reflection of MHD fast-mode waves propagating upward in an isothermal atmosphere. A new method of transforming the linearized wave equation into Klein-Gordon form is utilized to calculate a local cutoff (critical) frequency for these waves. This critical frequency determines the height in the atmosphere at which reflection dominates and above which wave propagation is effectively cut off. Comparison of our results to those previously obtained shows that earlier calculations of the critical frequency for MHD fast mode waves were done incorrectly. The results may be helpful in explaining the short-period end of the spectrum of the solar global p-mode oscillations. They may also be important in studies of wave propagation and wave trapping in highly magnetized stellar atmospheres.

Stark, B. A.↗

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

Active Solid State Dosimetry for Lunar EVA

The primary threat to astronauts from space radiation is high-energy charged particles, such as electrons, protons, alpha and heavier particles, originating from galactic cosmic radiation (GCR), solar particle events (SPEs) and trapped radiation belts in Earth orbit. There is also the added threat of secondary neutrons generated as the space radiation interacts with atmosphere, soil and structural materials.[1] For Lunar exploration missions, the habitats and transfer vehicles are expected to provide shielding from standard background radiation. Unfortunately, the Lunar Extravehicular Activity (EVA) suit is not expected to afford such shielding. Astronauts need to be aware of potentially hazardous conditions in their immediate area on EVA before a health and hardware risk arises. These conditions would include fluctuations of the local radiation field due to changes in the space radiation field and unknown variations in the local surface composition. Should undue exposure occur, knowledge of the dynamic intensity conditions during the exposure will allow more precise diagnostic assessment of the potential health risk to the exposed individual.[2]

Wrbanek, John D.↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

Ripple-induced neoclassical toroidal viscous torque in Augmented-First Plasma operation phase in ITER

Abstract A systematic calculation is performed on the ripple-induced neoclassical toroidal viscous (NTV) torque for new ITER scenarios designed for the Augmented-First Plasma (A-FP) operation phase with the full tungsten wall, where the plasma-wall gap is varied in view of mitigating the impact of tungsten wall-plasma interactions. The torque calculation includes drift kinetic response of the plasma thermal and energetic particles to the n = 18 ( n is the toroidal harmonic number) ripple field. For the plasma scenario with ~45 cm plasma-wall gap at the outboard mid-plane and considering the corrected ripple level of 0.17% by the ferritic steel inserts, the computed net NTV torque acting on the plasma column is in the sub-Nm level. However, with decreasing the plasma-wall gap, the computed net NTV torque can reach a level comparable to that produced by the neutral-beam momentum injection in ITER. Ripple correction by ferritic inserts reduces the net torque by a factor of 3.3 for all the three A-FP scenarios considered. The n ω d = l ω b (with ω d and ω b being the toroidal precession and bounce frequencies of trapped particles, respectively, and l an integer number) type of resonance-enhancement of the NTV torque, due to thermal particles, is found to be weak in ITER despite high- n of 18. The same also holds for the ITER 10 MA steady state scenario from the D-T operation phase, where the aforementioned resonance associated with fusion-born alphas is also included. The ripple-induced NTV torque is well below that produced by the resonant magnetic perturbation applied for controlling the type-I edge-localized mode in ITER.

ITER↗

Low Frequency Modulation of Extreme Temperature Regimes in a Changing Climate

The project examines long-term changes in extreme temperature episodes (ETE) associated with planetary climate modes (PCMs) in both the real atmospheric and climate model simulations. The focus is on cold air outbreaks (CAOs) and warm waves (WWs) occurring over the continental US during the past 60 winters. No significant long-term trends in either WWs or CAOs are observed over the US. The annual frequency of CAOs is affected by the (i) North Atlantic Oscillation (NAO) over the Southeast US and (ii) Pacific–North American (PNA) pattern over the Northwest US. WW frequency is influenced by the (i) NAO over the eastern US and (ii) combined influence of PNA, Pacific decadal oscillation (PDO), and ENSO over the southern US. The collective influence of PCMs accounts for as much as 50% of the regional variability in ETE frequency. During CAO (WW) events occurring over the southeast US, there are low (high) pressure anomalies at higher atmospheric levels over the southeast US with oppositely-signed pressure anomalies in the lower atmosphere over the central US. These patterns lead to anomalous northerly (for CAOs) or southerly (for WWs) flow into the southeast leading to cold or warm surface air temperature anomalies, respectively. One distinction is that CAOs involve substantial air mass transport while WW formation is more local in nature. The primary differences among event categories are in the origin and nature of the pressure anomaly features linked to ETE onset. In some cases, PCMs help to provide a favorable environment for event onset. Heat budget analyses indicate that latitudinal transport in the lower atmosphere is the main contributor to regional cooling during CAO onset. This is partly offset by adiabatic warming associated with subsiding air. Additional diagnoses reveal that this latitudinal transport is partly due to the remote physical influence of a shallow cold pool of air trapped along the east side of the Rocky Mountains. ETE and PCM behavior is also studied in (CMIP5) climate model simulations. Although the climate models considered are able to represent the overall behavior of ETEs, the frequency of WWs (CAOs) is too high (low) in many models. While all models qualitatively replicate the overall structure of the PNA pattern, a small minority of models fails to properly simulate the NAO pattern. Model shortcomings in representing the NAO and PNA patterns have important consequences for simulating associated regional variability in surface air temperature and storm track behavior. The influence of PCMs on ETEs is underestimated in most CMIP5 models. In particular, none of the models are able to accurately simulate observed linkages between ETEs and the PDO, due to a gross misrepresentation of the PDO pattern in most models. Our results indicate that predictions of future CAO and WW behavior are currently limited by the ability of climate models to accurately represent PCM characteristics. Our study also considers the behavior of PCMs known as annular modes. It is determined that north-south movements in the stratospheric jet stream (related to the Polar Annular Mode) result in long-lasting impacts upon surface weather conditions including regional air temperature anomalies. The structure and dynamics of the stratospheric northern annular mode (or SNAM, related to changes in the strength of the stratospheric jet stream) was studied in CMIP5 models. In models with poorly-resolved stratospheres, the amplitude of SNAM at stratospheric altitudes is typically too weak, consistent with weaker stratospheric jet variability. However, this distinction does not carry over to the associated tropospheric signature of SNAM. A regional analysis illustrates that most CMIP5 models (regardless of whether the stratosphere is well-resolved) have anomalously weak and eastward shifted (compared to observed SNAM events) storm track and sea level pressure anomaly patterns during SNAM events. Analyses of stratosphere–troposphere coupling reveal that large-scale wave activity in the stratosphere is anomalously weak in CMIP5 models having poorly-resolved stratospheres, due to correspondingly weak upward propagation of wave activity from the troposphere. This suggests that such models have anomalously weak stratosphere-troposphere coupling in association with SNAM events. Paradoxically, however, there is little apparent impact upon the attendant tropospheric variability.

54 ENVIRONMENTAL SCIENCES↗

Possible role of grain-boundary and dislocation structure for the magnetic-flux trapping behavior of niobium: A first-principles study

In this work, first-principles methods were used to understand magnetic flux trapping at vacancies, dislocations, and grain boundaries in high-purity superconducting niobium. Full-potential linear augmented plane-wave methods were applied in progressively greater complexity, starting at simple vacancies and extending to screw dislocations and tilt grain boundaries to analyze the effects of magnetic field on the superconducting state surrounding these defects. Density-functional theory calculations identified changes in electronic structure at the dislocation core and different types of symmetric tilt grain boundaries relative to bulk niobium. Electron redistribution enhanced nonparamagnetic effects and thus perturb superconductivity, resulting in local conditions suitable for premature flux penetration and subsequently flux pinning. Since the coherence length of superconducting niobium at 0 K is significantly larger than the lattice parameter, the effects of line and planar defects in niobium are predicted to be stronger for defect clusters than single defects in isolation, which is consistent with recent experimental observations. Controlling accumulation or depletion of charge at the defects, e.g., by segregation of an impurity atom, can mitigate these tendencies thus increasing the quality of superconducting niobium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dissipative superfluid hydrodynamics for the unitary Fermi gas

In this work we establish constraints on the temperature dependence of the shear viscosity η in the superfluid phase of a dilute Fermi gas in the unitary limit. Our results are based on analyzing experiments that measure the aspect ratio of a deformed cloud after release from an optical trap. Here, we discuss how to apply the two-fluid formalism to the unitary gas and provide a suitable parametrization of the equation of state. We show that in expansion experiments the difference between the normal and superfluid velocities remains small and can be treated as a perturbation. We find that expansion experiments favor a shear viscosity that decreases significantly in the superfluid regime. Using an exponential parametrization, we find η(T c /2T F ) ≲ 0.37[η(T c /T F )], where T c is the critical temperature and T F is the local Fermi temperature of the gas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pressure Ratio to Thermal Environments

A pressure ratio to thermal environments (PRatTlE.pl) program is a Perl language code that estimates heating at requested body point locations by scaling the heating at a reference location times a pressure ratio factor. The pressure ratio factor is the ratio of the local pressure at the reference point and the requested point from CFD (computational fluid dynamics) solutions. This innovation provides pressure ratio-based thermal environments in an automated and traceable method. Previously, the pressure ratio methodology was implemented via a Microsoft Excel spreadsheet and macro scripts. PRatTlE is able to calculate heating environments for 150 body points in less than two minutes. PRatTlE is coded in Perl programming language, is command-line-driven, and has been successfully executed on both the HP and Linux platforms. It supports multiple concurrent runs. PRatTlE contains error trapping and input file format verification, which allows clear visibility into the input data structure and intermediate calculations.

Lopez, Pedro↗