Sculplexity: Sculptures of Complexity using 3D printing

We show how to convert models of complex systems such as 2D cellular automata into a 3D printed object. Our method takes into account the limitations inherent to 3D printing processes and materials. Our approach automates the greater part of this task, bypassing the use of CAD software and the need for manual design. As a proof of concept, a physical object representing a modified forest fire model was successfully printed. Automated conversion methods similar to the ones developed here can be used to create objects for research, for demonstration and teaching, for outreach, or simply for aesthetic pleasure. As our outputs can be touched, they may be particularly useful for those with visual disabilities.Comment: Free access to article on European Physics Letter

Forecasting in the presence of recent structural change

We examine how to forecast after a recent break. We consider monitoring for change and then combining forecasts from models that do and do not use data before the change; and robust methods, namely rolling regressions, forecast averaging over different windows and exponentially weighted moving average (EWMA) forecasting. We derive analytical results for the performance of the robust methods relative to a full-sample recursive benchmark. For a location model subject to stochastic breaks the relative MSFE ranking is EWMA < rolling regression < forecast averaging. No clear ranking emerges under deterministic breaks. In Monte Carlo experiments forecast averaging improves performance in many cases with little penalty where there are small or infrequent changes. Similar results emerge when we examine a large number of UK and US macroeconomic series

Coupling iterated Kolmogorov diffusions

The Kolmogorov (1934) diffusion is the two-dimensional diffusion generated by real Brownian motion B and its time integral integral B d t. In this paper we construct successful co-adapted couplings for iterated Kolmogorov diffusions defined by adding iterated time integrals integral integral B d s d t,... as further components to the original Kolmogorov diffusion. A Laplace-transform argument shows it is not possible successfully to couple all iterated time integrals at once; however we give an explicit construction of a successful co-adapted coupling method for (B, integral B d t, integral integral B d s d t); and a more implicit construction of a successful co-adapted coupling method which works for finite sets of iterated time integrals

Ridge Fusion in Statistical Learning

We propose a penalized likelihood method to jointly estimate multiple precision matrices for use in quadratic discriminant analysis and model based clustering. A ridge penalty and a ridge fusion penalty are used to introduce shrinkage and promote similarity between precision matrix estimates. Block-wise coordinate descent is used for optimization, and validation likelihood is used for tuning parameter selection. Our method is applied in quadratic discriminant analysis and semi-supervised model based clustering.Comment: 24 pages and 9 tables, 3 figure

Is the dust-to-gas ratio constant in molecular clouds?

We perform numerical simulations of dusty, supersonic turbulence in molecular clouds. We model 0.1, 1 and 10 {\mu}m sized dust grains at an initial dust-to-gas mass ratio of 1:100, solving the equations of combined gas and dust dynamics where the dust is coupled to the gas through a drag term. We show that, for 0.1 and 1 {\mu}m grains, the dust-to-gas ratio deviates by typically 10-20% from the mean, since the stopping time of the dust due to gas drag is short compared to the dynamical time. Contrary to previous findings, we find no evidence for orders of magnitude fluctuation in the dust-to-gas ratio for 0.1 {\mu}m grains. Larger, 10 {\mu}m dust grains may have dust-to-gas ratios increased by up to an order of magnitude locally. Both small (0.1 {\mu}m) and large ($\gtrsim$ 1 {\mu}m) grains trace the large-scale morphology of the gas, however we find evidence for 'size-sorting' of grains, where turbulence preferentially concentrates larger grains into dense regions. Size-sorting may help to explain observations of 'coreshine' from dark clouds, and why extinction laws differ along lines of sight through molecular clouds in the Milky Way compared to the diffuse interstellar medium.Comment: 6 pages, 4 figures, accepted for publication in MNRAS Letters, videos available at https://www.youtube.com/channel/UC7J6IDzQklFzKV3c6pBqxU

Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics

We present a constrained formulation of Dedner et al's hyperbolic/parabolic divergence cleaning scheme for enforcing the \nabla\dot B = 0 constraint in Smoothed Particle Magnetohydrodynamics (SPMHD) simulations. The constraint we impose is that energy removed must either be conserved or dissipated, such that the scheme is guaranteed to decrease the overall magnetic energy. This is shown to require use of conjugate numerical operators for evaluating \nabla\dot B and \nabla{\psi} in the SPMHD cleaning equations. The resulting scheme is shown to be stable at density jumps and free boundaries, in contrast to an earlier implementation by Price & Monaghan (2005). Optimal values of the damping parameter are found to be {\sigma} = 0.2-0.3 in 2D and {\sigma} = 0.8-1.2 in 3D. With these parameters, our constrained Hamiltonian formulation is found to provide an effective means of enforcing the divergence constraint in SPMHD, typically maintaining average values of h |\nabla\dot B| / |B| to 0.1-1%, up to an order of magnitude better than artificial resistivity without the associated dissipation in the physical field. Furthermore, when applied to realistic, 3D simulations we find an improvement of up to two orders of magnitude in momentum conservation with a corresponding improvement in numerical stability at essentially zero additional computational expense.Comment: 28 pages, 25 figures, accepted to J. Comput. Phys. Movies at http://www.youtube.com/playlist?list=PL215D649FD0BDA466 v2: fixed inverted figs 1,4,6, and several color bar

Body shape preferences: Associations with rater body shape and sociosexuality

Copyright @ 2013 Price et al.This article has been made available through the Brunel Open Access Publishing Fund.There is accumulating evidence of condition-dependent mate choice in many species, that is, individual preferences varying in strength according to the condition of the chooser. In humans, for example, people with more attractive faces/bodies, and who are higher in sociosexuality, exhibit stronger preferences for attractive traits in opposite-sex faces/bodies. However, previous studies have tended to use only relatively simple, isolated measures of rater attractiveness. Here we use 3D body scanning technology to examine associations between strength of rater preferences for attractive traits in opposite-sex bodies, and raters' body shape, self-perceived attractiveness, and sociosexuality. For 118 raters and 80 stimuli models, we used a 3D scanner to extract body measurements associated with attractiveness (male waist-chest ratio [WCR], female waist-hip ratio [WHR], and volume-height index [VHI] in both sexes) and also measured rater self-perceived attractiveness and sociosexuality. As expected, WHR and VHI were important predictors of female body attractiveness, while WCR and VHI were important predictors of male body attractiveness. Results indicated that male rater sociosexuality scores were positively associated with strength of preference for attractive (low) VHI and attractive (low) WHR in female bodies. Moreover, male rater self-perceived attractiveness was positively associated with strength of preference for low VHI in female bodies. The only evidence of condition-dependent preferences in females was a positive association between attractive VHI in female raters and preferences for attractive (low) WCR in male bodies. No other significant associations were observed in either sex between aspects of rater body shape and strength of preferences for attractive opposite-sex body traits. These results suggest that among male raters, rater self-perceived attractiveness and sociosexuality are important predictors of preference strength for attractive opposite-sex body shapes, and that rater body traits -with the exception of VHI in female raters- may not be good predictors of these preferences in either sex.This study is funded by the Brunel University School of Social Sciences and made available through the Brunel Open Access Publishing Fund

Surface roughness detector Patent

Roughness detector for recording surface pattern of irregularitie