Relatively few studies have investigated the genetic population structure of sea anemones. This is particularly true for sea anemones that host some of the most iconic fishes on coral reefs, the anemonefishes. One of the main reasons for this knowledge gap is the lack of appropriate genetic markers. We developed and characterized a total of 47 novel polymorphic microsatellite markers for four host sea anemone species from the Indo-Pacific: Entacmaea quadricolor (n = 16 microsatellite markers), Heteractis magnifica (n = 8), Stichodactyla mertensii (n = 13), and Stichodactyla gigantea (n = 10). Here, we report genetic diversity statistics from two different sampling locations for each anemone species. Overall, we found that most markers were highly polymorphic. On average, we found a mean of seven alleles per locus. Observed and expected heterozygosities displayed high variation among loci, ranging from 0.033 to 0.980 and from 0.038 to 0.927, respectively. Only four loci showed deviations of Hardy–Weinberg equilibrium in both populations and were identified as having null alleles. Additionally, two pairs of loci were identified to be in linkage disequilibrium in only one population. Host anemones are highly sought after in the marine aquarium trade and are susceptible to thermal bleaching. Although most studies focus on their obligate symbionts (the anemonefish), genetic analyses of host sea anemones can expand our understanding of the biology, connectivity, and population structure of these organisms and potentially help develop conservation strategies that will aid both the host and its symbionts.

Direct numerical simulation is performed for flow past an isolated cylinder at Re=1,000. The corners of the cylinder are rounded at different radii, with the non-dimensional radius of curvature varying from R+=R/D=0.000 (square cylinder with sharp corners) to 0.500 (circular cylinder), in which R is the corner radius and D is the cylinder diameter. Our objective is to investigate the effect of the rounded corners on the development of the separated and transitional flow past the cylinder in terms of time-averaged statistics, time-dependent behavior, turbulent statistics and three-dimensional flow patterns. Numerical results reveal that the rounding of the corners significantly reduces the time-averaged drag and the force fluctuations. The wake flow downstream of the square cylinder recovers the slowest and has the largest wake width. However, the statistical quantities do not monotonically vary with the corner radius, but exhibit drastic variations between the cases of square cylinder and partially rounded cylinders, and between the latter and the circular cylinder. The free shear layer separated from the R+=0.125 cylinder is the most stable in which the first roll up of the wake vortex occurs furthest from the cylinder and results in the largest recirculation bubble, whose size reduces as R+ further increases. The coherent and incoherent Reynolds stresses are most pronounced in the near-wake close to the reattachment point, while also being noticeable in the shear layer for the square and R+=0.125 cylinders. The wake vortices translate in the streamwise direction with a convection velocity that is almost constant at approximately 80% of the incoming flow velocity. These vortices exhibit nearly the same trajectory for the rounded cylinders and are furthest away from the wake centerline for the square one. The flow past the square cylinder is strongly three-dimensional as indicated by the significant primary and secondary enstrophy, while it is dominated by the primary enstrophy (View the MathML source) for the rounded cylinders.

We investigate the performance of various discrete Hodge star operators for discrete exterior calculus (DEC) using circumcentric and barycentric dual meshes. The performance is evaluated through the DEC solution of Darcy and incompressible Navier–Stokes flows over surfaces. While the circumcentric Hodge operators may be favorable due to their diagonal structure, the barycentric (geometric) and the Galerkin Hodge operators have the advantage of admitting arbitrary simplicial meshes. Numerical experiments reveal that the barycentric and the Galerkin Hodge operators retain the numerical convergence order attained through the circumcentric (diagonal) Hodge operators. Furthermore, when the barycentric or the Galerkin Hodge operators are employed, a super-convergence behavior is observed for the incompressible flow solution over unstructured simplicial surface meshes generated by successive subdivision of coarser meshes. Insofar as the computational cost is concerned, the Darcy flow solutions exhibit a moderate increase in the solution time when using the barycentric or the Galerkin Hodge operators due to a modest decrease in the linear system sparsity. On the other hand, for the incompressible flow simulations, both the solution time and the linear system sparsity do not change for either the circumcentric or the barycentric and the Galerkin Hodge operators.

