Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup

Background: Sexual dimorphism of body size has been the subject of numerous studies, but few have examined sexual shape dimorphism (SShD) and its evolution. Allometry, the shape change associated with size variation, has been suggested to be a main component of SShD. Yet little is known about the relative importance of the allometric and non-allometric components for the evolution of SShD. Results: We investigated sexual dimorphism in wing shape in the nine species of the Drosophila melanogaster subgroup. We used geometric morphometrics to characterise wing shape and found significant SShD in all nine species. The amount of shape difference and the diversity of the shape changes evolved across the group. However, mapping the divergence of SShD onto the phylogeny of the Drosophila melanogaster subgroup indicated that there is little phylogenetic signal. Finally, allometry accounted for a substantial part of SShD, but did not explain the bulk of evolutionary divergence in SShD because allometry itself was found to be evolutionarily plastic. Conclusion: SShD in the Drosophila wing can evolve rapidly and therefore shows only weak phylogenetic structure. The variable contribution of allometric and non-allometric components to the evolutionary divergence of SShD and the evolutionary plasticity of allometry suggest that SShD and allometry are influenced by a complex interaction of processes

Thermonuclear explosions of rapidly rotating white dwarfs - I. Deflagrations

Context: Turbulent deflagrations of Chandrasekhar mass White Dwarfs are commonly used to model Type Ia Supernova explosions. In this context, rapid rotation of the progenitor star is plausible but has so far been neglected. Aims: The aim of this work is to explore the influence of rapid rotation on the deflagration scenario. Methods: We use three dimensional hydrodynamical simulations to model turbulent deflagrations ignited within a variety of rapidly rotating CO WDs obeying rotation laws suggested by accretion studies. Results: We find that rotation has a significant impact on the explosion. The flame develops a strong anisotropy with a preferred direction towards the stellar poles, leaving great amounts of unburnt matter along the equatorial plane. Conclusions: The large amount of unburnt matter is contrary to observed spectral features of SNe Ia. Thus, rapid rotation of the progenitor star and the deflagration scenario are incompatible in order to explain SNe Ia.Comment: 13 pages, 10 figures, accepted for publication by A&

Computer simulations of electrorheological fluids in the dipole-induced dipole model

We have employed the multiple image method to compute the interparticle force for a polydisperse electrorheological (ER) fluid in which the suspended particles can have various sizes and different permittivites. The point-dipole (PD) approximation being routinely adopted in computer simulation of ER fluids is shown to err considerably when the particles approach and finally touch due to multipolar interactions. The PD approximation becomes even worse when the dielectric contrast between the particles and the host medium is large. From the results, we show that the dipole-induced-dipole (DID) model yields very good agreements with the multiple image results for a wide range of dielectric contrasts and polydispersity. As an illustration, we have employed the DID model to simulate the athermal aggregation of particles in ER fluids both in uniaxial and rotating fields. We find that the aggregation time is significantly reduced. The DID model accounts for multipolar interaction partially and is simple to use in computer simulation of ER fluids.Comment: 22 pages, 7 figures, submitted to Phys. Rev.

ARBITRARY ORDER FINITE VOLUME WELL-BALANCED SCHEMES FOR THE EULER EQUATIONS WITH GRAVITY

This work presents arbitrary high order well balanced finite volume schemes for the Euler equations with a prescribed gravitational field. It is assumed that the desired equilibrium solution is known, and we construct a scheme which is exactly well balanced for that particular equilibrium. The scheme is based on high order reconstructions of the fluctuations from equilibrium of density, velocity, and pressure, and on a well-balanced integration of the source terms, while no assumptions are needed on the numerical flux, beside consistency. This technique also allows one to construct well-balanced methods for a class of moving equilibria. Several numerical tests demonstrate the performance of the scheme on different scenarios, from equilibrium solutions to nonsteady problems involving shocks. The numerical tests are carried out with methods up to fifth order in one dimension, and third order accuracy in two dimensions

AnAll Speed SecondOrder IMEXRelaxation Scheme for the Euler Equations

We present an implicit-explicit finite volume scheme for the Euler equations. We start from the non-dimensionalised Euler equations where we split the pressure in a slow and a fast acoustic part. We use a Suliciu type relaxation model which we split in an explicit part, solved using a Godunov-type scheme based on an approximate Riemann solver, and an implicit part where we solve an elliptic equation for the fast pressure. The relaxation source terms are treated projecting the solution on the equilibrium manifold. The proposed scheme is positivity preserving with respect to the density and internal energy and asymptotic preserving towards the incompressible Euler equations. For this first order scheme we give a second order extension which maintains the positivity property. We perform numerical experiments in 1D and 2D to show the applicability of the proposed splitting and give convergence results for the second order extension

Applications of the Gauss-Bonnet theorem to gravitational lensing

In this geometrical approach to gravitational lensing theory, we apply the Gauss-Bonnet theorem to the optical metric of a lens, modelled as a static, spherically symmetric, perfect non-relativistic fluid, in the weak deflection limit. We find that the focusing of the light rays emerges here as a topological effect, and we introduce a new method to calculate the deflection angle from the Gaussian curvature of the optical metric. As examples, the Schwarzschild lens, the Plummer sphere and the singular isothermal sphere are discussed within this framework.Comment: 10 pages, 1 figure, IoP styl

An Algorithm for constructing Hjelmslev planes

Projective Hjelmslev planes and Affine Hjelmselv planes are generalisations of projective planes and affine planes. We present an algorithm for constructing a projective Hjelmslev planes and affine Hjelsmelv planes using projective planes, affine planes and orthogonal arrays. We show that all 2-uniform projective Hjelmslev planes, and all 2-uniform affine Hjelsmelv planes can be constructed in this way. As a corollary it is shown that all 2-uniform Affine Hjelmselv planes are sub-geometries of 2-uniform projective Hjelmselv planes.Comment: 15 pages. Algebraic Design Theory and Hadamard matrices, 2014, Springer Proceedings in Mathematics & Statistics 13