Spectra of graphene nanoribbons with armchair and zigzag boundary conditions

We study the spectral properties of the two-dimensional Dirac operator on bounded domains together with the appropriate boundary conditions which provide a (continuous) model for graphene nanoribbons. These are of two types, namely, the so-called armchair and zigzag boundary conditions, depending on the line along which the material was cut. In the former case, we show that the spectrum behaves in what might be called a classical way; while in the latter, we prove the existence of a sequence of finite multiplicity eigenvalues converging to zero and which correspond to edge states.

On the similarity of Sturm-Liouville operators with non-Hermitian boundary conditions to self-adjoint and normal operators

We consider one-dimensional Schrödinger-type operators in a bounded interval with non-self-adjoint Robin-type boundary conditions. It is well known that such operators are generically conjugate to normal operators via a similarity transformation. Motivated by recent interests in quasi-Hermitian Hamiltonians in quantum mechanics, we study properties of the transformations and similar operators in detail. In the case of parity and time reversal boundary conditions, we establish closed integral-type formulae for the similarity transformations, derive a non-local self-adjoint operator similar to the Schrödinger operator and also find the associated “charge conjugation” operator, which plays the role of fundamental symmetry in a Krein-space reformulation of the problem.

Diagenesis of a light, tight-oil chert reservoir at the Ediacaran/Cambrian boundary, Sultanate of Oman

The Al Shomou Silicilyte Member (Athel Formation) in the South Oman Salt Basin shares many of the characteristics of a light, tight-oil (LTO) reservoir: it is a prolifi c source rock mature for light oil, it produces light oil from a very tight matrix and reservoir, and hydraulic fracking technology is required to produce the oil. What is intriguing about the Al Shomou Silicilyte, and different from other LTO reservoirs, is its position related to the Precambrian/Cambrian Boundary (PCB) and the fact that it is a ‘laminated chert‘ rather than a shale. In an integrated diagenetic study we applied microstructural analyses (SEM, BSE) combined with state-of-the-art stable isotope and trace element analysis of the silicilyte matrix and fractures. Fluid inclusion microthermometry was applied to record the salinity and minimum trapping temperatures. The microstructural investigations reveal a fi ne lamination of the silicilyte matrix with a mean lamina thickness of ca. 20 μm consisting of predominantly organic matter-rich and fi nely crystalline quartz-rich layers, respectively. Authigenic, micron-sized idiomorphic quartz crystals are the main matrix components of the silicilyte. Other diagenetic phases are pyrite, apatite, dolomite, magnesite and barite cements. Porosity values based on neutron density logs and core plug data indicate porosity in the silicilyte ranges from less than 2% to almost to 40%. The majority of the pore space in the silicilyte is related to (primary) inter-crystalline pores, with locally important oversized secondary pores. Pore casts of the silica matrix show that pores are extremely irregular in three dimensions, and are generally interconnected by a complex web or meshwork of fi ne elongate pore throats. Mercury injection capillary data are in line with the microstructural observations suggesting two populations of pore throats, with an effective average modal diameter of 0.4 μm. The acquired geochemical data support the interpretation that the primary source of the silica is the ambient seawater rather than hydrothermal or biogenic. A maximum temperature of ca. 45°C for the formation of microcrystalline quartz in the silicilyte is good evidence that the lithifi cation and crystallization of quartz occurred in the fi rst 5 Ma after deposition. Several phases of brittle fracturing and mineralization occurred in response to salt tectonics during burial. The sequences of fracture-fi lling mineral phases (dolomite - layered chalcedony – quartz – apatite - magnesite I+II - barite – halite) indicates a complex fl uid evolution after silicilyte lithifi cation. Primary, all-liquid fl uid inclusions in the fracturefi lling quartz are good evidence of growth beginning at low temperatures, i.e. ≤ 50ºC. Continuous precipitation during increasing temperature and burial is documented by primary two-phase fl uid inclusions in quartz cements that show brines at 50°C and fi rst hydrocarbons at ca. 70°C. The absolute timing of each mineral phase can be constrained based on U-Pb geochronometry, and basin modelling. Secondary fl uid inclusions in quartz, magnesite and barite indicate reactivation of the fracture system after peak burial temperature during the major cooling event, i.e. uplift, between 450 and 310 Ma. A number of fi rst-order trends in porosity and reservoir-quality distribution are observed which are strongly related to the diagenetic and fl uid history of the reservoir: the early in-situ generation of hydrocarbons and overpressure development arrests diagenesis and preserves matrix porosity. Chemical compaction by pressure dissolution in the fl ank areas could be a valid hypothesis to explain the porosity variations in the silicilitye slabs resulting in lower porosity and poorer connectivity on the fl anks of the reservoir. Most of the hydrocarbon storage and production comes from intervals characterized by Amthor et al. 114488 preserved micropores, not hydrocarbon storage in a fracture system. The absence of oil expulsion results in present-day high oil saturations. The main diagenetic modifi cations of the silicilyte occurred and were completed relatively early in its history, i.e. before 300 Ma. An instrumental factor for preserving matrix porosity is the diffi culty for a given slab to evacuate all the fl uids (water and hydrocarbons), or in other words, the very good sealing capacity of the salt embedding the slab.

FEM/FD Immersed Boundary FSI Simulations

Immersed boundary simulations have been under development for physiological flows, allowing for elegant handling of fluid-structure interaction modelling with large deformations due to retained domain-specific meshing. We couple a structural system in Lagrangian representation that is formulated in a weak form with a Navier-Stokes system discretized through a finite differences scheme. We build upon a proven highly scalable imcompressible flow solver that we extend to handle FSI. We aim at applying our method to investigating the hemodynamics of Aortic Valves. The code is going to be extended to conform to the new hybrid-node supercomputers.

