Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
20topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2023arXiv

Architecting Safer Autonomous Aviation Systems

The aviation literature gives relatively little guidance to practitioners about the specifics of architecting systems for safety, particularly the impact of architecture on allocating safety requirements, or the relative ease of system assurance resulting from system or subsystem level architectural choices. As an exemplar, this paper considers common architectural patterns used within traditional aviation systems and explores their safety and safety assurance implications when applied in the context of integrating artificial intelligence (AI) and machine learning (ML) based functionality. Considering safety as an architectural property, we discuss both the allocation of safety requirements and the architectural trade-offs involved early in the design lifecycle. This approach could be extended to other assured properties, similar to safety, such as security. We conclude with a discussion of the safety considerations that emerge in the context of candidate architectural patterns that have been proposed in the recent literature for enabling autonomy capabilities by integrating AI and ML. A recommendation is made for the generation of a property-driven architectural pattern catalogue.

preprint2022arXiv

Lacunarity Transition

Experiments investigating particles floating on a randomly stirred fluid show regions of very low density, which are not well understood. We introduce a simplified model for understanding sparsely occupied regions of the phase space of non-autonomous, chaotic dynamical systems, based upon an extension of the skinny bakers' map. We show how the distribution of the sizes of voids in the phase space can be mapped to the statistics of the running maximum of a Wiener process. We find that the model exhibits a lacunarity transition, which is characterised by regions of the phase space remaining empty as the number of trajectories is increased.

preprint2022arXiv

Mean square winding angle of Brownian motion around an impenetrable cylinder

An exact formula is derived, as an integral, for the mean square winding angle of Brownian motion (that is, diffusion) after time t, around an infinitely long impenetrable cylinder of radius a, having started at radius R(>a) from the axis. Strikingly, for the simpler problem with a=0, the mean square winding angle around a straight line, is long known to be instantly infinite however far away the starting point lies. the fractally small, fast, random walk steps of mathematical Brownian motion allow unbounded windings around the zero thickness of the straight line. A remedy if it is required, is to accord the line non-zero thickness, an impenetrable cylinder, as analysed here. The problem straight away reduces to a 2D one of winding around a disc in a plane since the axial component of the 3D Brownian motion is independent of the others. After deriving the exact mean square winding angle, the integral is evaluated in the limit of a narrow cylinder a<<R, highlighting the limits of short and long diffusion times addressed by previous approximate treatments.

preprint2022arXiv

Quantifying the lucky droplet model for rainfall

It is difficult to explain rainfall from ice-free clouds, because the timescale for the onset of rain showers is shorter than the mean time for collisions between microscopic water droplets. It has been suggested that raindrops are produced from very rare &#39; lucky&#39; droplets, which undergo a large number of collisions on a timescale which is short compared to the mean time for a the first collision. This work uses large deviation theory to develop estimates for the timescale for the onset of a rain shower, as a function of the collision rate coefficients. The growth history of the fast-growing droplets which do become raindrops is discussed. It is shown that their first few collisions are always approximately equally spaced in time, regardless of how the mean time for typical droplets varies as a function of the number of collisions.

preprint2020arXiv

Flooding dynamics of diffusive dispersion in a random potential

We discuss the combined effects of overdamped motion in a quenched random potential and diffusion, in one dimension, in the limit where the diffusion coefficient is small. Our analysis considers the statistics of the mean first-passage time $T(x)$ to reach position $x$, arising from different realisations of the random potential: specifically, we contrast the median $\bar T(x)$, which is an informative description of the typical course of the dispersion, with the expectation value $\langle T(x)\rangle$, which is dominated by rare events where there is an exceptionally high barrier to diffusion. We show that at relatively short times the median $\bar T(x)$ is explained by a &#39;flooding&#39; model, where $T(x)$ is predominantly determined by the highest barriers which is encountered before reaching position $x$. These highest barriers are quantified using methods of extreme value statistics.

preprint2019arXiv

Nests and Chains of Hofstadter Butterflies

The \lq Hofstadter butterfly&#39;, a plot of the spectrum of an electron in a two-dimensional periodic potential with a uniform magnetic field, contains subsets which resemble small, distorted images of the entire plot. We show how the sizes of these sub-images are determined, and calculate scaling factors describing their self-similar nesting, revealing an un-expected simplicity in the fractal structure of the spectrum. We also characterise semi-infinite chains of sub-images, showing one end of the chain is a result of gap closure, and the other end is at an accumulation point.

