The development of parallel implementation for a CFD based fire field model utilising conventional office based PCs

Parallel processing techniques have been used in the past to provide high performance computing resources for activities such as fire-field modelling. This has traditionally been achieved using specialized hardware and software, the expense of which would be difficult to justify for many fire engineering practices. In this article we demonstrate how typical office-based PCs attached to a Local Area Network has the potential to offer the benefits of parallel processing with minimal costs associated with the purchase of additional hardware or software. It was found that good speedups could be achieved on homogeneous networks of PCs, for example a problem composed of ~100,000 cells would run 9.3 times faster on a network of 12 800MHz PCs than on a single 800MHz PC. It was also found that a network of eight 3.2GHz Pentium 4 PCs would run 7.04 times faster than a single 3.2GHz Pentium computer. A dynamic load balancing scheme was also devised to allow the effective use of the software on heterogeneous PC networks. This scheme also ensured that the impact between the parallel processing task and other computer users on the network was minimized.

Droplet generation with modulated AC electromagnetic field at nozzle exit

A method of droplet generation based on applying a modulated AC high frequency magnetic field in the localized region of capillary breakup is considered as ans alternative to traditional methods for high temperature liquid melt droplet serial production by pressure variation. The method is based on a pseudo-spectral approximation with a coordinate transformation adaptin to the developing free surface. The electromagnetic field is recomputed continuously with the domain shape change. Practical application cases for liquid silicon droplets of 0.5 - 2 mm diameter are considered in detail.

Levitated liquid droplets in AC and DC magnetic field

Electromagnetic levitation of liquid metal droplets can be used to measure the properties of highly reactive liquid materials. Two independent numerical models, the commercial COMSOL and the spectral-collocation based free surface code SPHINX, have been applied to solve the transient electromagnetic, fluid flow and thermodynamic equations, which describe the levitated liquid motion and heating processes. The SPHINX model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the electromagnetic and gravity forces, temperature dependent surface tension, magnetically controlled turbulent momentum transport. The models are adapted to incorporate periodic laser heating at the top of the droplet, which is used to measure the thermal conductivity of the material. Novel effects in the levitated droplet of magnetically damped turbulence and nonlinear growth of velocities in high DC magnetic field are analysed.

Simulation of mackerel (Scomber scombrus) recruitment with an individual-based model and comparison with field data

An individual-based model (IBM) for the simulation of year-to-year survival during the early life-history stages of the north-east Atlantic stock of mackerel (Scomber scombrus) was developed within the EU funded Shelf-Edge Advection, Mortality and Recruitment (SEAMAR) programme. The IBM included transport, growth and survival and was used to track the passive movement of mackerel eggs, larvae and post-larvae and determine their distribution and abundance after approximately 2 months of drift. One of the main outputs from the IBM, namely distributions and numbers of surviving post-larvae, are compared with field data as recruit (age-0/age-1 juveniles) distribution and abundance for the years 1998, 1999 and 2000. The juvenile distributions show more inter-annual and spatial variability than the modelled distributions of survivors; this may be due to the restriction of using the same initial egg distribution for all 3 yr of simulation. The IBM simulations indicate two main recruitment areas for the north-east Atlantic stock of mackerel, these being Porcupine Bank and the south-eastern Bay of Biscay. These areas correspond to areas of high juvenile catches, although the juveniles generally have a more widespread distribution than the model simulations. The best agreement between modelled data and field data for distribution (juveniles and model survivors) is for the year 1998. The juvenile catches in different representative nursery areas are totalled to give a field abundance index (FAI). This index is compared with a model survivor index (MSI) which is calculated from the total of survivors for the whole spawning season. The MSI compares favourably with the FAI for 1998 and 1999 but not for 2000; in this year, juvenile catches dropped sharply compared with the previous years but there was no equivalent drop in modelled survivors.

Estimation of atmospheric nutrient inputs to the Atlantic Ocean from 50°N to 50°S based on large-scale field sampling: Iron and other dust-associated elements

Atmospheric inputs of mineral dust supply iron and other trace metals to the remote ocean and can influence the marine carbon cycle due to iron's role as a potentially limiting micronutrient. Dust generation, transport, and deposition are highly heterogeneous, and there are very few remote marine locations where dust concentrations and chemistry (e.g., iron solubility) are routinely monitored. Here we use aerosol and rainwater samples collected during 10 large-scale research cruises to estimate the atmospheric input of iron, aluminum, and manganese to four broad regions of the Atlantic Ocean over two 3 month periods for the years 2001–2005. We estimate total inputs of these metals to our study regions to be 4.2, 17, and 0.27 Gmol in April–June and 4.9, 14, and 0.19 Gmol in September–November, respectively. Inputs were highest in regions of high rainfall (the intertropical convergence zone and South Atlantic storm track), and rainfall contributed higher proportions of total input to wetter regions. By combining input estimates for total and soluble metals for these time periods, we calculated overall percentage solubilities for each metal that account for the contributions from both wet and dry depositions and the relative contributions from different aerosol types. Calculated solubilities were in the range 2.4%–9.1% for iron, 6.1%–15% for aluminum, and 54%–73% for manganese. We discuss sources of uncertainty in our estimates and compare our results to some recent estimates of atmospheric iron input to the Atlantic.

