Generalized quasiperiodic Rauzy tilings

We present a geometrical description of new canonical $d$-dimensional codimension one quasiperiodic tilings based on generalized Fibonacci sequences. These tilings are made up of rhombi in 2d and rhombohedra in 3d as the usual Penrose and icosahedral tilings. Thanks to a natural indexing of the sites according to their local environment, we easily write down, for any approximant, the sites coordinates, the connectivity matrix and we compute the structure factor.Comment: 11 pages, 3 EPS figures, final version with minor change

Easy plane baby skyrmions

The baby Skyrme model is studied with a novel choice of potential, $V=1/2 \phi_3^2$. This "easy plane" potential vanishes at the equator of the target two-sphere. Hence, in contrast to previously studied cases, the boundary value of the field breaks the residual SO(2) internal symmetry of the model. Consequently, even the unit charge skyrmion has only discrete symmetry and consists of a bound state of two half lumps. A model of long-range inter-skyrmion forces is developed wherein a unit skyrmion is pictured as a single scalar dipole inducing a massless scalar field tangential to the vacuum manifold. This model has the interesting feature that the two-skyrmion interaction energy depends only on the average orientation of the dipoles relative to the line joining them. Its qualitative predictions are confirmed by numerical simulations. Global energy minimizers of charges B=1,...,14,18,32 are found numerically. Up to charge B=6, the minimizers have 2B half lumps positioned at the vertices of a regular 2B-gon. For charges B >= 7, rectangular or distorted rectangular arrays of 2B half lumps are preferred, as close to square as possible.Comment: v3: replaced with journal version, one new reference, one deleted reference; 8 pages, 5 figures v2: fixed some typos and clarified the relationship with condensed matter systems 8 pages, 5 figure

Adaptive finite element method assisted by stochastic simulation of chemical systems

Stochastic models of chemical systems are often analysed by solving the corresponding\ud Fokker-Planck equation which is a drift-diffusion partial differential equation for the probability\ud distribution function. Efficient numerical solution of the Fokker-Planck equation requires adaptive mesh refinements. In this paper, we present a mesh refinement approach which makes use of a stochastic simulation of the underlying chemical system. By observing the stochastic trajectory for a relatively short amount of time, the areas of the state space with non-negligible probability density are identified. By refining the finite element mesh in these areas, and coarsening elsewhere, a suitable mesh is constructed and used for the computation of the probability density

Low zinc status and absorption exist in infants with jejunostomies or ileostomies which persists after intestinal repair.

There is very little data regarding trace mineral nutrition in infants with small intestinal ostomies. Here we evaluated 14 infants with jejunal or ileal ostomies to measure their zinc absorption and retention and biochemical zinc and copper status. Zinc absorption was measured using a dual-tracer stable isotope technique at two different time points when possible. The first study was conducted when the subject was receiving maximal tolerated feeds enterally while the ostomy remained in place. A second study was performed as soon as feasible after full feeds were achieved after intestinal repair. We found biochemical evidence of deficiencies of both zinc and copper in infants with small intestinal ostomies at both time points. Fractional zinc absorption with an ostomy in place was 10.9% ± 5.3%. After reanastamosis, fractional zinc absorption was 9.4% ± 5.7%. Net zinc balance was negative prior to reanastamosis. In conclusion, our data demonstrate that infants with a jejunostomy or ileostomy are at high risk for zinc and copper deficiency before and after intestinal reanastamosis. Additional supplementation, especially of zinc, should be considered during this time period

Improving Performance of Iterative Methods by Lossy Checkponting