preprint2019arXiv

Numerical Coalescence of Chaotic Trajectories

Pairs of numerically computed trajectories of a chaotic system may coalesce because of finite arithmetic precision. We analyse an example of this phenomenon, showing that it occurs surprisingly frequently. We argue that our model belongs to a universality class of chaotic systems where this numerical coincidence effect can be described by mapping it to a first-passage process. Our results are applicable to aggregation of small particles in random flows, as well as to numerical investigation of chaotic systems.

preprint2016arXiv

Non-Thermal Einstein Relations

We consider a particle moving with equation of motion $\dot x=f(t)$, where $f(t)$ is a random function with statistics which are independent of $x$ and $t$, with a finite drift velocity $v=\langle f\rangle$ and in the presence of a reflecting wall. Far away from the wall, translational invariance implies that the stationary probability distribution is $P(x)\sim \exp(αx)$. A classical example of a problem of this type is sedimentation equilibrium, where $α$ is determined by temperature. In this work we do not introduce a thermal reservoir and $α$ is determined from the equation of motion. We consider a general approach to determining $α$ which is not always in agreement with Einstein&#39;s relation between the mean velocity and the diffusion coefficient. We illustrate our results with a model inspired by the Boltzmann equation.

preprint2015arXiv

Analysis of the Correlation Dimension for Inertial Particles

We obtain an implicit equation for the correlation dimension which describes clustering of inertial particles in a complex flow onto a fractal measure. Our general equation involves a propagator of a nonlinear stochastic process in which the velocity gradient of the fluid appears as additive noise. When the long-time limit of the propagator is considered our equation reduces to an existing large-deviation formalism, from which it is difficult to extract concrete results. In the short-time limit, however, our equation reduces to a solvability condition on a partial differential equation. We show how this approach leads to a perturbative expansion of the correlation dimension, for which the coefficients can be obtained exactly and in principle to any order. We derive the perturbation series for the correlation dimension of inertial particles suspended in three-dimensional spatially smooth random flows with white-noise time correlations, obtaining the first 33 non-zero coefficients exactly.

preprint2015arXiv

Collisional Aggregation due to Turbulence

Collisions between particles suspended in a fluid play an important role in many physical processes. As an example, collisions of microscopic water droplets in clouds are a necessary step in the production of macroscopic raindrops. Collisions of dust grains are also conjectured to be important for planet formation in the gas surrounding young stars, and also to play a role in the dynamics of sand storms. In these processes, collisions are favoured by fast turbulent motions. Here we review recent advances in the understanding of collisional aggregation due to turbulence. We discuss the role of fractal clustering of particles, and caustic singularities of their velocities. We also discuss limitations of the Smoluchowski equation for modelling these processes. These advances lead to a semi-quantitative understanding on the influence of turbulence on collision rates, and point to deficiencies in the current understanding of rainfall and planet formation.

preprint2015arXiv

Stochastic Cascade Amplification of Fluctuations

We consider a dynamical system which has a stable attractor and which is perturbed by an additive noise. Under some quite typical conditions, the fluctuations from the attractor are intermittent and have a probability distribution with power-law tails. We show that this results from a stochastic cascade of amplification of fluctuations due to transient periods of instability. The exponent of the power-law is interpreted as a negative fractal dimension.

preprint2014arXiv

A test-tube model for rainfall

If the temperature of a cell containing two partially miscible liquids is changed very slowly, so that the miscibility is decreased, microscopic droplets nucleate, grow and migrate to the interface due to their buoyancy. The system may show an approximately periodic variation of the turbidity of the mixture, as the mean droplet size fluctuates. These precipitation events are analogous to rainfall from warm clouds. This paper considers a theoretical model for these experiments. After nucleation the initial growth is by Ostwald ripening, followed by a finite-time runaway growth of droplet sizes due to larger droplets sweeping up smaller ones. The model predicts that the period $Δt$ and the temperature sweep rate $ξ$ are related by $Δt\sim C ξ^{-3/7}$, and is in good agreement with experiments. The coefficient $C$ has a power-law divergence approaching the critical point of the miscibility transition: $C\sim (T-T_{\rm c})^{-η}$, and the critical exponent $η$ is determined.

preprint2014arXiv

Advection diffusion equation with absorbing boundary

