
We propose a general technique for extracting a larger set of stable
information from persistent homology computations than is currently done. The
persistent homology algorithm is usually viewed as a procedure which starts
with a filtered complex and ends with a persistence diagram. This procedure is
stable (at least to certain types of perturbations of the input). This
justifies the use of the diagram as a signature of the input, and the use of
features derived from it in statistics and machine learning. However, these
computations also produce other information of great interest to practitioners
that is unfortunately unstable. For example, each point in the diagram
corresponds to a simplex whose addition in the filtration results in the birth
of the corresponding persistent homology class, but this correspondence is
unstable. In addition, the persistence diagram is not stable with respect to
other procedures that are employed in practice, such as thresholding a point
cloud by density. We recast these problems as realvalued functions which are
discontinuous but measurable, and then observe that convolving such a function
with a suitable function produces a Lipschitz function. The resulting stable
function can be estimated by perturbing the input and averaging the output. We
illustrate this approach with a number of examples, including a stable
localization of a persistent homology generator from brain imaging data.

Combined freestream disturbance measurements and receptivity studies in
hypersonic wind tunnels were conducted by means of a slender wedge probe and
direct numerical simulation. The study comprises comparative tunnel noise
measurements at Mach 3, 6 and 7.4 in two Ludwieg tube facilities and a shock
tunnel. Surface pressure fluctuations were measured over a wide range of
frequencies and test conditions including harsh test environments not
accessible to measurement techniques such as pitot probes and hotwire
anemometry. Quantitative results of the tunnel noise are provided in frequency
ranges relevant for hypersonic boundary layer transition. In combination with
the experimental studies, direct numerical simulations of the leadingedge
receptivity to fast and slow acoustic waves were performed for the slender
wedge probe at conditions corresponding to the experimental freestream
conditions. The receptivity to fast acoustic waves was found to be
characterized by an early amplification of the induced fast mode. For slow
acoustic waves an initial decay was found close to the leading edge. At all
Mach numbers, and for all considered frequencies, the leadingedge receptivity
to fast acoustic waves was found to be higher than the receptivity to slow
acoustic waves. Further, the effect of inclination angles of the acoustic wave
with respect to the flow direction was investigated. The combined numerical and
experimental approach in the present study confirmed the previous suggestion
that the slow acoustic wave is the dominant acoustic mode in noisy hypersonic
wind tunnels.

The Astropy Project (http://astropy.org) is, in its own words, "a community
effort to develop a single core package for Astronomy in Python and foster
interoperability between Python astronomy packages." For five years this
project has been managed, written, and operated as a grassroots,
selforganized, almost entirely volunteer effort while the software is used by
the majority of the astronomical community. Despite this, the project has
always been and remains to this day effectively unfunded. Further, contributors
receive little or no formal recognition for creating and supporting what is now
critical software. This paper explores the problem in detail, outlines possible
solutions to correct this, and presents a few suggestions on how to address the
sustainability of general purpose astronomical software.

A general analysis of the hydrodynamic limit of multirelaxation time lattice
Boltzmann models is presented. We examine multirelaxation time BGK collision
operators that are constructed similarly to those for the MRT case, however,
without explicitly moving into a moment space representation. The corresponding
'moments' are derived as left eigenvectors of said collision operator in
velocity space. Consequently we can, in a representation independent of the
chosen base velocity set, generate the conservation equations. We find a
significant degree of freedom in the choice of the collision matrix and the
associated basis which leaves the collision operator invariant. Therefore we
can explain why MRT implementations in the literature reproduce identical
hydrodynamics despite being based on different orthogonalization relations.

We analyze the Lattice Boltzmann method for the simulation of fluctuating
hydrodynamics by Adhikari et al. [Europhys. Lett. 71, 473 (2005)] and find that
it shows excellent agreement with theory even for small wavelengths as long as
a stationary system is considered. This is in contrast to other finite
difference and older lattice Boltzmann implementations that show convergence
only in the limit of large wavelengths. In particular cross correlators vanish
to less than 0.5%. For larger mean velocities, however, Galilean invariance
violations manifest themselves through errors of a magnitude similar to those
of the earlier implementations.