A class of exact Navier-Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

We introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.

Stabilization of subcritical bypass transition in the leading-edge boundary layer by suction

This work investigates the subcritical spatial transition in the swept Hiemenz boundary layer by means of direct numerical simulations (DNS). A pair of steady co-rotating vortices located at the attachment line is enforced as a primary disturbance leading to streaks which are stable. A small secondary, time-dependent disturbance interacts with these streaks such that instability and breakdown to turbulence may occur. The instability only occurs for a certain band of secondary disturbance frequencies. Positive secondary instability growth rates could be observed for Reynolds numbers as low as , whereas the linear critical Reynolds number is. Uniform wall suction is shown to stabilise this transition mechanism, analogously to results from linear stability theory. The effects of suction on the formation of primary streaks and on the secondary growth rate are decoupled. For streaks of different suction whose amplitude is held constant by adjusting the Reynolds number, the suction is shown to increase the growth rate of the secondary instability. The stabilising influence of wall suction consists in decreasing the streak amplitude only. Depending on the Reynolds number and the suction strength, breakdown may either occur locally and may be convected along the far-field streamlines, or occur globally and cover broad regions in the downstream direction.

A deep crust–mantle boundary in the asteroid 4 Vesta

The asteroid 4 Vesta was recently found to have two large impact craters near its south pole, exposing subsurface material. Modelling suggested that surface material in the northern hemisphere of Vesta came from a depth of about 20 kilometres, whereas the exposed southern material comes from a depth of 60 to 100 kilometres. Large amounts of olivine from the mantle were not seen, suggesting that the outer 100 kilometres or so is mainly igneous crust. Here we analyse the data on Vesta and conclude that the crust–mantle boundary (or Moho) is deeper than 80 kilometres.

Topology of the ambient boundary and the convergence of causal curves

We discuss the topological nature of the boundary spacetime, the conformal infinity of the ambient cosmological metric. Due to the existence of a homothetic group, the bounding spacetime must be equipped not with the usual Euclidean metric topology but with the Zeeman fine topology. This then places severe constraints to the convergence of a sequence of causal curves and the extraction of a limit curve, and also to our ability to formulate conditions for singularity formation.

Kelvin-Helmholtz instabilities at the magnetic cavity boundary of comet 67P/Churyumov-Gerasimenko

We investigate the plasma environment of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. Rosetta will rendezvous with the comet in 2014 at almost 3.5 AU and follow it all the way to and past perihelion at 1.3 AU. During its journey towards the inner solar system the comet's environment will significantly change. The interaction of the solar wind with a well developed neutral coma leads to the formation of an upstream bow shock and, closer to the comet, the inner shock separating the solar wind, with cometary pick-up ions mass-loaded, from the inner cometary ions which are dragged outward through abundant collisions and charge exchange with the expanding neutral gas. As a consequence the interplanetary magnetic field is prevented from penetrating the innermost region of the comet, the so-called magnetic cavity. We use our magnetohydrodynamics model BATSRUS (Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme) to simulate the solar wind - comet interaction. The model includes photoionization, ion-electron recombination, and charge exchange. Under certain conditions our model predicts an unstable plasma flow at the inner shock. We show that the plasma shear flow around the magnetic cavity can lead to Kelvin-Helmholtz instabilities. We investigate the onset of this phenomenon with change of heliocentric distance and furthermore show that a previously stable magnetic cavity boundary can become unstable when the neutral gas is predominately released from the dayside of the comet.

Study of mineral dust entrainment in the planetary boundary layer by lidar depolarisation technique

Measurements on 27 June 2011 were performed over the Southern Iberian Peninsula at Granada EARLINET station, using active and passive remote sensing and airborne and surface in-situ data in order to study the entrainment processes between aerosols in the free troposphere and those in the planetary boundary layer (PBL). To this aim the temporal evolution of the lidar depolarisation, backscatter-related Angström exponent and potential temperature profiles were used in combination with the PBL contribution to the aerosol optical depth (AOD). Our results show that the mineral dust entrainment in the PBL was caused by the convective processes which ‘trapped’ the lofted mineral dust layer, distributing the mineral dust particles within the PBL. The temporal evolution of ground-based in-situ data evidenced the impact of this process at surface level. Finally, the amount of mineral dust in the atmospheric column available to be dispersed into the PBL was estimated by means of POLIPHON (Polarizing Lidar Photometer Networking). The dust mass concentration derived from POLIPHON was compared with the coarse-mode mass concentration retrieved with airborne in-situ measurements. Comparison shows differences below 50 µg/m³ (30% relative difference) indicating a relative good agreement between both techniques.

Comparison of Mixed and Kinematic Uniform Boundary Conditions in Homogenized Elasticity of Femoral Trabecular Bone Using Microfinite Element Analyses

Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Micro-finite element (microFE) analysis allows computing the apparent elastic properties of trabecular bone biopsies, but the results depend on the type of applied boundary conditions (BCs). In this study, 167 femoral trabecular cubic biopsies with a side length of 5.3 mm were analyzed using microFE analysis to compare their stiffness systematically with kinematic uniform boundary conditions (KUBCs) and periodicity-compatible mixed uniform boundary conditions (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset-Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length. Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset-Curnier model predicted apparent elastic constants successfully in both cases with adjusted coefficients of determination of 0.986 for KUBCs and 0.975 for PMUBCs. The role of these boundary conditions in finite element analyses of whole bones and bone-implant systems will need to be investigated in future work.