We investigate the inception of cavitation and resulting structures when a sphere collides with a solid surface covered with a layer of non-Newtonian liquid having a kinematic viscosity of up to (Formula presented.) cSt. We show the existence of shear-stress-induced cavitation during sphere approach towards the base wall (i.e. the pressurization stage) in ultra-viscous films using a synchronized dual-view high-speed imaging system. For the experimental parameters employed, liquids having viscoelastic properties of (Formula presented.) are shown to enable sphere rebound without any prior contact with the solid wall. Cavitation by depressurization (i.e. during rebound) in such non-contact cases is observed to onset after a noticeable delay from when the minimum gap distance is reached. Also, the cavities created originate from remnant bubbles, being the remains of the primary bubble entrapment formed by the lubrication pressure of the air during film entry. Cases where physical contact occurs (contact cases) in 10 000 cSt (Formula presented.) cSt films produce cavities attached to the base wall, which extend into an hourglass shape. In contrast, strikingly different structures occur in the most viscous liquids due to the disproportionality in radial expansion and longitudinal extension along the cavity length. Horizontal shear rates calculated using particle image velocimetry (PIV) measurements show the apparent fluid viscosity to vary substantially as the sphere approaches and rebounds away from the base wall. A theoretical model based on the lubrication assumption is solved for the squeeze flow in the regime identified for shear-induced cavity events, to investigate the criterion for cavity inception in further detail. © 2016 Cambridge University Press

High-surface-area π-conjugated polymeric networks have the potential to lend outstanding capacitance to supercapacitors because of the pronounced faradaic processes that take place across the dense intimate interface between active material and electrolytes. In this report, we describe how benzo[1,2-b:3,4-b’:5,6-b’’]trithiophene (BTT) and tris-EDOT-benzo[1,2-b:3,4-b’:5,6-b’’]trithiophene (TEBTT) can serve as 2D (trivalent) building blocks in the development of electropolymerized hierarchical π-conjugated frameworks with particularly high areal capacitance. In comparing electropolymerized networks of BTT, TEBTT, and their copolymers with EDOT, we show that P(TEBTT/EDOT)-based frameworks can achieve higher areal capacitance (e.g., as high as 443.8 mF cm-2 at 1 mA cm-2) than those achieved by their respective homopolymers (PTEBTT and PEDOT) in the same experimental conditions of electrodeposition (PTEBTT: 271.1 mF cm-2 (at 1 mA cm-2) and PEDOT: 12.1 mF cm-2 (at 1 mA cm-2)). For example, P(TEBTT/EDOT)-based frameworks synthesized in a 1:1 monomer-to-comonomer ratio show a ca. 35x capacitance improvement over PEDOT. The high areal capacitance measured for P(TEBTT/EDOT) copolymers can be explained by the open, highly porous hierarchical morphologies formed during the electropolymerization step. With >70% capacitance retention over 1,000 cycles (up to 89% achieved), both PTEBTT- and P(TEBTT/EDOT)-based frameworks are resilient to repeated electrochemical cycling and can be considered promising systems for high life cycle capacitive electrode applications.

Closely related species can provide valuable insights into evolutionary processes through comparison of their ecology, geographic distribution and the history recorded in their genomes. In the Indo-Pacific, many reef fishes are divided into sister species that come into secondary contact at biogeographic borders, most prominently where Indian Ocean and Pacific Ocean faunas meet. It is unclear whether hybridization in this contact zone represents incomplete speciation, secondary contact, an evolutionary dead-end (for hybrids) or some combination of the above. To address these issues, we conducted comprehensive surveys of two widely-distributed surgeonfish species, Acanthurus leucosternon (N = 141) and A. nigricans (N = 412), with mtDNA cytochrome b sequences and ten microsatellite loci. These surgeonfishes are found primarily in the Indian and Pacific Oceans, respectively, but overlap at the Christmas and Cocos-Keeling Islands hybrid zone in the eastern Indian Ocean. We also sampled the two other Pacific members of this species complex, A. achilles (N = 54) and A. japonicus (N = 49), which are known to hybridize with A. nigricans where their ranges overlap. Our results indicate separation between the four species that range from the recent Pleistocene to late Pliocene (235,000 to 2.25 million years ago). The Pacific A. achilles is the most divergent (and possibly ancestral) species with mtDNA dcorr ≈ 0.04, whereas the other two Pacific species (A. japonicus and A. nigricans) are distinguishable only at a population or subspecies level (ΦST = 0.6533, P < 0.001). Little population structure was observed within species, with evidence of recent population expansion across all four geographic ranges. We detected sharing of mtDNA haplotypes between species and extensive hybridization based on microsatellites, consistent with later generation hybrids but also the effects of allele homoplasy. Despite extensive introgression, 98% of specimens had concordance between mtDNA lineage and species identification based on external morphology, indicating that species integrity may not be eroding. The A. nigricans complex demonstrates a range of outcomes from incomplete speciation to secondary contact to decreasing hybridization with increasing evolutionary depth.

