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[Interfacial load transfer mechanisms in carbon nanotube-polymer nanocomposites]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20170705?rss=1
Carbon nanotubes (CNTs) are highly promising for strength reinforcement in polymer nanocomposites, but conflicting interfacial properties have been reported by single nanotube pull-out experiments. Here, we report the interfacial load transfer mechanisms during pull-out of CNTs from PMMA matrices, using massively- parallel molecular dynamics simulations. We show that the pull-out forces associated with non-bonded interactions between CNT and PMMA are generally small, and are weakly-dependent on the embedment length of the nanotube. These pull-out forces do not significantly increase with the presence of Stone Wales or vacancy defects along the nanotube. In contrast, low-density distribution of cross-links along the CNT-PMMA interface increases the pull-out forces by an order of magnitude. At each cross-linked site, mechanical unfolding and pull-out of single or pair polymer chain(s) attached to the individual cross-link bonds result in substantial interfacial strengthening and toughening, while contributing to interfacial slip between CNT and PMMA. Our interfacial shear-slip model shows that the interfacial loads are evenly-distributed among the finite number of cross-link bonds at low cross-link densities or for nanotubes with short embedment lengths. At higher cross-link densities or for nanotubes with longer embedment lengths, a no-slip zone now develops where shear-lag effects become important. Implications of these results, in the context of recent nanotube pull-out experiments, are discussed.
]]>2018-08-08T00:05:15-07:00info:doi/10.1098/rspa.2017.0705hwp:master-id:royprsa;rspa.2017.07052018-08-08Research articles47422162017070520170705<![CDATA[Radon signals in geological (natural) geogas and in a simultaneous enhanced confined mode simulation experiment]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20170787?rss=1
An enhanced confined mode (ECM) radon simulation experiment, tested in the laboratory in Jerusalem, was relocated to a subsurface geophysical observatory located 400 km apart, at a depth of 150 m and with a stable temperature. Five gamma sensors are placed around the ECM canister and lead shielding minimizes the influence of natural local gamma radiation. Simultaneous measurement of the geological radon and from radon in the ECM system indicates that the temporal variation of gamma radiation from radon in the ECM system contains annual, multi-day and daily signals, that correspond to signals in the local geological radon. This implies that a common external driver influences the radiation pattern of the geological radon and from radon inside the ECM canister. Once activated at BGO the typical variation pattern of the experimental system occurring in the laboratory changed to that occurring at the observatory. This is interpreted to indicate that the overall style of the temporal patterns of radiation from radon is site dependent. The outcome of this investigation conforms and further substantiates the recent suggestion that a component in solar radiation is driving the annual and daily periodic components in the variation of radon. New geophysical research potential is indicated.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2017.0787hwp:master-id:royprsa;rspa.2017.07872018-08-01Research articles47422162017078720170787<![CDATA[Analytical estimation of non-local deformation-mediated magneto-electric coupling in soft composites]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20170803?rss=1
For a long time, the search for magneto-electric materials concentrated on multi-ferroics and hard-matter composites. By contrast, rather recently the exploitation of strain-mediated magneto-electric (ME) coupling in soft composites was proposed. The basic idea behind this approach is to combine the magneto- and electro-mechanical responses of composites consisting of a soft matrix carrying magnetic inclusions. Despite that such composites are straightforward to manufacture and have cheap constituents, they did not gain much attention up to now. In this contribution, we demonstrate that ME coupling induced by finite deformations could be of significant magnitude. Our approach relies on shape effects as a special non-local phenomenon in magneto- and electro-elasticity. Based on that we characterize an up to now overlooked ME coupling mechanism which purely relies on these shape effects in soft-matter-based magnetic and electric media. While soft magnetic media are commonly realized as composites, the coupling effect to be highlighted exists independently of the origin of a body's magnetic and electric properties. We show that the magnitude of the effect is indeed significant and, among ellipsoidal bodies, most pronounced for those of spherical to moderately prolate shape. Finite-element simulations are performed to assess the quality of the analytical predictions.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2017.0803hwp:master-id:royprsa;rspa.2017.08032018-08-01Research articles47422162017080320170803<![CDATA[Small-amplitude static periodic patterns at a fluid-ferrofluid interface]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180038?rss=1
We establish the existence of static doubly periodic patterns (in particular rolls, squares and hexagons) on the free surface of a ferrofluid near onset of the Rosensweig instability, assuming a general (nonlinear) magnetization law. A novel formulation of the ferrohydrostatic equations in terms of Dirichlet–Neumann operators for nonlinear elliptic boundary-value problems is presented. We demonstrate the analyticity of these operators in suitable function spaces and solve the ferrohydrostatic problem using an analytic version of Crandall–Rabinowitz local bifurcation theory. Criteria are derived for the bifurcations to be sub-, super- or transcritical with respect to a dimensionless physical parameter.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2018.0038hwp:master-id:royprsa;rspa.2018.00382018-08-01Research articles47422162018003820180038<![CDATA[Wrinkling in engineering fabrics: a comparison between two different comprehensive modelling approaches]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180063?rss=1
We consider two ‘comprehensive’ modelling approaches for engineering fabrics. We distinguish the two approaches using the terms ‘semi-discrete’ and ‘continuum’, reflecting their natures. We demonstrate a fitting procedure, used to identify the constitutive parameters of the continuum model from predictions of the semi-discrete model, the parameters of which are in turn fitted to experimental data. We, then, check the effectiveness of the continuum model by verifying the correspondence between semi-discrete and continuum model predictions using test cases not previously used in the identification process. Predictions of both modelling approaches are compared against full-field experimental kinematic data, obtained using stereoscopic digital image correlation techniques, and also with measured force data. Being a reduced order model and being implemented in an implicit rather than an explicit finite-element code, the continuum model requires significantly less computational power than the semi-discrete model and could therefore be used to more efficiently explore the mechanical response of engineering fabrics.
]]>2018-08-08T00:05:15-07:00info:doi/10.1098/rspa.2018.0063hwp:master-id:royprsa;rspa.2018.00632018-08-08Research articles47422162018006320180063<![CDATA[Comment on 'The cutting of metals by plastic buckling by Udupa et al.]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180170?rss=1
In a recent paper on machining annealed copper at a low cutting speed, and at an uncut chip thickness one-tenth of the mean grain size of the copper, Udupa et al. (Proc. R. Soc A473, 20160863, doi:10.1098/rspa.2016.0863) report chip thicknesses larger than 10 times the uncut thickness and then a new mode of chip formation. Plastic bulging occurs in the surface of the copper ahead of the tool, leading to chip formation by a series of folds. The strain in the chip is less than that expected in a chip formed by shear according to long-standing classical theory. The authors suggest that the foundations of that theory need to be re-examined. In response, continuum mechanics numerical simulations presented here show a continuous transition from the classical condition towards that observed by Udupa et al. as the ratio of chip thickness to uncut thickness increases above approximately 7. Bulging is obtained by introducing (approximately) material heterogeneity to the simulations at a grain size scale but whether such heterogeneity is essential for the bulging flows remains an open question.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2018.0170hwp:master-id:royprsa;rspa.2018.01702018-08-01Comment47422162018017020180170<![CDATA[Dispersive shock waves governed by the Whitham equation and their stability]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180278?rss=1
Dispersive shock waves (DSWs), also termed undular bores in fluid mechanics, governed by the non-local Whitham equation are studied in order to investigate short wavelength effects that lead to peaked and cusped waves within the DSW. This is done by combining the weak nonlinearity of the Korteweg–de Vries equation with full linear dispersion relations. The dispersion relations considered are those for surface gravity waves, the intermediate long wave equation and a model dispersion relation introduced by Whitham to investigate the 120° peaked Stokes wave of highest amplitude. A dispersive shock fitting method is used to find the leading (solitary wave) and trailing (linear wave) edges of the DSW. This method is found to produce results in excellent agreement with numerical solutions up until the lead solitary wave of the DSW reaches its highest amplitude. Numerical solutions show that the DSWs for the water wave and Whitham peaking kernels become modulationally unstable and evolve into multi-phase wavetrains after a critical amplitude which is just below the DSW of maximum amplitude.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2018.0278hwp:master-id:royprsa;rspa.2018.02782018-08-01Research articles47422162018027820180278<![CDATA[From arteries to boreholes: transient response of a poroelastic cylinder to fluid injection]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180284?rss=1
The radially outward flow of fluid through a porous medium occurs in many practical problems, from transport across vascular walls to the pressurization of boreholes in the subsurface. When the driving pressure is non-negligible relative to the stiffness of the solid structure, the poromechanical coupling between the fluid and the solid can control both the steady state and the transient mechanics of the system. Very large pressures or very soft materials lead to large deformations of the solid skeleton, which introduce kinematic and constitutive nonlinearity that can have a non-trivial impact on these mechanics. Here, we study the transient response of a poroelastic cylinder to sudden fluid injection. We consider the impacts of kinematic and constitutive nonlinearity, both separately and in combination, and we highlight the central role of driving method in the evolution of the response. We show that the various facets of nonlinearity may either accelerate or decelerate the transient response relative to linear poroelasticity, depending on the boundary conditions and the initial geometry, and that an imposed fluid pressure leads to a much faster response than an imposed fluid flux.
]]>2018-08-01T00:05:26-07:00info:doi/10.1098/rspa.2018.0284hwp:master-id:royprsa;rspa.2018.02842018-08-01Research articles47422162018028420180284<![CDATA[On the geometrically exact low-order modelling of a flexible beam: formulation and numerical tests]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180423?rss=1
This paper proposes a low-order geometrically exact flexible beam formulation based on the utilization of generic beam shape functions to approximate distributed kinematic properties of the deformed structure. The proposed nonlinear beam shapes approach is in contrast to the majority of geometrically nonlinear treatments in the literature in which element-based—and hence high-order—discretizations are adopted. The kinematic quantities approximated specifically pertain to shear and extensional gradients as well as local orientation parameters based on an arbitrary set of globally referenced attitude parameters. In developing the dynamic equations of motion, an Euler angle parametrization is selected as it is found to yield fast computational performance. The resulting dynamic formulation is closed using an example shape function set satisfying the single generic kinematic constraint. The formulation is demonstrated via its application to the modelling of a series of static and dynamic test cases of both simple and non-prismatic structures; the simulated results are verified using MSC Nastran and an element-based intrinsic beam formulation. Through these examples, it is shown that the nonlinear beam shapes approach is able to accurately capture the beam behaviour with a very minimal number of system states.
]]>2018-08-08T00:05:15-07:00info:doi/10.1098/rspa.2018.0423hwp:master-id:royprsa;rspa.2018.04232018-08-08Research articles47422162018042320180423<![CDATA[Designing heterogeneous catalysts]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2216/20180514?rss=1
2018-08-15T00:05:20-07:00info:doi/10.1098/rspa.2018.0514hwp:master-id:royprsa;rspa.2018.05142018-08-15Special feature47422162018051420180514<![CDATA[The underwater resonant airbag: a new wave energy converter]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170192?rss=1
The time-domain simulation follows the heaving of the conical float in waves and calculates the bag shape, ballast motion, adiabatic air pressure and the flow through the turbine. There are two independent oscillators, the float with its resonance and the bag/ballast with its resonance. The coupling of the two oscillators gives rise to a wide band response with two peaks in the capture width each reaching the theoretical /2. In this new wave energy converter, apart from the turbine, there are no mechanical moving parts, no joints nor pistons, no end stops nor sliding seals, no flaps nor one-way valves. The expected life of the airtight flexible bag remains to be determined, but potential manufacturers are optimistic.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2017.0192hwp:master-id:royprsa;rspa.2017.01922018-07-18Research articles47422152017019220170192<![CDATA[Shock waves in media with negative electromagnetic susceptibilities]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170522?rss=1
The Fock–Kemmer theory of precursor shock wave speeds in classical field theory is extended to the case of media with negative electromagnetic susceptibilities.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2017.0522hwp:master-id:royprsa;rspa.2017.05222018-07-18Research articles47422152017052220170522<![CDATA[Large numbers of explanatory variables: a probabilistic assessment]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170631?rss=1
Recently, Cox and Battey (2017 Proc. Natl Acad. Sci. USA114, 8592–8595 (doi:10.1073/pnas.1703764114)) outlined a procedure for regression analysis when there are a small number of study individuals and a large number of potential explanatory variables, but relatively few of the latter have a real effect. The present paper reports more formal statistical properties. The results are intended primarily to guide the choice of key tuning parameters.
]]>2018-07-04T00:09:48-07:00info:doi/10.1098/rspa.2017.0631hwp:master-id:royprsa;rspa.2017.06312018-07-04Research articles47422152017063120170631<![CDATA[Effects of anisotropy and regime of diffusion on the measurement of lattice diffusion coefficient of hydrogen in metals]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170677?rss=1
The lattice diffusion coefficient of hydrogen in metals is commonly measured by permeation tests. Such tests assume isotropic diffusivity with a set-up which measures permeation flux only along one direction. The measured values of the lattice diffusion coefficient are strongly influenced by the trapping of hydrogen at microstructural defects. These factors lead to highly inaccurate determination of the lattice diffusion coefficient, more so in an anisotropic medium. In this work, we present a three-dimensional (3D) diffusion equation in non-dimensional form for an anisotropic medium with source and sink terms which account for detrapping and trapping of hydrogen. The concentration of hydrogen at lattice and trap sites is assumed to be in a local equilibrium. An initial boundary value problem of the permeation test is formulated and the governing partial differential equation is implemented in a 3D finite-element code. The influence of anisotropic diffusivity on the measurement of lattice diffusion coefficient is shown by numerical simulations. Asymptotic analysis of the governing equation reveals that the lattice diffusion coefficient can only be measured in certain regimes when performing permeation tests at varying temperatures. The nonlinear behaviour of Arrhenius plots of diffusion coefficient versus inverse of temperature due to trapping is analytically and numerically predicted.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2017.0677hwp:master-id:royprsa;rspa.2017.06772018-07-18Research articles47422152017067720170677<![CDATA[The n-component nonlinear Schrödinger equations: dark-bright mixed N- and high-order solitons and breathers, and dynamics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170688?rss=1
The general n-component nonlinear Schrödinger equations are systematically investigated with the aid of the Darboux transformation method and its extension. Firstly, we explore the condition of the existence for dark–bright mixed soliton solutions and derive an explicit formula of dark–bright mixed multi-soliton solutions in terms of the determinant. Secondly, we present the formula of dark–bright mixed high-order semi-rational solitons, and predict their general Nth-order wave structures. Thirdly, we investigate the wing-shaped structures of breather. Finally, we perform the numerical simulations for some representative solitons to study their dynamical behaviours.
]]>2018-07-11T00:09:41-07:00info:doi/10.1098/rspa.2017.0688hwp:master-id:royprsa;rspa.2017.06882018-07-11Research articles47422152017068820170688<![CDATA[Determination of macro-scale soil properties from pore scale structures: image-based modelling of poroelastic structures]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170745?rss=1
We show how a combination of X-ray computed tomography (X-CT) and image-based modelling can be used to calculate the effect of moisture content and compaction on the macroscopic structural properties of soil. Our method is based on the equations derived in Daly & Roose (2018 Proc. R. Soc. A474, 20170141. (doi:10.1098/rspa.2017.0141)), which we have extended so they can be directly applied to the segmented images obtained from X-CT. We assume that the soils are composed of air-filled pore space, solid mineral grains and a mixed phase composed of both clay particles and water. We considered three different initial soil treatments, composed of two different compaction levels and two different moisture contents. We found that the effective properties of the soils were unaffected by compaction over the range tested in this paper. However, changing the moisture content significantly altered the hydraulic and mechanical properties of the soils. A key strength of this method is that it enables the optimization or even design of soils composed from different constituents, with specific mechanical and hydraulic properties.
]]>2018-07-11T00:09:41-07:00info:doi/10.1098/rspa.2017.0745hwp:master-id:royprsa;rspa.2017.07452018-07-11Research articles47422152017074520170745<![CDATA[We-T classification of diesel fuel droplet impact regimes]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170759?rss=1
A combined experimental and computational investigation of micrometric diesel droplets impacting on a heated aluminium substrate is presented. Dual view high-speed imaging has been employed to visualize the evolution of the impact process at various conditions. The parameters investigated include wall-surface temperature ranging from 140 to 400°C, impact Weber and Reynolds numbers of 19–490 and 141–827, respectively, and ambient pressure of 1 and 2 bar. Six possible post-impact regimes were identified, termed as Stick, Splash, Partial-Rebound, Rebound, Breakup-Rebound and Breakup-Stick, and plotted on the We-T map. Additionally, the temporal variation of the apparent dynamic contact angle and spreading factor have been determined as a function of the impact Weber number and surface temperature. Numerical simulations have also been performed using a two-phase flow model with interface capturing, phase-change and variable physical properties. Increased surface temperature resulted to increased maximum spreading diameter and induced quicker and stronger recoiling behaviour, mostly attributed to the change of liquid viscosity.
]]>2018-07-04T00:09:48-07:00info:doi/10.1098/rspa.2017.0759hwp:master-id:royprsa;rspa.2017.07592018-07-04Research articles47422152017075920170759<![CDATA[Effect of surface bending and stress on the transmission of line force to an elastic substrate]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170775?rss=1
For a broad class of soft materials their surface stress can strongly influence mechanical behaviour. For example, a line force applied to the surface of an elastic substrate is locally supported by surface stress over an elasto-capillary length lc (surface stress/elastic modulus). Surface stress regularizes the otherwise highly singular stress and strain fields. However, surface such as lipid bilayer interfaces can also resist deformation by bending. This has not been studied either by experiments or theories. We analyse a theoretical model of the response of a half-space to a line force when the surface carries both a stress and resistance to bending. We find that surface bending further regularizes the singular fields. The local stress field near the line load can be separated into three regions. Region 1 occupies distances from the line load smaller than an elasto-capillary bending length lb (bending stiffness/elastic modulus to the 1/3 power) where surface bending dominates and the elastic stress and strains are continuous. Region 2 occupies intermediate distances between lb and lc ( > lb) where surface stress dominates. At distances larger than lc we retrieve the classical elasticity solution. The size of region 2 depends on =lc/lb and vanishes for small lc.
]]>2018-07-11T00:09:41-07:00info:doi/10.1098/rspa.2017.0775hwp:master-id:royprsa;rspa.2017.07752018-07-11Research articles47422152017077520170775<![CDATA[A generalized model for dynamic deformation and failure of clamped beam under impact loading]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170799?rss=1
A mathematical model representing the dynamic behaviour (both plastic deformation and fracture) of a clamped beam under impact loading is developed. Indentation at the impact point, transverse deformation, formation and propagation of plastic hinge, arrest of plastic hinge leading to plastic work concentration and finally failure are the physical processes which constitute the basis of the derived model. The effect of imperfection is also incorporated. Imperfection is considered in the form of a -shaped notch located at the impact point, at supports or at both. The distinct feature of the present formulation is that it accommodates different possible deformation and failure modes in a single model. Final plastic deformation, time histories of different field variables and failure modes as predicted from the derived model are found to be in good agreement with the corresponding experimental and numerical results. This model provides a quick understanding of the dynamic behaviour of beam under impact and also the effect of various underlying parameters which may be useful for forming design provisions for impact-resistant structures.
]]>2018-07-04T00:18:17-07:00info:doi/10.1098/rspa.2017.0799hwp:master-id:royprsa;rspa.2017.07992018-07-04Research articles47422152017079920170799<![CDATA[An elliptical liquid inclusion in an infinite elastic plane]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170813?rss=1
Beyond recent related literature, which focused on spherical incompressible liquid inclusions, the present work studies an elliptical compressible liquid inclusion in an infinite elastic plane under static remote mechanical loading. Here, it is assumed that the change of pressure inside the liquid inclusion is linearly related to the change of inclusion volume with the bulk modulus of the liquid as the proportionality coefficient. Also, the role of the liquid surface tension on the solid–liquid interface is examined especially when the size of the liquid inclusion is comparable to or smaller than the elastocapillary length. Our results show that both the surface tension and the change of liquid pressure have a significant effect on reducing the stress concentration factor at the endpoints of an elliptical liquid inclusion. In addition, the pressure change inside the liquid inclusion is studied when a uniaxial remote stress is applied perpendicular or parallel to the major axis of the elliptical liquid inclusion. In particular, the effective plane-strain Young's modulus of a solid–liquid composite containing circular liquid inclusions predicted by the present model is linearly related to the volume fraction of the liquid inclusions, in reasonable agreement with existing experimental data.
]]>2018-07-04T00:09:48-07:00info:doi/10.1098/rspa.2017.0813hwp:master-id:royprsa;rspa.2017.08132018-07-04Research articles47422152017081320170813<![CDATA[Elastic waves in periodically heterogeneous two-dimensional media: locally periodic and anti-periodic modes]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20170908?rss=1
Propagation of anti-plane waves through a discrete square lattice and through a continuous fibrous medium is studied. In the long-wave limit, for periodically heterogeneous structures the solution can be periodic or anti-periodic across the unit cell. It is shown that combining periodicity and anti-periodicity conditions in different directions of the translational symmetry allows one to detect different types of modes that do not arise in the purely periodic case. Such modes may be interpreted as counterparts of non-classical waves appearing in phenomenological theories. Dispersion diagrams of the discrete square lattice are evaluated in a closed analytical from. Dispersion properties of the fibrous medium are determined using Floquet–Bloch theory and Fourier series approximations. Influence of a viscous damping is taken into account.
]]>2018-07-25T00:05:18-07:00info:doi/10.1098/rspa.2017.0908hwp:master-id:royprsa;rspa.2017.09082018-07-25Research articles47422152017090820170908<![CDATA[Negative and positive stiffness in auxetic magneto-mechanical metamaterials]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180003?rss=1
This work discusses the concept of allowing the control of the stiffness of a particular class of re-entrant auxetic magneto-mechanical metamaterials through the introduction of magnets to the system. It is shown, through experimental testing backed up by a theoretical model, that the appropriate insertion of magnets in such a system will alter its stiffness, possibly even making it exhibit ‘negative stiffness’. This leads to a completely different behaviour of the structure in terms of stability. It is also reported that the investigated mechanical metamaterials may exhibit both negative stiffness and negative Poisson's ratio at the same time. Moreover, it is shown that the effect which magnets have on the stiffness of the system may be fine-tuned upon replacing magnets with electromagnets. Such systems have the potential to be used in a wide range of practical applications such as vibration damping devices where achieving a negative stiffness is of fundamental importance.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2018.0003hwp:master-id:royprsa;rspa.2018.00032018-07-18Research articles47422152018000320180003<![CDATA[Extraordinary interactions between light and matter determined by anomalous weak values]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180030?rss=1
Some predictions regarding pre- and post-selected states are far-reaching, thereby requiring validation with standard quantum measurements in addition to the customary weak measurements used so far, as well as other advanced techniques. We go further pursuing this goal, proposing two thought experiments which incorporate novel yet feasible validation methods of unconventional light-matter interactions. An excited atom traverses a Mach–Zehnder interferometer (MZI) under a special combination of pre- and post-selection. In the first experiment, photons emitted by the superposed atom, after being hit by two laser beams, are individually counted. Despite the interaction having definitely taken place, as revealed by the atom becoming ground, the numbers of photons emitted from each arm of the MZI are predicted, at the ensemble level, to be different from those expected with standard stimulated emission. In the second experiment, the atom spontaneously emits a photon while still in the MZI. This photon later serves as a strong measurement of the atom's energy upon hitting a photographic plate. The experiment is repeated to enable an interference effect of the emitted photons. Interestingly, the latter gives the appearance that the photons have been emitted by the atom from a position much farther from the two MZI arms L and R, as if in a ‘phantom arm’ R'. Nevertheless, their time of arrival is similar to that of photons coming from L and R. These experiments also emphasize the key role of anomalous weak values in determining light–matter interactions. In fact, they present a straightforward realization of an entity we term counter-particles, namely pre- and post-selected states acting as if they have negative physical variables such as mass and energy. The novel verification methods we suggest for testing these predictions resemble weak measurements in some aspects, yet result from definite atomic transitions verified by the detected photons.
]]>2018-07-04T00:09:48-07:00info:doi/10.1098/rspa.2018.0030hwp:master-id:royprsa;rspa.2018.00302018-07-04Research articles47422152018003020180030<![CDATA[Shapes of river networks]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180081?rss=1
River network scaling laws describe how their shape varies with their size. However, the regional variation of this size-dependence remains poorly understood. Here we show that river network scaling laws vary systematically with the climatic aridity index. We find that arid basins do not change their proportions with size, while humid basins do. To explore why, we study an aspect ratio L_{}/L_{||} between basin width L_{} and basin length L_{||}. We find that the aspect ratio exhibits a dependence on climate and argue that this can be understood as a structural consequence of the confluence angle. We then find that, in humid basins, the aspect ratio decreases with basin size, which we attribute to a common hydrogeological hierarchy. Our results offer an explanation of the variability in network scaling exponents and suggest that the absence of self-similarity in humid basins can be understood as a morphological expression of subsurface processes.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2018.0081hwp:master-id:royprsa;rspa.2018.00812018-07-18Research articles47422152018008120180081<![CDATA[Bayesian inference of spreading processes on networks]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180129?rss=1
Infectious diseases are studied to understand their spreading mechanisms, to evaluate control strategies and to predict the risk and course of future outbreaks. Because people only interact with few other individuals, and the structure of these interactions influence spreading processes, the pairwise relationships between individuals can be usefully represented by a network. Although the underlying transmission processes are different, the network approach can be used to study the spread of pathogens in a contact network or the spread of rumours in a social network. We study simulated simple and complex epidemics on synthetic networks and on two empirical networks, a social/contact network in an Indian village and an online social network. Our goal is to learn simultaneously the spreading process parameters and the first infected node, given a fixed network structure and the observed state of nodes at several time points. Our inference scheme is based on approximate Bayesian computation, a likelihood-free inference technique. Our method is agnostic about the network topology and the spreading process. It generally performs well and, somewhat counter-intuitively, the inference problem appears to be easier on more heterogeneous network topologies, which enhances its future applicability to real-world settings where few networks have homogeneous topologies.
]]>2018-07-18T00:09:44-07:00info:doi/10.1098/rspa.2018.0129hwp:master-id:royprsa;rspa.2018.01292018-07-18Research articles47422152018012920180129<![CDATA[Interfacial waveforms in chiral lattices with gyroscopic spinners]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180132?rss=1
We demonstrate a new method of achieving topologically protected states in an elastic hexagonal system of trusses by attaching gyroscopic spinners, which bring chirality to the system. Dispersive features of this medium are investigated in detail, and it is shown that one can manipulate the locations of stop-bands and Dirac points by tuning the parameters of the spinners. We show that, in the proximity of such points, uni-directional interfacial waveforms can be created in an inhomogeneous lattice and the direction of such waveforms can be controlled. The effect of inserting additional soft internal links into the system, which is thus transformed into a heterogeneous triangular lattice, is also investigated, as the hexagonal lattice represents the limit case of the heterogeneous triangular lattice with soft links. This work introduces a new perspective in the design of periodic media possessing non-trivial topological features.
]]>2018-07-25T00:05:18-07:00info:doi/10.1098/rspa.2018.0132hwp:master-id:royprsa;rspa.2018.01322018-07-25Research articles47422152018013220180132<![CDATA[Three-dimensional visualization and quantification of the fracture mechanisms in sparse fibre networks using multiscale X-ray microtomography]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180175?rss=1
The structural changes that are induced by the initiation and the propagation of a crack in a low-density paper (LDP) were studied using single edge-notched fracture tests that were imaged under an optical microscope or in laboratory or synchrotron X-ray microtomographs. The two-dimensional optical images were used to analyse the links between the mesoscale structural variations of LDP and the crack path. Medium-resolution X-ray three-dimensional images were used to analyse the variations in the thickness and local porosity of samples as well as their displacement field that were induced by the LDP fracture. High-resolution three-dimensional images showed that these mesostructural variations were accompanied by complex fibre and bond deformation mechanisms that were, for the first time, in situ imaged. These mechanisms occurred in the fracture process zone that developed ahead of the crack tip before the crack path became distinct and visible. They were at the origin of the aforementioned thickness variations that developed more particularly along the crack path. They eventually led to fibre–fibre bond detachment phenomena and crack propagation through the fibrous network. These results can be used to enhance the current structural and mechanical models for the prediction of the fracture behaviour of papers.
]]>2018-07-11T00:09:41-07:00info:doi/10.1098/rspa.2018.0175hwp:master-id:royprsa;rspa.2018.01752018-07-11Research articles47422152018017520180175<![CDATA[Quantifying plasma immersion ion implantation of insulating surfaces in a dielectric barrier discharge: how to control the dose]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180263?rss=1
The plasma physics of dielectric barrier discharges (DBD) for carrying out ion implantation in insulators is investigated. A hollow cathode DBD excited by high-voltage pulses is suitable for ion bombardment of the surfaces of insulating tubing, porous material, particles and sheets. Plasma immersion ion implantation of insulating surfaces is useful for many applications in medicine and engineering. The ion bombardment of glass is useful for cleaning and surface modification. The ion implantation of polymers creates radicals that are able to bind molecules to their surfaces for applications in medical procedures and diagnostics. A wire diagnostic probe and optical emission spectroscopy are used for experimental work. A theory based on mutual capacitance is developed to convert data from the probe to give implanted charge as a function of pressure, voltage and pulse duration. We find the operating conditions that allow for charge to be implanted and those that achieve the highest implanted charge.
]]>2018-07-25T00:05:18-07:00info:doi/10.1098/rspa.2018.0263hwp:master-id:royprsa;rspa.2018.02632018-07-25Research articles47422152018026320180263<![CDATA[Providing sustainable catalytic solutions for a rapidly changing world]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2215/20180414?rss=1
2018-07-25T00:05:18-07:00info:doi/10.1098/rspa.2018.0414hwp:master-id:royprsa;rspa.2018.04142018-07-25Editorial47422152018041420180414<![CDATA[Water-limited vegetated ecosystems driven by stochastic rainfall: feedbacks and bimodality]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170649?rss=1
In arid or semi-arid ecosystems, water availability is one of the primary controls on vegetation growth. When subsurface water resources are unavailable, the vegetation growth is dictated by the rainfall, and the random nature of the rainfall arrivals and quantities induces a probability distribution of soil moisture and vegetation biomass via the coupled dynamic equations of biomass balance and water balance. We have previously obtained an exact solution for these distributions under certain conditions, and shown that the mapping of rainfall variability to observed biomass variability can be successfully applied to a field site. Here, we expand upon our earlier theoretical work to show how the dynamics can give rise to more complicated, bimodal (and multimodal) structures in the biomass distribution when positive feedbacks between growth and water availability are included. We also derive some new analytical results for the crossing properties of this system, which enable us to determine on what time scale the effects of these feedbacks will be felt, and, relatedly, how long the system will take to cross between different modes.
]]>2018-06-20T00:09:52-07:00info:doi/10.1098/rspa.2017.0649hwp:master-id:royprsa;rspa.2017.06492018-06-20Research articles47422142017064920170649<![CDATA[Nematic liquid crystals on curved surfaces: a thin film limit]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170686?rss=1
We consider a thin film limit of a Landau–de Gennes Q-tensor model. In the limiting process, we observe a continuous transition where the normal and tangential parts of the Q-tensor decouple and various intrinsic and extrinsic contributions emerge. The main properties of the thin film model, like uniaxiality and parameter phase space, are preserved in the limiting process. For the derived surface Landau–de Gennes model, we consider an L^{2}-gradient flow. The resulting tensor-valued surface partial differential equation is numerically solved to demonstrate realizations of the tight coupling of elastic and bulk free energy with geometric properties.
]]>2018-06-20T00:09:52-07:00info:doi/10.1098/rspa.2017.0686hwp:master-id:royprsa;rspa.2017.06862018-06-20Research articles47422142017068620170686<![CDATA[Exact eigenstates of a nanometric paraboloidal emitter and field emission quantities]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170692?rss=1
The progress in field emission theory from its initial Fowler–Nordheim form is centred on the transmission coefficient. For the supply (of electrons) function one still uses the constant value due to a supply of plane-waves states. However, for emitting tips of apex radius of 1–5 nm this is highly questionable. To address this issue, we have solved the Schrödinger equation in a sharp paraboloidally shaped quantum box. The Schrödinger equation is separable in the rotationally parabolic coordinate system and we hence obtain the exact eigenstates of the system. Significant differences from the usual Cartesian geometry are obtained. (1) Both the normally incident and parallel electron fluxes are functions of the angle to the emitter axis and affect the emission angle. (2) The WKB approximation fails for this system. (3) The eigenfunctions of the nanoemitter form a continuum only in one dimension while complete discretization occurs in the other two directions. (4) The parallel electron velocity vanishes at the apex which may explain the recent spot-size measurements in near-field scanning electron microscopy. (5) Competing effects are found as the tip radius decreases to 1 nm: The electric field increases but the total supply function decreases so that possibly an optimum radius exists.
]]>2018-06-13T00:09:52-07:00info:doi/10.1098/rspa.2017.0692hwp:master-id:royprsa;rspa.2017.06922018-06-13Research articles47422142017069220170692<![CDATA[A no-go theorem for theories that decohere to quantum mechanics]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170732?rss=1
To date, there has been no experimental evidence that invalidates quantum theory. Yet it may only be an effective description of the world, in the same way that classical physics is an effective description of the quantum world. We ask whether there exists an operationally defined theory superseding quantum theory, but which reduces to it via a decoherence-like mechanism. We prove that no such post-quantum theory exists if it is demanded that it satisfy two natural physical principles: causality and purification. Causality formalizes the statement that information propagates from present to future, and purification that each state of incomplete information arises in an essentially unique way due to lack of information about an environment. Hence, our result can be viewed either as evidence that the fundamental theory of Nature is quantum or as showing in a rigorous manner that any post-quantum theory must abandon causality, purification or both.
]]>2018-06-27T00:09:55-07:00info:doi/10.1098/rspa.2017.0732hwp:master-id:royprsa;rspa.2017.07322018-06-27Research articles47422142017073220170732<![CDATA[Surface theory in discrete projective differential geometry. I. A canonical frame and an integrable discrete Demoulin system]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170770?rss=1
We present the first steps of a procedure which discretizes surface theory in classical projective differential geometry in such a manner that underlying integrable structure is preserved. We propose a canonical frame in terms of which the associated projective Gauss-Weingarten and Gauss-Mainardi-Codazzi equations adopt compact forms. Based on a scaling symmetry which injects a parameter into the linear Gauss-Weingarten equations, we set down an algebraic classification scheme of discrete projective minimal surfaces which turns out to admit a geometric counterpart formulated in terms of discrete notions of Lie quadrics and their envelopes. In the case of discrete Demoulin surfaces, we derive a Bäcklund transformation for the underlying discrete Demoulin system and show how the latter may be formulated as a two-component generalization of the integrable discrete Tzitzéica equation which has originally been derived in a different context. At the geometric level, this connection leads to the retrieval of the standard discretization of affine spheres in affine differential geometry.
]]>2018-06-13T00:17:57-07:00info:doi/10.1098/rspa.2017.0770hwp:master-id:royprsa;rspa.2017.07702018-06-13Research articles47422142017077020170770<![CDATA[How realistic are painted lightnings? Quantitative comparison of the morphology of painted and real lightnings: a psychophysical approach]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170859?rss=1
Inspired by the pioneer work of the nineteenth century photographer, William Nicholson Jennings, we studied quantitatively how realistic painted lightnings are. In order to answer this question, we examined 100 paintings and 400 photographs of lightnings. We used our software package to process and evaluate the morphology of lightnings. Three morphological parameters of the main lightning branch were analysed: (i) number of branches N_{b}, (ii) relative length r, and (iii) number of local maxima (peaks) N_{p} of the turning angle distribution. We concluded: (i) Painted lightnings differ from real ones in N_{b} and N_{p}. (ii) The r-values of painted and real lightnings vary in the same range. (iii) 67 and 22% of the studied painted and real lightnings were non-bifurcating (N_{b} = 1, meaning only the main branch), the maximum of N_{b} of painted and real lightnings is 11 and 51, respectively, and painted bifurcating lightnings possess mostly 2–4 branches, while real lightnings have mostly 2–10 branches. To understand these findings, we performed two psychophysical experiments with 10 test persons, whose task was to guess N_{b} on photographs of real lightnings which were flashed for short time periods t = 0.5, 0.75 and 1 s (characteristic to lightnings) on a monitor. We obtained that (i) test persons can estimate the number of lightning branches quite correctly if N_{b} ≤ 11. (ii) If N_{b} > 11, its value is strongly underestimated with exponentially increasing difference between the real and estimated numbers. (iii) The estimation is independent of the flashing period t of lightning photos/pictures. (iv) The estimation is more accurate, if skeletonized lightning pictures are flashed, rather than real lightning photos. These findings explain why artists usually illustrate lightnings with branches not larger than 11.
]]>2018-06-06T00:09:54-07:00info:doi/10.1098/rspa.2017.0859hwp:master-id:royprsa;rspa.2017.08592018-06-06Research47422142017085920170859<![CDATA[Force appropriation of nonlinear structures]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20170880?rss=1
Nonlinear normal modes (NNMs) are widely used as a tool for developing mathematical models of nonlinear structures and understanding their dynamics. NNMs can be identified experimentally through a phase quadrature condition between the system response and the applied excitation. This paper demonstrates that this commonly used quadrature condition can give results that are significantly different from the true NNM, in particular, when the excitation applied to the system is limited to one input force, as is frequently used in practice. The system studied is a clamped–clamped cross-beam with two closely spaced modes. This paper shows that the regions where the quadrature condition is (in)accurate can be qualitatively captured by analysing transfer of energy between the modes of the system, leading to a discussion of the appropriate number of input forces and their locations across the structure.
]]>2018-06-13T00:09:52-07:00info:doi/10.1098/rspa.2017.0880hwp:master-id:royprsa;rspa.2017.08802018-06-13Research articles47422142017088020170880<![CDATA[An asymptotic higher-order theory for rectangular beams]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180001?rss=1
A direct asymptotic integration of the full three-dimensional problem of elasticity is employed to derive a consistent governing equation for a beam with the rectangular cross section. The governing equation is consistent in the sense that it has the same long-wave low-frequency behaviour as the exact solution of the original three-dimensional problem. Performance of the new beam equation is illustrated by comparing its predictions against the results of direct finite-element computations. Limiting behaviours for beams with large (and small) aspect ratios, which can be established using classical plate theories, are recovered from the new governing equation to illustrate its consistency and also to illustrate the importance of using plate theories with the correctly refined boundary conditions. The implications for the correct choice of the shear correction factor in Timoshenko's beam theory are also discussed.
]]>2018-06-13T00:09:52-07:00info:doi/10.1098/rspa.2018.0001hwp:master-id:royprsa;rspa.2018.00012018-06-13Research articles47422142018000120180001<![CDATA[An algorithm to explore entanglement in small systems]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180023?rss=1
A quantum state’s entanglement across a bipartite cut can be quantified with entanglement entropy or, more generally, Schmidt norms. Using only Schmidt decompositions, we present a simple iterative algorithm to maximize Schmidt norms. Depending on the choice of norm, the optimizing states maximize or minimize entanglement, possibly across several bipartite cuts at the same time and possibly only among states in a specified subspace. Recognizing that convergence but not success is certain, we use the algorithm to explore topics ranging from fermionic reduced density matrices and varieties of pure quantum states to absolutely maximally entangled states and minimal output entropy of channels.
]]>2018-06-13T00:17:57-07:00info:doi/10.1098/rspa.2018.0023hwp:master-id:royprsa;rspa.2018.00232018-06-13Research articles47422142018002320180023<![CDATA[Critical curvature localization in graphene. I. Quantum-flexoelectricity effect]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180054?rss=1
Here, we report the discovery of a new, curvature-localizing, subcritical buckling mode that produces shallow-kink corrugation in multi-layer graphene. Our density functional theory (DFT) analysis reveals the mode configuration—an approximately 2 nm wide boundary layer of highly localized curvature that connects two regions of uniformly but oppositely sheared stacks of flat atomic sheets. The kink angle between the two regions is limited to a few degrees, ensuring elastic deformation. By contrast, a purely mechanical model of sandwich structures shows progressive supercritical curvature localization spread over a 50–100 nm wide boundary layer. Our effective-locality model of electromechanics reveals that coupling between atomic-layer curvature and electric-charge polarization, i.e. quantum flexoelectricity, leads to emergence of a boundary layer in which curvature is focused primarily within a 0.86 nm fixed band width. Both DFT and the model analyses show focused distributions of curvature and polarization exhibiting oscillating decay within the approximately 2 nm wide boundary layer. The results show that dipole–dipole interaction lowers the potential energy with such a distribution. Furthermore, this model predicts peak-polarization density approximately 0.12 e^{-} nm^{–1} for 3° tilt angle. This high polarization concentration can be controlled by macroscopic deformation and is expected to be useful in studies of selective graphene-surface functionalization for various applications.
]]>2018-06-27T00:37:07-07:00info:doi/10.1098/rspa.2018.0054hwp:master-id:royprsa;rspa.2018.00542018-06-27Research articles47422142018005420180054<![CDATA[Tension-dependent transverse buckles and wrinkles in twisted elastic sheets]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180062?rss=1
We investigate with experiments the twist-induced transverse buckling instabilities of an elastic sheet of length L, width W and thickness t, that is clamped at two opposite ends while held under a tension T. Above a critical tension T_{} and critical twist angle _{tr}, we find that the sheet buckles with a mode number n≥1 transverse to the axis of twist. Three distinct buckling regimes characterized as clamp-dominated, bendable and stiff are identified, by introducing a bendability length L_{B} and a clamp length L_{C}(<L_{B}). In the stiff regime (L>L_{B}), we find that mode n=1 develops above _{tr}_{S}~(t/W)T^{–1/2}, independent of L. In the bendable regime L_{C}<L<L_{B}, n=1 as well as n>1 occur above trB~t/LT–1/4. Here, we find the wavelength B~LtT–1/4, when n>1. These scalings agree with those derived from a covariant form of the Föppl-von Kármán equations, however, we find that the n=1 mode also occurs over a surprisingly large range of L in the bendable regime. Finally, in the clamp-dominated regime (L<L_{C}), we find that _{tr} is higher compared to _{B} due to additional stiffening induced by the clamped boundary conditions.
]]>2018-06-13T00:09:52-07:00info:doi/10.1098/rspa.2018.0062hwp:master-id:royprsa;rspa.2018.00622018-06-13Research articles47422142018006220180062<![CDATA[A damage initiation criterion for a class of viscoelastic solids]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180064?rss=1
We extend the methodology introduced for the initiation of damage within the context of a class of elastic solids to a class of viscoelastic solids (Alagappan et al. 2016 Proc. R. Soc. Lond. A: Math. Phys. Eng. Sci.472, 20160231. (doi:10.1098/rspa.2016.0231)). In a departure from studies on damage that consider the body to be homogeneous, with initiation of damage being decided by parameters that are based on a quantity such as the strain, that requires information concerning a special reference configuration, or the use of ad hoc parameters that have no physically meaningful origins, in this study we use a physically relevant parameter that is completely determined in the current deformed state of the body to predict the initiation of damage. Damage is initiated due to the inhomogeneity of the body wherein certain regions in the body are unable to withstand the stresses, strains, etc. The specific inhomogeneity that is considered is the variation of the density in the body. We consider damage within the context of the deformation of two representative viscoelastic solids, a generalization of a model proposed by Gent (1996 Rubber Chemistry and Technology69, 59–61. (doi:10.5254/1.3538357)) for polymeric solids and a generalization of the Kelvin–Voigt model. We find that the criterion leads to results that are in keeping with the experiments of Gent & Lindley (1959 Proc. R. Soc. Lond. A: Math. Phys. Eng. Sci.249, 195–205. (doi:10.1098/rspa.1959.0016)).
]]>2018-06-27T00:09:55-07:00info:doi/10.1098/rspa.2018.0064hwp:master-id:royprsa;rspa.2018.00642018-06-27Research articles47422142018006420180064<![CDATA[The nature of Earth's correlation wavefield: late coda of large earthquakes]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180082?rss=1
The seismic correlation wavefield constructed from the stacked cross-correlograms of the late coda of earthquake signals at stations across the globe provides a wealth of observed pulses as a function of inter-station distance. The interval from 3 to 10 h after the onset of major earthquakes is employed for the period range from 15 to 50 s. The observations can be well matched by synthetic seismograms for a radially stratified Earth. Many of the correlation phases have similar time behaviour to those in the regular wavefield, but others have no correspondence. All such correlation phases can be explained by the interaction of arrivals with a common slowness at the each of the stations being correlated. Using a generalized ray description of the seismic wavefield, the time-distance behaviour of these correlation phases arises from differences in accumulated phase on different propagation paths through the Earth. Distinct arrivals emerge from the correlation field when there are many ways in which combinations of seismic phases can arise with the same difference in propagation legs. The constituents of the late coda are dominated by steeply travelling waves, and in consequence features associated with multiple passages through the whole Earth emerge distinctly, such as high-order multiples of PKIKP.
]]>2018-06-13T00:09:52-07:00info:doi/10.1098/rspa.2018.0082hwp:master-id:royprsa;rspa.2018.00822018-06-13Research articles47422142018008220180082<![CDATA[Composite wave models for elastic plates]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180103?rss=1
The long-term challenge of formulating an asymptotically motivated wave theory for elastic plates is addressed. Composite two-dimensional models merging the leading or higher-order parabolic equations for plate bending and the hyperbolic equation for the Rayleigh surface wave are constructed. Analysis of numerical examples shows that the proposed approach is robust not only at low- and high-frequency limits but also over the intermediate frequency range.
]]>2018-06-20T00:09:53-07:00info:doi/10.1098/rspa.2018.0103hwp:master-id:royprsa;rspa.2018.01032018-06-20Research articles47422142018010320180103<![CDATA[Bond graph modelling of the cardiac action potential: implications for drift and non-unique steady states]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2214/20180106?rss=1
Mathematical models of cardiac action potentials have become increasingly important in the study of heart disease and pharmacology, but concerns linger over their robustness during long periods of simulation, in particular due to issues such as model drift and non-unique steady states. Previous studies have linked these to violation of conservation laws, but only explored those issues with respect to charge conservation in specific models. Here, we propose a general and systematic method of identifying conservation laws hidden in models of cardiac electrophysiology by using bond graphs, and develop a bond graph model of the cardiac action potential to study long-term behaviour. Bond graphs provide an explicit energy-based framework for modelling physical systems, which makes them well suited for examining conservation within electrophysiological models. We find that the charge conservation laws derived in previous studies are examples of the more general concept of a ‘conserved moiety’. Conserved moieties explain model drift and non-unique steady states, generalizing the results from previous studies. The bond graph approach provides a rigorous method to check for drift and non-unique steady states in a wide range of cardiac action potential models, and can be extended to examine behaviours of other excitable systems.
]]>2018-06-27T00:09:55-07:00info:doi/10.1098/rspa.2018.0106hwp:master-id:royprsa;rspa.2018.01062018-06-27Research articles47422142018010620180106<![CDATA[Determining erosion rates in Allchar (Macedonia) to revive the lorandite neutrino experiment]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170470?rss=1
^{205}Tl in the lorandite (TiAsS_{2}) mine of Allchar (Majdan, FYR Macedonia) is transformed to ^{205}Pb by cosmic ray reactions with muons and neutrinos. At depths of more than 300 m, muogenic production would be sufficiently low for the 4.3 Ma old lorandite deposit to be used as a natural neutrino detector. Unfortunately, the Allchar deposit currently sits at a depth of only 120 m below the surface, apparently making the lorandite experiment technically infeasible. We here present 25 erosion rate estimates for the Allchar area using in situ produced cosmogenic ^{36}Cl in carbonates and ^{10}Be in alluvial quartz. The new measurements suggest long-term erosion rates of 100–120 m Ma^{–1} in the silicate lithologies that are found at the higher elevations of the Majdanksa River valley, and 200–280 m Ma^{–1} in the underlying marbles and dolomites. These values indicate that the lorandite deposit has spent most of its existence at depths of more than 400 m, sufficient for the neutrinogenic ^{205}Pb component to dominate the muon contribution. Our results suggest that this unique particle physics experiment is theoretically feasible and merits further development.
]]>2018-05-16T00:26:05-07:00info:doi/10.1098/rspa.2017.0470hwp:master-id:royprsa;rspa.2017.04702018-05-16Research articles47422132017047020170470<![CDATA[A mathematical model for fitting and predicting relaxation modulus and simulating viscoelastic responses]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170540?rss=1
We propose a mathematical model for relaxation modulus and its numerical solution. The model formula is extended from sigmoidal function considering nonlinear strain hardening. Its physical meaning can be interpreted by a macroscale elastic network-viscous medium model with only five model parameters in a simpler format than the molecular-chain-based polymer models to represent general solid materials. We also developed a finite-element (FE) framework and robust numerical algorithm to implement this model for simulating responses under both static and dynamic loadings. We validated the model through both experimental data and numerical simulations on a variety of materials including asphalt concrete, polymer, spider silk, hydrogel, agar and bone. By satisfying the second law of thermodynamics in the form of Calusius–Duhem inequality, the model is able to simulate creep and sinusoidal deformation as well as energy dissipation. Compared to the Prony series, the widely used model with a large number of model parameters, the proposed model has improved accuracy in fitting experimental data and prediction stability outside of the experimental range with competitive numerical stability and computation speed. We also present simulation results of nonlinear stress–strain relationships of spider silk and hydrogels, and dynamic responses of a multilayer structure.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2017.0540hwp:master-id:royprsa;rspa.2017.05402018-05-16Research articles47422132017054020170540<![CDATA[Power packet transferability via symbol propagation matrix]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170552?rss=1
A power packet is a unit of electric power composed of a power pulse and an information tag. In Shannon’s information theory, messages are represented by symbol sequences in a digitized manner. Referring to this formulation, we define symbols in power packetization as a minimum unit of power transferred by a tagged pulse. Here, power is digitized and quantized. In this paper, we consider packetized power in networks for a finite duration, giving symbols and their energies to the networks. A network structure is defined using a graph whose nodes represent routers, sources and destinations. First, we introduce the concept of a symbol propagation matrix (SPM) in which symbols are transferred at links during unit times. Packetized power is described as a network flow in a spatio-temporal structure. Then, we study the problem of selecting an SPM in terms of transferability, that is, the possibility to represent given energies at sources and destinations during the finite duration. To select an SPM, we consider a network flow problem of packetized power. The problem is formulated as an M-convex submodular flow problem which is a solvable generalization of the minimum cost flow problem. Finally, through examples, we verify that this formulation provides reasonable packetized power.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2017.0552hwp:master-id:royprsa;rspa.2017.05522018-05-16Research articles47422132017055220170552<![CDATA[Multi-vortex crystal lattices in Bose-Einstein condensates with a rotating trap]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170553?rss=1
We consider vortex dynamics in the context of Bose–Einstein condensates (BECs) with a rotating trap, with or without anisotropy. Starting with the Gross–Pitaevskii (GP) partial differential equation (PDE), we derive a novel reduced system of ordinary differential equations (ODEs) that describes stable configurations of multiple co-rotating vortices (vortex crystals). This description is found to be quite accurate quantitatively especially in the case of multiple vortices. In the limit of many vortices, BECs are known to form vortex crystal structures, whereby vortices tend to arrange themselves in a hexagonal-like spatial configuration. Using our asymptotic reduction, we derive the effective vortex crystal density and its radius. We also obtain an asymptotic estimate for the maximum number of vortices as a function of rotation rate. We extend considerations to the anisotropic trap case, confirming that a pair of vortices lying on the long (short) axis is linearly stable (unstable), corroborating the ODE reduction results with full PDE simulations. We then further investigate the many-vortex limit in the case of strong anisotropic potential. In this limit, the vortices tend to align themselves along the long axis, and we compute the effective one-dimensional vortex density, as well as the maximum admissible number of vortices. Detailed numerical simulations of the GP equation are used to confirm our analytical predictions.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2017.0553hwp:master-id:royprsa;rspa.2017.05532018-05-16Research articles47422132017055320170553<![CDATA[Revealing new dynamical patterns in a reaction-diffusion model with cyclic competition via a novel computational framework]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170608?rss=1
Understanding how patterns and travelling waves form in chemical and biological reaction–diffusion models is an area which has been widely researched, yet is still experiencing fast development. Surprisingly enough, we still do not have a clear understanding about all possible types of dynamical regimes in classical reaction–diffusion models, such as Lotka–Volterra competition models with spatial dependence. In this study, we demonstrate some new types of wave propagation and pattern formation in a classical three species cyclic competition model with spatial diffusion, which have been so far missed in the literature. These new patterns are characterized by a high regularity in space, but are different from patterns previously known to exist in reaction–diffusion models, and may have important applications in improving our understanding of biological pattern formation and invasion theory. Finding these new patterns is made technically possible by using an automatic adaptive finite element method driven by a novel a posteriori error estimate which is proved to provide a reliable bound for the error of the numerical method. We demonstrate how this numerical framework allows us to easily explore the dynamical patterns in both two and three spatial dimensions.
]]>2018-05-23T00:10:03-07:00info:doi/10.1098/rspa.2017.0608hwp:master-id:royprsa;rspa.2017.06082018-05-23Research articles47422132017060820170608<![CDATA[Time periodic solutions of one-dimensional forced Kirchhoff equations with x-dependent coefficients]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170620?rss=1
In this paper, we consider the one-dimensional Kirchhoff equation with x-dependent coefficients under Dirichlet boundary conditions, which models the forced vibrations of a clamped inhomogeneous string in the presence of a time periodic external forcing with period 2/ and amplitude . By using the Nash–Moser iteration technique, we obtain the existence, regularity and local uniqueness of time periodic solutions with period 2/ and order . Such results hold for parameters (,) in a positive measure Cantor set that has asymptotically full measure as the amplitude goes to zero.
]]>2018-05-02T00:10:07-07:00info:doi/10.1098/rspa.2017.0620hwp:master-id:royprsa;rspa.2017.06202018-05-02Research articles47422132017062020170620<![CDATA[On the dynamic homogenization of periodic media: Willis approach versus two-scale paradigm]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170638?rss=1
When considering an effective, i.e. homogenized description of waves in periodic media that transcends the usual quasi-static approximation, there are generally two schools of thought: (i) the two-scale approach that is prevalent in mathematics and (ii) the Willis’ homogenization framework that has been gaining popularity in engineering and physical sciences. Notwithstanding a mounting body of literature on the two competing paradigms, a clear understanding of their relationship is still lacking. In this study, we deploy an effective impedance of the scalar wave equation as a lens for comparison and establish a low-frequency, long-wavelength dispersive expansion of the Willis’ effective model, including terms up to the second order. Despite the intuitive expectation that such obtained effective impedance coincides with its two-scale counterpart, we find that the two descriptions differ by a modulation factor which is, up to the second order, expressible as a polynomial in frequency and wavenumber. We track down this inconsistency to the fact that the two-scale expansion is commonly restricted to the free-wave solutions and thus fails to account for the body source term which, as it turns out, must also be homogenized—by the reciprocal of the featured modulation factor. In the analysis, we also (i) reformulate for generality the Willis’ effective description in terms of the eigenfunction approach, and (ii) obtain the corresponding modulation factor for dipole body sources, which may be relevant to some recent efforts to manipulate waves in metamaterials.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2017.0638hwp:master-id:royprsa;rspa.2017.06382018-05-09Research articles47422132017063820170638<![CDATA[The correct and unusual coordinate transformation rules for electromagnetic quadrupoles]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170652?rss=1
Despite being studied for over a century, the use of quadrupoles have been limited to Cartesian coordinates in flat space–time due to the incorrect transformation rules used to define them. Here the correct transformation rules are derived, which are particularly unusual as they involve second derivatives of the coordinate transformation and an integral. Transformations involving integrals have not been seen before. This is significantly different from the familiar transformation rules for a dipole, where the components transform as tensors. It enables quadrupoles to be correctly defined in general relativity and to prescribe the equations of motion for a quadrupole in a coordinate system adapted to its motion and then transform them to the laboratory coordinates. An example is given of another unusual feature: a quadrupole which is free of dipole terms in polar coordinates has dipole terms in Cartesian coordinates. It is shown that dipoles, electric dipoles, quadrupoles and electric quadrupoles can be defined without reference to a metric and in a coordinates-free manner. This is particularly useful given their complicated coordinate transformation.
]]>2018-05-09T02:05:01-07:00info:doi/10.1098/rspa.2017.0652hwp:master-id:royprsa;rspa.2017.06522018-05-09Research articles47422132017065220170652<![CDATA[Stochastic modelling of urban structure]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170700?rss=1
The building of mathematical and computer models of cities has a long history. The core elements are models of flows (spatial interaction) and the dynamics of structural evolution. In this article, we develop a stochastic model of urban structure to formally account for uncertainty arising from less predictable events. Standard practice has been to calibrate the spatial interaction models independently and to explore the dynamics through simulation. We present two significant results that will be transformative for both elements. First, we represent the structural variables through a single potential function and develop stochastic differential equations to model the evolution. Second, we show that the parameters of the spatial interaction model can be estimated from the structure alone, independently of flow data, using the Bayesian inferential framework. The posterior distribution is doubly intractable and poses significant computational challenges that we overcome using Markov chain Monte Carlo methods. We demonstrate our methodology with a case study on the London, UK, retail system.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2017.0700hwp:master-id:royprsa;rspa.2017.07002018-05-09Research articles47422132017070020170700<![CDATA[Quantum stopwatch: how to store time in a quantum memory]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170773?rss=1
Quantum mechanics imposes a fundamental trade-off between the accuracy of time measurements and the size of the systems used as clocks. When the measurements of different time intervals are combined, the errors due to the finite clock size accumulate, resulting in an overall inaccuracy that grows with the complexity of the set-up. Here, we introduce a method that, in principle, eludes the accumulation of errors by coherently transferring information from a quantum clock to a quantum memory of the smallest possible size. Our method could be used to measure the total duration of a sequence of events with enhanced accuracy, and to reduce the amount of quantum communication needed to stabilize clocks in a quantum network.
]]>2018-05-23T00:10:03-07:00info:doi/10.1098/rspa.2017.0773hwp:master-id:royprsa;rspa.2017.07732018-05-23Research articles47422132017077320170773<![CDATA[Stability transitions of an axially moving string subjected to a distributed follower force]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170779?rss=1
The transverse vibrations of an axially moving string that is subjected to a distributed follower force are examined here. This model provides an insight into the complex dynamics of seemingly simpler systems such as silicon wafer cutting using wire saws, and aerial or marine towing, where a relatively long flexible structure is dragged through fluid. The equation of motion is derived and it includes the axial variation in the tension that arises due to acceleration and the follower force. As the exact analytical solution of this equation is difficult to determine, the approximate closed-form modal solution of a non-travelling counterpart of the system is obtained using the asymptotic technique, which is then used as a basis to obtain the numerical solution for the axially moving string. The effect of the follower force and viscous dissipation on the eigenstructure of the system is investigated. Mathematical operations such as the Hermite form and the Routh–Hurwitz criterion are applied to the characteristic polynomial to investigate the dynamic behaviour of these modes. The semi-analytical approach presented explains the ‘mathematical’ instability (in the absence of damping) that arises when both axial transport and follower force are simultaneously present. An unusual transition of the dynamic behaviour from the stable to the overdamped and then directly to the unstable regime is observed.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2017.0779hwp:master-id:royprsa;rspa.2017.07792018-05-09Research articles47422132017077920170779<![CDATA[Simultaneity of centres in Zq-equivariant systems]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170811?rss=1
We study the simultaneous existence of centres for two families of planar Zq-equivariant systems. First, we give a short review about Zq-equivariant systems. Next, we present the necessary and sufficient conditions for the simultaneous existence of centres for a Z2-equivariant cubic system and for a Z2-equivariant quintic system.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2017.0811hwp:master-id:royprsa;rspa.2017.08112018-05-09Research articles47422132017081120170811<![CDATA[Data-driven forecasting of high-dimensional chaotic systems with long short-term memory networks]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170844?rss=1
We introduce a data-driven forecasting method for high-dimensional chaotic systems using long short-term memory (LSTM) recurrent neural networks. The proposed LSTM neural networks perform inference of high-dimensional dynamical systems in their reduced order space and are shown to be an effective set of nonlinear approximators of their attractor. We demonstrate the forecasting performance of the LSTM and compare it with Gaussian processes (GPs) in time series obtained from the Lorenz 96 system, the Kuramoto–Sivashinsky equation and a prototype climate model. The LSTM networks outperform the GPs in short-term forecasting accuracy in all applications considered. A hybrid architecture, extending the LSTM with a mean stochastic model (MSM–LSTM), is proposed to ensure convergence to the invariant measure. This novel hybrid method is fully data-driven and extends the forecasting capabilities of LSTM networks.
]]>2018-05-23T00:10:03-07:00info:doi/10.1098/rspa.2017.0844hwp:master-id:royprsa;rspa.2017.08442018-05-23Research articles47422132017084420170844<![CDATA[Surface instability of elastic half-spaces by using the energy method]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170854?rss=1
Finding the complete set of stability conditions of an elastic half-space has been an open problem ever since Biot (Biot 1963 Appl. Sci. Res.12, 168–182 (doi:10.1007/BF03184638)) first studied the surface instability of half-spaces by seeking solutions of the incremental equilibrium equations. Towards solving this problem, a method based on the energy stability criterion is developed in the present work. A variational problem of minimizing the elastic energy associated with a half-space is formulated. The second variation condition is derived and is converted to an eigenvalue problem. For a half-space of neo-Hookean materials, the eigenvalue problem is solved, which leads to complete descriptions of stability and instability regions in the deformation space.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2017.0854hwp:master-id:royprsa;rspa.2017.08542018-05-09Research articles47422132017085420170854<![CDATA[Static bistability of spherical caps]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20170910?rss=1
Depending on its geometry, a spherical shell may exist in one of two stable states without the application of any external force: there are two ‘self-equilibrated’ states, one natural and the other inside out (or ‘everted’). Though this is familiar from everyday life—an umbrella is remarkably stable, yet a contact lens can be easily turned inside out—the precise shell geometries for which bistability is possible are not known. Here, we use experiments and finite-element simulations to determine the threshold between bistability and monostability for shells of different solid angle. We compare these results with the prediction from shallow shell theory, showing that, when appropriately modified, this offers a very good account of bistability even for relatively deep shells. We then investigate the robustness of this bistability against pointwise indentation. We find that indentation provides a continuous route for transition between the two states for shells whose geometry makes them close to the threshold. However, for thinner shells, indentation leads to asymmetrical buckling before snap-through, while also making these shells more ‘robust’ to snap-through. Our work sheds new light on the robustness of the ‘mirror buckling’ symmetry of spherical shell caps.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2017.0910hwp:master-id:royprsa;rspa.2017.09102018-05-16Research articles47422132017091020170910<![CDATA[Horizontally isotropic bidispersive thermal convection]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180018?rss=1
A bidispersive porous material is one which has usual pores but additionally contains a system of micro pores due to cracks or fissures in the solid skeleton. We present general equations for thermal convection in a bidispersive porous medium when the permeabilities, interaction coefficient and thermal conductivity are anisotropic but symmetric tensors. In this case, we show exchange of stabilities holds and fluid movement will commence via stationary convection, and additionally we show the global nonlinear stability threshold is the same as the linear instability one. Attention is then focused on the case where the interaction coefficient and thermal conductivity are isotropic, and the permeability is isotropic in the horizontal directions, although the permeability in the vertical direction is different. The nonlinear stability threshold is calculated in this case and numerical results are presented and discussed in detail.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2018.0018hwp:master-id:royprsa;rspa.2018.00182018-05-16Research articles47422132018001820180018<![CDATA[Trapped modes and reflectionless modes as eigenfunctions of the same spectral problem]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180050?rss=1
We consider the reflection–transmission problem in a waveguide with obstacle. At certain frequencies, for some incident waves, intensity is perfectly transmitted and the reflected field decays exponentially at infinity. In this work, we show that such reflectionless modes can be characterized as eigenfunctions of an original non-self-adjoint spectral problem. To select ingoing waves on one side of the obstacle and outgoing waves on the other side, we use complex scalings (or perfectly matched layers) with imaginary parts of different signs. We prove that the real eigenvalues of the obtained spectrum correspond either to trapped modes (or bound states in the continuum) or to reflectionless modes. Interestingly, complex eigenvalues also contain useful information on weak reflection cases. When the geometry has certain symmetries, the new spectral problem enters the class of PT-symmetric problems.
]]>2018-05-09T00:15:38-07:00info:doi/10.1098/rspa.2018.0050hwp:master-id:royprsa;rspa.2018.00502018-05-09Research articles47422132018005020180050<![CDATA[New variational and multisymplectic formulations of the Euler-Poincare equation on the Virasoro-Bott group using the inverse map]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180052?rss=1
We derive a new variational principle, leading to a new momentum map and a new multisymplectic formulation for a family of Euler–Poincaré equations defined on the Virasoro–Bott group, by using the inverse map (also called ‘back-to-labels’ map). This family contains as special cases the well-known Korteweg–de Vries, Camassa–Holm and Hunter–Saxton soliton equations. In the conclusion section, we sketch opportunities for future work that would apply the new Clebsch momentum map with 2-cocycles derived here to investigate a new type of interplay among nonlinearity, dispersion and noise.
]]>2018-05-09T00:10:05-07:00info:doi/10.1098/rspa.2018.0052hwp:master-id:royprsa;rspa.2018.00522018-05-09Research articles47422132018005220180052<![CDATA[Explicit backbone curves from spectral submanifolds of forced-damped nonlinear mechanical systems]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180083?rss=1
Spectral submanifolds (SSMs) have recently been shown to provide exact and unique reduced-order models for nonlinear unforced mechanical vibrations. Here, we extend these results to periodically or quasi-periodically forced mechanical systems, obtaining analytic expressions for forced responses and backbone curves on modal (i.e. two dimensional) time-dependent SSMs. A judicious choice of the parametrization of these SSMs allows us to simplify the reduced dynamics considerably. We demonstrate our analytical formulae on three numerical examples and compare them to results obtained from available normal-form methods.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2018.0083hwp:master-id:royprsa;rspa.2018.00832018-05-16Research articles47422132018008320180083<![CDATA[Energy extraction from vortex-induced vibrations using period-1 rotation of an autoparametric pendulum]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180086?rss=1
We propose and analyse the feasibility of extracting energy from vortex-induced vibrations using rotating motion of an attached pendulum. The resulting autoparametric pendulum system is studied primarily to understand the effect of pendulum motion on the performance of the harvester which is typically ignored to result in a simple parametric pendulum. We find that rotating motions are possible only for small values of the pendulum mass when compared with the effective mass of the vibrating structure. However, the pendulum motion reduces the basin of attraction as well as the range of system parameters corresponding to the existence of rotary solutions. This significantly alters the harvester performance. By contrast, the evolution of the pendulum coordinates (angular position and velocity) remains largely unaffected by this interaction. Hence, for the purpose of design of controllers to robustly initiate/maintain rotation from arbitrary disturbances, the simplification to a parametric pendulum is reasonable while for the design of the harvester, this exercise is completely unsatisfactory.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2018.0086hwp:master-id:royprsa;rspa.2018.00862018-05-16Research articles47422132018008620180086<![CDATA[Periodic behaviour in ground-level environmental radioactivity: fingerprints of solar activity?]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180109?rss=1
Long-term data on ground-level environmental radioactivity are analysed, these referring to daily records of the total β activity in air and absorbed dose rate between 1976 and 2016 in Zagreb, Croatia. A Fourier analysis of these two datasets reveals a periodic behaviour with numerous periods ranging from approximately four months to roughly 21 years. Ninety per cent of the periods agree remarkably well with known periodicities in solar activity throughout solar cycles 21–24, which suggests that the observed modulations were most probably caused by variations in the cosmogenic radionuclide production in the upper atmosphere. Hence, one can extract the fine structure of solar activity from environmental radioactivity monitoring data, which implies that it could be possible to supplement the existing information on the exposure of Earth to cosmic rays by analysing results from radioactivity monitoring stations all over the world. These datasets cover Earth densely due to the existence of a global radioecological network and are complementary to those resulting from dedicated cosmic ray detection experiments. Because of that, such analyses might possibly lead to new findings.
]]>2018-05-30T00:09:59-07:00info:doi/10.1098/rspa.2018.0109hwp:master-id:royprsa;rspa.2018.01092018-05-30Research articles47422132018010920180109<![CDATA[Could hydrodynamic Rossby waves explain the westward drift?]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180119?rss=1
A novel theory for the origin of the westward drift of the Earth’s magnetic field is proposed, based upon the propagation of hydrodynamic Rossby waves in the liquid outer core. These waves have the obscure property that their crests always progress eastwards—but, for a certain subset, energy can nevertheless be transmitted westwards. In fact, this subset corresponds to sheet-like flow structures, extended in both the axial and radial directions, which are likely to be preferentially excited by convective upwellings in the Earth’s rapidly rotating outer core. To enable their analysis, the quasi-geostrophic (QG) approximation is employed, which assumes horizontal motions to be independent of distance along the rotation axis, yet accounts for variations in the container height (i.e. the slope of the core–mantle boundary). By projecting the momentum equation onto flows of a QG form, a general equation governing their evolution is derived, which is then adapted for the treatment of two initial value problems—in both Cartesian and spherical geometries—which demonstrate the preference for westward energy propagation by the waves in question. The merits of this mechanism as an explanation for westward drift are discussed.
]]>2018-05-16T00:10:11-07:00info:doi/10.1098/rspa.2018.0119hwp:master-id:royprsa;rspa.2018.01192018-05-16Research articles47422132018011920180119<![CDATA[Avoiding/inducing dynamic buckling in a thermomechanically coupled plate: a local and global analysis of slow/fast response]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180206?rss=1
The nonlinear response of a reduced model of an orthotropic single-layered plate with thermomechanical coupling is investigated in the presence of thermal excitations, in addition to mechanical ones. Different issues are addressed via accurate and extended local and global analyses. (i) Assessing the possible occurrence, disappearance or modification of mechanical buckling as a result of thermal aspects; (ii) exploiting global dynamics to unveil the effects of coupling; (iii) highlighting the crucial role played by the slow thermal transient evolution in modifying the fast steady mechanical response; (iv) framing the influence of coupling and underlining the need to use a thermomechanical model to grasp the actual plate dynamics; and (v) getting hints of technical interest as to the outcome robustness with respect to variations in the external/internal thermal parameters.
]]>2018-05-30T00:09:59-07:00info:doi/10.1098/rspa.2018.0206hwp:master-id:royprsa;rspa.2018.02062018-05-30Research articles47422132018020620180206<![CDATA[Correction to 'Antiplane wave scattering from a cylindrical cavity in pre-stressed nonlinear elastic media]]>
http://rspa.royalsocietypublishing.org/cgi/content/short/474/2213/20180268?rss=1
2018-05-23T00:10:03-07:00info:doi/10.1098/rspa.2018.0268hwp:master-id:royprsa;rspa.2018.02682018-05-23Correction47422132018026820180268