Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences recent issues
http://rspa.royalsocietypublishing.org
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences RSS feed -- recent issues1471-2946Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences1364-5021<![CDATA[Ducks in space: from nonlinear absolute instability to noise-sustained structures in a pattern-forming system]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170018?rss=1
A subcritical pattern-forming system with nonlinear advection in a bounded domain is recast as a slow–fast system in space and studied using a combination of geometric singular perturbation theory and numerical continuation. Two types of solutions describing the possible location of stationary fronts are identified, whose origin is traced to the onset of convective and absolute instability when the system is unbounded. The former are present only for non-zero upstream boundary conditions and provide a quantitative understanding of noise-sustained structures in systems of this type. The latter correspond to the onset of a global mode and are present even with zero upstream boundary conditions. The role of canard trajectories in the nonlinear transition between these states is clarified and the stability properties of the resulting spatial structures are determined. Front location in the convective regime is highly sensitive to the upstream boundary condition, and its dependence on this boundary condition is studied using a combination of numerical continuation and Monte Carlo simulations of the partial differential equation. Statistical properties of the system subjected to random or stochastic boundary conditions at the inlet are interpreted using the deterministic slow–fast spatial dynamical system.
]]>2017-11-08T00:05:22-08:00info:doi/10.1098/rspa.2017.0018hwp:master-id:royprsa;rspa.2017.00182017-11-08Research articles47322072017001820170018<![CDATA[Three-dimensional quasi-periodic shifted Green function throughout the spectrum, including Wood anomalies]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170242?rss=1
This work, part II in a series, presents an efficient method for evaluation of wave scattering by doubly periodic diffraction gratings at or near what are commonly called ‘Wood anomaly frequencies’. At these frequencies, there is a grazing Rayleigh wave, and the quasi-periodic Green function ceases to exist. We present a modification of the Green function by adding two types of terms to its lattice sum. The first type are transversely shifted Green functions with coefficients that annihilate the growth in the original lattice sum and yield algebraic convergence. The second type are quasi-periodic plane wave solutions of the Helmholtz equation which reinstate certain necessary grazing modes without leading to blow-up at Wood anomalies. Using the new quasi-periodic Green function, we establish, for the first time, that the Dirichlet problem of scattering by a smooth doubly periodic scattering surface at a Wood frequency is uniquely solvable. We also present an efficient high-order numerical method based on this new Green function for scattering by doubly periodic surfaces at and around Wood frequencies. We believe this is the first solver able to handle Wood frequencies for doubly periodic scattering problems in three dimensions. We demonstrate the method by applying it to acoustic scattering.
]]>2017-11-01T02:09:49-07:00info:doi/10.1098/rspa.2017.0242hwp:master-id:royprsa;rspa.2017.02422017-11-01Research articles47322072017024220170242<![CDATA[A post Gurney quantum mechanical perspective on the electrolysis of water: ion neutralization in solution]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170371?rss=1
Electron fluxes crossing the interface between a metallic conductor and an aqueous environment are important in many fields; hydrogen production, environmental scanning tunnelling microscopy, scanning electrochemical microscopy being some of them. Gurney (Gurney 1931 Proc. R. Soc. Lond.134, 137 (doi:10.1098/rspa.1931.0187)) provided in 1931 a scheme for tunnelling during electrolysis and outlined conditions for it to occur. We measure the low-voltage current flows between gold electrodes in pure water and use the time-dependent behaviour at voltage switch-on and switch-off to evaluate the relative contribution to the steady current arising from tunnelling of electrons between the electrodes and ions in solution and from the neutralization of ions adsorbed onto the electrode surface. We ascribe the larger current contribution to quantum tunnelling of electrons to and from ions in solution near the electrodes. We refine Gurney's barrier scheme to include solvated electron states and quantify energy differences using updated information. We show that Gurney's conditions would prevent the current flow at low voltages we observe but outline how the ideas of Marcus (Marcus 1956 J. Chem. Phys.24, 966–978 (doi:10.1063/1.1742723)) concerning solvation fluctuations enable the condition to be relaxed. We derive an average barrier tunnelling model and a multiple pathways tunnelling model and compare predictions with measurements of the steady-state current–voltage relation. The tunnelling barrier was found to be wide and low in agreement with other experimental studies. Applications as a biosensing mechanism are discussed that exploit the fast tunnelling pathways along molecules in solution.
]]>2017-11-08T00:05:22-08:00info:doi/10.1098/rspa.2017.0371hwp:master-id:royprsa;rspa.2017.03712017-11-08Research articles47322072017037120170371<![CDATA[Force direction patterns promote whole body stability even in hip-flexed walking, but not upper body stability in human upright walking]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170404?rss=1
Directing the ground reaction forces to a focal point above the centre of mass of the whole body promotes whole body stability in human and animal gaits similar to a physical pendulum. Here we show that this is the case in human hip-flexed walking as well. For all upper body orientations (upright, 25°, 50°, maximum), the focal point was well above the centre of mass of the whole body, suggesting its general relevance for walking. Deviations of the forces' lines of action from the focal point increased with upper body inclination from 25 to 43 mm root mean square deviation (RMSD). With respect to the upper body in upright gait, the resulting force also passed near a focal point (17 mm RMSD between the net forces' lines of action and focal point), but this point was 18 cm below its centre of mass. While this behaviour mimics an unstable inverted pendulum, it leads to resulting torques of alternating sign in accordance with periodic upper body motion and probably provides for low metabolic cost of upright gait by keeping hip torques small. Stabilization of the upper body is a consequence of other mechanisms, e.g. hip reflexes or muscle preflexes.
]]>2017-11-08T00:05:22-08:00info:doi/10.1098/rspa.2017.0404hwp:master-id:royprsa;rspa.2017.04042017-11-08Research articles47322072017040420170404<![CDATA[Age distribution dynamics with stochastic jumps in mortality]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170451?rss=1
While deterministic age distribution models have been extensively studied and applied in various disciplines, little work has been devoted to understanding the role of stochasticity in birth and mortality terms. In this paper, we analyse a stochastic M’Kendrick–von Foerster equation in which jumps in mortality represent intense losses of population due to external events. We present explicit solutions for the probability density functions of the age distribution and the total population and for the temporal dynamics of their moments. We also derive the dynamics of the mean age of the population and its harmonic mean. The framework is then used to calculate the age distribution of salt in the soil root zone, where the accumulation of salt by atmospheric deposition is counteracted by plant uptake and by jump losses due to percolation events.
]]>2017-11-01T02:09:49-07:00info:doi/10.1098/rspa.2017.0451hwp:master-id:royprsa;rspa.2017.04512017-11-01Research articles47322072017045120170451<![CDATA[Swelling-induced wrinkling in layered gel beams]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170454?rss=1
Gels are widely employed in smart mechanical devices and biomedical applications. Swelling-induced bending actuation can be obtained by means of a simple bilayer gel beam. We show that this system can also exhibit wrinkling patterns of potential interest for structural morphing and sensing. We study swelling-induced wrinkling at the extrados of a bilayer gel beam with the softer layer on top. The bent configuration at finite strain is recovered first and, starting from it, a linear perturbation analysis is performed. We delimit the zone corresponding to wrinkling modes in a parameter plane encompassing a mechanical stiffness ratio and a geometric top layer to total height ratio. Interestingly, we observe that surface instability precedes and envelopes wrinkling modes of finite wavelength. Finally, we discuss the effect of changes in stiffness and of the Flory–Huggins parameters on the size of the wrinkling domain.
]]>2017-11-01T02:09:49-07:00info:doi/10.1098/rspa.2017.0454hwp:master-id:royprsa;rspa.2017.04542017-11-01Research articles47322072017045420170454<![CDATA[On the buckling of an elastic holey column]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170477?rss=1
We report the results of a numerical and theoretical study of buckling in elastic columns containing a line of holes. Buckling is a common failure mode of elastic columns under compression, found over scales ranging from metres in buildings and aircraft to tens of nanometers in DNA. This failure usually occurs through lateral buckling, described for slender columns by Euler’s theory. When the column is perforated with a regular line of holes, a new buckling mode arises, in which adjacent holes collapse in orthogonal directions. In this paper, we firstly elucidate how this alternate hole buckling mode coexists and interacts with classical Euler buckling modes, using finite-element numerical calculations with bifurcation tracking. We show how the preferred buckling mode is selected by the geometry, and discuss the roles of localized (hole-scale) and global (column-scale) buckling. Secondly, we develop a novel predictive model for the buckling of columns perforated with large holes. This model is derived without arbitrary fitting parameters, and quantitatively predicts the critical strain for buckling. We extend the model to sheets perforated with a regular array of circular holes and use it to provide quantitative predictions of their buckling.
]]>2017-11-15T00:05:25-08:00info:doi/10.1098/rspa.2017.0477hwp:master-id:royprsa;rspa.2017.04772017-11-15Research articles47322072017047720170477<![CDATA[Bidispersive vertical convection]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170481?rss=1
A bidispersive porous material is one which has usual pores but additionally contains a system of micro pores. We consider a fluid-saturated bidispersive porous medium in the vertical layer x(–1/2,1/2) with gravity in the –z (downward) direction. The walls of the layer are maintained at different constant temperatures. A suitable Rayleigh number is defined and we derive a global stability threshold below which no instability may arise. We additionally show that the porous layer is stable for all Rayleigh numbers provided the initial temperature gradient is bounded in a precise sense.
]]>2017-11-01T02:09:49-07:00info:doi/10.1098/rspa.2017.0481hwp:master-id:royprsa;rspa.2017.04812017-11-01Research articles47322072017048120170481<![CDATA[Analysing causal structures with entropy]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170483?rss=1
A central question for causal inference is to decide whether a set of correlations fits a given causal structure. In general, this decision problem is computationally infeasible and hence several approaches have emerged that look for certificates of compatibility. Here, we review several such approaches based on entropy. We bring together the key aspects of these entropic techniques with unified terminology, filling several gaps and establishing new connections, all illustrated with examples. We consider cases where unobserved causes are classical, quantum and post-quantum, and discuss what entropic analyses tell us about the difference. This difference has applications to quantum cryptography, where it can be crucial to eliminate the possibility of classical causes. We discuss the achievements and limitations of the entropic approach in comparison to other techniques and point out the main open problems.
]]>2017-11-01T03:11:33-07:00info:doi/10.1098/rspa.2017.0483hwp:master-id:royprsa;rspa.2017.04832017-11-01Review articles47322072017048320170483<![CDATA[Multi-scale mechanics of granular solids from grain-resolved X-ray measurements]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170491?rss=1
This work discusses an experimental technique for studying the mechanics of three-dimensional (3D) granular solids. The approach combines 3D X-ray diffraction and X-ray computed tomography to measure grain-resolved strains, kinematics and contact fabric in the bulk of a granular solid, from which continuum strains, grain stresses, interparticle forces and coarse-grained elasto-plastic moduli can be determined. We demonstrate the experimental approach and analysis of selected results on a sample of 1099 stiff, frictional grains undergoing multiple uniaxial compression cycles. We investigate the inter-particle force network, elasto-plastic moduli and associated length scales, reversibility of mechanical responses during cyclic loading, the statistics of microscopic responses and microstructure–property relationships. This work serves to highlight both the fundamental insight into granular mechanics that is furnished by combined X-ray measurements and describes future directions in the field of granular materials that can be pursued with such approaches.
]]>2017-11-01T02:09:49-07:00info:doi/10.1098/rspa.2017.0491hwp:master-id:royprsa;rspa.2017.04912017-11-01Research articles47322072017049120170491<![CDATA[On the analysis of the double Hopf bifurcation in machining processes via centre manifold reduction]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170502?rss=1
The single-degree-of-freedom model of orthogonal cutting is investigated to study machine tool vibrations in the vicinity of a double Hopf bifurcation point. Centre manifold reduction and normal form calculations are performed to investigate the long-term dynamics of the cutting process. The normal form of the four-dimensional centre subsystem is derived analytically, and the possible topologies in the infinite-dimensional phase space of the system are revealed. It is shown that bistable parameter regions exist where unstable periodic and, in certain cases, unstable quasi-periodic motions coexist with the equilibrium. Taking into account the non-smoothness caused by loss of contact between the tool and the workpiece, the boundary of the bistable region is also derived analytically. The results are verified by numerical continuation. The possibility of (transient) chaotic motions in the global non-smooth dynamics is shown.
]]>2017-11-15T00:05:25-08:00info:doi/10.1098/rspa.2017.0502hwp:master-id:royprsa;rspa.2017.05022017-11-15Research articles47322072017050220170502<![CDATA[Advection and Taylor-Aris dispersion in rivulet flow]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170524?rss=1
Motivated by the need for a better understanding of the transport of solutes in microfluidic flows with free surfaces, the advection and dispersion of a passive solute in steady unidirectional flow of a thin uniform rivulet on an inclined planar substrate driven by gravity and/or a uniform longitudinal surface shear stress are analysed. Firstly, we describe the short-time advection of both an initially semi-infinite and an initially finite slug of solute of uniform concentration. Secondly, we describe the long-time Taylor–Aris dispersion of an initially finite slug of solute. In particular, we obtain the general expression for the effective diffusivity for Taylor–Aris dispersion in such a rivulet, and discuss in detail its different interpretations in the special case of a rivulet on a vertical substrate.
]]>2017-11-15T00:05:27-08:00info:doi/10.1098/rspa.2017.0524hwp:master-id:royprsa;rspa.2017.05242017-11-15Research articles47322072017052420170524<![CDATA[Symmetric rearrangement of groundwater-fed streams]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170539?rss=1
Streams shape landscapes through headward growth and lateral migration. When these streams are primarily fed by groundwater, recent work suggests that their tips advance to maximize the symmetry of the local Laplacian field associated with groundwater flow. We explore the extent to which such forcing is responsible for the lateral migration of streams by studying two features of groundwater-fed streams in Bristol, Florida: their confluence angle near junctions and their curvature. First, we find that, while streams asymptotically form a 72° angle near their tips, they simultaneously exhibit a wide 120° confluence angle within approximately 10 m of their junctions. We show that this wide angle maximizes the symmetry of the groundwater field near the junction. Second, we argue that streams migrate laterally within valleys and present a new spectral analysis method to relate planform curvature to the surrounding groundwater field. Our results suggest that streams migrate laterally in response to fluxes from the surrounding groundwater table, providing evidence of a new mechanism that complements Laplacian growth at their tips.
]]>2017-11-08T00:59:42-08:00info:doi/10.1098/rspa.2017.0539hwp:master-id:royprsa;rspa.2017.05392017-11-08Research articles47322072017053920170539<![CDATA[Lagrangian averaging with geodesic mean]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170558?rss=1
This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler-α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.
]]>2017-11-15T00:43:18-08:00info:doi/10.1098/rspa.2017.0558hwp:master-id:royprsa;rspa.2017.05582017-11-15Research articles47322072017055820170558<![CDATA[Nonlinear mode interaction in equal-leg angle struts susceptible to cellular buckling]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170583?rss=1
A variational model that describes the interactive buckling of a thin-walled equal-leg angle strut under pure axial compression is presented. A formulation combining the Rayleigh–Ritz method and continuous displacement functions is used to derive a system of differential and integral equilibrium equations for the structural component. Solving the equations using numerical continuation reveals progressive cellular buckling (or snaking) arising from the nonlinear interaction between the weak-axis flexural buckling mode and the strong-axis flexural–torsional buckling mode for the first time—the resulting behaviour being highly unstable. Physical experiments conducted on 10 cold-formed steel specimens are presented and the results show good agreement with the variational model.
]]>2017-11-01T07:43:36-07:00info:doi/10.1098/rspa.2017.0583hwp:master-id:royprsa;rspa.2017.05832017-11-01Research articles47322072017058320170583<![CDATA[Continuum modelling of pantographic sheets for out-of-plane bifurcation and vibrational analysis]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2207/20170636?rss=1
A nonlinear two-dimensional (2D) continuum with a latent internal structure is introduced as a coarse model of a plane network of beams which, in turn, is assumed as a model of a pantographic structure made up by two families of equispaced beams, superimposed and connected by pivots. The deformation measures of the beams of the network and that of the 2D body are introduced and the former are expressed in terms of the latter by making some kinematical assumptions. The expressions for the strain and kinetic energy densities of the network are then introduced and given in terms of the kinematic quantities of the 2D continuum. To account for the modelling abilities of the 2D continuum in the linear range, the eigenmode and eigenfrequencies of a given specimen are determined. The buckling and post-buckling behaviour of the same specimen, subjected to two different loading conditions are analysed as tests in the nonlinear range. The problems have been solved numerically by means of the COMSOL Multiphysics finite element software.
]]>2017-11-15T03:22:36-08:00info:doi/10.1098/rspa.2017.0636hwp:master-id:royprsa;rspa.2017.06362017-11-15Research articles47322072017063620170636<![CDATA[Critical ignition conditions in exothermically reacting systems: first-order reactions]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170145?rss=1
In this paper, the comparative analysis of the thermal explosion (TE) critical conditions on the planes temperature–conversion degree and temperature–time was conducted. It was established that the ignition criteria are almost identical only at relatively small values of Todes parameter. Otherwise, the results of critical conditions analysis on the plane temperature–conversion degree may be wrong. The asymptotic method of critical conditions calculation for the first-order reactions was proposed (taking into account the reactant consumption). The degeneration conditions of TE were determined. The calculation of critical conditions for specific first-order reaction was made. The comparison of the analytical results obtained with the results of numerical calculations and experimental data showed that they are in good agreement.
]]>2017-10-25T00:05:19-07:00info:doi/10.1098/rspa.2017.0145hwp:master-id:royprsa;rspa.2017.01452017-10-25Research articles47322062017014520170145<![CDATA[Extended lubrication theory: improved estimates of flow in channels with variable geometry]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170234?rss=1
Lubrication theory is broadly applicable to the flow characterization of thin fluid films and the motion of particles near surfaces. We offer an extension to lubrication theory by starting with Stokes equations and considering higher-order terms in a systematic perturbation expansion to describe the fluid flow in a channel with features of a modest aspect ratio. Experimental results qualitatively confirm the higher-order analytical solutions, while numerical results are in very good agreement with the higher-order analytical results. We show that the extended lubrication theory is a robust tool for an accurate estimate of pressure drop in channels with shape changes on the order of the channel height, accounting for both smooth and sharp changes in geometry.
]]>2017-10-04T00:05:17-07:00info:doi/10.1098/rspa.2017.0234hwp:master-id:royprsa;rspa.2017.02342017-10-04Research articles47322062017023420170234<![CDATA[Modelling wave-induced sea ice break-up in the marginal ice zone]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170258?rss=1
A model of ice floe break-up under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave break-up events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a break-up criterion. A closed-feedback loop algorithm is devised, which (i) solves the wave-scattering problem for a given FSD under time-harmonic plane wave forcing, (ii) computes the stress field in all the floes, (iii) fractures the floes satisfying the break-up criterion, and (iv) generates an updated FSD, initializing the geometry for the next iteration of the loop. The FSD after 50 break-up events is unimodal and near normal, or bimodal, suggesting waves alone do not govern the power law observed in some field studies. Multiple scattering is found to enhance break-up for long waves and thin ice, but to reduce break-up for short waves and thick ice. A break-up front marches forward in the latter regime, as wave-induced fracture weakens the ice cover, allowing waves to travel deeper into the MIZ.
]]>2017-10-04T00:05:17-07:00info:doi/10.1098/rspa.2017.0258hwp:master-id:royprsa;rspa.2017.02582017-10-04Research articles47322062017025820170258<![CDATA[Comparing the constructions of Goldberg, Fuller, Caspar, Klug and Coxeter, and a general approach to local symmetry-preserving operations]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170267?rss=1
The use of operations on polyhedra possibly dates back to the ancient Greeks, who were the first to describe the Archimedean solids, which can be constructed from the Platonic solids by local symmetry-preserving operations (e.g. truncation) on the solid. By contrast, the results of decorations of polyhedra, e.g. by Islamic artists and by Escher, have been interpreted as decorated polyhedra—and not as new and different polyhedra. Only by interpreting decorations as combinatorial operations does it become clear how closely these two approaches are connected. In this article, we first sketch and compare the operations of Goldberg, Fuller, Caspar & Klug and Coxeter to construct polyhedra with icosahedral symmetry, where all faces are pentagons or hexagons and all vertices have three neighbours. We point out and correct an error in Goldberg’s construction. In addition, we transform the term symmetry-preserving into an exact requirement. This goal, symmetry-preserving, could also be obtained by taking global properties into account, e.g. the symmetry group itself, so we make precise the terms local and operation. As a result, we can generalize Goldberg’s approach to a systematic one that uses chamber operations to encompass all local symmetry-preserving operations on polyhedra.
]]>2017-10-25T00:05:19-07:00info:doi/10.1098/rspa.2017.0267hwp:master-id:royprsa;rspa.2017.02672017-10-25Research articles47322062017026720170267<![CDATA[On three-dimensional misorientation spaces]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170274?rss=1
Determining the local orientation of crystals in engineering and geological materials has become routine with the advent of modern crystallographic mapping techniques. These techniques enable many thousands of orientation measurements to be made, directing attention towards how such orientation data are best studied. Here, we provide a guide to the visualization of misorientation data in three-dimensional vector spaces, reduced by crystal symmetry, to reveal crystallographic orientation relationships. Domains for all point group symmetries are presented and an analysis methodology is developed and applied to identify crystallographic relationships, indicated by clusters in the misorientation space, in examples from materials science and geology. This analysis aids the determination of active deformation mechanisms and evaluation of cluster centres and spread enables more accurate description of transformation processes supporting arguments regarding provenance.
]]>2017-10-25T00:05:19-07:00info:doi/10.1098/rspa.2017.0274hwp:master-id:royprsa;rspa.2017.02742017-10-25Research articles47322062017027420170274<![CDATA[A Fokker-Planck reaction model for the epitaxial growth and shape transition of quantum dots]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170308?rss=1
We construct a Fokker–Planck reaction (FPR) model to investigate the dynamics of the coupled epitaxial growth and shape transition process of an array of quantum dots (QDs). The FPR model is based on a coupled system of Fokker–Planck equations wherein the distribution of each island type is governed by its own Fokker–Planck equation for growth, with reaction terms describing the shape transitions between islands of different types including asymmetric shapes. The reaction terms for the shape transitions depend on the island size and are determined from explicit calculations of the lowest barrier pathway for each shape transition. This mean-field model enables us to consider the kinetics of asymmetric shape transitions and study the evolution of island shape distributions during the coupled growth and transition process. Asymmetric metastable shapes play a crucial role in the dynamics, with asymmetric QDs comprising up to 10% of the population, and with up to 100% of the shape transitions passing through asymmetric shapes. Moreover, we find that the characteristic multimodal distribution of pyramid/dome QD coarsening can be eliminated at sufficiently high temperature and deposition rate.
]]>2017-10-18T00:05:20-07:00info:doi/10.1098/rspa.2017.0308hwp:master-id:royprsa;rspa.2017.03082017-10-18Research articles47322062017030820170308<![CDATA[Indentation of a floating elastic sheet: geometry versus applied tension]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170335?rss=1
The localized loading of an elastic sheet floating on a liquid bath occurs at scales from a frog sitting on a lily pad to a volcano supported by the Earth’s tectonic plates. The load is supported by a combination of the stresses within the sheet (which may include applied tensions from, for example, surface tension) and the hydrostatic pressure in the liquid. At the same time, the sheet deforms, and may wrinkle, because of the load. We study this problem in terms of the (relatively weak) applied tension and the indentation depth. For small indentation depths, we find that the force–indentation curve is linear with a stiffness that we characterize in terms of the applied tension and bending stiffness of the sheet. At larger indentations, the force–indentation curve becomes nonlinear and the sheet is subject to a wrinkling instability. We study this wrinkling instability close to the buckling threshold and calculate both the number of wrinkles at onset and the indentation depth at onset, comparing our theoretical results with experiments. Finally, we contrast our results with those previously reported for very thin, highly bendable membranes.
]]>2017-10-11T00:28:31-07:00info:doi/10.1098/rspa.2017.0335hwp:master-id:royprsa;rspa.2017.03352017-10-11Research articles47322062017033520170335<![CDATA[Nonlinear mechanics of non-rigid origami: an efficient computational approach]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170348?rss=1
Origami-inspired designs possess attractive applications to science and engineering (e.g. deployable, self-assembling, adaptable systems). The special geometric arrangement of panels and creases gives rise to unique mechanical properties of origami, such as reconfigurability, making origami designs well suited for tunable structures. Although often being ignored, origami structures exhibit additional soft modes beyond rigid folding due to the flexibility of thin sheets that further influence their behaviour. Actual behaviour of origami structures usually involves significant geometric nonlinearity, which amplifies the influence of additional soft modes. To investigate the nonlinear mechanics of origami structures with deformable panels, we present a structural engineering approach for simulating the nonlinear response of non-rigid origami structures. In this paper, we propose a fully nonlinear, displacement-based implicit formulation for performing static/quasi-static analyses of non-rigid origami structures based on ‘bar-and-hinge’ models. The formulation itself leads to an efficient and robust numerical implementation. Agreement between real models and numerical simulations demonstrates the ability of the proposed approach to capture key features of origami behaviour.
]]>2017-10-11T00:28:31-07:00info:doi/10.1098/rspa.2017.0348hwp:master-id:royprsa;rspa.2017.03482017-10-11Research articles47322062017034820170348<![CDATA[Kinematic dynamos in spheroidal geometries]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170432?rss=1
The kinematic dynamo problem is solved numerically for a spheroidal conducting fluid of possibly large aspect ratio with an insulating exterior. The solution method uses solenoidal representations of the magnetic field and the velocity by spheroidal toroidal and poloidal fields in a non-orthogonal coordinate system. Scaling of coordinates and fields to a spherical geometry leads to a modified form of the kinematic dynamo problem with a geometric anisotropic diffusion and an anisotropic current-free condition in the exterior, which is solved explicitly. The scaling allows the use of well-developed spherical harmonic techniques in angle. Dynamo solutions are found for three axisymmetric flows in oblate spheroids with semi-axis ratios 1≤a/c≤25. For larger aspect ratios strong magnetic fields may occur in any region of the spheroid, depending on the flow, but the external fields for all three flows are weak and concentrated near the axis or periphery of the spheroid.
]]>2017-10-04T00:05:17-07:00info:doi/10.1098/rspa.2017.0432hwp:master-id:royprsa;rspa.2017.04322017-10-04Research articles47322062017043220170432<![CDATA[Nonlinear graphene plasmonics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170433?rss=1
The rapid development of graphene has opened up exciting new fields in graphene plasmonics and nonlinear optics. Graphene's unique two-dimensional band structure provides extraordinary linear and nonlinear optical properties, which have led to extreme optical confinement in graphene plasmonics and ultrahigh nonlinear optical coefficients, respectively. The synergy between graphene's linear and nonlinear optical properties gave rise to nonlinear graphene plasmonics, which greatly augments graphene-based nonlinear device performance beyond a billion-fold. This nascent field of research will eventually find far-reaching revolutionary technological applications that require device miniaturization, low power consumption and a broad range of operating wavelengths approaching the far-infrared, such as optical computing, medical instrumentation and security applications.
]]>2017-10-25T00:05:19-07:00info:doi/10.1098/rspa.2017.0433hwp:master-id:royprsa;rspa.2017.04332017-10-25Review articles47322062017043320170433<![CDATA[The motion of a vortex on a closed surface of constant negative curvature]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170447?rss=1
The purpose of this work is to present an algorithm to determine the motion of a single hydrodynamic vortex on a closed surface of constant curvature and of genus greater than one. The algorithm is based on a relation between the Laplace–Beltrami Green function and the heat kernel. The algorithm is used to compute the motion of a vortex on the Bolza surface. This is the first determination of the orbits of a vortex on a closed surface of genus greater than one. The numerical results show that all the 46 vortex equilibria can be explicitly computed using the symmetries of the Bolza surface. Some of these equilibria allow for the construction of the first two examples of infinite vortex crystals on the hyperbolic disc. The following theorem is proved: ‘a Weierstrass point of a hyperellitic surface of constant curvature is always a vortex equilibrium’.
]]>2017-10-25T00:05:19-07:00info:doi/10.1098/rspa.2017.0447hwp:master-id:royprsa;rspa.2017.04472017-10-25Research articles47322062017044720170447<![CDATA[Twistor theory at fifty: from contour integrals to twistor strings]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170530?rss=1
We review aspects of twistor theory, its aims and achievements spanning the last five decades. In the twistor approach, space–time is secondary with events being derived objects that correspond to compact holomorphic curves in a complex threefold—the twistor space. After giving an elementary construction of this space, we demonstrate how solutions to linear and nonlinear equations of mathematical physics—anti-self-duality equations on Yang–Mills or conformal curvature—can be encoded into twistor cohomology. These twistor correspondences yield explicit examples of Yang–Mills and gravitational instantons, which we review. They also underlie the twistor approach to integrability: the solitonic systems arise as symmetry reductions of anti-self-dual (ASD) Yang–Mills equations, and Einstein–Weyl dispersionless systems are reductions of ASD conformal equations. We then review the holomorphic string theories in twistor and ambitwistor spaces, and explain how these theories give rise to remarkable new formulae for the computation of quantum scattering amplitudes. Finally, we discuss the Newtonian limit of twistor theory and its possible role in Penrose’s proposal for a role of gravity in quantum collapse of a wave function.
]]>2017-10-11T00:28:31-07:00info:doi/10.1098/rspa.2017.0530hwp:master-id:royprsa;rspa.2017.05302017-10-11Review articles47322062017053020170530<![CDATA[A bioinspired study on the compressive resistance of helicoidal fibre structures]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2206/20170538?rss=1
Helicoidal fibre structures are widely observed in natural materials. In this paper, an integrated experimental and analytical approach was used to investigate the compressive resistance of helicoidal fibre structures. First, helicoidal fibre-reinforced composites were created using three-dimensionally printed helicoids and polymeric matrices, including plain, ring-reinforced and helix-reinforced helicoids. Then, load–displacement curves under monotonic compression tests were collected to measure the compressive strengths of helicoidal fibre composites. Fractographic characterization was performed using an X-ray microtomographer and scanning electron microscope, through which crack propagations in helicoidal structures were illustrated. Finally, mathematical modelling was performed to reveal the essential fibre architectures in the compressive resistance of helicoidal fibre structures. This work reveals that fibre–matrix ratios, helix pitch angles and interlayer rotary angles are critical to the compressive resistance of helicoidal structures.
]]>2017-10-12T02:18:05-07:00info:doi/10.1098/rspa.2017.0538hwp:master-id:royprsa;rspa.2017.05382017-10-12Research articles47322062017053820170538<![CDATA[Density of convex intersections and applications]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20160919?rss=1
In this paper, we address density properties of intersections of convex sets in several function spaces. Using the concept of -convergence, it is shown in a general framework, how these density issues naturally arise from the regularization, discretization or dualization of constrained optimization problems and from perturbed variational inequalities. A variety of density results (and counterexamples) for pointwise constraints in Sobolev spaces are presented and the corresponding regularity requirements on the upper bound are identified. The results are further discussed in the context of finite-element discretizations of sets associated with convex constraints. Finally, two applications are provided, which include elasto-plasticity and image restoration problems.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2016.0919hwp:master-id:royprsa;rspa.2016.09192017-09-20Research articles47322052016091920160919<![CDATA[On the Keller-Rubinow model for Liesegang ring formation]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170128?rss=1
We study the model of Keller & Rubinow (Keller & Rubinow 1981 J. Chem. Phys74, 5000–5007. (doi:10.1063/1.441752)) describing the formation of Liesegang rings due to Ostwald's supersaturation mechanism. Keller and Rubinow provided an approximate solution both for the growth and equilibration of the first band, and also for the formation of secondary bands, based on a presumed asymptotic limit. However, they did not provide a parametric basis for the assumptions in their solution, nor did they provide any numerical corroboration, particularly of the secondary band formation. Here, we provide a different asymptotic solution, based on a specific parametric limit, and we show that the growth and subsequent cessation of the first band can be explained. We also show that the model is unable to explain the formation of finite width secondary bands, and we confirm this result by numerical computation. We conclude that the model is not fully posed, lacking a transition variable which can describe the hysteretic switch across the nucleation threshold.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2017.0128hwp:master-id:royprsa;rspa.2017.01282017-09-20Research articles47322052017012820170128<![CDATA[Bitwise efficiency in chaotic models]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170144?rss=1
Motivated by the increasing energy consumption of supercomputing for weather and climate simulations, we introduce a framework for investigating the bit-level information efficiency of chaotic models. In comparison with previous explorations of inexactness in climate modelling, the proposed and tested information metric has three specific advantages: (i) it requires only a single high-precision time series; (ii) information does not grow indefinitely for decreasing time step; and (iii) information is more sensitive to the dynamics and uncertainties of the model rather than to the implementation details. We demonstrate the notion of bit-level information efficiency in two of Edward Lorenz’s prototypical chaotic models: Lorenz 1963 (L63) and Lorenz 1996 (L96). Although L63 is typically integrated in 64-bit ‘double’ floating point precision, we show that only 16 bits have significant information content, given an initial condition uncertainty of approximately 1% of the size of the attractor. This result is sensitive to the size of the uncertainty but not to the time step of the model. We then apply the metric to the L96 model and find that a 16-bit scaled integer model would suffice given the uncertainty of the unresolved sub-grid-scale dynamics. We then show that, by dedicating computational resources to spatial resolution rather than numeric precision in a field programmable gate array (FPGA), we see up to 28.6% improvement in forecast accuracy, an approximately fivefold reduction in the number of logical computing elements required and an approximately 10-fold reduction in energy consumed by the FPGA, for the L96 model.
]]>2017-09-06T00:09:07-07:00info:doi/10.1098/rspa.2017.0144hwp:master-id:royprsa;rspa.2017.01442017-09-06Research articles47322052017014420170144<![CDATA[The problem of camouflaging via mirror reflections]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170147?rss=1
This work is related to billiards and their applications in geometric optics. It is known that perfectly invisible bodies with mirror surface do not exist. It is, therefore, natural to search for bodies that are, in a sense, close to invisible. We introduce a visibility index of a body measuring the mean angle of deviation of incident light rays, and derive a lower estimate for this index. This estimate is a function of the body’s volume and of the minimal radius of a ball containing the body. This result is far from being final and opens a possibility for further research.
]]>2017-09-26T07:50:32-07:00info:doi/10.1098/rspa.2017.0147hwp:master-id:royprsa;rspa.2017.01472017-09-26Research articles47322052017014720170147<![CDATA[Fully localized post-buckling states of cylindrical shells under axial compression]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170177?rss=1
We compute nonlinear force equilibrium solutions for a clamped thin cylindrical shell under axial compression. The equilibrium solutions are dynamically unstable and located on the stability boundary of the unbuckled state. A fully localized single dimple deformation is identified as the edge state—the attractor for the dynamics restricted to the stability boundary. Under variation of the axial load, the single dimple undergoes homoclinic snaking in the azimuthal direction, creating states with multiple dimples arranged around the central circumference. Once the circumference is completely filled with a ring of dimples, snaking in the axial direction leads to further growth of the dimple pattern. These fully nonlinear solutions embedded in the stability boundary of the unbuckled state constitute critical shape deformations. The solutions may thus be a step towards explaining when the buckling and subsequent collapse of an axially loaded cylinder shell is triggered.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0177hwp:master-id:royprsa;rspa.2017.01772017-09-13Research articles47322052017017720170177<![CDATA[Meso-scale defect evaluation of selective laser melting using spatially resolved acoustic spectroscopy]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170194?rss=1
Developments in additive manufacturing technology are serving to expand the potential applications. Critical developments are required in the supporting areas of measurement and in process inspection to achieve this. CM247LC is a nickel superalloy that is of interest for use in aerospace and civil power plants. However, it is difficult to process via selective laser melting (SLM) as it suffers from cracking during rapid cooling and solidification. This limits the viability of CM247LC parts created using SLM. To quantify part integrity, spatially resolved acoustic spectroscopy (SRAS) has been identified as a viable non-destructive evaluation technique. In this study, a combination of optical microscopy and SRAS was used to identify and classify the surface defects present in SLM-produced parts. By analysing the datasets and scan trajectories, it is possible to correlate morphological information with process parameters. Image processing was used to quantify porosity and cracking for bulk density measurement. Analysis of surface acoustic wave data showed that an error in manufacture in the form of an overscan occurred. Comparing areas affected by overscan with a bulk material, a change in defect density from 1.17% in the bulk material to 5.32% in the overscan regions was observed, highlighting the need to reduce overscan areas in manufacture.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0194hwp:master-id:royprsa;rspa.2017.01942017-09-13Research articles47322052017019420170194<![CDATA[Elastic field of a spherical inclusion with non-uniform eigenfields in second strain gradient elasticity]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170254?rss=1
In this paper, within the complete form of Mindlin’s second strain gradient theory, the elastic field of an isolated spherical inclusion embedded in an infinitely extended homogeneous isotropic medium due to a non-uniform distribution of eigenfields is determined. These eigenfields, in addition to eigenstrain, comprise eigen double and eigen triple strains. After the derivation of a closed-form expression for Green’s function associated with the problem, two different cases of non-uniform distribution of the eigenfields are considered as follows: (i) radial distribution, i.e. the distributions of the eigenfields are functions of only the radial distance of points from the centre of inclusion, and (ii) polynomial distribution, i.e. the distributions of the eigenfields are polynomial functions in the Cartesian coordinates of points. While the obtained solution for the elastic field of the latter case takes the form of an infinite series, the solution to the former case is represented in a closed form. Moreover, Eshelby’s tensors associated with the two mentioned cases are obtained.
]]>2017-09-06T00:09:07-07:00info:doi/10.1098/rspa.2017.0254hwp:master-id:royprsa;rspa.2017.02542017-09-06Research articles47322052017025420170254<![CDATA[The steady aerodynamics of aerofoils with porosity gradients]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170266?rss=1
This theoretical study determines the aerodynamic loads on an aerofoil with a prescribed porosity distribution in a steady incompressible flow. A Darcy porosity condition on the aerofoil surface furnishes a Fredholm integral equation for the pressure distribution, which is solved exactly and generally as a Riemann–Hilbert problem provided that the porosity distribution is Hölder-continuous. The Hölder condition includes as a subset any continuously differentiable porosity distributions that may be of practical interest. This formal restriction on the analysis is examined by a class of differentiable porosity distributions that approach a piecewise, discontinuous function in a certain parametric limit. The Hölder-continuous solution is verified in this limit against analytical results for partially porous aerofoils in the literature. Finally, a comparison made between the new theoretical predictions and experimental measurements of SD7003 aerofoils presented in the literature. Results from this analysis may be integrated into a theoretical framework to optimize turbulence noise suppression with minimal impact to aerodynamic performance.
]]>2017-09-26T07:50:32-07:00info:doi/10.1098/rspa.2017.0266hwp:master-id:royprsa;rspa.2017.02662017-09-26Research articles47322052017026620170266<![CDATA[Rate-based structural health monitoring using permanently installed sensors]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170270?rss=1
Permanently installed sensors are becoming increasingly ubiquitous, facilitating very frequent in situ measurements and consequently improved monitoring of ‘trends’ in the observed system behaviour. It is proposed that this newly available data may be used to provide prior warning and forecasting of critical events, particularly system failure. Numerous damage mechanisms are examples of positive feedback; they are ‘self-accelerating’ with an increasing rate of damage towards failure. The positive feedback leads to a common time-response behaviour which may be described by an empirical relation allowing prediction of the time to criticality. This study focuses on Structural Health Monitoring of engineering components; failure times are projected well in advance of failure for fatigue, creep crack growth and volumetric creep damage experiments. The proposed methodology provides a widely applicable framework for using newly available near-continuous data from permanently installed sensors to predict time until failure in a range of application areas including engineering, geophysics and medicine.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0270hwp:master-id:royprsa;rspa.2017.02702017-09-13Research articles47322052017027020170270<![CDATA[Numerical studies of cavitation erosion on an elastic-plastic material caused by shock-induced bubble collapse]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170315?rss=1
We present a study of shock-induced collapse of single bubbles near/attached to an elastic–plastic solid using the free-Lagrange method, which forms the latest part of our shock-induced collapse studies. We simulated the collapse of 40 μm radius single bubbles near/attached to rigid and aluminium walls by a 60 MPa lithotripter shock for various scenarios based on bubble–wall separations, and the collapse of a 255 μm radius bubble attached to aluminium foil with a 65 MPa lithotripter shock. The coupling of the multi-phases, compressibility, axisymmetric geometry and elastic–plastic material model within a single solver has enabled us to examine the impingement of high-speed liquid jets from the shock-induced collapsing bubbles, which imposes an extreme compression in the aluminium that leads to pitting and plastic deformation. For certain scenarios, instead of the high-speed jet, a radially inwards flow along the aluminium surface contracts the bubble to produce a ‘mushroom shape’. This work provides methods for quantifying which parameters (e.g. bubble sizes and separations from the solid) might promote or inhibit erosion on solid surfaces.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0315hwp:master-id:royprsa;rspa.2017.03152017-09-13Research articles47322052017031520170315<![CDATA[Singular diffusionless limits of double-diffusive instabilities in magnetohydrodynamics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170344?rss=1
We study local instabilities of a differentially rotating viscous flow of electrically conducting incompressible fluid subject to an external azimuthal magnetic field. In the presence of the magnetic field, the hydrodynamically stable flow can demonstrate non-axisymmetric azimuthal magnetorotational instability (AMRI) both in the diffusionless case and in the double-diffusive case with viscous and ohmic dissipation. Performing stability analysis of amplitude transport equations of short-wavelength approximation, we find that the threshold of the diffusionless AMRI via the Hamilton–Hopf bifurcation is a singular limit of the thresholds of the viscous and resistive AMRI corresponding to the dissipative Hopf bifurcation and manifests itself as the Whitney umbrella singular point. A smooth transition between the two types of instabilities is possible only if the magnetic Prandtl number is equal to unity, Pm=1. At a fixed Pm=1, the threshold of the double-diffusive AMRI is displaced by finite distance in the parameter space with respect to the diffusionless case even in the zero dissipation limit. The complete neutral stability surface contains three Whitney umbrella singular points and two mutually orthogonal intervals of self-intersection. At these singularities, the double-diffusive system reduces to a marginally stable system which is either Hamiltonian or parity–time-symmetric.
]]>2017-09-13T00:57:51-07:00info:doi/10.1098/rspa.2017.0344hwp:master-id:royprsa;rspa.2017.03442017-09-13Research articles47322052017034420170344<![CDATA[Irreversible particle motion in surfactant-laden interfaces due to pressure-dependent surface viscosity]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170346?rss=1
The surface shear viscosity of an insoluble surfactant monolayer often depends strongly on its surface pressure. Here, we show that a particle moving within a bounded monolayer breaks the kinematic reversibility of low-Reynolds-number flows. The Lorentz reciprocal theorem allows such irreversibilities to be computed without solving the full nonlinear equations, giving the leading-order contribution of surface pressure-dependent surface viscosity. In particular, we show that a disc translating or rotating near an interfacial boundary experiences a force in the direction perpendicular to that boundary. In unbounded monolayers, coupled modes of motion can also lead to non-intuitive trajectories, which we illustrate using an interfacial analogue of the Magnus effect. This perturbative approach can be extended to more complex geometries, and to two-dimensional suspensions more generally.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2017.0346hwp:master-id:royprsa;rspa.2017.03462017-09-20Research articles47322052017034620170346<![CDATA[Stability of the shear-thinning boundary-layer flow over a flat inclined plate]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170350?rss=1
In this study, we consider the boundary-layer flow of an inelastic non-Newtonian fluid over an inclined flat plate. Using two popular generalized Newtonian models, we determine base flow profiles and associated linear stability results for a range shear-thinning fluids. In addition to neutral stability curves, we also present results concerning the linear growth of the Tollmien–Schlichting waves as they propagate downstream. Furthermore, to gain an insight into the underlying physical mechanisms affecting the destabilization of the disturbances, an integral energy equation is derived and energy calculations are presented. Results from all three analyses suggest that the effect of shear-thinning will act to stabilize the boundary-layer flow. Consequently, it can be argued that the addition of shear-thinning agents could act as a passive control mechanism for flows of this nature.
]]>2017-09-06T00:09:07-07:00info:doi/10.1098/rspa.2017.0350hwp:master-id:royprsa;rspa.2017.03502017-09-06Research articles47322052017035020170350<![CDATA[Uncertainties in the 2004 Sumatra-Andaman source through nonlinear stochastic inversion of tsunami waves]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170353?rss=1
Numerical inversions for earthquake source parameters from tsunami wave data usually incorporate subjective elements to stabilize the search. In addition, noisy and possibly insufficient data result in instability and non-uniqueness in most deterministic inversions, which are barely acknowledged. Here, we employ the satellite altimetry data for the 2004 Sumatra–Andaman tsunami event to invert the source parameters. We also include kinematic parameters that improve the description of tsunami generation and propagation, especially near the source. Using a finite fault model that represents the extent of rupture and the geometry of the trench, we perform a new type of nonlinear joint inversion of the slips, rupture velocities and rise times with minimal a priori constraints. Despite persistently good waveform fits, large uncertainties in the joint parameter distribution constitute a remarkable feature of the inversion. These uncertainties suggest that objective inversion strategies should incorporate more sophisticated physical models of seabed deformation in order to significantly improve the performance of early warning systems.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2017.0353hwp:master-id:royprsa;rspa.2017.03532017-09-20Research articles47322052017035320170353<![CDATA[Accuracy of the hypothetical sky-polarimetric Viking navigation versus sky conditions: revealing solar elevations and cloudinesses favourable for this navigation method]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170358?rss=1
According to Thorkild Ramskou's theory proposed in 1967, under overcast and foggy skies, Viking seafarers might have used skylight polarization analysed with special crystals called sunstones to determine the position of the invisible Sun. After finding the occluded Sun with sunstones, its elevation angle had to be measured and its shadow had to be projected onto the horizontal surface of a sun compass. According to Ramskou's theory, these sunstones might have been birefringent calcite or dichroic cordierite or tourmaline crystals working as polarizers. It has frequently been claimed that this method might have been suitable for navigation even in cloudy weather. This hypothesis has been accepted and frequently cited for decades without any experimental support. In this work, we determined the accuracy of this hypothetical sky-polarimetric Viking navigation for 1080 different sky situations characterized by solar elevation and cloudiness , the sky polarization patterns of which were measured by full-sky imaging polarimetry. We used the earlier measured uncertainty functions of the navigation steps 1, 2 and 3 for calcite, cordierite and tourmaline sunstone crystals, respectively, and the newly measured uncertainty function of step 4 presented here. As a result, we revealed the meteorological conditions under which Vikings could have used this hypothetical navigation method. We determined the solar elevations at which the navigation uncertainties are minimal at summer solstice and spring equinox for all three sunstone types. On average, calcite sunstone ensures a more accurate sky-polarimetric navigation than tourmaline and cordierite. However, in some special cases (generally at 35° ≤ ≤ 40°, 1 okta ≤ ≤ 6 oktas for summer solstice, and at 20° ≤ ≤ 25°, 0 okta ≤ ≤ 4 oktas for spring equinox), the use of tourmaline and cordierite results in smaller navigation uncertainties than that of calcite. Generally, under clear or less cloudy skies, the sky-polarimetric navigation is more accurate, but at low solar elevations its accuracy remains relatively large even at high cloudiness. For a given , the absolute value of averaged peak North uncertainties dramatically decreases with increasing until the sign (±) change of these uncertainties. For a given , this absolute value can either decrease or increase with increasing . The most advantageous sky situations for this navigation method are at summer solstice when the solar elevation and cloudiness are 35° ≤ ≤ 40° and 2 oktas ≤ ≤ 3 oktas.
]]>2017-09-13T00:57:51-07:00info:doi/10.1098/rspa.2017.0358hwp:master-id:royprsa;rspa.2017.03582017-09-13Research articles47322052017035820170358<![CDATA[On the topology of the Lorenz system]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170374?rss=1
We present a new paradigm for three-dimensional chaos, and specifically for the Lorenz equations. The main difficulty in these equations and for a generic flow in dimension 3 is the existence of singularities. We show how to use knot theory as a way to remove the singularities. Specifically, we claim: (i) for certain parameters, the Lorenz system has an invariant one-dimensional curve, which is a trefoil knot. The knot is a union of invariant manifolds of the singular points. (ii) The flow is topologically equivalent to an Anosov flow on the complement of this curve, and moreover to a geodesic flow. (iii) When varying the parameters, the system exhibits topological phase transitions, i.e. for special parameter values, it will be topologically equivalent to an Anosov flow on a knot complement. Different knots appear for different parameter values and each knot controls the dynamics at nearby parameters.
]]>2017-09-20T00:09:30-07:00info:doi/10.1098/rspa.2017.0374hwp:master-id:royprsa;rspa.2017.03742017-09-20Research articles47322052017037420170374<![CDATA[Toroidal drop under electric field: arbitrary drop-to-ambient fluid viscosity ratio]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170379?rss=1
In the absence of external forces, a liquid toroidal drop freely suspended in another fluid shrinks towards its centre. It is shown that if the two phases are slightly conducting viscous incompressible fluids with the drop-to-ambient fluid ratios of electric conductivities, dielectric constants and viscosities to be 1/R, Q and , respectively, then the toroidal drop with volume 4/3 and having major radius can become almost stationary when subjected to a uniform electric field aligned with the drop’s axis of symmetry. In this case, Q and electric capillary number Ca_{E} that defines the ratio of electric stress to surface tension, are functions of R, and and are found analytically. Those functions are relatively insensitive to , and for ≥3, they admit simple approximations, which coincide with those obtained recently for =1. Streamlines inside and outside the toroidal drop for the same R and but different are qualitatively similar.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2017.0379hwp:master-id:royprsa;rspa.2017.03792017-09-20Research articles47322052017037920170379<![CDATA[A priori estimation of memory effects in reduced-order models of nonlinear systems using the Mori-Zwanzig formalism]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170385?rss=1
Reduced models of nonlinear dynamical systems require closure, or the modelling of the unresolved modes. The Mori–Zwanzig procedure can be used to derive formally closed evolution equations for the resolved physics. In these equations, the unclosed terms are recast as a memory integral involving the time history of the resolved variables. While this procedure does not reduce the complexity of the original system, these equations can serve as a mathematically consistent basis to develop closures based on memory approximations. In this scenario, knowledge of the memory kernel is paramount in assessing the validity of a memory approximation. Unravelling the memory kernel requires solving the orthogonal dynamics, which is a high-dimensional partial differential equation that is intractable, in general. A method to estimate the memory kernel a priori, using full-order solution snapshots, is proposed. The key idea is to solve a pseudo orthogonal dynamics equation, which has a convenient Liouville form, instead. This ersatz arises from the assumption that the semi-group of the orthogonal dynamics is a composition operator for one observable. The method is exact for linear systems. Numerical results on the Burgers and Kuramoto–Sivashinsky equations demonstrate that the proposed technique can provide valuable information about the memory kernel.
]]>2017-09-26T07:50:32-07:00info:doi/10.1098/rspa.2017.0385hwp:master-id:royprsa;rspa.2017.03852017-09-26Research articles47322052017038520170385<![CDATA[On the differing growth mechanisms of black-smoker and Lost City-type hydrothermal vents]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170387?rss=1
Black smokers and Lost City-type springs are varieties of hydrothermal vents on the ocean floors that emit hot, acidic water and cool, alkaline water, respectively. While both produce precipitation structures as the issuing fluid encounters oceanic water, Lost City-type hydrothermal vents in particular have been implicated in the origin of life on the Earth. We present a parallel-velocity flow model for the radius and flow rate of a cylindrical jet of fluid that forms the template for the growth of a tube precipitated about itself and we compare the solution with previous laboratory experimental results from growth of silicate chemical gardens. We show that when the growth of the solid structure is determined by thermal diffusion, fluid flow is slow at the solid–liquid contact. However, in the case of chemical diffusive transport, the fluid jet effectively drags the liquid in the pores of the solid precipitate. These findings suggest a continuum in the diffusive growth rate of hydrothermal vent structures, where Lost City-type hydrothermal vents favour contact between the vent fluid and the external seawater. We explore the implications for the road to life.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0387hwp:master-id:royprsa;rspa.2017.03872017-09-13Research articles47322052017038720170387<![CDATA[Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170388?rss=1
In Holm (Holm 2015 Proc. R. Soc. A471, 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow.
]]>2017-09-20T01:05:12-07:00info:doi/10.1098/rspa.2017.0388hwp:master-id:royprsa;rspa.2017.03882017-09-20Research articles47322052017038820170388<![CDATA[Can a quantum state over time resemble a quantum state at a single time?]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170395?rss=1
The standard formalism of quantum theory treats space and time in fundamentally different ways. In particular, a composite system at a given time is represented by a joint state, but the formalism does not prescribe a joint state for a composite of systems at different times. If there were a way of defining such a joint state, this would potentially permit a more even-handed treatment of space and time, and would strengthen the existing analogy between quantum states and classical probability distributions. Under the assumption that the joint state over time is an operator on the tensor product of single-time Hilbert spaces, we analyse various proposals for such a joint state, including one due to Leifer and Spekkens, one due to Fitzsimons, Jones and Vedral, and another based on discrete Wigner functions. Finding various problems with each, we identify five criteria for a quantum joint state over time to satisfy if it is to play a role similar to the standard joint state for a composite system: that it is a Hermitian operator on the tensor product of the single-time Hilbert spaces; that it represents probabilistic mixing appropriately; that it has the appropriate classical limit; that it has the appropriate single-time marginals; that composing over multiple time steps is associative. We show that no construction satisfies all these requirements. If Hermiticity is dropped, then there is an essentially unique construction that satisfies the remaining four criteria.
]]>2017-09-20T00:09:30-07:00info:doi/10.1098/rspa.2017.0395hwp:master-id:royprsa;rspa.2017.03952017-09-20Research articles47322052017039520170395<![CDATA[Coordinate representation for non-Hermitian position and momentum operators]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170434?rss=1
In this paper, we undertake an analysis of the eigenstates of two non-self-adjoint operators q^ and p^ similar, in a suitable sense, to the self-adjoint position and momentum operators q^0 and p^0 usually adopted in ordinary quantum mechanics. In particular, we discuss conditions for these eigenstates to be biorthogonal distributions, and we discuss a few of their properties. We illustrate our results with two examples, one in which the similarity map between the self-adjoint and the non-self-adjoint is bounded, with bounded inverse, and the other in which this is not true. We also briefly propose an alternative strategy to deal with q^ and p^, based on the so-called quasi *-algebras.
]]>2017-09-13T00:09:15-07:00info:doi/10.1098/rspa.2017.0434hwp:master-id:royprsa;rspa.2017.04342017-09-13Research articles47322052017043420170434<![CDATA[Cautionary tales on air-quality improvement in Beijing]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170457?rss=1
The official air-quality statistic reported that Beijing had a 9.9% decline in the annual concentration of PM_{2.5} in 2016. While this statistic offered some relief for the inhabitants of the capital, we present several analyses based on Beijing's PM_{2.5} data of the past 4 years at 36 monitoring sites along with meteorological data of the past 7 years. The analyses reveal the air pollution situation in 2016 was not as rosy as the 9.9% decline would convey, and improvement if any was rather uncertain. The paper also provides an assessment on the city's PM_{2.5} situation in the past 4 years.
]]>2017-09-20T00:09:30-07:00info:doi/10.1098/rspa.2017.0457hwp:master-id:royprsa;rspa.2017.04572017-09-20Research articles47322052017045720170457<![CDATA[Elastodynamic image forces on screw dislocations in the presence of phase boundaries]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170484?rss=1
The elastodynamic image forces acting on straight screw dislocations in the presence of planar phase boundaries are derived. Two separate dislocations are studied: (i) the injected, non-moving screw dislocation and (ii) the injected (or pre-existing), generally non-uniformly moving screw dislocation. The image forces are derived for both the case of a rigid surface and of a planar interface between two homogeneous, isotropic phases. The case of a rigid interface is shown to be solvable employing Head's image dislocation construction. The case of the elastodynamic image force due to an interface is solved by deriving the reflected wave's contribution to the global solution across the interface. This entails obtaining the fundamental solution (Green's function) for a point unit force via Cagniard's method, and then applying the convolution theorem for a screw dislocation modelled as a force distribution. Complete, explicit formulae are provided when available. It is shown that the elastodynamic image forces are generally affected by retardation effects, and that those acting on the moving dislocations display a dynamic magnification that exceed the attraction (or repulsion) predicted in classical elastostatic calculations.
]]>2017-09-20T00:09:30-07:00info:doi/10.1098/rspa.2017.0484hwp:master-id:royprsa;rspa.2017.04842017-09-20Research articles47322052017048420170484<![CDATA[Correction to 'On a new class of electroelastic bodies. I]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2205/20170523xx?rss=1
2017-09-28T23:25:40-07:00info:doi/10.1098/rspa.2017.0523hwp:master-id:royprsa;rspa.2017.05232017-09-28Corrections473220520170523xx20170523xx<![CDATA[Whitham modulation theory for the Kadomtsev- Petviashvili equation]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20160695?rss=1
The genus-1 Kadomtsev–Petviashvili (KP)-Whitham system is derived for both variants of the KP equation; namely the KPI and KPII equations. The basic properties of the KP-Whitham system, including symmetries, exact reductions and its possible complete integrability, together with the appropriate generalization of the one-dimensional Riemann problem for the Korteweg–de Vries equation are discussed. Finally, the KP-Whitham system is used to study the linear stability properties of the genus-1 solutions of the KPI and KPII equations; it is shown that all genus-1 solutions of KPI are linearly unstable, while all genus-1 solutions of KPII are linearly stable within the context of Whitham theory.
]]>2017-08-02T00:08:23-07:00info:doi/10.1098/rspa.2016.0695hwp:master-id:royprsa;rspa.2016.06952017-08-02Research articles47322042016069520160695<![CDATA[A variation on the Donsker-Varadhan inequality for the principal eigenvalue]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20160877?rss=1
The purpose of this short paper is to give a variation on the classical Donsker–Varadhan inequality, which bounds the first eigenvalue of a second-order elliptic operator on a bounded domain by the largest mean first exit time of the associated drift–diffusion process via 1≥1supxExc.Instead of looking at the mean of the first exit time, we study quantiles: let dp,:->R≥0 be the smallest time t such that the likelihood of exiting within that time is p, then 1≥log(1/p)supxdp,(x).Moreover, as p->0, this lower bound converges to _{1}.
]]>2017-08-23T00:08:55-07:00info:doi/10.1098/rspa.2016.0877hwp:master-id:royprsa;rspa.2016.08772017-08-23Research articles47322042016087720160877<![CDATA[Recursive modular modelling methodology for lumped-parameter dynamic systems]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20160891?rss=1
This paper proposes a novel approach to the modelling of lumped-parameter dynamic systems, based on representing them by hierarchies of mathematical models of increasing complexity instead of a single (complex) model. Exploring the multilevel modularity that these systems typically exhibit, a general recursive modelling methodology is proposed, in order to conciliate the use of the already existing modelling techniques. The general algorithm is based on a fundamental theorem that states the conditions for computing projection operators recursively. Three procedures for these computations are discussed: orthonormalization, use of orthogonal complements and use of generalized inverses. The novel methodology is also applied for the development of a recursive algorithm based on the Udwadia–Kalaba equation, which proves to be identical to the one of a Kalman filter for estimating the state of a static process, given a sequence of noiseless measurements representing the constraints that must be satisfied by the system.
]]>2017-08-30T00:09:15-07:00info:doi/10.1098/rspa.2016.0891hwp:master-id:royprsa;rspa.2016.08912017-08-30Research articles47322042016089120160891<![CDATA[Rebound mechanics of micrometre-scale, spherical particles in high-velocity impacts]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20160936?rss=1
The impact mechanics of micrometre-scale metal particles with flat metal surfaces is investigated for high-velocity impacts ranging from 50 m s^{–1} to more than 1 km s^{–1}, where impact causes predominantly plastic deformation. A material model that includes high strain rate and temperature effects on the yield stress, heat generation due to plasticity, material damage due to excessive plastic strain and heat transfer is used in the numerical analysis. The coefficient of restitution e is predicted by the classical work using elastic–plastic deformation analysis with quasi-static impact mechanics to be proportional to Vi–1/4 and Vi–1/2 for the low and moderate impact velocities that span the ranges of 0–10 and 10–100 m s^{–1}, respectively. In the elastic–plastic and fully plastic deformation regimes the particle rebound is attributed to the elastic spring-back that initiates at the particle–substrate interface. At higher impact velocities (0.1–1 km s^{–1}) e is shown to be proportional to approximately Vi–1. In this deeply plastic deformation regime various deformation modes that depend on plastic flow of the material including the time lag between the rebound instances of the top and bottom points of particle and the lateral spreading of the particle are identified. In this deformation regime, the elastic spring-back initiates subsurface, in the substrate.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2016.0936hwp:master-id:royprsa;rspa.2016.09362017-08-16Research articles47322042016093620160936<![CDATA[Model selection for dynamical systems via sparse regression and information criteria]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170009?rss=1
We develop an algorithm for model selection which allows for the consideration of a combinatorially large number of candidate models governing a dynamical system. The innovation circumvents a disadvantage of standard model selection which typically limits the number of candidate models considered due to the intractability of computing information criteria. Using a recently developed sparse identification of nonlinear dynamics algorithm, the sub-selection of candidate models near the Pareto frontier allows feasible computation of Akaike information criteria (AIC) or Bayes information criteria scores for the remaining candidate models. The information criteria hierarchically ranks the most informative models, enabling the automatic and principled selection of the model with the strongest support in relation to the time-series data. Specifically, we show that AIC scores place each candidate model in the strong support, weak support or no support category. The method correctly recovers several canonical dynamical systems, including a susceptible-exposed-infectious-recovered disease model, Burgers’ equation and the Lorenz equations, identifying the correct dynamical system as the only candidate model with strong support.
]]>2017-08-30T00:09:15-07:00info:doi/10.1098/rspa.2017.0009hwp:master-id:royprsa;rspa.2017.00092017-08-30Research articles47322042017000920170009<![CDATA[Organic molecule fluorescence as an experimental test-bed for quantum jumps in thermodynamics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170099?rss=1
We demonstrate with an experiment how molecules are a natural test bed for probing fundamental quantum thermodynamics. Single-molecule spectroscopy has undergone transformative change in the past decade with the advent of techniques permitting individual molecules to be distinguished and probed. We demonstrate that the quantum Jarzynski equality for heat is satisfied in this set-up by considering the time-resolved emission spectrum of organic molecules as arising from quantum jumps between states. This relates the heat dissipated into the environment to the free energy difference between the initial and final state. We demonstrate also how utilizing the quantum Jarzynski equality allows for the detection of energy shifts within a molecule, beyond the relative shift.
]]>2017-08-30T00:09:15-07:00info:doi/10.1098/rspa.2017.0099hwp:master-id:royprsa;rspa.2017.00992017-08-30Research articles47322042017009920170099<![CDATA[Initial boundary-value problem for the spherically symmetric Einstein equations with fluids with tangential pressure]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170113?rss=1
We prove that, for a given spherically symmetric fluid distribution with tangential pressure on an initial space-like hypersurface with a time-like boundary, there exists a unique, local in time solution to the Einstein equations in a neighbourhood of the boundary. As an application, we consider a particular elastic fluid interior matched to a vacuum exterior.
]]>2017-08-09T00:08:43-07:00info:doi/10.1098/rspa.2017.0113hwp:master-id:royprsa;rspa.2017.01132017-08-09Research articles47322042017011320170113<![CDATA[Designing the optimal bit: balancing energetic cost, speed and reliability]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170117?rss=1
We consider the challenge of operating a reliable bit that can be rapidly erased. We find that both erasing and reliability times are non-monotonic in the underlying friction, leading to a trade-off between erasing speed and bit reliability. Fast erasure is possible at the expense of low reliability at moderate friction, and high reliability comes at the expense of slow erasure in the underdamped and overdamped limits. Within a given class of bit parameters and control strategies, we define ‘optimal’ designs of bits that meet the desired reliability and erasing time requirements with the lowest operational work cost. We find that optimal designs always saturate the bound on the erasing time requirement, but can exceed the required reliability time if critically damped. The non-trivial geometry of the reliability and erasing time scales allows us to exclude large regions of parameter space as suboptimal. We find that optimal designs are either critically damped or close to critical damping under the erasing procedure.
]]>2017-08-23T00:08:55-07:00info:doi/10.1098/rspa.2017.0117hwp:master-id:royprsa;rspa.2017.01172017-08-23Research articles47322042017011720170117<![CDATA[Characterization of columnar inertial modes in rapidly rotating spheres and spheroids]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170181?rss=1
We consider fluid-filled spheres and spheroidal containers of eccentricity in rapid rotation, as a proxy for the interior dynamics of stars and planets. The fluid motion is assumed to be quasi-geostrophic (QG): horizontal motions are invariant parallel to the rotation axis z, a characteristic which is handled by use of a stream function formulation which additionally enforces mass conservation and non-penetration at the boundary. By linearizing about a quiescent background state, we investigate a variety of methods to study the QG inviscid inertial wave modes which are compared with fully three-dimensional (3D) calculations. We consider the recently proposed weak formulation of the inviscid system valid in spheroids of arbitrary eccentricity, to which we present novel closed-form polynomial solutions. Our modal solutions accurately represent, in both spatial structure and frequency, the most z-invariant of the inertial wave modes in a spheroid, and constitute a simple basis set for the analysis of rotationally dominated fluids. We further show that these new solutions are more accurate than those of the classical axial-vorticity equation, which is independent of and thus fails to properly encode the container geometry. We also consider the effects of viscosity for the cases of both no-slip and stress-free boundary conditions for a spherical container. Calculations performed under the columnar approximation are compared with 3D solutions and excellent agreement has been found despite fundamental differences in the two formulations.
]]>2017-08-09T00:08:43-07:00info:doi/10.1098/rspa.2017.0181hwp:master-id:royprsa;rspa.2017.01812017-08-09Research articles47322042017018120170181<![CDATA[Effective potentials in nonlinear polycrystals and quadrature formulae]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170213?rss=1
This study presents a family of estimates for effective potentials in nonlinear polycrystals. Noting that these potentials are given as averages, several quadrature formulae are investigated to express these integrals of nonlinear functions of local fields in terms of the moments of these fields. Two of these quadrature formulae reduce to known schemes, including a recent proposition (Ponte Castañeda 2015 Proc. R. Soc. A471, 20150665 (doi:10.1098/rspa.2015.0665)) obtained by completely different means. Other formulae are also reviewed that make use of statistical information on the fields beyond their first and second moments. These quadrature formulae are applied to the estimation of effective potentials in polycrystals governed by two potentials, by means of a reduced-order model proposed by the authors (non-uniform transformation field analysis). It is shown how the quadrature formulae improve on the tangent second-order approximation in porous crystals at high stress triaxiality. It is found that, in order to retrieve a satisfactory accuracy for highly nonlinear porous crystals under high stress triaxiality, a quadrature formula of higher order is required.
]]>2017-08-30T00:09:15-07:00info:doi/10.1098/rspa.2017.0213hwp:master-id:royprsa;rspa.2017.02132017-08-30Research articles47322042017021320170213<![CDATA[Numerical simulations of drumlin formation]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170220?rss=1
We summarize the present form of the instability theory for drumlin formation, which describes the coupled subglacial flow of ice, water and sediment. This model has evolved over the last 20 years, and is now at the point where it can predict instabilities corresponding to ribbed moraine, drumlins and mega-scale glacial lineations, but efforts to provide numerical solutions of the model have been limited. The present summary adds some slight nuances to previously published versions of the theory, notably concerning the constitutive description of the subglacial water film and its flow. A new numerical method is devised to solve the model, and we show that it can be solved for realistic values of most of the parameters, with the exception of that corresponding to the water film thickness. We show that evolved bedforms can be three-dimensional and of the correct sizes, and we explore to some extent the variation of the solutions with the model’s parameters.
]]>2017-08-30T00:09:15-07:00info:doi/10.1098/rspa.2017.0220hwp:master-id:royprsa;rspa.2017.02202017-08-30Research articles47322042017022020170220<![CDATA[Temperature profile in a liquid-vapour interface near the critical point]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170229?rss=1
Thanks to an expansion with respect to densities of energy, mass and entropy, we discuss the concept of thermocapillary fluid for inhomogeneous fluids. The non-convex state law valid for homogeneous fluids is modified by adding terms taking account of the gradients of these densities. This seems more realistic than Cahn and Hilliard’s model which uses a density expansion in mass-density gradient only. Indeed, through liquid–vapour interfaces, realistic potentials in molecular theories show that entropy density and temperature do not vary with the mass density as it would do in bulk phases. In this paper, we prove using a rescaling process near the critical point, that liquid–vapour interfaces behave essentially in the same way as in Cahn and Hilliard’s model.
]]>2017-08-09T00:08:43-07:00info:doi/10.1098/rspa.2017.0229hwp:master-id:royprsa;rspa.2017.02292017-08-09Research articles47322042017022920170229<![CDATA[Flexural edge waves generated by steady-state propagation of a loaded rectilinear crack in an elastically supported thin plate]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170265?rss=1
The problem of a rectilinear crack propagating at constant speed in an elastically supported thin plate and acted upon by an equally moving load is considered. The full-field solution is obtained and the spotlight is set on flexural edge wave generation. Below the critical speed for the appearance of travelling waves, a threshold speed is met which marks the transformation of decaying edge waves into edge waves propagating along the crack and dying away from it. Yet, besides these, and for any propagation speed, a pair of localized edge waves, which rapidly decay behind the crack tip, is also shown to exist. These waves are characterized by a novel dispersion relation and fade off from the crack line in an oscillatory manner, whence they play an important role in the far field behaviour. Dynamic stress intensity factors are obtained and, for speed close to the critical speed, they show a resonant behaviour which expresses the most efficient way to channel external work into the crack. Indeed, this behaviour is justified through energy considerations regarding the work of the applied load and the energy release rate. Results might be useful in a wide array of applications, ranging from fracturing and machining to acoustic emission and defect detection.
]]>2017-08-30T00:09:14-07:00info:doi/10.1098/rspa.2017.0265hwp:master-id:royprsa;rspa.2017.02652017-08-30Research articles47322042017026520170265<![CDATA[The structure of turbulence in a rapid tidal flow]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170295?rss=1
The structure of turbulence in a rapid tidal flow is characterized through new observations of fundamental statistical properties at a site in the UK which has a simple geometry and sedate surface wave action. The mean flow at the Sound of Islay exceeded 2.5 m s^{–1} and the turbulent boundary layer occupied the majority of the water column, with an approximately logarithmic mean velocity profile identifiable close to the seabed. The anisotropic ratios, spectral scales and higher-order statistics of the turbulence generally agree well with values reported for two-dimensional open channels in the laboratory and other tidal channels, therefore providing further support for the application of universal models. The results of the study can assist in developing numerical models of turbulence in rapid tidal flows such as those proposed for tidal energy generation.
]]>2017-08-23T00:08:55-07:00info:doi/10.1098/rspa.2017.0295hwp:master-id:royprsa;rspa.2017.02952017-08-23Research articles47322042017029520170295<![CDATA[Nonlinear electroelasticity: material properties, continuum theory and applications]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170311?rss=1
In the last few years, it has been recognized that the large deformation capacity of elastomeric materials that are sensitive to electric fields can be harnessed for use in transducer devices such as actuators and sensors. This has led to the reassessment of the mathematical theory that is needed for the description of the electromechanical (in particular, electroelastic) interactions for purposes of material characterization and prediction. After a review of the key experiments concerned with determining the nature of the electromechanical interactions and a discussion of the range of applications to devices, we provide a short account of the history of developments in the nonlinear theory. This is followed by a succinct modern treatment of electroelastic theory, including the governing equations and constitutive laws needed for both material characterization and the analysis of general electroelastic coupling problems. For illustration, the theory is then applied to two simple representative boundary-value problems that are relevant to the geometries of activation devices; in particular, (a) a rectangular plate and (b) a circular cylindrical tube, in each case with compliant electrodes on the major surfaces and a potential difference between them. In (a), an electric field is generated normal to the major surfaces and in (b), a radial electric field is present. This is followed by a short section in which other problems addressed on the basis of the general theory are described briefly.
]]>2017-08-02T00:08:23-07:00info:doi/10.1098/rspa.2017.0311hwp:master-id:royprsa;rspa.2017.03112017-08-02Review articles47322042017031120170311<![CDATA[Mittag-Leffler synchronization of delayed fractional-order bidirectional associative memory neural networks with discontinuous activations: state feedback control and impulsive control schemes]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170322?rss=1
This paper is concerned with the drive–response synchronization for a class of fractional-order bidirectional associative memory neural networks with time delays, as well as in the presence of discontinuous activation functions. The global existence of solution under the framework of Filippov for such networks is firstly obtained based on the fixed-point theorem for condensing map. Then the state feedback and impulsive controllers are, respectively, designed to ensure the Mittag-Leffler synchronization of these neural networks and two new synchronization criteria are obtained, which are expressed in terms of a fractional comparison principle and Razumikhin techniques. Numerical simulations are presented to validate the proposed methodologies.
]]>2017-08-02T00:08:23-07:00info:doi/10.1098/rspa.2017.0322hwp:master-id:royprsa;rspa.2017.03222017-08-02Research articles47322042017032220170322<![CDATA[Elastic deformation of twinned microstructures]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170330?rss=1
Many crystalline materials exhibit twinned microstructures, where well-defined orientation relationships define the special symmetry between different, elastically anisotropic twin variants. When such twins are subjected to external loading, additional internal stresses necessarily occur at the twin boundaries in order to maintain compatibility. These compatibility stresses are constant inside each variant in repeating stacks of twins and considerably affect the local stress state. In this paper, we use anisotropic linear elasticity to derive general analytical solutions for compatibility stresses in a stack of twin variants in arbitrary materials, for arbitrary variant volume fractions and twin types, subjected to arbitrary applied stresses. By considering two examples, growth twins in electrodeposited Cu and B19' martensite twins in the shape memory alloy NiTi, we further demonstrate that compatibility stresses can considerably alter the preferred slip systems for dislocation plasticity as well as the effective macroscopic behaviour of twinned microstructures.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2017.0330hwp:master-id:royprsa;rspa.2017.03302017-08-16Research articles47322042017033020170330<![CDATA[Adaptive compliant structures for flow regulation]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170334?rss=1
This paper introduces conceptual design principles for a novel class of adaptive structures that provide both flow regulation and control. While of general applicability, these design principles, which revolve around the idea of using the instabilities and elastically nonlinear behaviour of post-buckled panels, are exemplified through a case study: the design of a shape-adaptive air inlet. The inlet comprises a deformable post-buckled member that changes shape depending on the pressure field applied by the surrounding fluid, thereby regulating the inlet aperture. By tailoring the stress field in the post-buckled state and the geometry of the initial, stress-free configuration, the deformable section can snap through to close or open the inlet completely. Owing to its inherent ability to change shape in response to external stimuli—i.e. the aerodynamic loads imposed by different operating conditions—the inlet does not have to rely on linkages and mechanisms for actuation, unlike conventional flow-controlling devices.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2017.0334hwp:master-id:royprsa;rspa.2017.03342017-08-16Research articles47322042017033420170334<![CDATA[A tale of two nested elastic rings]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170340?rss=1
Elastic rods in contact provide a rich paradigm for understanding shape and deformation in interacting elastic bodies. Here, we consider the problem of determining the static solutions of two nested elastic rings in the plane. If the inner ring is longer than the outer ring, it will buckle creating a space between the two rings. This deformation can be further influenced by either adhesion between the rings or if pressure is applied externally or internally. We obtain an exact solution of this problem when both rings are assumed inextensible and unshearable. Through a variational formulation of the problem, we identify the boundary conditions at the contact point and use the Kirchhoff analogy to give exact solutions of the problems in terms of elliptic functions. The role of both adhesion and pressure is explored.
]]>2017-08-09T00:08:43-07:00info:doi/10.1098/rspa.2017.0340hwp:master-id:royprsa;rspa.2017.03402017-08-09Research articles47322042017034020170340<![CDATA[Spectral derivation of the classic laws of wall-bounded turbulent flows]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170354?rss=1
We show that the classic laws of the mean-velocity profiles (MVPs) of wall-bounded turbulent flows—the ‘law of the wall,’ the ‘defect law’ and the ‘log law’—can be predicated on a sufficient condition with no manifest ties to the MVPs, namely that viscosity and finite turbulent domains have a depressive effect on the spectrum of turbulent energy. We also show that this sufficient condition is consistent with empirical data on the spectrum and may be deemed a general property of the energetics of wall turbulence. Our findings shed new light on the physical origin of the classic laws and their immediate offshoot, Prandtl’s theory of turbulent friction.
]]>2017-08-09T00:08:43-07:00info:doi/10.1098/rspa.2017.0354hwp:master-id:royprsa;rspa.2017.03542017-08-09Research articles47322042017035420170354<![CDATA[Multi-parameter actuation of a neutrally stable shell: a flexible gear-less motor]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170364?rss=1
We have designed and tested experimentally a morphing structure consisting of a neutrally stable thin cylindrical shell driven by a multi-parameter piezoelectric actuation. The shell is obtained by plastically deforming an initially flat copper disc, so as to induce large isotropic and almost uniform inelastic curvatures. Following the plastic deformation, in a perfectly isotropic system, the shell is theoretically neutrally stable, having a continuous set of stable cylindrical shapes corresponding to the rotation of the axis of maximal curvature. Small imperfections render the actual structure bistable, giving preferred orientations. A three-parameter piezoelectric actuation, exerted through micro-fibre-composite actuators, allows us to add a small perturbation to the plastic inelastic curvature and to control the direction of maximal curvature. This actuation law is designed through a geometrical analogy based on a fully nonlinear inextensible uniform-curvature shell model. We report on the fabrication, identification and experimental testing of a prototype and demonstrate the effectiveness of the piezoelectric actuators in controlling its shape. The resulting motion is an apparent rotation of the shell, controlled by the voltages as in a ‘gear-less motor’, which is, in reality, a precession of the axis of principal curvature.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2017.0364hwp:master-id:royprsa;rspa.2017.03642017-08-16Research articles47322042017036420170364<![CDATA[Origin of the onset of meandering of a straight river]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170376?rss=1
In this paper, to explore the origin of the onset of meandering of a straight river, we, first, analyse the linear stability of a straight river. We discover that the natural perturbation modes of a straight river maintain an equilibrium state by confining themselves to an onset wavenumber band that is dependent on the flow regimes, aspect ratio, relative roughness number and Shields number. Then, we put forward a phenomenological description of the onset of meandering of a straight river. Its mechanism is governed by turbulent flow, with counter-rotation of neighbouring large-scale or macro-turbulent eddies in succession to generate the processes of alternating erosion and deposition of sediment grains of the riverbed. This concept is explained by a theorem (universal scaling law) stemming from the phenomenology of a turbulent energy cascade to provide a quantitative insight into the criterion for the onset of meandering of a straight river. It is revealed from this universal scaling law that, at the onset of meandering of a river, the longitudinal riverbed slope is a unique function of the river width, flow discharge and sediment grain size. This unique functional relationship is corroborated by the data obtained from the measurements in natural and model rivers.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2017.0376hwp:master-id:royprsa;rspa.2017.03762017-08-16Research articles47322042017037620170376<![CDATA[Some applications of mathematics in golf]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170392?rss=1
At its core, like many other sports, golf is a game of integers. The minimization of the number of strokes played is generally what determines the winner, whether each of these are associated with the shortest of putts or the longest of drives. The outcomes of these shots are influenced by very slight changes, but hopefully in a deterministic sense. Understanding the mechanics of golf necessitates the development of models and this is coupled more often than not to the use of statistics. In essence, the individual aspects of the sport can be modelled adequately via fairly simplistic models, but the presence of a human at one end of the kinematic chain has a significant impact on the variability of the entire process. In this paper, we will review some of the ways that mathematics has been used to develop the understanding of the physical processes involved in the sport, including some of the analysis which is exploited within the Equipment Rules. We will also discuss some of the future challenges.
]]>2017-08-02T00:08:23-07:00info:doi/10.1098/rspa.2017.0392hwp:master-id:royprsa;rspa.2017.03922017-08-02Special feature47322042017039220170392<![CDATA[Unsteady solute dispersion in small blood vessels using a two-phase Casson model]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/473/2204/20170427?rss=1
This study explores the transport of a solute in an unsteady blood flow in small arteries with and without absorption at the wall. The Casson fluid model is suitable for blood flow in small vessels. Owing to the aggregation of red cells in the central region of the small vessels, a two-phase model is considered in this investigation. Using the generalized dispersion model (Sankarasubramanian & Gill 1973 Proc. R. Soc. Lond. A333, 115–132. (doi:10.1098/rspa.1973.0051)), the convection, dispersion and mean concentration of the solute are analysed at all times in small arteries of different radii. The effects of the yield stress, wall absorption, the amplitude of the fluctuating pressure gradient component, the peripheral layer thickness, the Womersley frequency parameter, the Schmidt number and the Peclet number on the dispersion process are discussed. A comparative study of solute dispersion among single- and two-phase fluid models is presented. For small vessels, a significant difference between these models is observed during the solute dispersion; however, this difference becomes insignificant for large vessels. The mean concentration of solute reduces with increasing radius of the vessels. The present investigation is more realistic for understanding the transportation process of drugs in blood flow in small arteries using the non-Newtonian fluid model.
]]>2017-08-16T00:08:59-07:00info:doi/10.1098/rspa.2017.0427hwp:master-id:royprsa;rspa.2017.04272017-08-16Research articles47322042017042720170427