Tree-code simulations of proton acceleration from laser-irradiated wire targets

Recent experiments using Terawatt lasers to accelerate protons deposited on thin wire targets are modeled with a new type of gridless plasma simulation code. In contrast to conventional mesh-based methods, this technique offers a unique capability in emulating the complex geometry and open-ended boundary conditions characteristic of contemporary experimental conditions. Comparisons of ion acceleration are made between the tree code and standard particle-in-cell simulations for a typical collisionless

Radiochromic film spectroscopy of laser-accelerated proton beams using the FLUKA code and dosimetry traceable to primary standards

A new approach to spectroscopy of laser induced proton beams using radiochromic film (RCF) is presented. This approach allows primary standards of absorbed dose-to-water as used in radiotherapy to be transferred to the calibration of GafChromic HD-810 and EBT in a 29 MeV proton beam from the Birmingham cyclotron. These films were then irradiated in a common stack configuration using the TARANIS Nd:Glass multi-terawatt laser at Queens University Belfast, which can accelerate protons to 10-12 MeV, and a depth-dose curve was measured from a collimated beam. Previous work characterizing the relative effectiveness (RE) of GafChromic film as a function of energy was implemented into Monte Carlo depth-dose curves using FLUKA. A Bragg peak (BP) "library" for proton energies 0-15 MeV was generated, both with and without the RE function. These depth-response curves were iteratively summed in a FORTRAN routine to solve for the measured RCF depth-dose using a simple direct search algorithm. By comparing resultant spectra with both BP libraries, it was found that the effect of including the RE function accounted for an increase in the total number of protons by about 50%. To account for the energy loss due to a 20 mu m aluminum filter in front of the film stack, FLUKA was used to create a matrix containing the energy loss transformations for each individual energy bin. Multiplication by the pseudo-inverse of this matrix resulted in "up-shifting" protons to higher energies. Applying this correction to two laser shots gave further increases in the total number of protons, N of 31% and 56%. Failure to consider the relative response of RCF to lower proton energies and neglecting energy losses in a stack filter foil can potentially lead to significant underestimates of the total number of protons in RCF spectroscopy of the low energy protons produced by laser ablation of thin targets.

nero- a post-maximum supernova radiation transport code

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a non-local thermodynamic equilibrium (NLTE) treatment of radiation transport has rarely been studied. In this paper, we present a new SN radiation transport code, nero, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare nero to the radiation transport code of Jerkstrand, Fransson & Kozma and to the nebular code of Mazzali et al. All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, nero is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' Type Ia supernova (SN Ia) 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post-maximum epochs. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

Accelerating code on multi-cores with FastFlow

FastFlow is a programming framework specifically targeting cache-coherent shared-memory multi-cores. It is implemented as a stack of C++ template libraries built on top of lock-free (and memory fence free) synchronization mechanisms. Its philosophy is to combine programmability with performance. In this paper a new FastFlow programming methodology aimed at supporting parallelization of existing sequential code via offloading onto a dynamically created software accelerator is presented. The new methodology has been validated using a set of simple micro-benchmarks and some real applications. © 2011 Springer-Verlag.

Blown Jet Vortex Generator Control of a Separated Diffuser Flow

The application of blown jet vortex generators to control flow separation in a diffuser with an opening angle of 10° has been studied using the computational fluid dynamics (CFD) code Fluent 6™. Experimental data is available for the uncontrolled flow in the diffuser. The section of the duct upstream of the diffuser has a height H equal to 15 mm; its length and breadth are 101H and 41H respectively; the diffuser has an expansion ratio of 4.7:1. Fully developed flow is achieved upstream of the diffuser. Pipes of diameters equal to 1.5%, 2.5% and 5% of H were considered; pitch angle was constant at 45° and yaw angle was fixed at 60°; velocity ratio was varied from 1.7 to 8.0; both co-rotating and counter-rotating arrays were studied. The best results were obtained with a counter-rotating array of generators with a hole diameter of 5% of H and a velocity ratio of 3.7.

Impact of Wind Generation Mix on Transient Stability for an Interconnected Power System

The rapid growth of wind generation in many European countries is pushing power systems into
uncharted territory. As additional wind generators are installed, the changing generation mix may
impact on power system stability. This paper adopts the New England 39 bus system as a test
system for transient stability analysis. Thermal generator models are based on a likely future plant
mix for existing systems, while varying capacities of fixed-speed induction generators (FSIG) and
doubly-fed induction generators (DFIG) are considered. The main emphasis here has been placed
on the impact of wind technology mix on inter-area oscillations following transient grid
disturbances. In addition, both rotor angle stability and transient voltage stability are examined, and
results are compared with current grid code requirements and standards. Results have shown that
FSIGs can reduce tie-line oscillations and improve damping following a transient disturbance, but
they also cause voltage stability and rotor angle stability problems at high wind penetrations. In
contrast, DFIGs can improve both voltage and rotor angle stability, but their power output
noticeably oscillates during disturbances.

Making heads or tails of probability:an experiment with random generators

The equiprobability bias is a tendency for individuals to think of probabilistic events as 'equiprobable' by nature, and to judge outcomes that occur with different probabilities as equally likely. The equiprobability bias has been repeatedly found to be related to formal education in statistics, and it is claimed to be based on a misunderstanding of the concept of randomness.

A spectral synthesis code for rapid modelling of supernovae

We present TARDIS-an open-source code for rapid spectral modelling of supernovae (SNe). Our goal is to develop a tool that is sufficiently fast to allow exploration of the complex parameter spaces of models for SN ejecta. This can be used to analyse the growing number of highquality SN spectra being obtained by transient surveys. The code uses Monte Carlo methods to obtain a self-consistent description of the plasma state and to compute a synthetic spectrum. It has a modular design to facilitate the implementation of a range of physical approximations that can be compared to assess both accuracy and computational expediency. This will allow users to choose a level of sophistication appropriate for their application. Here, we describe the operation of the code and make comparisons with alternative radiative transfer codes of differing levels of complexity (SYN++, PYTHON and ARTIS). We then explore the consequence of adopting simple prescriptions for the calculation of atomic excitation, focusing on four species of relevance to Type Ia SN spectra-Si II, SII, MgII and Ca II. We also investigate the influence of three methods for treating line interactions on our synthetic spectra and the need for accurate radiative rate estimates in our scheme.

A new method for transmission loss allocation considering the circulating currents between generators

A new method is presented for transmission loss allocation based on the separation of transmission loss caused by load and the loss due to circulating currents between generators. The theoretical basis for and derivation of the loss formulae are presented using simple systems. The concept is then extended to a general power system using the Ybus model. Details of the application of the proposed method to a typical power system are presented along with results from the IEEE 30 bus test system. The results from both the small system and the standard IEEE test system demonstrate the validity of the proposed method.