Iterative methods are commonly used approaches to solve large, sparse linear systems, which are fundamental operations for many modern scientific simulations. When the large-scale iterative methods are running with a large number of ranks in parallel, they have to checkpoint the dynamic variables periodically in case of unavoidable fail-stop errors, requiring fast I/O systems and large storage space. To this end, significantly reducing the checkpointing overhead is critical to improving the overall performance of iterative methods. Our contribution is fourfold. (1) We propose a novel lossy checkpointing scheme that can significantly improve the checkpointing performance of iterative methods by leveraging lossy compressors. (2) We formulate a lossy checkpointing performance model and derive theoretically an upper bound for the extra number of iterations caused by the distortion of data in lossy checkpoints, in order to guarantee the performance improvement under the lossy checkpointing scheme. (3) We analyze the impact of lossy checkpointing (i.e., extra number of iterations caused by lossy checkpointing files) for multiple types of iterative methods. (4)We evaluate the lossy checkpointing scheme with optimal checkpointing intervals on a high-performance computing environment with 2,048 cores, using a well-known scientific computation package PETSc and a state-of-the-art checkpoint/restart toolkit. Experiments show that our optimized lossy checkpointing scheme can significantly reduce the fault tolerance overhead for iterative methods by 23%~70% compared with traditional checkpointing and 20%~58% compared with lossless-compressed checkpointing, in the presence of system failures.Comment: 14 pages, 10 figures, HPDC'1

Black Hole-Neutron Star Binaries in General Relativity: Quasiequilibrium Formulation

We present a new numerical method for the construction of quasiequilibrium models of black hole-neutron star binaries. We solve the constraint equations of general relativity, decomposed in the conformal thin-sandwich formalism, together with the Euler equation for the neutron star matter. We take the system to be stationary in a corotating frame and thereby assume the presence of a helical Killing vector. We solve these coupled equations in the background metric of a Kerr-Schild black hole, which accounts for the neutron star's black hole companion. In this paper we adopt a polytropic equation of state for the neutron star matter and assume large black hole--to--neutron star mass ratios. These simplifications allow us to focus on the construction of quasiequilibrium neutron star models in the presence of strong-field, black hole companions. We summarize the results of several code tests, compare with Newtonian models, and locate the onset of tidal disruption in a fully relativistic framework.Comment: 17 pages, 7 figures; added discussion, tables; PRD in pres

Comparison of experiments on bio-inspired hover kinematics with the unsteady vortex model and CFD

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106510/1/AIAA2013-66.pd

How universal is the fractional-quantum-Hall edge Luttinger liquid?

This article reports on our microscopic investigations of the edge of the fractional quantum Hall state at filling factor $\nu=1/3$. We show that the interaction dependence of the wave function is well described in an approximation that includes mixing with higher composite-fermion Landau levels in the lowest order. We then proceed to calculate the equal time edge Green function, which provides evidence that the Luttinger exponent characterizing the decay of the Green function at long distances is interaction dependent. The relevance of this result to tunneling experiments is discussed.Comment: 5 page

Full sphere hydrodynamic and dynamo benchmarks

Convection in planetary cores can generate fluid flow and magnetic fields, and a number of sophisticated codes exist to simulate the dynamic behaviour of such systems. We report on the first community activity to compare numerical results of computer codes designed to calculate fluid flow within a whole sphere. The flows are incompressible and rapidly rotating and the forcing of the flow is either due to thermal convection or due to moving boundaries. All problems defined have solutions that allow easy comparison, since they are either steady, slowly drifting or perfectly periodic. The first two benchmarks are defined based on uniform internal heating within the sphere under the Boussinesq approximation with boundary conditions that are uniform in temperature and stress-free for the flow. Benchmark 1 is purely hydrodynamic, and has a drifting solution. Benchmark 2 is a magnetohydrodynamic benchmark that can generate oscillatory, purely periodic, flows and magnetic fields. In contrast, Benchmark 3 is a hydrodynamic rotating bubble benchmark using no slip boundary conditions that has a stationary solution. Results from a variety of types of code are reported, including codes that are fully spectral (based on spherical harmonic expansions in angular coordinates and polynomial expansions in radius), mixed spectral and finite difference, finite volume, finite element and also a mixed Fourier–finite element code. There is good agreement between codes. It is found that in Benchmarks 1 and 2, the approximation of a whole sphere problem by a domain that is a spherical shell (a sphere possessing an inner core) does not represent an adequate approximation to the system, since the results differ from whole sphere results