We present observations from an experimental study of the ablation of a free liquid surface promoted by a focused laser pulse, causing a rapid discharge of liquid in the form of a very thin conical-shaped sheet. In order to capture the dynamics, we employ a state-of-the-art ultra-high-speed video camera capable of capturing events at (Formula presented.) fps with shutter speeds down to 20 ns, whereby we were able to capture not only the ejecta sheet, but also the shock wave, emerging at speeds of up to 1.75 km/s, which is thus found to be hypersonic (Mach 5). Experiments were performed at a range of ambient pressures in order to study the effect of air drag on the evolution of the sheet, which was always observed to dome over, even at pressures as low as 3.8 kPa. At reduced pressures, the extended sheet evolution led to the formation of interference fringe patterns from which, by comparison with the opening speed of rupture, we were able to determine the ejecta thickness. © 2016, Springer-Verlag Berlin Heidelberg.

The regal angelfish (Pygoplites diacanthus; family Pomacanthidae) occupies reefs from the Red Sea to the central Pacific, with an Indian Ocean/Rea Sea color morph distinct from a Pacific Ocean morph. To assess population differentiation and evaluate the possibility of cryptic evolutionary partitions in this monotypic genus, we surveyed mtDNA cytochrome b and two nuclear introns (S7 and RAG2) in 547 individuals from 15 locations. Phylogeographic analyses revealed four mtDNA lineages (d = 0.006 – 0.015) corresponding to the Pacific Ocean, the Red Sea, and two admixed lineages in the Indian Ocean, a pattern consistent with known biogeographical barriers. Christmas Island in the eastern Indian Ocean had both Indian and Pacific lineages. Both S7 and RAG2 showed strong population-level differentiation between the Red Sea, Indian Ocean, and Pacific Ocean (ΦST = 0.066 – 0.512). The only consistent population sub-structure within these three regions was at the Society Islands (French Polynesia), where surrounding oceanographic conditions may reinforce isolation. Coalescence analyses indicate the Pacific (1.7 Ma) as the oldest extant lineage followed by the Red Sea lineage (1.4 Ma). Results from a median-joining network suggest radiations of two lineages from the Red Sea that currently occupy the Indian Ocean (0.7 – 0.9 Ma). Persistence of a Red Sea lineage through Pleistocene glacial cycles suggests a long-term refuge in this region. The affiliation of Pacific and Red Sea populations, apparent in cytochrome b and S7 (but equivocal in RAG2) raises the hypthosis that the Indian Ocean was recolonized from the Red Sea, possibly more than once. Assessing the genetic architecture of this widespread monotypic genus reveals cryptic evolutionary diversity that merits subspecific recognition.

We present new observations from an experimental investigation of the classical problem of the crown splash and sealing phenomena observed during the impact of spheres onto quiescent liquid pools. In the experiments, a 6 m tall vacuum chamber was used to provide the required ambient conditions from atmospheric pressure down to of an atmosphere, whilst high-speed videography was exploited to focus primarily on the above-surface crown formation and ensuing dynamics, paying particular attention to the moments just prior to the surface seal. In doing so, we have observed a buckling-type azimuthal instability of the crown. This instability is characterised by vertical striations along the crown, between which thin films form that are more susceptible to the air flow and thus are drawn into the closing cavity, where they atomize to form a fine spray within the cavity. To elucidate to the primary mechanisms and forces at play, we varied the sphere diameter, liquid properties and ambient pressure. Furthermore, a comparison between the entry of room-temperature spheres, where the contact line pins around the equator, and Leidenfrost spheres (i.e. an immersed superheated sphere encompassed by a vapour layer), where there is no contact line, indicates that the buckling instability appears in all crown sealing events, but is intensified by the presence of a pinned contact line. © 2016 Cambridge University Press.

We perform BiGlobal linear stability analysis on flow past a NACA0012 airfoil at 16° angle of attack and Reynolds number ranging from 400 to 1000. The steady-state two-dimensional base flows are computed using a well-tested finite difference code in combination with the selective frequency damping method. The base flow is characterized by two asymmetric recirculation bubbles downstream of the airfoil whose streamwise extent and the maximum reverse flow velocity increase with the Reynolds number. The stability analysis of the flow past the airfoil is carried out under very small spanwise wavenumber β = 10−4 to approximate the two-dimensional perturbation, and medium and large spanwise wavenumbers (β = 1–8) to account for the three-dimensional perturbation. Numerical results reveal that under small spanwise wavenumber, there are at most two oscillatory unstable modes corresponding to the near wake and far wake instabilities; the growth rate and frequency of the perturbation agree well with the two-dimensional direct numerical simulation results under all Reynolds numbers. For a larger spanwise wavenumber β = 1, there is only one oscillatory unstable mode associated with the wake instability at Re = 400 and 600, while at Re = 800 and 1000 there are two oscillatory unstable modes for the near wake and far wake instabilities, and one stationary unstable mode for the monotonically growing perturbation within the recirculation bubble via the centrifugal instability mechanism. All the unstable modes are weakened or even suppressed as the spanwise wavenumber further increases, among which the stationary mode persists until β = 4.