SquidLab—A user-friendly program for background subtraction and fitting of magnetization data

We present an open-source program free to download for academic use with a full user-friendly graphical interface for performing flexible and robust background subtraction and dipole fitting on magnetization data. For magnetic samples with small moment sizes or sample environments with large or asymmetric magnetic backgrounds, it can become necessary to separate background and sample contributions to each measured raw voltage measurement before fitting the dipole signal to extract magnetic moments. Originally designed for use with pressure cells on a Quantum Design MPMS3 SQUID magnetometer, SquidLab is a modular object-oriented platform implemented in Matlab with a range of importers for different widely available magnetometer systems (including MPMS, MPMS-XL, MPMS-IQuantum, MPMS3, and S700X models) and has been tested with a broad variety of background and signal types. The software allows background subtraction of baseline signals, signal preprocessing, and performing fits to dipole data using Levenberg–Marquardt non-linear least squares or a singular value decomposition linear algebra algorithm that excels at picking out noisy or weak dipole signals. A plugin system allows users to easily extend the built-in functionality with their own importers, processes, or fitting algorithms. SquidLab can be downloaded, under Academic License, from the University of Warwick depository (wrap.warwick.ac.uk/129665)

Goal-oriented a posteriori error estimation for the travel time functional in porous media flows

In this article we consider the a posteriori error estimation and adaptive mesh refinement for the numerical approximation of the travel time functional arising in porous media flows. The key application of this work is in the safety assessment of radioactive waste facilities; in this setting, the travel time functional measures the time taken for a non-sorbing radioactive solute, transported by groundwater, to travel from a potential site deep underground to the biosphere. To ensure the computability of the travel time functional, we employ a mixed formulation of Darcy's law and conservation of mass, together with Raviart-Thomas H(div) conforming finite elements. The proposed a posteriori error bound is derived based on a variant of the standard Dual-Weighted-Residual approximation, which takes into account the lack of smoothness of the underlying functional of interest. The proposed adaptive refinement strategy is tested on both a simple academic test case and a problem based on the geological units found at the Sellafield site in the UK

Direct imaging of correlated defect nanodomains in a metal-organic framework

Defect engineering can enhance key properties of metal-organic frameworks (MOFs). Tailoring the distribution of defects, for example in correlated nanodomains, requires characterization across length scales. However, a critical nanoscale characterization gap has emerged between the bulk diffraction techniques used to detect defect nanodomains and the sub-nanometer imaging used to observe individual defects. Here, we demonstrate that the emerging technique of scanning electron diffraction (SED) can bridge this gap uniquely enabling both nanoscale crystallographic analysis and the lowdose formation of multiple diffraction contrast images for defect analysis in MOFs. We directly image defect nanodomains in the MOF UiO-66(Hf) over an area of ca. 1 000 nm and with a spatial resolution ca. 5 nm to reveal domain morphology and distribution. Based on these observations, we suggest possible crystal growth processes underpinning synthetic control of defect nanodomains. We also identify likely dislocations and small angle grain boundaries, illustrating that SED could be a key technique in developing the potential for engineering the distribution of defects, or “microstructure”, in functional MOF design

Collimated Outflow Formation via Binary Stars. 3-D Simulations of AGB Wind and Disk Wind Interactions

We present three-dimensional hydrodynamic simulations of the interaction of a slow wind from an asymptotic giant branch(AGB) star and a jet blown by an orbiting companion. The jet or "Collimated Fast Wind" is assumed to originate from an accretion disk which forms via Bondi accretion of the AGB wind or Roche lobe overflow. We present two distinct regimes in the wind-jet interaction determined by the ratio of the AGB wind to jet momentum flux. Our results show that when the wind momentum flux overwhelms the flux in the jet a more dis-ordered outflow outflow results with the jet assuming a corkscrew pattern and multiple shock structures driven into the AGB wind. In the opposite regime the jet dominates and will drive a highly collimated narrow waisted outflow. We compare our results with scenarios described by Soker & Rappaport (2000) and extrapolate the structures observed in PNe and Symbiotic stars.Comment: 22 pages, 8 figures, submitted to Ap