We consider a spatially homogeneous advection-diffusion equation in which the diffusion tensor and drift velocity are time-independent, but otherwise general. We derive asymptotic expressions, valid at large distances from a steady point source, for the flux onto a completely permeable boundary and onto an absorbing boundary. The absorbing case is treated by making a source of antiparticles at the boundary. In both cases there is an exponential decay as the distance from the source increases; we find that the exponent is the same for both boundary conditions.

preprint2014arXiv

Collision rate for suspensions at large Stokes numbers - comparing Navier-Stokes and synthetic turbulence

The use of simplified models of turbulent flows provides an appealing possibility to study the collision rate of turbulent suspensions, especially in conditions relevant to astrophysics, which require large time scale separations. To check the validity of such approaches, we used a direct numerical simulation (DNS) velocity field, which satisfies the Navier-Stokes equations (although it neglects the effect of the suspended particles on the flow field), and a kinematic simulation (KS) velocity field, which is a random field designed so that its statistics are in accord with the Kolmogorov theory for fully-developed turbulence. In the limit where the effects of particle inertia (characterised by the Stokes number) are negligible, the collision rates from the two approaches agree. As the Stokes number St increases, however, we show that the DNS collision rate exceeds the KS collision rate by orders of magnitude. We propose an explanation for this phenomenon and explore its consequences. We discuss the collision rate $R$ for particles in high Reynolds number flows at large Stokes number, and present evidence that $R\propto \sqrt{\rm St}$.

preprint2014arXiv

Convective Ripening and Initiation of Rainfall

This paper discusses the evolution of the droplet size distribution for a liquid-in-gas aerosol contained in a Rayleigh-Bénard cell. It introduces a non-collisional model for broadening the droplet size distribution, termed &#39;convective ripening&#39;. The paper also considers the initiation of rainfall from ice-free cumulus clouds. It is argued that while collisional mechanisms cannot explain the production of rain from clouds with water droplet diameters of $20\ μ{\rm m}$, the non-collisional convective ripening mechanism gives a much faster route to increasing the size of the small fraction of droplets that grow into raindrops.

preprint2014arXiv

Triangular Constellations in Fractal Measures

The local structure of a fractal set is described by its dimension $D$, which is the exponent of a power-law relating the mass ${\cal N}$ in a ball to its radius $ε$: ${\cal N}\sim ε^D$. It is desirable to characterise the {\em shapes} of constellations of points sampling a fractal measure, as well as their masses. The simplest example is the distribution of shapes of triangles formed by triplets of points, which we investigate for fractals generated by chaotic dynamical systems. The most significant parameter describing the triangle shape is the ratio $z$ of its area to the radius of gyration squared. We show that the probability density of $z$ has a phase transition: $P(z)$ is independent of $ε$ and approximately uniform below a critical flow compressibility $β_{\rm c}$, but for $β>β_{\rm c}$ it is described by two power laws: $P(z)\sim z^{α_1}$ when $1\gg z\gg z_{\rm c}(ε)$, and $P(z)\sim z^{α_2}$ when $z\ll z_{\rm c}(ε)$.

preprint2013arXiv

A model for the shapes of advected triangles

Three particles floating on a fluid surface define a triangle. The aim of this paper is to characterise the shape of the triangle, defined by two of its angles, as the three vertices are subject to a complex or turbulent motion. We consider a simple class of models for this process, involving a combination of a random strain of the fluid and Brownian motion of the particles. Following D. G. Kendall, we map the space of triangles to a sphere, whose equator corresponds to degenerate triangles with colinear vertices, with equilaterals at the poles. We map our model to a diffusion process on the surface of the sphere and find an exact solution for the shape distribution. Whereas the action of the random strain tends to make the shape of the triangles infinitely elongated, in the presence of a Brownian diffusion of the vertices, the model has an equilibrium distribution of shapes. We determine here exactly this shape distribution in the simple case where the increments of the strain are diffusive.

preprint2013arXiv

Multiple collisions in turbulent flows

