An open source, parallel DSMC code for rarefied gas flows in arbitrary geometries

This paper presents the results of validation of an open source Direct Simulation Monte Carlo (DSMC) code for general application to rarefied gas flows. The new DSMC code, called dsmcFoam, has been written within the framework of the open source C++ CFD toolbox OpenFOAM. The main features of dsmcFoam code include the capability to perform both steady and transient solutions, to model arbitrary 2D/3D geometries, and unlimited parallel processing. Test cases have been selected to cover a wide range of benchmark examples from 1D to 3D. These include relaxation to equilibrium, 2D flow over a flat plate and a cylinder, and 3D supersonic flows over complex geometries. In all cases, dsmcFoam shows very good agreement with data provided by both analytical solutions and other contemporary DSMC codes

Efficient History Matching of a High Dimensional Individual-Based HIV Transmission Model

History matching is a model (pre-)calibration method that has been applied to computer models from a wide range of scientific disciplines. In this work we apply history matching to an individual-based epidemiological model of HIV that has 96 input and 50 output parameters, a model of much larger scale than others that have been calibrated before using this or similar methods. Apart from demonstrating that history matching can analyze models of this complexity, a central contribution of this work is that the history match is carried out using linear regression, a statistical tool that is elementary and easier to implement than the Gaussian process--based emulators that have previously been used. Furthermore, we address a practical difficulty with history matching, namely, the sampling of tiny, nonimplausible spaces, by introducing a sampling algorithm adjusted to the specific needs of this method. The effectiveness and simplicity of the history matching method presented here shows that it is a useful tool for the calibration of computationally expensive, high dimensional, individual-based models

Bayesian history matching of complex infectious disease models using emulation: A tutorial and a case study on HIV in Uganda

Advances in scientific computing have allowed the development of complex models that are being routinely applied to problems in disease epidemiology, public health and decision making. The utility of these models depends in part on how well they can reproduce empirical data. However, fitting such models to real world data is greatly hindered both by large numbers of input and output parameters, and by long run times, such that many modelling studies lack a formal calibration methodology. We present a novel method that has the potential to improve the calibration of complex infectious disease models (hereafter called simulators). We present this in the form of a tutorial and a case study where we history match a dynamic, event-driven, individual-based stochastic HIV simulator, using extensive demographic, behavioural and epidemiological data available from Uganda. The tutorial describes history matching and emulation. History matching is an iterative procedure that reduces the simulator's input space by identifying and discarding areas that are unlikely to provide a good match to the empirical data. History matching relies on the computational efficiency of a Bayesian representation of the simulator, known as an emulator. Emulators mimic the simulator's behaviour, but are often several orders of magnitude faster to evaluate. In the case study, we use a 22 input simulator, fitting its 18 outputs simultaneously. After 9 iterations of history matching, a non-implausible region of the simulator input space was identified that was times smaller than the original input space. Simulator evaluations made within this region were found to have a 65% probability of fitting all 18 outputs. History matching and emulation are useful additions to the toolbox of infectious disease modellers. Further research is required to explicitly address the stochastic nature of the simulator as well as to account for correlations between outputs

Houdini: a remote mobile platform for tank waste retrieval tasks

RedZone has developed Houdini{trademark}, a folding frame vehicle for work in waste storage tanks and other confined-access areas. Houdini is a tethered, hydraulically-powered platform that folds to fit through small openings. Once deployed, the vehicle unfolds to provide a substantial work platform for the deployment of a wide variety of tools. The Houdini system will perform wheel removal, waste retrieval, waste mobilization, waste size reduction, and other tank waste retrieval and decommissioning tasks. Within the DOE Complex, 332 underground storage tanks have been used to process and store radioactive and chemical mixed waste generated from weapon materials production. The ultimate goal of the program is to develop and commercialize the Houdini system for broad application throughout the DOE Complex

Effect of nearest neighbor repulsion on the low frequency phase diagram of a quarter-filled Hubbard-Holstein chain

We have studied the influence of nearest-neighbor (NN) repulsion on the low frequency phase diagram of a quarter-filled Hubbard-Holstein chain. The NN repulsion term induces the apparition of two new long range ordered phases (one $4k_F$ CDW for positive $U_{eff} = U-2g^2/\omega$ and one $2k_F$ CDW for negative $U_{eff}$) that did not exist in the V=0 phase diagram. These results are put into perspective with the newly observed charge ordered phases in organic conductors and an interpretation of their origin in terms of electron-molecular vibration coupling is suggested.Comment: 10 pages, 10 figure

The James Clerk Maxwell Telescope Spectral Legacy Survey

Original article can be found at: http://www.journals.uchicago.edu/loi/pasp Copyright University of Chicago Press / AAS. DOI: 10.1086/511161Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes that can reveal the dynamics, physics, chemistry, and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the James Clerk Maxwell Telescope (JCMT) Board of Directors. Beginning in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 and 373 GHz) toward a sample of five sources. Our intended targets are a low-mass core (NGC 1333 IRAS 4), three high-mass cores spanning a range of star-forming environments and evolutionary states (W49, AFGL 2591, and IRAS 20126), and a photodissociation region (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS (Heterodyne Array Receiver Programme B/Auto- Correlation Spectrometer and Imaging System) to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments and to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral data cubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.Peer reviewe

Entangled-Photon Generation from Parametric Down-Conversion in Media with Inhomogeneous Nonlinearity

We develop and experimentally verify a theory of Type-II spontaneous parametric down-conversion (SPDC) in media with inhomogeneous distributions of second-order nonlinearity. As a special case, we explore interference effects from SPDC generated in a cascade of two bulk crystals separated by an air gap. The polarization quantum-interference pattern is found to vary strongly with the spacing between the two crystals. This is found to be a cooperative effect due to two mechanisms: the chromatic dispersion of the medium separating the crystals and spatiotemporal effects which arise from the inclusion of transverse wave vectors. These effects provide two concomitant avenues for controlling the quantum state generated in SPDC. We expect these results to be of interest for the development of quantum technologies and the generation of SPDC in periodically varying nonlinear materials.Comment: submitted to Physical Review

Optimal imaging of remote bodies using quantum detectors

We implement a general imaging method by measuring the complex degree of coherence using linear optics and photon number resolving detectors. In the absence of collective or entanglement-assisted measurements, our method is optimal over a large range of practically relevant values of the complex degree of coherence. We measure the size and position of a small distant source of pseudothermal light, and show that our method outperforms the traditional imaging method by an order of magnitude in precision. Finally, we show that a lack of photon-number resolution in the detectors has only a modest detrimental effect on measurement precision and simulate imaging using the new and traditional methods with an array of detectors, showing that the new method improves both image clarity and contrast