Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak localization”

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 1,081 records · Page 60

Evaluating some computer enhancement algorithms that improve the visibility of cometary morphology

The observed morphology of cometary comae is determined by ejection circumstances and the interaction of the ejected material with the local environment. Anisotropic emission can provide useful information on such things as orientation of the nucleus, location of active areas on the nucleus, and the formation of ion structure near the nucleus. However, discrete coma features are usually diffuse, of low amplitude, and superimposed on a steep intensity gradient radial to the nucleus. To improve the visibility of these features, a variety of digital enhancement algorithms were employed with varying degrees of success. They usually produce some degree of spatial filtering, and are chosen to optimize visibility of certain detail. Since information in the image is altered, it is important to understand the effects of parameter selection and processing artifacts can have on subsequent interpretation. Using the criteria that the ideal algorithm must enhance low contrast features while not introducing misleading artifacts (or features that cannot be seen in the stretched, unprocessed image), the suitability of various algorithms that aid cometary studies were assessed. The strong and weak points of each are identified in the context of maintaining positional integrity of features at the expense of photometric information.

Larson, S. M.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(p)(k), where A(p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for negative values of A(p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. it is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the long-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, S. B.↗

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

Computer experiments on the structure and dynamics of spiral galaxies

The evolution of an initially balanced rotating disk of stars with an initial velocity dispersion given by Toomre's local criterion was investigated by means of a computer model for isolated disks of stars. It was found that the disk is unstable against very large scale modes. A stable axisymmetric disk with a velocity dispersion much larger than that given by Toomre's criterion was generated. The final mass distribution for the disk gives a high density central core and a disk population of stars that is closely approximated by an exponential variation. Various methods and rates of cooling the hot axisymmetric disks were investigated. It was found that the cooling resulted in the development of two-arm spiral structures which persisted as long as cooling continued. An experiment was performed to induce spiral structure in a galaxy by means of the close passage of a companion galaxy. Parameters similar to those expected for M51 and its companion were used. It was found that because of the high velocity dispersion of the disturbed disk galaxy, only a weak two-arm spiral structure appeared. The evolution of a uniformly rotating disk galaxy which is a stationary solution of the collisionless Boltzmann equation was investigated for various values of the initial rms velocity dispersion. It was found that the disk becomes stable at a value of the velocity dispersion predicted by theory.

Hohl, F.↗

Electronic density of states of a U ( 1 ) quantum spin liquid with spinon Fermi surface. I. Orbital magnetic field effects

A quantum spin liquid with a spinon Fermi surface is an exotic insulator that hosts neutral Fermi surfaces inside the insulating gap. In an external magnetic field, it has been pointed out that the neutral Fermi surfaces are Landau quantized to form Landau levels due to the coupling to the induced emergent gauge magnetic field. In this paper, we calculate the electronic density of states (as observed in tunneling experiments) of the quantum spin liquid in an orbital magnetic field. Here, we find that the Landau levels from the neutral Fermi surfaces give rise to a set of steps emerging at the upper and lower Hubbard band edges. Each of the Hubbard band edge steps further develop into a band edge resonance peak when a weak gauge binding arising from the gauge field fluctuations is taken into account. Importantly, each Hubbard band edge step and its resulting resonance peak in the weak gauge binding are found to have a corresponding Landau level from the neutral Fermi surfaces, so the Hubbard band edge steps and the band edge resonance peaks provide signatures to the unique feature that characterizes the Landau quantization of the in-gap neutral Fermi surfaces in the spin liquid. We further consider the strong gauge-binding regime where the band edge resonance peaks move into the Mott gap and develop into true in-gap bound states. In the strong gauge binding regime, we solve the Landau-level spectrum of the in-gap bound states in an orbital magnetic field. For the in-gap bound state with a Mexican-hat-like band dispersion, we find that the envelop energy to have a state excited from the bound state Landau levels decreases quadratically with the magnetic field. The quadratic decrease behavior of the envelop energy is consistent with the intuition that applying magnetic field localizes the states and energetically promotes the in-gap bound states formation. Finally, we discuss the connection of our results to the electronic density of states spectra measured in the layered 1T-TaS 2 . We point out that a quantum spin liquid with a quasibound state in the upper Hubbard band can give the density of states spectra similar to the one measured in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spectroscopic Determination of Stellar Parameters and Oxygen Abundances for Hyades/Field G–K Dwarfs

It has been suggested that Fe abundances of K dwarfs derived from Fe i and Fe ii lines show considerable discrepancies, and oxygen abundances determined from high-excitation O i 7771–5 triplet lines are appreciably overestimated (the problem becoming more serious toward lower T {sub eff}), which, however, has not yet been widely confirmed. With the aim of clarifying this issue, we spectroscopically determined the atmospheric parameters of 148 G–K dwarfs (Hyades cluster stars and field stars) by assuming the classical Fe i/Fe ii ionization equilibrium as usual, and determined their oxygen abundances by applying the non-local thermal equilibrium spectrum fitting analysis to O i 7771–5 lines. It turned out that the resulting parameters did not show any significant inconsistency with those determined by other methods (for example, the mean differences in T {sub eff} and logg from the well-determined solutions of Hyades dwarfs are mostly ≲100 K and ≲0.1 dex). Likewise, the oxygen abundances of Hyades stars are around [O/H] ∼ +0.2 dex (consistent with the metallicity of this cluster) without exhibiting any systematic T {sub eff}-dependence. Accordingly, we conclude that parameters can be spectroscopically evaluated to a sufficient precision in the conventional manner (based on the Saha–Boltzmann equation for Fe i/Fe ii) and oxygen abundances can be reliably determined from the O i 7771–5 triplet for K dwarfs as far as stars of T {sub eff} ≳ 4500 K are concerned. We suspect that previously reported strongly T {sub eff}-dependent discrepancies may have stemmed mainly from overestimation of weak-line strengths and/or improper T {sub eff} scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Ground Based Studies of Thermocapillary Flows in Levitated Drops

Ground-based experiments together with analytical studies are presently being conducted for levitated drops. Both acoustic and electrostatic techniques are being employed to achieve levitation of drops in a gaseous environment. The scientific effort is principally on the thermal and the fluid phenomena associated with the local heating of levitated drops, both at 1-g and at low-g. In particular, the thermocapillary flow associated with local spot heating is being studied. Fairly stable acoustic levitation of drops has been achieved with some exceptions when random rotational motion of the drop persists. The flow visualization has been carried out by light scattering from smoke particles for the exterior flow and fluorescent tracer particles in the drop. The results indicate a lack of axial symmetry in the internal flow even though the apparatus and the heating are symmetric. The theoretical studies for the past year have included fundamental analyses of acoustically levitated spherical drops. The flow associated with a particle near the velocity antinode is being investigated by the singular perturbation technique. As a first step towards understanding the effect of the particle displacement from the antinode, the flow field about the node has been calculated for the first time. The effect of the acoustic field on the interior of a liquid drop has also been investigated. The results predict that the internal flow field is very weak.

Sadhal, Satwindar Singh↗

Missile impact craters (White Sands Missile Range, New Mexico) and applications to lunar research: Contributions to astrogeology

Craters in natural materials at White Sands Missile Range, N. Mex., were produced by the impact of high-velocity to hypervelocity missiles traveling along oblique trajectories with kinetic energies between 2.1 and 81 × 1014 ergs. The oblique impacts produce craters 2 to 10 m across with morphologies and ejecta that are bilaterally symmetrical with respect to the plane of the missile trajectory. Rims are high and the amount of ejecta large in down-trajectory and lateral directions, whereas rims are low to nonexistent and ejecta thin to absent up-trajectory. Symmetry development and modifications of the symmetry are a function of target material, local topography, and angle of impact. Seven mappable units can be recognized in and around the craters. Three of these are ejecta: thick ejecta near the crater, thin to discontinuous ejecta at greater distances, and scattered ejecta at the greatest distances to the limit of throwout. These ejecta units may be absent on the up-trajectory side; if present, they are rarely as thick or continuous as on other sides of the crater. Three units are target materials: undeformed target material exposed in local patches through thin to discontinuous ejecta and everywhere between the fragments of scattered ejecta, tilted and broken target material exposed in upper crater walls, and shattered and fractured target material exposed on the up-trajectory crater wall. The seventh unit is slope material composed of talus and fallback within the crater. Development, character, and exposure of these units varies chiefly with the target material. Ejecta from the craters is chiefly broken but relatively undeformed target material that may range in size from very fine grained debris to large blocks. Where the target is porous, significant amounts of the ejecta are composed of sheared and compressed fragments, some coated with dark layers of mixed projectile pieces, powder, and fused metal mixed with crushed target material. For layered targets, the original stratigraphic sequence is crudely preserved and in inverted order in thick ejecta. Secondary impact craters are produced by the impact of ejected fragments when the surrounding surface materials are sufficiently weak. A wide variety of secondary impact crater relations may result. Secondary craters nearest the primary crater have blocks in them that are larger than or the same size as the crater they produced. Farther from the primary crater, the fragments are generally smaller than the secondary crater and are ejected from it. Excavation of four craters revealed a mixed breccia beneath the crater floor composed of missile pieces, sheared and compressed target material, and crushed debris. Banded, disaggregated target material and nonmixed breccia surrounded the mixed breccia, and these breccias were surrounded by a zone of conjugate fractures. Beneath the ejecta on the lateral and down-trajectory crater flanks, the target materials were tilted upward and broken. Up-trajectory, open fractures and downward displacement occurred in two of the craters. No displacement was observed for the other two. Beneath the down-trajectory rims of craters with distinct layering, overturned synclines were observed. Missile breakup and behavior during cratering are a function of target and missile properties. Missile breakup depends on missile velocity and is most extensive at high velocities, where the missile is fragmented, powdered, and partly fused. Burial of missile or its fragmented, powdered, and fused remains is greatest for porous targets and least for dense cohesive targets. For very porous targets, camouflet structures containing the fragmented missile may form. Least squares fit to the data on craters in dry to moist targets indicate V(a) = 10^(-11.433)E(p)^(1.205) where V(a) is the volume of the apparent crater and E(p) is the kinetic energy of the missile. This equation is consistent with expectations of the equations relating apparent depth and radius to kinetic energy. Extrapolation of displaced masses and kinetic energies for laboratory impacts with sand and rock converge near 10^(15) to 10^(16) ergs, where the extrapolations are near the data on missile impact craters, corrected for impact angle. Displaced masses of craters produced by missile impacts and by chemical explosives with small scaled depths of burial are about the same when the kinetic energies of the missiles (corrected for angle of impact) are equal to the TNT equivalent energy of the explosive. The problem of equivalent scaled depth of burst for an impact crater is complicated and not entirely resolved, however. Both missile impact craters and chemical explosive craters in water-saturated targets are larger than their counterparts in dry to moist materials. Data collected during the study of missile impact craters have helped resolve a number of problems in lunar research: (1) the soillike nature of lunar surface materials was predicted, (2) sizes of craters produced by artificial impacts were correctly predicted, (3) certain features imaged by Surveyor were found to be analogous to features associated with missile impact craters, {4) missile impacts were used in support of the Apollo passive seismic experiment, (5) craters seen in Apollo orbital photographs were found to be similar to some missile impact craters, (6) missile impact craters supplied data on sample collection and crater phenomenology used in training astronauts, and (7) some returned lunar samples are similar to coated, sheared, and compressed fragments ejected from missile impact craters.

H. J. Moore↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. It is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the large-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, Stephen B.↗

Seasonal Evolution and Interannual Variability of the Local Solar Energy Absorbed by the Arctic Sea Ice-Ocean System

The melt season of the Arctic sea ice cover is greatly affected by the partitioning of the incident solar radiation between reflection to the atmosphere and absorption in the ice and ocean. This partitioning exhibits a strong seasonal cycle and significant interannual variability. Data in the period 1998, 2000-2004 were analyzed in this study. Observations made during the 1997-1998 SHEBA (Surface HEat Budget of the Arctic Ocean) field experiment showed a strong seasonal dependence of the partitioning, dominated by a five-phase albedo evolution. QuikSCAT scatterometer data from the SHEBA region in 1999-2004 were used to further investigate solar partitioning in summer. The time series of scatterometer data were used to determine the onset of melt and the beginning of freezeup. This information was combined with SSM/I-derived ice concentration, TOVS-based estimates of incident solar irradiance, and SHEBA results to estimate the amount of solar energy absorbed in the ice-ocean system for these years. The average total solar energy absorbed in the ice-ocean system from April through September was 900 MJ m(sup -2). There was considerable interannual variability, with a range of 826 to 1044 MJ m(sup -2). The total amount of solar energy absorbed by the ice and ocean was strongly related to the date of melt onset, but only weakly related to the total duration of the melt season or the onset of freezeup. The timing of melt onset is significant because the incident solar energy is large and a change at this time propagates through the entire melt season, affecting the albedo every day throughout melt and freezeup.

Perovich, Donald K.↗

HostSub_GP: Precise Galaxy Background Subtraction in Transient Long-slit Spectroscopy with Gaussian Processes

We present a novel host galaxy subtraction technique in long-slit spectroscopy for extragalactic transients. Unlike classic methods which generally estimate the background using simple interpolation of local galaxy flux in the 2D spectrum, our approach leverages multi-band archival images of the host galaxies to model the background emission from the galaxy in the 2D spectrum. Such imaging encodes the wavelength-dependent galaxy profile along the slit, and is readily accessible through wide-field imaging surveys. We construct a smooth prior for the 2D galaxy profile with a Gaussian process (GP) based on these reference images, and use another GP to model the correlated deviations from the prior in the observed spectrum. This enables accurate inference of the galaxy flux blended with the transient. On synthetic long-slit data of a spiral galaxy extracted from a Multi Unit Spectroscopic Explorer hyper-spectral cube, the GP method remains robust as long as the host galaxy is spatially resolved and consistently outperforms classic methods. We apply the method to archival Keck spectra of two real transients, SN 2019eix and AT 2019qiz, to further demonstrate how the method uniquely recovers weak spectral features amid strong galaxy contamination, enabling refined constraints on the properties of both transients. We have released the software implementation, HostSub_GP, a scalable toolkit that leverages JAX, with an MIT license.

79 ASTRONOMY AND ASTROPHYSICS↗

Mirror dark sector solution of the Hubble tension with time-varying fine-structure constant

Here we explore a model introduced by Cyr-Racine et al. that resolves the Hubble tension by invoking a “mirror world” dark sector with energy density a fixed fraction of the “ordinary” sector of Λ+ cold dark matter. Although it reconciles cosmic microwave background and large-scale structure observations with local measurements of the Hubble constant, H 0 , the model requires a value of the primordial helium mass fraction, Y p = 0.170 ± 0.025, that is discrepant with observations and with the predictions of big bang nucleosynthesis. We consider a variant of the model with standard helium mass fraction but with the value of the electromagnetic fine-structure constant, α, slightly different during photon decoupling from its present value. If α at that epoch is lower than its current value by Δα ≃ -2 x 10 -5 , then we can achieve the same Hubble tension resolution as in the work by Cyr-Racine et al. but with consistent helium abundance. As an example of such time evolution of α, we consider a toy model of an ultralight scalar field, with mass m < 4 x 10 -29 eV, coupled to electromagnetism, which evolves after photon decoupling at redshift z ≃ 10 3 , and that with appropriate coupling appears to be consistent with late-time constraints on α variation and the weak equivalence principle.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion

The classical Landau/Levich models of liquid propellant combustion, which serve as seminal examples of hydrodynamic instability in reactive systems, have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and/or temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity, surface tension and viscosity on the hydrodynamic instability of the propagating liquid/gas interface. In particular, a composite asymptotic expression, spanning three distinguished wavenumber regimes, is derived for both cellular and pulsating hydrodynamic neutral stability boundaries A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. For the case of cellular (Landau) instability, the results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that cellular hydrodynamic instability in this context is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(l) wavenumber disturbances. It is also demonstrated that, in the large wavenumber regime, surface tension and both liquid and gas viscosity all produce comparable stabilizing effects in the large-wavenumber regime, thereby providing significant modifications to previous analyses of Landau instability in which one or more of these effects were neglected. In contrast, the pulsating hydrodynamic stability boundary is found to be insensitive to gravitational and surface-tension effects, but is more sensitive to the effects of liquid viscosity, which is a significant stabilizing effect for O(l) and higher wavenumbers. Liquid-propellant combustion is predicted to be stable (i.e., steady and planar) only for a range of negative pressure sensitivities that lie between the two types of hydrodynamic stability boundaries.

Margolis, Stephen B.↗

Theory of the phase transition in random unitary circuits with measurements

Here, we present a theory of the entanglement transition tuned by measurement strength in qudit chains evolved by random unitary circuits and subject to either weak or random projective measurements. The transition can be understood as a nonanalytic change in the amount of information extracted by the measurements about the initial state of the system, quantified by the Fisher information. To compute the von Neumann entanglement entropy $\textit{S}$ and the Fisher information $\mathcal{F}$, we apply a replica method based on a sequence of quantities $\tilde{S}^{(n)}$ and $\mathcal{F}^{(n)}$ that depend on the $\textit{n}$th moments of density matrices and reduce to $\textit{S}$ and $\mathcal{F}$ in the limit $\textit{n}$ → 1. These quantities with $\textit{n}$ ≥ 2 are mapped to free energies of a classical spin model with $\textit{n}$! internal states in two dimensions with specific boundary conditions. In particular, $\tilde{S}^{(n)}$ is the excess free energy of a domain wall terminating on the top boundary, and $\mathcal{F}^{(n)}$ is related to the magnetization on the bottom boundary. Phase transitions occur as the spin models undergo ordering transitions in the bulk. Taking the limit of large local Hilbert space dimension $\textit{q}$ followed by the replica limit $\textit{n}$ → 1, we obtain the critical measurement probability $p_c$ = 1/2 and identify the transition as a bond percolation in the 2D square lattice in this limit. Finally, we show there is no phase transition if the measurements are allowed in an arbitrary nonlocal basis, thereby highlighting the relation between the phase transition and information scrambling. We establish an explicit connection between the entanglement phase transition and the purification dynamics of a mixed state evolution and discuss implications of our results to experimental observations of the transition and simulability of quantum dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonlinear dynamics, avalanches, and noise for driven Wigner crystals

Here we consider the driven dynamics of Wigner crystals interacting with random disorder. Using numerical simulations, we find a rich variety of transport phenomena as a function of charge density, drive, and pinning strength. For weak pinning, the system forms a defect-free crystal that depins elastically. When the pinning is stronger, a pinned glass phase appears that depins into a filamentary flow state, transitions at higher drives into a disordered flow phase, and finally forms a moving smectic. Within the filamentary flow phase, the conduction curves can show switching dynamics as well as negative differential conductivity in which switching events cause some flow channels to be blocked. The velocity noise in the filamentary flow regime exhibits narrow band characteristics due to the one-dimensional nature of the motion, while the moving smectic has narrow band velocity noise with a washboard frequency. In the disordered flow state, the noise power reaches a peak value and the noise has a 1/$\mathcal{f}$ character. Our transport results are consistent with recent experimental transport studies in systems where Wigner crystal states are believed to be occurring. Below the conduction threshold, we find that avalanches with a power-law size distribution appear when there are sudden local rearrangements of charges from one pinned configuration to another.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Detection of radio continuum emission from Procyon

We have detected the F5 IV-V star Procyon as a weak and variable 3.6 cm radio continuum source using the VLA. The inferred radio luminosity is similar to, though some-what higher than, the X-band luminosity of the active and flaring sun. The 33 micro-Jy flux density level at which we detected Procyon on four of five occasions is close to the 36 micro-Jy radio flux density expected from a model in which the radio emission consists of two components: optically thick 'stellar disk' emission with a 3.6 cm brightness temperature of 20,000 K that is 50 percent larger than the solar value, and optically thin coronal emission with an emission measure the same as that indicated by Einstein and EXOSAT X-ray flux measurements in 1981 and 1983. The maximum mass-loss rate of a warm stellar wind is less than 2 x 10 exp -11 solar mass/yr. An elevated flux density of 115 micro-Jy observed on a single occasion provides circumstantial evidence for the existence of highly localized magnetic fields on the surface of Procyon.

Drake, Stephen A.↗

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Joining and Assembly of Bulk Metallic Glass Composites Through Capacitive Discharge

Bulk metallic glasses (BMGs), a class of amorphous metals defined as having a thickness greater than 1 mm, are being broadly investigated by NASA for use in spacecraft hardware. Their unique properties, attained from their non-crystalline structure, motivate several game-changing aerospace applications. BMGs have low melting temperatures so they can be cheaply and repeatedly cast into complex net shapes, such as mirrors or electronic casings. They are extremely strong and wear-resistant, which motivates their use in gears and bearings. Amorphous metal coatings are hard, corrosion-resistant, and have high reflectivity. BMG composites, reinforced with soft second phases, can be fabricated into energy-absorbing cellular panels for orbital debris shielding. One limitation of BMG materials is their inability to be welded, bonded, brazed, or fastened in a convenient method to form larger structures. Cellular structures (which can be classified as trusses, foams, honeycombs, egg boxes, etc.) are useful for many NASA, commercial, and military aerospace applications, including low-density paneling and shields. Although conventional cellular structures exhibit high specific strength, their porous structures make them challenging to fabricate. In particular, metal cellular structures are extremely difficult to fabricate due to their high processing temperatures. Aluminum honeycomb sandwich panels, for example, are used widely as spacecraft shields due to their low density and ease of fabrication, but suffer from low strength. A desirable metal cellular structure is one with high strength, combined with low density and simple fabrication. The thermoplastic joining process described here allows for the fabrication of monolithic BMG truss-like structures that are 90% porous and have no heat-affected zone, weld, bond, or braze. This is accomplished by welding the nodes of stacked BMG composite panels using a localized capacitor discharge, forming a single monolithic structure. This removes many complicated and costly fabrication steps. Moreover, the cellular structures detailed in this work are among the highest- strength and most energy-absorbent materials known. This implies that a fabricated structure made from these materials would have unequaled mechanical properties compared to other metal foams or trusses. The process works by taking advantage of the electrical properties of the matrix material in the metal-matrix composite, which in this case is a metallic glass. Due to the random nanoscale arrangement of atoms (without any grain boundaries), the matrix glass exhibits a near-constant electrical resistivity as a function of temperature. By placing the composite panels between two copper electrode plates and discharging a capacitor, the entire matrix of the panel can be heated to approximately 700 C in 10 milliseconds, which is above the alloy s solidus but below the liquidus. By designing the geometry of the panels into the shape of an egg box, the electrical discharge localizes only in the tips of each pyramidal cell. By applying a forging load during discharge, the nodes of the panels can be fused together into a single piece, which then dissipates heat through radiation back into a glassy state. This means that two panels can be metallurgically fused into one panel with no heat-affected zone, creating a seamless connection between panels. During the process, the soft metal particles (dendrites) that are uniformly distributed in the glassy matrix to increase the toughness are completely unaffected by the thermoplastic joining. The novelty is that a truss (or foam-like) structure can be formed with excellent energy- absorbing capabilities without the need for machining. The technique allows for large-scale fabrication of panels, well-suited for spacecraft shields or military vehicle door panels. Crystalline metal cellular structures cannot be fabricated using the thermoplastic joining technique described here. If metal panels were te assembled into a cellular structure, they would either have to be welded, brazed, bonded, or fastened together, creating a weak spot in the structure at each connection. Welded parts require a welding material to be added to the joint and exhibit a soft and weak heat-affected zone. Brazing and bonding do not form a metallurgical joint and thus exhibit low strengths, especially when the panels are pulled apart and fasteners require high-stress-concentration holes to be drilled. No equivalent rapid heating method exists for assembling metal panels together into cellular structures, and thus, those parts must be foamed, machined, or investment cast if they are to form a monolithic structure. If the crystalline panels were to be joined using capacitive discharge, as with a spot welder, their bond would be very weak, and the panels would have to be extremely thin. In contrast, the strength of joined BMG parts has been demonstrated to have strength comparable to the parent material. This technique opens up the possibility of using large-scale BMG hardware in spacecraft, military, or commercial applications.

Hofmann, Douglas C.↗