Consistent increase in dimethyl sulphide (DMS) in response to high CO2 in five shipboard bioassays from contrasting NW European waters

The ubiquitous marine trace gas dimethyl sulphide (DMS) comprises the greatest natural source of sulphur to the atmosphere and is a key player in atmospheric chemistry and climate. We explore the short term response of DMS and its algal precursor dimethyl sulphoniopropionate (DMSP) production and cycling to elevated carbon dioxide (CO2) and ocean acidification (OA) in five highly replicated 96 h shipboard bioassay experiments from contrasting sites in NW European shelf waters. In general, the response to OA throughout this region showed little variation, despite encompassing a range of biological and biogeochemical conditions. We observed consistent and marked increases in DMS concentrations relative to ambient controls, and decreases in DMSP concentrations. Quantification of rates of specific DMSP synthesis by phytoplankton and bacterial DMS gross production/consumption suggest algal processes dominated the CO2 response, likely due to a physiological response manifested as increases in direct cellular exudation of DMS and/or DMSP lyase enzyme activities. The variables and rates we report increase our understanding of the processes behind the response to OA. This could provide the opportunity to improve upon mesocosm-derived empirical modelling relationships, and move towards a mechanistic approach for predicting future DMS concentrations.

Consistent increase in dimethyl sulfide (DMS) in response to high CO2 in five shipboard bioassays from contrasting NW European waters

The ubiquitous marine trace gas dimethyl sulfide (DMS) comprises the greatest natural source of sulfur to the atmosphere and is a key player in atmospheric chemistry and climate. We explore the short-term response of DMS production and cycling and that of its algal precursor dimethyl sulfoniopropionate (DMSP) to elevated carbon dioxide (CO₂) and ocean acidification (OA) in five 96 h shipboard bioassay experiments. Experiments were performed in June and July 2011, using water collected from contrasting sites in NW European waters (Outer Hebrides, Irish Sea, Bay of Biscay, North Sea). Concentrations of DMS and DMSP, alongside rates of DMSP synthesis and DMS production and consumption, were determined during all experiments for ambient CO₂ and three high-CO₂ treatments (550, 750, 1000 μatm). In general, the response to OA throughout this region showed little variation, despite encompassing a range of biological and biogeochemical conditions. We observed consistent and marked increases in DMS concentrations relative to ambient controls (110% (28–223%) at 550 μatm, 153% (56–295%) at 750 μatm and 225% (79–413%) at 1000 μatm), and decreases in DMSP concentrations (28% (18–40%) at 550 μatm, 44% (18–64%) at 750 μatm and 52% (24–72%) at 1000 μatm). Significant decreases in DMSP synthesis rate constants (μDMSP, d⁻¹) and DMSP production rates (nmol d⁻¹) were observed in two experiments (7–90% decrease), whilst the response under high CO₂ from the remaining experiments was generally indistinguishable from ambient controls. Rates of bacterial DMS gross consumption and production gave weak and inconsistent responses to high CO₂. The variables and rates we report increase our understanding of the processes behind the response to OA. This could provide the opportunity to improve upon mesocosm-derived empirical modelling relationships and to move towards a mechanistic approach for predicting future DMS concentrations.

Genetic maternity and paternity in a local population of armadillos assessed by microsatellite DNA markers and field data.

Genetic data from polymorphic microsatellite loci were employed to estimate paternity and maternity in a local population of nine-banded armadillos (Dasypus novemcinctus) in northern Florida. The parentage assessments took advantage of maximum likelihood procedures developed expressly for situations when individuals of neither gender can be excluded a priori as candidate parents. The molecular data for 290 individuals, interpreted alone and in conjunction with detailed biological and spatial information for the population, demonstrate high exclusion probabilities and reasonably strong likelihoods of genetic parentage assignment in many cases; low mean probabilities of successful reproductive contribution to the local population by individual armadillo adults in a given year; and statistically significant microspatial associations of parents and their offspring. Results suggest that molecular assays of highly polymorphic genetic systems can add considerable power to assessments of biological parentage in natural populations even when neither parent is otherwise known.

A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.

Physicochemical and mineralogical characteristics of soil/saprolite cores from a field research site, Tennessee

