Accuracy of FEM and BEM for electromagnetic flux calculations in high power applications

The development and design of electric high power devices with electromagnetic computer-aided engineering (EM-CAE) software such as the Finite Element Method (FEM) and Boundary Element Method (BEM) has been widely adopted. This paper presents the analysis of a Fault Current Limiter (FCL), which acts as a high-voltage surge protector for power grids. A prototype FCL was built. The magnetic flux in the core and the resulting electromagnetic forces in the winding of the FCL were analyzed using both FEM and BEM. An experiment on the prototype was conducted in a laboratory. The data obtained from the experiment is compared to the numerical solutions to determine the suitability and accuracy of the two methods.

The development of a Monte Carlo system to verify radiotherapy treatment dose calculations

Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the “gold standard” for predicting dose deposition in the patient. In this study, software has been developed that enables the transfer of treatment plan information from the treatment planning system to a Monte Carlo dose calculation engine. A database of commissioned linear accelerator models (Elekta Precise and Varian 2100CD at various energies) has been developed using the EGSnrc/BEAMnrc Monte Carlo suite. Planned beam descriptions and CT images can be exported from the treatment planning system using the DICOM framework. The information in these files is combined with an appropriate linear accelerator model to allow the accurate calculation of the radiation field incident on a modelled patient geometry. The Monte Carlo dose calculation results are combined according to the monitor units specified in the exported plan. The result is a 3D dose distribution that could be used to verify treatment planning system calculations. The software, MCDTK (Monte Carlo Dicom ToolKit), has been developed in the Java programming language and produces BEAMnrc and DOSXYZnrc input files, ready for submission on a high-performance computing cluster. The code has been tested with the Eclipse (Varian Medical Systems), Oncentra MasterPlan (Nucletron B.V.) and Pinnacle3 (Philips Medical Systems) planning systems. In this study the software was validated against measurements in homogenous and heterogeneous phantoms. Monte Carlo models are commissioned through comparison with quality assurance measurements made using a large square field incident on a homogenous volume of water. This study aims to provide a valuable confirmation that Monte Carlo calculations match experimental measurements for complex fields and heterogeneous media.

C-BN Single-Walled Nanotubes from Hybrid Connection of BN/C Nanoribbons: Prediction by ab initio Density Functional Calculations

We demonstrated for the first time by ab initio density functional calculation and molecular dynamics simulation that C0.5(BN)0.5 armchair single-walled nanotubes (NT) are gapless semiconductors and can be spontaneously formed via the hybrid connection of graphene/BN Nanoribbons (GNR/BNNR) at room temperature. The direct synthesis of armchair C0.5(BN)0.5 via the hybrid connection of GNR/BNNR is predicted to be both thermodynamically and dynamically stable. Such novel armchair C0.5(BN)0.5 NTs possess enhanced conductance as that observed in GNRs. Additionally, the zigzag C0.5(BN)0.5 SWNTs are narrow band gap semiconductors, which may have potential application for light emission. In light of recent experimental progress and the enhanced degree of control in the synthesis of GNRs and BNNR, our results highlight an interesting avenue for synthesizing a novel specific type of C0.5(BN)0.5 nanotube (gapless or narrow direct gap semiconductor), with potentially important applications in BNC-based nanodevices.

Multifunctional porous graphene for nanoelectronics and hydrogen storage : new properties revealed by first principle calculations