In turbulent suspensions, collision rates determine how rapidly particles coalesce or react with each other. To determine the collision rate, many numerical studies rely on the &#39;Ghost Collision Approximation&#39; (GCA), which simply records how often pairs of point particles come within a threshold distance. In many applications, the suspended particles stick (or in the case of liquid droplets, coalesce) upon collision, and it is the frequency of first contact which is of interest. If a pair of &#39;ghost&#39; particles undergoes multiple collisions, the GCA may overestimate the true collision rate. Here, using fully resolved Direct Numerical Simulations of turbulent flows at moderate Reynolds number (R_λ= 130), we investigate the prevalence and properties of multiple collisions. We demonstrate that the GCA leads to a systematic overestimate of the collision rate, which is of the order of 15% when the particle inertia is small, and slowly decreases when inertia increases. We investigate the probability P(N) for a given pair of ghost particles colliding N times. We find P(N)=βα^N for N_>1, where αand βare coefficients which depend upon the particle inertia. This result is used to explain the discrepancy between the GCA and the true collision rates. We also investigate the statistics of the times that ghost particles remain in contact. We show that the probability density function of the contact time is different for the first collision. The difference is explained by the effect of caustics in the phase space of the suspended particles.

preprint2013arXiv

Prevalence of the sling effect for enhancing collision rates in turbulent suspensions

Turbulence facilitates collisions between particles suspended in a turbulent flow. Two effects have been proposed which can enhance the collision rate at high turbulence intensities: &#39;preferential concentration&#39; (a clustering phenomenon) and the &#39;sling effect&#39; (arising from the formation of caustic folds in the phase-space of the suspended particles). We have determined numerically the collision rate of small heavy particles as a function of their size and densities. The dependence on particle densities reveals that the enhancement by turbulence of the collision rate of particles with significant inertia is due almost entirely to the sling effect.

preprint2012arXiv

Ostwald ripening and the kinetics of rain initiation

A central problem in cloud physics is understanding the kinetics of the growth of water droplets. It is believed that there exists a &#39;condensation-coalescence bottleneck&#39; in the growth of intermediate size droplets. Here the Lifshitz-Slezov theory of Ostwald ripening is applied to the kinetics of the growth of rain drops. The analysis shows that kinetic barriers to rain initiation are not significant. It is also shown that Ostwald ripening can greatly enhance the apparent collision efficiency of droplets falling under gravity.

preprint2012arXiv

Super-diffusion in optical realizations of Anderson localization

We discuss the dynamics of particles in one dimension in potentials that are random both in space and in time. The results are applied to recent optics experiments on Anderson localization, in which the transverse spreading of a beam is suppressed by random fluctuations in the refractive index. If the refractive index fluctuates along the direction of the paraxial propagation of the beam, the localization is destroyed. We analyze this broken localization, in terms of the spectral decomposition of the potential. When the potential has a discrete spectrum, the spread is controlled by the overlap of Chirikov resonances in phase space. As the number of Fourier components is increased, the resonances merge into a continuum, which is described by a Fokker-Planck equation. We express the diffusion coefficient in terms of the spectral intensity of the potential. For a general class of potentials that are commonly used in optics, the solutions of the Fokker-Planck equation exhibit anomalous diffusion in phase space, implying that when Anderson localization is broken by temporal fluctuations of the potential, the result is transport at a rate similar to a ballistic one or even faster. For a class of potentials which arise in some existing realizations of Anderson localization atypical behavior is found.

preprint2012arXiv

Universal anomalous diffusion of weakly damped particles

We show that anomalous diffusion arises in two different models for the motion of randomly forced and weakly damped particles: one is a generalisation of the Ornstein-Uhlenbeck process with a random force which depends on position as well as time, the other is a generalisation of the Chandrasekhar-Rosenbluth model of stellar dynamics, encompassing non-Coulombic potentials. We show that both models exhibit anomalous diffusion of position $x$ and momentum $p$ with the same exponents: $<x^2> \sim C_x t^2$ and $<p^2> \sim C_p t^{2/5}$. We are able to determine the prefactors $C_x$, $C_p$ analytically.

preprint2011arXiv

Spherical Ornstein-Uhlenbeck processes

The paper considers random motion of a point on the surface of a sphere, in the case where the angular velocity is determined by an Ornstein-Uhlenbeck process. The solution is fully characterized by only one dimensionless number, the persistence angle, which is the typical angle of rotation of the object during the correlation time of the angular velocity. We first show that the two-dimensional case is exactly solvable. When the persistence angle is large, a series for the correlation function has the surprising property that its sum varies much more slowly than any of its individual terms. In three dimensions we obtain asymptotic forms for the correlation function, in the limits where the persistence angle is very small and very large. The latter case exhibits a complicated transient, followed by a much slower exponential decay. The decay rate is determined by the solution of a radial Schrodinger equation in which the angular momentum quantum number takes an irrational value, namely j=(sqrt{17}-1)/2. Possible applications of the model to objects tumbling in a turbulent environment are discussed.

preprint2010arXiv