Equilibrium states and chaos in an oscillating double-well potential

We investigate numerically parametrically driven coupled nonlinear Schrödinger equations modeling the dynamics of coupled wave fields in a periodically oscillating double-well potential. The equations describe, among other things, two coupled periodically curved optical waveguides with Kerr nonlinearity or Bose-Einstein condensates in a double-well potential that is shaken horizontally and periodically in time. In particular, we study the persistence of equilibrium states of the undriven system due to the presence of the parametric drive. Using numerical continuations of periodic orbits and calculating the corresponding Floquet multipliers, we find that the drive can (de)stabilize a continuation of an equilibrium state indicated by the change in the (in)stability of the orbit, showing that parametric drives can provide a powerful control to nonlinear (optical- or matter-wave-) field tunneling. We also discuss the appearance of chaotic regions reported in previous studies that is due to destabilization of a periodic orbit. Analytical approximations based on an averaging method are presented. Using perturbation theory, the influence of the drive on the symmetry-breaking bifurcation point is analyzed. © 2014 American Physical Society

A Dimension-Adaptive Multi-Index Monte Carlo Method Applied to a Model of a Heat Exchanger

We present an adaptive version of the Multi-Index Monte Carlo method, introduced by Haji-Ali, Nobile and Tempone (2016), for simulating PDEs with coefficients that are random fields. A classical technique for sampling from these random fields is the Karhunen-Lo\`eve expansion. Our adaptive algorithm is based on the adaptive algorithm used in sparse grid cubature as introduced by Gerstner and Griebel (2003), and automatically chooses the number of terms needed in this expansion, as well as the required spatial discretizations of the PDE model. We apply the method to a simplified model of a heat exchanger with random insulator material, where the stochastic characteristics are modeled as a lognormal random field, and we show consistent computational savings

Numerical and Experimental Investigation of Circulation in Short Cylinders

In preparation for an experimental study of magnetorotational instability (MRI) in liquid metal, we explore Couette flows having height comparable to the gap between cylinders, centrifugally stable rotation, and high Reynolds number. Experiments in water are compared with numerical simulations. Simulations show that endcaps corotating with the outer cylinder drive a strong poloidal circulation that redistributes angular momentum. Predicted azimuthal flow profiles agree well with experimental measurements. Spin-down times scale with Reynolds number as expected for laminar Ekman circulation; extrapolation from two-dimensional simulations at $Re\le 3200$ agrees remarkably well with experiment at $Re\sim 10^6$. This suggests that turbulence does not dominate the effective viscosity. Further detailed numerical studies reveal a strong radially inward flow near both endcaps. After turning vertically along the inner cylinder, these flows converge at the midplane and depart the boundary in a radial jet. To minimize this circulation in the MRI experiment, endcaps consisting of multiple, differentially rotating rings are proposed. Simulations predict that an adequate approximation to the ideal Couette profile can be obtained with a few rings

Can digital discussion support tools provide cost-effective options for agricultural extension services?

Agricultural extension that delivers timely, targeted, and cost-effective support to farmers will help ensure the sustainability and adaptive capacity of agriculture, enhancing both food security and environmental security. Leveraging advances in agriclimate science and adult education, innovative digital technologies offer significant new opportunities to engage with farmers and to support decision making. In this study, animated video clips (machinimas), developed using the Second LifeTM virtual world gaming platform, model conversations around climate risk and critical on-farm decisions in the Australian sugarcane farming industry. Early evaluation indicates that this is an engaging format that promotes discussion by leveraging farmers’ natural modes of information gathering and social learning. Comparison with conventional extension practices indicates that these discussion support tools may be a cost effective addition to existing approaches. The format’s flexibility means machinimas are readily updated with new information and customized to meet the needs of different farmer groups. Rapid growth in digital access globally and the scalability of such approaches promise greater equity of access to high-value information, critical to better risk management decision making, at minimal cost, for millions of farmers

Observational Implications of Precessing Protostellar Discs and Jets

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.Comment: MNRAS, in press. 10 pages. Also available at http://www.ast.cam.ac.uk/~mbat