Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

Turbulent spots in boundary layers in a free-piston shock-tunnel flow

In this paper, experiments to detect turbulent spots in the transitional boundary layers, formed on a flat plate in a free-piston shock tunnel how, are reported. Experiments indicate that thin-film heat-transfer gauges are suitable for identifying turbulent-spot activity and can be used to identify parameters such as the convection rate of spots and the intermittency of turbulence.

A compact shock-assisted free-piston driver for impulse facilities

A free-piston driver that employs entropy-raising shock processes with diaphragm rupture has been constructed, which promises significant theoretical advantages over isentropic compression. Results from a range of conditions with helium and argon driver gases are reported. Significant performance gains were achieved in some test cases. Heat losses are shown to have a strong effect on driver processes. Measurements compare well with predictions from a quasi-one-dimensional numerical code.

A generalized Markov sampler

A recent development of the Markov chain Monte Carlo (MCMC) technique is the emergence of MCMC samplers that allow transitions between different models. Such samplers make possible a range of computational tasks involving models, including model selection, model evaluation, model averaging and hypothesis testing. An example of this type of sampler is the reversible jump MCMC sampler, which is a generalization of the Metropolis-Hastings algorithm. Here, we present a new MCMC sampler of this type. The new sampler is a generalization of the Gibbs sampler, but somewhat surprisingly, it also turns out to encompass as particular cases all of the well-known MCMC samplers, including those of Metropolis, Barker, and Hastings. Moreover, the new sampler generalizes the reversible jump MCMC. It therefore appears to be a very general framework for MCMC sampling. This paper describes the new sampler and illustrates its use in three applications in Computational Biology, specifically determination of consensus sequences, phylogenetic inference and delineation of isochores via multiple change-point analysis.

Segmenting eukaryotic genomes with the generalized Gibbs sampler

Eukaryotic genomes display segmental patterns of variation in various properties, including GC content and degree of evolutionary conservation. DNA segmentation algorithms are aimed at identifying statistically significant boundaries between such segments. Such algorithms may provide a means of discovering new classes of functional elements in eukaryotic genomes. This paper presents a model and an algorithm for Bayesian DNA segmentation and considers the feasibility of using it to segment whole eukaryotic genomes. The algorithm is tested on a range of simulated and real DNA sequences, and the following conclusions are drawn. Firstly, the algorithm correctly identifies non-segmented sequence, and can thus be used to reject the null hypothesis of uniformity in the property of interest. Secondly, estimates of the number and locations of change-points produced by the algorithm are robust to variations in algorithm parameters and initial starting conditions and correspond to real features in the data. Thirdly, the algorithm is successfully used to segment human chromosome 1 according to GC content, thus demonstrating the feasibility of Bayesian segmentation of eukaryotic genomes. The software described in this paper is available from the author's website (www.uq.edu.au/similar to uqjkeith/) or upon request to the author.

Emission and holographic interferometry measurements in a superorbital expansion tube

Free-piston-driven expansion tubes are capable of generating flaw conditions over a wide range of enthalpies ranging from orbital up to superorbital velocities. Initial optical measurements aimed at investigating the flow in such a facility are presented. Emission studies were used to identify impurities in the how and to investigate spectral regions that are accessible by optical techniques. At moderate enthalpies, it was found that significant radiation resulted from metallic contaminants. At high enthalpies, the spectrum consisted of a number of atomic lines together with a broadband background component indicative of the presence of electrons. The presence of this radiation may limit the applicability of optical techniques that require spectral regions free from the influence of atomic transitions or background radiation. Emission spectroscopy (through Stark broadened hydrogen lines) and two-wavelength holographic interferometry were used to measure the electron number density behind a bow shock on a blunt body at conditions where significant ionization was observed. They yielded average concentrations of (3 +/- 1) x 10(17) cm(-3) from the emission measurements and (3.8 +/- 0.6) x 10(17) cm(-3) from the interferometry.

Fuel stagnation temperature effects on mixing with supersonic combustion flows

An experimental study of the effect of fuel stagnation temperature on mixing in a supersonic hydrogen-air flame is described, The combustor consisted of a constant-area rectangular duct with a centrally located fuel-injection strut that spanned the width. A high-enthalpy stream of air was supplied by a free-piston shock tunnel, and heated hydrogen fuel, supplied by a gun-tunnel, was injected into the freestream as a coflowing planar jet. The freestream total enthalpies were 5.6, 6.5, and 9 MJ/kg, and fuel stagnation temperatures were 300, 450, and 700 K, Raising the fuel stagnation temperature increased the fuel velocity to be near that of the airstream and resulted in a decrease in the mixing rate, Even as the fuel and air velocities became equal, significant mixing still occurred because of a large difference in density, Increasing the freestream enthalpy reduced the difference between the initial air temperature and the adiabatic flame temperature, which in turn reduced the heat addition, and subsequently, the amount of pressure rise in the duct.

Boundary-layer transition measurements in hypervelocity flows in a shock tunnel

Experiments to investigate the transition process in hypervelocity boundary layers were performed in the T4 free-piston shock tunnel. An array of thin-film heat-transfer gauges was used to detect the location and extent of the transitional region on a 1500 mm long x 120 turn wide flat plate, which formed one of the walls of a duct. The experiments were performed in a Mach 6 flow of air with 6- and 12-MJ/kg nozzle-supply enthalpies at unit Reynolds numbers ranging from 1.6 x 10(6) to 4.9 x 10(6) m(-1). The results show that the characteristics typical of transition taking place through the initiation, growth, and merger of turbulent spots are evident in the heat-transfer signals. A 2-mm-high excrescence located 440 turn from the leading edge was found to be capable of generating a turbulent wedge within an otherwise laminar boundary layer at a unit Reynolds number of 2.6 x 10(6) m(-1) at the 6-MJ/kg condition. A tripping strip, located 100 mm from the leading edge and consisting of a line 37 teeth of 2 rum height equally spaced and spanning the test surface, was also found to be capable of advancing the transition location at the same condition and at the higher enthalpy condition.

Insights into collisional magmatism from isotopic fingerprints of melting reactions

Piston-cylinder experiments in the granite system demonstrate that a variety of isotopically distinct melts can arise from progressive melting of a single source. The relation between the isotopic composition of Sr and the stoichiometry of the observed melting reactions suggests that isotopic signatures of anatectic magmas can be used to infer melting reactions in natural systems. Our results also indicate that distinct episodes of dehydration and fluid-fluxed melting of a single, metapelitic source region may have contributed to the bimodal geochemistry of crustally derived leucogranites of the Himalayan orogen.