Site characterization is an essential initial step in determining the feasibility of remedial alternatives at hazardous waste sites. Physicochemical and mineralogical characterization of U-contaminated soils in deeply weathered saprolite at Area 2 of the DOE Field Research Center (FRC) site, Oak Ridge, TN, was accomplished to examine the feasibility of bioremediation. Concentrations of U in soil–saprolite (up to 291 mg kg–1 in oxalate-extractable Uo) were closely related to low pH (ca. 4–5), high effective cation exchange capacity without Ca (64.7–83.2 cmolc kg–1), amorphous Mn content (up to 9910 mg kg–1), and the decreased presence of relative clay mineral contents in the bulk samples (i.e., illite 2.5–12 wt. %, average 32 wt. %). The pH of the fill material ranged from 7.0 to 10.5, whereas the pH of the saprolite ranged from 4.5 to 8. Uranium concentration was highest (about 300 mg kg–1) at around 6 m below land surface near the saprolite–fill interface. The pH of ground water at Area 2 tended to be between 6 and 7 with U concentrations of about 0.9 to 1.7 mg L–1. These site specific characteristics of Area 2, which has lower U and nitrate contamination levels and more neutral ground water pH compared with FRC Areas 1 and 3 (ca. 5.5 and

HI observations of the high-velocity cloud in the direction of M 92

We present wide-field neutral hydrogen (H I) Lovell telescope multibeam, and Dominion Radio Astrophysical Observatory Hi synthesis observations, of the high velocity cloud (HVC) located in the general direction of the globular cluster M92. This cloud is part of the larger Complex C and lies at velocities between similar to -80 and -130 km s(-1) in the Local Standard of Rest. The Lovell telescope observations, of resolution 12 arcmin spatially and 3.0 km s(-1) in velocity, fully sampling a 3.1 degrees x 12.6 degrees RA-Dec grid, have found that this part of HVC Complex C comprises two main condensations, lying approximately north-south in declination, separated by similar to2 degrees and being parallel to the Galactic plane. At this resolution, peak values of the brightness temperature and Hi column density of similar to1.4 K and similar to5 x 10(19) cm(-2) are determined, with relatively high values of the full width half maximum velocity (FWHM) of similar to 22 km s(-1) being observed, equivalent to a gas kinetic temperature, in the absence of turbulence and geometric effects of similar to 10 000 K. Each of these properties, as well as the sizes of the clouds, are similar in the two components. The DRAO observations, towards the Northern HVC condensation, are the first high-resolution Hi spectra of Complex C. When smoothed to a resolution of 3 arcmin, they identify several Hi intensity peaks with column densities in the range 4-7 x 10(19) cm(-2). Further smoothing of these data to 6 arcmin resolution tentatively indicates that parts of the HVC consist of two velocity components, of similar brightness temperature, separated by similar to7 km s(-1) in velocity, and with FWHM velocity widths of similar to5-7 km s(-1). No IRAS 60 or 100 micron flux is associated with the M92 HVC. Cloud properties are briefly discussed and compared to previous observations of HVCs.

Strong parallel magnetic field effects on the hydrogen molecular ion

Equilibrium distances, binding energies and dissociation energies for the ground and low-lying states of the hydrogen molecular ion in a strong magnetic field parallel to the internuclear axis are calculated and refined, by using the two- dimensional pseudospectral method. High-precision results are presented for the binding energies over a wider field regime than already given in the literature (Kravchenko and Liberman 1997 Phys. Rev. A 55 2701). The present work removes a long- standing discrepancy for the R-eq value in the 1sigma(u) state at a field strength of 1.0 x 10(6) T. The dissociation energies of the antibonding 1pi(g) state induced by magnetic fields are determined accurately. We have also observed that the antibonding 1pi(g) potential energy curve develops a minimum if the field is sufficiently strong. Some unreliable results in the literature are pointed out and discussed. A way to efficiently treat vibrational processes and coupling between the nuclear and the electronic motions in magnetic fields is also suggested within a three-dimensional pseudospectral scheme.

Interpolation formula between very low and intermediate-to-high damping Kramers escape rates for single-domain ferromagnetic particles

It is shown that the Mel'nikov-Meshkov formalism for bridging the very low damping (VLD) and intermediate-to-high damping (IHD) Kramers escape rates as a function of the dissipation parameter for mechanical particles may be extended to the rotational Brownian motion of magnetic dipole moments of single-domain ferromagnetic particles in nonaxially symmetric potentials of the magnetocrystalline anisotropy so that both regimes of damping, occur. The procedure is illustrated by considering the particular nonaxially symmetric problem of superparamagnetic particles possessing uniaxial anisotropy subject to an external uniform field applied at an angle to the easy axis of magnetization. Here the Mel'nikov-Meshkov treatment is found to be in good agreement with an exact calculation of the smallest eigenvalue of Brown's Fokker-Planck equation, provided the external field is large enough to ensure significant departure from axial symmetry, so that the VLD and IHD formulas for escape rates of magnetic dipoles for nonaxially symmetric potentials are valid.