A multi-dimensional, adiabatic, hydrodynamics code for studying tidal excitation

We have developed a parallel, simple, and fast hydrodynamics code for multi-dimensional, self-gravitating, adiabatic flows. Our primary motivation is the study of the non-linear evolution of white dwarf oscillations excited via tidal resonances, typically over hundreds of stellar dynamical times. Consequently, we require long term stability, low diffusivity, and high algorithmic efficiency. An explicit, Eulerian, finite-difference scheme on a regular Cartesian grid fulfills these requirements. It provides uniform resolution throughout the flow, as well as simplifying the computation of the self-gravitational potential, which is done via spectral methods. In this paper, we describe the numerical scheme and present the results of some diagnostic problems. We also demonstrate the stability of a cold white dwarf in three dimensions over hundreds of dynamical times. Finally, we compare the results of the numerical scheme to the linear theory of adiabatic oscillations, finding numerical quality factors on the order of 6000, and excellent agreement with the oscillation frequency obtained by the linear analysis.Comment: 11 pages, 10 figure

Secondary circulations in the bottom boundary layer over sedimentary furrows

Secondary circulation is known to be an important feature of many atmospheric and laboratory boundary layers. The presence of streamwise, helical, counterrotating vortices is documented here for the first time in the bottom boundary layer of a large natural body of water. Simultaneous vertical profiles of velocity and temperature were recorded on either side of a sedimentary furrow on the floor of Lake Superior (depth = 100 m) in November 1985. Flow roughness length zo was estimated to be 0.3 cm. Friction velocity u* and turbulent boundary layer thickness LW were estimated for each profile allowing for stratification effects. Typically, LW ≈ 10 m. Thermal stratification near the lake bed was an important constraint to boundary layer development; bottom mixed layers were absent in most cases. Mean flow toward the furrow at lower levels within the boundary layer (z/LW\u3c0.37) and mean flow away from the furrow at higher levels (0.37 \u3c z/LW \u3c 1.65) were observed for near‐bottom speeds greater than 6 cm/s when mean flow direction was within 25° of the furrow direction. This implies helical vortex pair circulations with upward motion over the furrow, consistent with earlier hypotheses. Cross‐stream (secondary) flows were typically 5% of the free stream (primary flow) speed. Streamwise vorticity in the range 1.5 m ≤ z ≤ 5 m was estimated to be 2.1×10−3 s−1

It's a wonderful tail: the mass loss history of Mira

Recent observations of the Mira AB binary system have revealed a surrounding arc-like structure and a stream of material stretching 2 degrees away in opposition to the arc. The alignment of the proper motion vector and the arc-like structure shows the structures to be a bow shock and accompanying tail. We have successfully hydrodynamically modelled the bow shock and tail as the interaction between the asymptotic giant branch (AGB) wind launched from Mira A and the surrounding interstellar medium. Our simulations show that the wake behind the bow shock is turbulent: this forms periodic density variations in the tail similar to those observed. We investigate the possiblity of mass-loss variations, but find that these have limited effect on the tail structure. The tail is estimated to be approximately 450,000 years old, and is moving with a velocity close to that of Mira itself. We suggest that the duration of the high mass-loss phase on the AGB may have been underestimated. Finally, both the tail curvature and the rebrightening at large distance can be qualitatively understood if Mira recently entered the Local Bubble. This is estimated to have occured 17 pc downstream from its current location.Comment: 12 pages, 3 colour figures, accepted by ApJ Part II (Letters

Dispersive wave runup on non-uniform shores

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto a non-uniform beach are presented along with great lines of the employed numerical method (see D. Dutykh et al. (2011) for more details).Comment: 8 pages, 6 figures, 18 references. This preprint is submitted to FVCA6 conference proceedings. Other author papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh

Observation and analysis of in situ carbonaceous matter in Nakhla: part II

Analysis of in situ carbonaceous matter in the Nakhla SNC meteorite has been carried out using a variety of techniques. Laser raman data shows the carbonaceous matter to be highly complex and static mass spectrometry has shown it to have an isotopic composition of '18 to '20' C

Athena: A New Code for Astrophysical MHD

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.Comment: 61 pages, 36 figures. accepted by ApJ

Enhanced Humoral Immune Response After COVID-19 Vaccination in Elderly Kidney Transplant Recipients on Everolimus Versus Mycophenolate Mofetil-containing Immunosuppressive Regimens

BACKGROUND: Elderly kidney transplant recipients (KTRs) represent almost one third of the total kidney transplant population. These patients have a very high coronavirus disease 2019 (COVID-19)-related mortality, whereas their response to COVID-19 vaccination is impaired. Finding ways to improve the COVID-19 vaccination response in this vulnerable population is of uttermost importance. METHODS: In the OPTIMIZE trial, we randomly assign elderly KTRs to an immunosuppressive regimen with standard-exposure calcineurin inhibitor (CNI), mycophenolate mofetil, and prednisolone or an adapted regimen with low dose CNI, everolimus, and prednisolone. In this substudy, we measured the humoral response after 2 (N = 32) and 3 (N = 22) COVID-19 mRNA vaccinations and the cellular response (N = 15) after 2 vaccinations. RESULTS: . The seroconversion rates of elderly KTRs on a standard immunosuppressive regimen were only 13% and 38% after 2 and 3 vaccinations, respectively, whereas the response rates of KTRs on the everolimus regimen were significantly higher at 56% (P = 0.009) and 100% (P = 0.006). Levels of severe acute respiratory syndrome coronaVirus 2 IgG antibodies were significantly higher at both time points in the everolimus group (P = 0.004 and P < 0.001). There were no differences in cellular response after vaccination. CONCLUSION: . An immunosuppressive regimen without mycophenolate mofetil, a lower CNI dose, and usage of everolimus is associated with a higher humoral response rate after COVID-19 vaccination in elderly KTRs after transplantation. This encouraging finding should be investigated in larger cohorts, including transplant recipients of all ages

Unsteady transonic cascade fow with in-passage shock wave

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76409/1/AIAA-25179-720.pd

Numerical simulations of the Accretion-Ejection Instability in magnetised accretion disks

The Accretion-Ejection Instability (AEI) described by Tagger & Pellat (1999) is explored numerically using a global 2d model of the inner region of a magnetised accretion disk. The disk is initially currentless but threaded by a vertical magnetic field created by external currents, and frozen in the flow. In agreement with the theory a spiral instability, similar in many ways to those observed in self-gravitating disks, develops when the magnetic field is, within a factor of a few, at equipartition with the disk thermal pressure. Perturbations in the flow build up currents and create a perturbed magnetic field within the disk. The present non-linear simulations give good evidence that such an instability can occur in the inner region of accretion disks, and generate accretion of gas and vertical magnetic flux toward the central object, if the equilibrium radial profiles of density and magnetic flux exceed a critical threshold.Comment: single tar file with GIF figure