The Numerical Simulation of Radiative Shocks I: The elimination of numerical shock instabilities using a localized oscillation filter

We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are exacerbated by thermal instability when cooling is present, forming post--shock `stripes'. These are orthogonal to the classic post--shock 'ringing' fluctuations. The resulting post--shock `striping' substantially modifies the flow. We discuss three different methods to resolve this problem. These include (1) a method based on artificial viscosity; (2) grid--jittering and (3) a new localized oscillation filter that acts on specific grid cells in the shock front. These methods are tested using a radiative wall shock problem with an embedded shear layer. The artificial viscosity method is unsatisfactory since, while it does reduce post--shock ringing, it does not eliminate the stripes and the excessive shock broadening renders the calculation of cooling inaccurate, resulting in an incorrect shock location. Grid--jittering effectively counteracts striping. However, elsewhere on the grid, the shear layer is unphysically diffused and this is highlighted in an extreme case. The oscillation filter method removes stripes and permits other high velocity gradient regions of the flow to evolve in a physically acceptable manner. It also has the advantage of only acting on a small fraction of the cells in a two or three dimensional simulation and does not significantly impair performance.Comment: 20 pages, 6 figures, revised version submitted to ApJ Supplement Serie

Exact calculation of the ground-state dynamical spin correlation function of a S=1/2 antiferromagnetic Heisenberg chain with free spinons

We calculate the exact dynamical magnetic structure factor S(Q,E) in the ground state of a one-dimensional S=1/2 antiferromagnet with gapless free S=1/2 spinon excitations, the Haldane-Shastry model with inverse-square exchange, which is in the same low-energy universality class as Bethe's nearest-neighbor exchange model. Only two-spinon excited states contribute, and S(Q,E) is found to be a very simple integral over these states.Comment: 11 pages, LaTeX, RevTeX 3.0, cond-mat/930903

Transport Properties of a One-Dimensional Two-Component Quantum Liquid with Hyperbolic Interactions

We present an investigation of the sinh-cosh (SC) interaction model with twisted boundary conditions. We argue that, when unlike particles repel, the SC model may be usefully viewed as a Heisenberg-Ising fluid with moving Heisenberg-Ising spins. We derive the Luttinger liquid relation for the stiffness and the susceptibility, both from conformal arguments, and directly from the integral equations. Finally, we investigate the opening and closing of the ground state gaps for both SC and Heisenberg-Ising models, as the interaction strength is varied.Comment: 10 REVTeX pages + 4 uuencoded figures, UoU-002029

Meltwater Intrusions Reveal Mechanisms for Rapid Submarine Melt at a Tidewater Glacier

Submarine melting has been implicated as a driver of glacier retreat and sea level rise, but to date melting has been difficult to observe and quantify. As a result, melt rates have been estimated from parameterizations that are largely unconstrained by observations, particularly at the near-vertical termini of tidewater glaciers. With standard coefficients, these melt parameterizations predict that ambient melting (the melt away from subglacial discharge outlets) is negligible compared to discharge-driven melting for typical tidewater glaciers. Here, we present new data from LeConte Glacier, Alaska, that challenges this paradigm. Using autonomous kayaks, we observe ambient meltwater intrusions that are ubiquitous within 400 m of the terminus, and we provide the first characterization of their properties, structure, and distribution. Our results suggest that ambient melt rates are substantially higher (×100) than standard theory predicts and that ambient melting is a significant part of the total submarine melt flux. We explore modifications to the prevalent melt parameterization to provide a path forward for improved modeling of ocean-glacier interactions.This work was funded by NSF OPP Grants 1503910, 1504191, 1504288, and 1504521 and National Geographic Grant CP4-171R-17. Additionally, this research was supported by the NOAA Climate and Global Change Postdoctoral Fellowship Program, administered by UCAR’s Cooperative Programs for the Advancement of Earth System Science (CPAESS) under award #NA18NWS4620043B. These observations would not be possible without the skilled engineering team who developed the autonomous kayaks—including Jasmine Nahorniak, June Marion, Nick McComb, Anthony Grana, and Corwin Perren—and also the Captain and crew of the M/V Amber Anne. We thank Donald Slater and an anonymous reviewer for valuable feedback that improved this manuscript. Data availability: All of the oceanographic data collected by ship and kayak have been archived with the National Centers for Environmental Information (Accession 0189574, https://accession.nodc.noaa.gov/ 0189574). The glacier data have been archived at the Arctic Data Center (https://doi.org/10.18739/A22G44).Ye

Short-timescale Fluctuations in the Difference Light Curves of QSO 0957+561A,B: Microlensing or Noise?

From optical R band data of the double quasar QSO 0957+561A,B, we made two new difference light curves (about 330 days of overlap between the time-shifted light curve for the A image and the magnitude-shifted light curve for the B image). We observed noisy behaviours around the zero line and no short-timescale events (with a duration of months), where the term event refers to a prominent feature that may be due to microlensing or another source of variability. Only one event lasting two weeks and rising - 33 mmag was found. Measured constraints on the possible microlensing variability can be used to obtain information on the granularity of the dark matter in the main lensing galaxy and the size of the source. In addition, one can also test the ability of the observational noise to cause the rms averages and the local features of the difference signals. We focused on this last issue. The combined photometries were related to a process consisting of an intrinsic signal plus a Gaussian observational noise. The intrinsic signal has been assumed to be either a smooth function (polynomial) or a smooth function plus a stationary noise process or a correlated stationary process. Using these three pictures without microlensing, we derived some models totally consistent with the observations. We finally discussed the sensitivity of our telescope (at Teide Observatory) to several classes of microlensing variability.Comment: MNRAS, in press (LaTeX, 14 pages, 22 eps figures