The lack of an obvious “band gap” is a formidable hurdle for making a nanotransistor from graphene. Here, we use density functional calculations to demonstrate for the first time that porosity such as evidenced in recently synthesized porous graphene (http://www.sciencedaily.com/releases/2009/11/091120084337.htm) opens a band gap. The size of the band gap (3.2 eV) is comparable to most popular photocatalytic titania and graphitic C3N4 materials. In addition, the adsorption of hydrogen on Li-decorated porous graphene is much stronger than that in regular Li-doped graphene due to the natural separation of Li cations, leading to a potential hydrogen storage gravimetric capacity of 12 wt %. In light of the most recent experimental progress on controlled synthesis, these results uncover new potential for the practical application of porous graphene in nanoelectronics and clean energy.

Formation of single-walled carbon nanotube via the interaction of graphene nanoribbons : ab initio density functional calculations

The interaction of bare graphene nanoribbons (GNRs) was investigated by ab initio density functional theory calculations with both the local density approximation (LDA) and the generalized gradient approximation (GGA). Remarkably, two bare 8-GNRs with zigzag-shaped edges are predicted to form an (8, 8) armchair single-wall carbon nanotube (SWCNT) without any obvious activation barrier. The formation of a (10, 0) zigzag SWCNT from two bare 10-GNRs with armchair-shaped edges has activation barriers of 0.23 and 0.61 eV for using the LDA and the revised PBE exchange correlation functional, respectively, Our results suggest a possible route to control the growth of specific types SWCNT via the interaction of GNRs.

Role of charge in destabilizing AlH4 and BH4 complex anions for hydrogen storage applications : Ab initio density functional calculations

NaAlH4 and LiBH4 are potential candidate materials for mobile hydrogen storage applications, yet they have the drawback of being highly stable and desorbing hydrogen only at elevated temperatures. In this letter, ab initio density functional theory calculations reveal how the stabilities of the AlH4 and BH4 complex anions will be affected by reducing net anionic charge. Tetrahedral AlH4 and BH4 complexes are found to be distorted with the decrease of negative charge. One H-H distance becomes smaller and the charge density will overlap between them at a small anion charge. The activation energies to release of H2 from AlH4 and BH4 complexes are thus greatly decreased. We demonstrate that point defects such as neutral Na vacancies or substitution of a Na atom with Ti on the NaAlH4(001) surface can potentially cause strong distortion of neighboring AlH4 complexes and even induce spontaneous dehydrogenation. Our results help to rationalize the conjecture that the suppression of charge transfer to AlH4 and BH4 anion as a consequence of surface defects should be very effective for improving the recycling performance of H2 in NaAlH4 and LiBH4. The understanding gained here will aid in the rational design and development of hydrogen storage materials based on these two systems.

Comparison of functionals for metal hexaboride band structure calculations

Density functional calculations of the electronic band structure for superconducting and semi-conducting metal hexaborides are compared using a consistent suite of assumptions and with emphasis on the physical implications of computed models. Spin polarization enhances mathematical accuracy of the functional approximations and adds significant physical meaning to model interpretation. For YB6 and LaB6, differences in alpha and beta projections occur near the Fermi energy. These differences are pronounced for superconducting hexaborides but do not occur for other metal hexaborides.

A comparison of surface doses for very small field size x-ray beams : Monte Carlo calculations and radiochromic film measurements

Stereotactic radiosurgery treatments involve the delivery of very high doses for a small number of fractions. To date, there is limited data in terms of the skin dose for the very small field sizes used in these treatments. In this work, we determine relative surface doses for small size circular collimators as used in stereotactic radiosurgery treatments. Monte Carlo calculations were performed using the BEAMnrc code with a model of the Novalis 15 Trilogy linear accelerator and the BrainLab circular collimators. The surface doses were calculated at the ICRU skin dose depth of 70 m all using the 6 MV SRS x-ray beam. The calculated surface doses varied between 15 – 12% with decreasing values as the field size increased from 4 to 30 mm. In comparison, surface doses were measured using Gafchromic EBT3 film positioned at the surface of a Virtual Water phantom. The absolute agreement between calculated and measured surface doses was better than 2.5% which is well within the 20 uncertainties of the Monte Carlo calculations and the film measurements. Based on these results, we have shown that the Gafchromic EBT3 film is suitable for surface dose estimates in very small size fields as used in SRS.

Spin polarized DFT calculations of metal borides

The metal borides, in particular the diborides and hexaborides, contain stoichiometric forms that include insulators, semiconductors and superconductors. In addition, their end-member structures have high symmetry and two atoms although, in general, substitution(s) of multi-valent ions into the metal site occurs consistent with Vegard’s law. These characteristics allow for fundamental comparison of important physical properties such as superconductivity and insulation within a relatively simple structure type. Our early work1,2 has demonstrated this for the hexaborides and this work compares similar attributes across a broader suite of boride structures. In all cases, theoretical calculations are referenced to structures determined via high resolution neutron or X-ray diffraction experiments.

Soil N2O emissions under N2-fixing legumes and N-fertilised canola: A reappraisal of emissions factor calculations

Introducing nitrogen (N)-fixing legumes into cereal-based crop rotations reduces synthetic fertiliser-N use and may mitigate soil emissions of nitrous oxide (N2O). Current IPCC calculations assume 100% of legume biomass N as the anthropogenic N input and use 1% of this as an emission factor (EF)—the percentage of input N emitted as N2O. However, legumes also utilise soil inorganic N, so legume-fixed N is typically less than 100% of legume biomass N. In two field experiments, we measured soil N2O emissions from a black Vertosol in sub-tropical Australia for 12 months after sowing of chickpea (Cicer arietinum L.), canola (Brassica napus L.), faba bean (Vicia faba L.), and field pea (Pisum sativum L.). Cumulative N2O emissions from N-fertilised canola (624 g N2O-N ha−1) greatly exceeded those from chickpea (127 g N2O-N ha−1) in Experiment 1. Similarly, N2O emitted from canola (385 g N2O-N ha−1) in Experiment 2 was significantly greater than chickpea (166 g N2O-N ha−1), faba bean (166 g N2O-N ha−1) or field pea (135 g N2O-N ha−1). Highest losses from canola were recorded during the growing season, whereas 75% of the annual N2O losses from the legumes occurred post-harvest. Legume N2-fixation provided 37–43% (chickpea), 54% (field pea) and 64% (faba bean) of total plant biomass N. Using only fixed-N inputs, we calculated EFs for chickpea (0.13–0.31%), field pea (0.18%) and faba bean (0.04%) that were significantly less than N-fertilised canola (0.48–0.78%) (P < 0.05), suggesting legume-fixed N is a less emissive form of N input to the soil than fertiliser N. Inputs of legume-fixed N should be more accurately quantified to properly gauge the potential for legumes to mitigate soil N2O emissions. EF’s from legume crops need to be revised and should include a factor for the proportion of the legume’s N derived from the atmosphere.

Electronic Structure of Vacancy Ordered Spinels, $GaMo_{4}S_{8}$ and $GaV_{4}S_{8}$, from ab Initio Calculations

We report ab initio calculations for the band dispersions and total as well as partial densities of states for vacancy ordered, clustered spinels, GaMo4S8 and GaV4S8. Results are presented for the high temperature cubic phase for both compounds. Additionally, we discuss results of similar calculations for GaMo4S8 in an idealized cubic structure, as well as the nonmagnetic and the ferromagnetic states of the low temperature rhombohedral structure. Comparison of these results allows us to discuss the unusual aspects of the electronic structure of this interesting class of compounds, and provide estimates of the crystal-field and exchange splitting strengths.

Flow behavior of long chain plasticizing liquids: Segment size calculations in carboxylate esters

This is the first comprehensive report on the calculation of segment size, which signifies the asic unit of flow in long chain plasticizing liquids, by a novel multi-pronged approach. Unlike,low molecular weight liquids and high polymer melts these complex long chain liquids encompasses the least understood domain of the liquid state. In the present work the flow behaviour of carboxylate ester (300-900 Da) has been explained through segmental motion taking into account the independence of molecular weight region. The segment size have been calculated by various methods based on satistical thermodynamics, molecular dynamics and group additivity nd their merits analysed.

An overview of land application of pig effluent-P using soil P chemistry and mass balance calculations.

Attention is directed at land application of piggery effluent (containing urine, faeces, water, and wasted feed) as a potential source of water resource contamination with phosphorus (P). This paper summarises P-related properties of soil from 0-0.05 m depth at 11 piggery effluent application sites, in order to explore the impact that effluent application has had on the potential for run-off transport of P. The sites investigated were situated on Alfisol, Mollisol, Vertisol, and Spodosol soils in areas that received effluent for 1.5-30 years (estimated effluent-P applications of 100-310000 kg P/ha in total). Total (PT), bicarbonate extractable (PB), and soluble P forms were determined for the soil (0-0.05 m) at paired effluent and no-effluent sites, as well as texture, oxalate-extractable Fe and Al, organic carbon, and pH. All forms of soil P at 0-0.05 m depth increased with effluent application (PB at effluent sites was 1.7-15 times that at no-effluent sites) at 10 of the 11 sites. Increases in PB were strongly related to net P applications (regression analysis of log values for 7 sites with complete data sets: 82.6 % of variance accounted for, p <0.01). Effluent irrigation tended to increase the proportion of soil PT in dilute CaCl2-extractable forms (PTC: effluent average 2.0 %; no-effluent average 0.6%). The proportion of PTC in non-molybdate reactive forms (centrifuged supernatant) decreased (no-effluent average, 46.4 %; effluent average, 13.7 %). Anaerobic lagoon effluent did not reliably acidify soil, since no consistent relationship was observed for pH with effluent application. Soil organic carbon was increased in most of the effluent areas relative to the no-effluent areas. The four effluent areas where organic carbon was reduced had undergone intensive cultivation and cropping. Current effluent management at many of the piggeries failed to maximise the potential for waste P recapture. Ten of the case-study effluent application areas have received effluent-P in excess of crop uptake. While this may not represent a significant risk of leaching where sorption retains P, it has increased the risk of transport of P by run-off. Where such sites are close to surface water, run-off P loads should be managed.