Hot-Carrier Model for an Anomalous Exponent of Photoconduction

Experiments often show that the photoconductance $σ$ of a semiconductor system and the light intensity $I$ are related by $σ\sim I^γ$. Conventional theories give a satisfactory explanation for $γ=1$ or $γ=1/2$, but anamalous exponents close to $γ=3/4$ are often observed. This paper argues that there is a universal anomalous regime for which $γ=3/4$ (or $γ=2/3$ in two-dimensional samples), resulting from the kinetics of electron-hole recombination being controlled by Coulombic attraction. Because the local electric fields are extremely high, the theory appeals to the &#39;hot-carrier&#39; model for transport in semiconductors.

preprint2009arXiv

Correlation Dimension of Inertial Particles in Random Flows

We obtain an implicit equation for the correlation dimension of dynamical systems in terms of an integral over a propagator. We illustrate the utility of this approach by evaluating the correlation dimension for inertial particles suspended in a random flow. In the limit where the correlation time of the flow field approaches zero, taking the short-time limit of the propagator enables the correlation dimension to be determined from the solution of a partial differential equation. We develop the solution as a power series in a dimensionless parameter which represents the strength of inertial effects.

preprint2009arXiv

Perturbation theory for a stochastic process with Ornstein-Uhlenbeck noise

The Ornstein-Uhlenbeck process may be used to generate a noise signal with a finite correlation time. If a one-dimensional stochastic process is driven by such a noise source, it may be analysed by solving a Fokker-Planck equation in two dimensions. In the case of motion in the vicinity of an attractive fixed point, it is shown how the solution of this equation can be developed as a power series. The coefficients are determined exactly by using algebraic properties of a system of annihilation and creation operators.

preprint2009arXiv

Poincare Indices of Rheoscopic Visualisations

Suspensions of small anisotropic particles, termed &#39;rheoscopic fluids&#39;, are used for flow visualisation. By illuminating the fluid with light of three different colours, it is possible to determine Poincare indices for vector fields formed by the longest axis of the particles. Because this vector field is non-oriented, half-integer Poincare indices are possible, and are observed experimentally. An exact solution for the direction vector appears to preclude the existence of topological singularities. However, we show that upon averaging over the random initial orientations of particles, singularities with half-integer Poincare index appear. We describe their normal forms.

preprint2008arXiv

Fingerprints of Random Flows?

We consider the patterns formed by small rod-like objects advected by a random flow in two dimensions. An exact solution indicates that their direction field is non-singular. However, we find from simulations that the direction field of the rods does appear to exhibit singularities. First, ` scar lines&#39; emerge where the rods abruptly change direction by $π$. Later, these scar lines become so narrow that they ` heal over&#39; and disappear, but their ends remain as point singularities, which are of the same type as those seen in fingerprints. We give a theoretical explanation for these observations.

preprint2008arXiv

Planet Formation by Concurrent Collapse

After reviewing the difficulties faced by the conventional theory of planet formation (based upon the aggregation of microscopic dust particles), we describe an alternative hypothesis. We propose that planets form by gravitational collapse at the same time as the star about which they orbit. This &#39;concurrent collapse&#39; hypothesis avoids theoretical difficulties associated with the conventional model and suggests satisfying explanations for various poorly understood phenomena. We introduce new explanations for FU Orionis outbursts seen in young stars, the discovery of exoplanets with eccentric orbits and the existence of small rocky objects such as chondrules in the solar system.

preprint2000arXiv

Absorption of Energy at a Metallic Surface due to a Normal Electric Field

The effect of an oscillating electric field normal to a metallic surface may be described by an effective potential. This induced potential is calculated using semiclassical variants of the random phase approximation (RPA). Results are obtained for both ballistic and diffusive electron motion, and for two and three dimensional systems. The potential induced within the surface causes absorption of energy. The results are applied to the absorption of radiation by small metal spheres and discs. They improve upon an earlier treatment which used the Thomas-Fermi approximation for the effective potential.

preprint1996arXiv

Distribution of Oscillator Strengths for Recombination of Localised Excitons in Two Dimensions

We investigate the distribution of oscillator strengths for the recombination of excitons in a two dimensional sample, trapped in local minima of the confinement potential: the results are derived from a statistical topographic model of the potential. The predicted distribution of oscillator strengths is very different from the Porter-Thomas disribution which usually characterises disordered systems, and is notable for the fact that small oscillator strengths are extremely rare.